CRISPRi screen for enhancing heterologous α-amylase yield in Bacillus subtilis.

Yield improvements in cell factories can potentially be obtained by fine-tuning the regulatory mechanisms for gene candidates. In pursuit of such candidates, we performed RNA-sequencing of two α-amylase producing Bacillus strains and predict hundreds of putative novel non-coding transcribed regions. Surprisingly, we found among hundreds of non-coding and structured RNA candidates that non-coding genomic regions are proportionally undergoing the highest changes in expression during fermentation. Since these classes of RNA are also understudied, we targeted the corresponding genomic regions with CRIPSRi knockdown to test for any potential impact on the yield. From differentially expression analysis, we selected 53 non-coding candidates. Although CRISPRi knockdowns target both the sense and the antisense strand, the CRISPRi experiment cannot link causes for yield changes to the sense or antisense disruption. Nevertheless, we observed on several instances with strong changes in enzyme yield. The knockdown targeting the genomic region for a putative antisense RNA of the 3' UTR of the skfA-skfH operon led to a 21% increase in yield. In contrast, the knockdown targeting the genomic regions of putative antisense RNAs of the cytochrome c oxidase subunit 1 (ctaD), the sigma factor sigH, and the uncharacterized gene yhfT decreased yields by 31 to 43%.

The impact of PrsA over-expression on the Bacillus subtilis transcriptome during fed-batch fermentation of alpha-amylase production.

The production of the alpha-amylase (AMY) enzyme in Bacillus subtilis at a high rate leads to the accumulation of unfolded AMY, which causes secretion stress. The over-expression of the PrsA chaperone aids enzyme folding and reduces stress. To identify affected pathways and potential mechanisms involved in the reduced growth, we analyzed the transcriptomic differences during fed-batch fermentation between a PrsA over-expressing strain and control in a time-series RNA-seq experiment. We observe transcription in 542 unannotated regions, of which 234 had significant changes in expression levels between the samples. Moreover, 1,791 protein-coding sequences, 80 non-coding genes, and 20 riboswitches overlapping UTR regions of coding genes had significant changes in expression. We identified putatively regulated biological processes via gene-set over-representation analysis of the differentially expressed genes; overall, the analysis suggests that the PrsA over-expression affects ATP biosynthesis activity, amino acid metabolism, and cell wall stability. The investigation of the protein interaction network points to a potential impact on cell motility signaling. We discuss the impact of these highlighted mechanisms for reducing secretion stress or detrimental aspects of PrsA over-expression during AMY production.

The Bacillaceae-1 RNA motif comprises two distinct classes.

Non-coding RNAs are key regulatory players in bacteria. Many computationally predicted non-coding RNAs, however, lack functional associations. An example is the Bacillaceae-1 RNA motif, whose Rfam model consists of two hairpin loops. We find the motif conserved in nine of 13 non-pathogenic strains of the genus Bacillus but only in one pathogenic strain. To elucidate functional characteristics, we studied 118 hits of the Rfam model in 11 Bacillus spp. and found two distinct classes based on the ensemble diversity of their RNA secondary structure and the genomic context concerning the ribosomal RNA (rRNA) cluster. Forty hits are associated with the rRNA cluster, of which all 19 hits upstream flanking of 16S rRNA have a reverse complementary structure of low structural diversity. Fifty-two hits have large ensemble diversity, of which 38 are located between two coding genes. For eight hits in Bacillus subtilis, we investigated public expression data under various conditions and observed either the forward or the reverse complementary motif expressed. Five hits are associated with the rRNA cluster. Four of them are located upstream of the 16S rRNA and are not transcriptionally active, but instead, their reverse complements with low structural diversity are expressed together with the rRNA cluster. The three other hits are located between two coding genes in non-conserved genomic loci. Two of them are independently expressed from their surrounding genes and are structurally diverse. In summary, we found that Bacillaceae-1 RNA motifs upstream flanking of ribosomal RNA clusters tend to have one stable structure with the reverse complementary motif expressed in B. subtilis. In contrast, a subgroup of intergenic motifs has the thermodynamic potential for structural switches.

BSGatlas: a unified Bacillus subtilis genome and transcriptome annotation atlas with enhanced information access.

A large part of our current understanding of gene regulation in Gram-positive bacteria is based on Bacillus subtilis, as it is one of the most well studied bacterial model systems. The rapid growth in data concerning its molecular and genomic biology is distributed across multiple annotation resources. Consequently, the interpretation of data from further B. subtilis experiments becomes increasingly challenging in both low- and large-scale analyses. Additionally, B. subtilis annotation of structured RNA and non-coding RNA (ncRNA), as well as the operon structure, is still lagging behind the annotation of the coding sequences. To address these challenges, we created the B. subtilis genome atlas, BSGatlas, which integrates and unifies multiple existing annotation resources. Compared to any of the individual resources, the BSGatlas contains twice as many ncRNAs, while improving the positional annotation for 70 % of the ncRNAs. Furthermore, we combined known transcription start and termination sites with lists of known co-transcribed gene sets to create a comprehensive transcript map. The combination with transcription start/termination site annotations resulted in 717 new sets of co-transcribed genes and 5335 untranslated regions (UTRs). In comparison to existing resources, the number of 5' and 3' UTRs increased nearly fivefold, and the number of internal UTRs doubled. The transcript map is organized in 2266 operons, which provides transcriptional annotation for 92 % of all genes in the genome compared to the at most 82 % by previous resources. We predicted an off-target-aware genome-wide library of CRISPR-Cas9 guide RNAs, which we also linked to polycistronic operons. We provide the BSGatlas in multiple forms: as a website (https://rth.dk/resources/bsgatlas/), an annotation hub for display in the UCSC genome browser, supplementary tables and standardized GFF3 format, which can be used in large scale -omics studies. By complementing existing resources, the BSGatlas supports analyses of the B. subtilis genome and its molecular biology with respect to not only non-coding genes but also genome-wide transcriptional relationships of all genes.

Ribosomal dimerization factor YfiA is the major protein synthesized after abrupt glucose depletion in Lactococcus lactis.

We analysed the response of the model bacterium Lactococcus lactis to abrupt depletion of glucose after several generations of exponential growth. Glucose depletion resulted in a drastic drop in the energy charge accompanied by an extremely low GTP level and an almost total arrest of protein synthesis. Strikingly, the cell prioritized the continued synthesis of a few proteins, of which the ribosomal dimerization factor YfiA was the most highly expressed. Transcriptome analysis showed no immediate decrease in total mRNA levels despite the lowered nucleotide pools and only marginally increased levels of the yfiA transcript. Severe up-regulation of genes in the FruR, CcpA, ArgR and AhrC regulons were consistent with a downshift in carbon and energy source. Based upon the results, we suggest that transcription proceeded long enough to record the transcriptome changes from activation of the FruR, CcpA, ArgR and AhrC regulons, while protein synthesis stopped due to an extremely low GTP concentration emerging a few minutes after glucose depletion. The yfiA deletion mutant exhibited a longer lag phase upon replenishment of glucose and a faster death rate after prolonged starvation supporting that YfiA-mediated ribosomal dimerization is important for keeping long-term starved cells viable and competent for growth initiation.

Global transcriptional analysis of Bacillus licheniformis reveals an overlap between heat shock and iron limitation stimulon.

In this study, we characterized the heat shock stimulon of the important industrial microorganism Bacillus licheniformis using DNA microarrays. While sharing a high degree of homology with the closely related model organism Bacillus subtilis, the heat shock stimulon of B. licheniformis exhibited several novel and unexpected features. Most notably, heat shock in B. licheniformis resulted in decreased amounts of mRNA from the ytrABCEF operon, encoding a putative acetoin uptake system, and stimulated the transcription of purine biosynthesis and iron uptake genes. Unexpectedly, deletion of the ytrEF genes did not affect acetoin uptake, but increased heat sensitivity. To investigate the connection between heat stress and iron uptake further, we analyzed the iron limitation response of B. licheniformis by DNA microarrays and concluded that the response mostly involves the genes related to iron uptake and metabolism, while the only heat shock gene affected by iron limitation was clpE. We also attempted to delete the fur gene (encoding the ferric uptake repressor), but unexpectedly found it to be essential in B. licheniformis. Using the fluorescent protein-encoding reporter gene under control of the dhb promoter, which responded to both heat shock and iron-starvation, we confirmed the overlap between these responses.

Resolvase-like recombination performed by the TP901-1 integrase.

The site-specific recombination system of temperate lactococcal bacteriophage TP901-1 is unusual in several respects. First, the integrase belongs to the family of extended resolvases rather than to the lambda integrase family and second, in the presence of this integrase, a 56 bp attP fragment is sufficient for efficient recombination with the chromosomal attB site in the host Lactococcus lactis subsp. cremoris MG1363. In the present work, this attB site was analysed and a 43 bp attB region was found to be the smallest fragment able to participate fully in recombination. In vitro studies showed that the TP901-1 integrase binds this 43 bp attB fragment, the 56 bp attP and a larger attP fragment with equal affinity. Mutational analysis of the 5 bp common core region (TCAAT) showed that the TC dinucleotide is essential for recombination, but not for binding of the integrase, whereas none of the last three bases are important for recombination. When a number of attL sites, obtained by recombination between an attB site containing a mutation in this TC dinucleotide and a wild-type attP site, were sequenced, a mix of sites with the wild-type or the mutated sequence was obtained. These results are consistent with the hypothesis that the TC dinucleotide constitutes the TP901-1 overlap region. A 2 bp overlap region has been observed in recombination reactions catalysed by all other members of the resolvase/invertase family tested so far. By selecting for attB sites with a decreased ability to participate in recombination, two bases located outside the core region of attB were shown to be involved in the in vitro binding of the TP901-1 integrase.