The role of mRNA structure in bacterial translational regulation

Machine learning-based classification reveals distinct clusters of non-coding genomic allelic variations associated with Erm-mediated antibiotic resistance

The erythromycin resistance RNA methyltransferase (erm) confers cross-resistance to all therapeutically important macrolides, lincosamides, and streptogramins (MLS phenotype). The expression of erm is often induced by the macrolide-mediated ribosome stalling in the upstream co-transcribed leader sequence, thereby triggering a conformational switch of the intergenic RNA hairpins to allow the translational initiation of erm. We investigated the evolutionary emergence of the upstream erm regulatory elements and the impact of allelic variation on erm expression and the MLS phenotype. Through systematic profiling of the upstream regulatory sequences across all known erm operons, we observed that specific erm subfamilies, such as ermB and ermC, have independently evolved distinct configurations of small upstream ORFs and palindromic repeats. A population-wide genomic analysis of the upstream ermB regions revealed substantial non-random allelic variation at numerous positions. Utilizing machine learning-based classification coupled with RNA structure modeling, we found that many alleles cooperatively influence the stability of alternative RNA hairpin structures formed by the palindromic repeats, which, in turn, affects the inducibility of ermB expression and MLS phenotypes. Subsequent experimental validation of 11 randomly selected variants demonstrated an impressive 91% accuracy in predicting MLS phenotypes. Furthermore, we uncovered a mixed distribution of MLS-sensitive and MLS-resistant ermB loci within the evolutionary tree, indicating repeated and independent evolution of MLS resistance. Taken together, this study not only elucidates the evolutionary processes driving the emergence and development of MLS resistance but also highlights the potential of using non-coding genomic allele data to predict antibiotic resistance phenotypes.

IMPORTANCE

Antibiotic resistance (AR) poses a global health threat as the efficacy of available antibiotics has rapidly eroded due to the widespread transmission of AR genes. Using Erm-dependent MLS resistance as a model, this study highlights the significance of non-coding genomic allelic variations. Through a comprehensive analysis of upstream regulatory elements within the erm family, we elucidated the evolutionary emergence and development of AR mechanisms. Leveraging population-wide machine learning (ML)-based genomic analysis, we transformed substantial non-random allelic variations into discernible clusters of elements, enabling precise prediction of MLS phenotypes from non-coding regions. These findings offer deeper insight into AR evolution and demonstrate the potential of harnessing non-coding genomic allele data for accurately predicting AR phenotypes.

Suppression of the Escherichia coli rnpA49 conditionally lethal phenotype by different compensatory mutations

RNase P is an essential enzyme found across all domains of life that is responsible for the 5'-end maturation of precursor tRNAs. For decades, numerous studies have sought to elucidate the mechanisms and biochemistry governing RNase P function. However, much remains unknown about the regulation of RNase P expression, the turnover and degradation of the enzyme, and the mechanisms underlying the phenotypes and complementation of specific RNase P mutations, especially in the model bacterium, Escherichia coli In E. coli, the temperature-sensitive (ts) rnpA49 mutation in the protein subunit of RNase P has arguably been one of the most well-studied mutations for examining the enzyme's activity in vivo. Here, we report for the first time naturally occurring temperature-resistant suppressor mutations of E. coli strains carrying the rnpA49 allele. We find that rnpA49 strains can partially compensate the ts defect via gene amplifications of either RNase P subunit (rnpA49 or rnpB) or by the acquisition of loss-of-function mutations in Lon protease or RNase R. Our results agree with previous plasmid overexpression and gene deletion complementation studies, and importantly suggest the involvement of Lon protease in the degradation and/or regulatory pathway(s) of the mutant protein subunit of RNase P. This work offers novel insights into the behavior and complementation of the rnpA49 allele in vivo and provides direction for follow-up studies regarding RNase P regulation and turnover in E. coli.

Illuminating the oral microbiome and its host interactions: tools and approaches for molecular microbiology studies

Advancements in DNA sequencing technologies within the last decade have stimulated an unprecedented interest in the human microbiome, largely due the broad diversity of human diseases found to correlate with microbiome dysbiosis. As a direct consequence of these studies, a vast number of understudied and uncharacterized microbes have been identified as potential drivers of mucosal health and disease. The looming challenge in the field is to transition these observations into defined molecular mechanistic studies of symbiosis and dysbiosis. In order to meet this challenge, many of these newly identified microbes will need to be adapted for use in experimental models. Consequently, this review presents a comprehensive overview of the molecular microbiology tools and techniques that have played crucial roles in genetic studies of the bacteria found within the human oral microbiota. Here, we will use specific examples from the oral microbiome literature to illustrate the biology supporting these techniques, why they are needed in the field, and how such technologies have been implemented. It is hoped that this information can serve as a useful reference guide to help catalyze molecular microbiology studies of the many new understudied and uncharacterized species identified at different mucosal sites in the body.

Analysis of key genes in Mycobacterium ulcerans reveals conserved RNA structural motifs and regions with apparent pressure to remain unstructured

Buruli Ulcer is a neglected tropical disease that results in disfiguring and dangerous lesions in affected persons across a wide geographic area, including much of West Africa. The causative agent of Buruli Ulcer is Mycobacterium ulcerans, a relative of the bacterium that causes tuberculosis and leprosy. Few therapeutic options exist for the treatment of this disease beyond antibiotics in the early stages, which are frequently ineffective, and surgical removal in the later stage. In this study we analyze six genes in Mycobacterium ulcerans that have high potential of therapeutic targeting. We focus our analysis on a combined in silico and comparative sequence study of potential RNA secondary structure across these genes. The result of this work was the comprehensive local RNA structural landscape across each of these significant genes. This revealed multiple sites of ordered and evolved RNA structure interspersed between sequences that either have no bias for structure or, indeed, appear to be ordered to be unstructured and (potentially) accessible. In addition to providing data that could be of interest to basic biology, our results provide guides for efforts aimed at targeting this pathogen at the RNA level. We explore this latter possibility through the in silico analysis of antisense oligonucleotides that could potentially be used to target pathogen RNA.

The transcriptomic landscape of Magnetospirillum gryphiswaldense during magnetosome biomineralization

BACKGROUND

One of the most complex prokaryotic organelles are magnetosomes, which are formed by magnetotactic bacteria as sensors for navigation in the Earth's magnetic field. In the alphaproteobacterium Magnetospirillum gryphiswaldense magnetosomes consist of chains of magnetite crystals (Fe₃O₄) that under microoxic to anoxic conditions are biomineralized within membrane vesicles. To form such an intricate structure, the transcription of > 30 specific structural genes clustered within the genomic magnetosome island (MAI) has to be coordinated with the expression of an as-yet unknown number of auxiliary genes encoding several generic metabolic functions. However, their global regulation and transcriptional organization in response to anoxic conditions most favorable for magnetite biomineralization are still unclear.

RESULTS

Here, we compared transcriptional profiles of anaerobically grown magnetosome forming cells with those in which magnetosome biosynthesis has been suppressed by aerobic condition. Using whole transcriptome shotgun sequencing, we found that transcription of about 300 of the > 4300 genes was significantly enhanced during magnetosome formation. About 40 of the top upregulated genes are directly or indirectly linked to aerobic and anaerobic respiration (denitrification) or unknown functions. The mam and mms gene clusters, specifically controlling magnetosome biosynthesis, were highly transcribed, but constitutively expressed irrespective of the growth condition. By Cappable-sequencing, we show that the transcriptional complexity of both the MAI and the entire genome decreased under anaerobic conditions optimal for magnetosome formation. In addition, predominant promoter structures were highly similar to sigma factor σ⁷⁰ dependent promoters in other Alphaproteobacteria.

CONCLUSIONS

Our transcriptome-wide analysis revealed that magnetite biomineralization relies on a complex interplay between generic metabolic processes such as aerobic and anaerobic respiration, cellular redox control, and the biosynthesis of specific magnetosome structures. In addition, we provide insights into global regulatory features that have remained uncharacterized in the widely studied model organism M. gryphiswaldense, including a comprehensive dataset of newly annotated transcription start sites and genome-wide operon detection as a community resource (GEO Series accession number GSE197098).

Simultaneous ribosome profiling of hundreds of microbes from the human microbiome

Ribosome profiling enables sequencing of ribosome-bound fragments of RNA, revealing which transcripts are being translated as well as the position of ribosomes along mRNAs. Although ribosome profiling has been applied to cultured bacterial isolates, its application to uncultured, mixed communities has been challenging. We present MetaRibo-Seq, a protocol that enables the application of ribosome profiling directly to the human fecal microbiome. MetaRibo-Seq is a benchmarked method that includes several modifications to existing ribosome profiling protocols, specifically addressing challenges involving fecal sample storage, purity and input requirements. We also provide a computational workflow to quality control and trim reads, de novo assemble a reference metagenome with metagenomic reads, align MetaRibo-Seq reads to the reference, and assess MetaRibo-Seq library quality ( https://github.com/bhattlab/bhattlab_workflows/tree/master/metariboseq ). This MetaRibo-Seq protocol enables researchers in standard molecular biology laboratories to study translation in the fecal microbiome in ~5 d.

RelA binding of mRNAs modulates translation or sRNA-mRNA basepairing depending on the position of the GGAG site

Previously, we reported that RelA protein facilitates Hfq-mediated mRNA-sRNA regulation by binding sRNAs carrying a Shine-Dalgarno-like GGAG sequence. In turn, sRNA-Hfq monomers are stabilized, enabling the attachment of more Hfq subunits to form a functional hexamer. Here, using CLIP-seq, we present a global analysis of RelA-bound RNAs showing that RelA interacts with sRNAs as well as with mRNAs carrying a GGAG motif. RelA binding of mRNAs carrying GGAG at position -7 relative to the initiation codon (AUG) inhibits translation by interfering with the binding of 30S ribosomes. The extent of inhibition depends on the distance of GGAG relative to the AUG, as shortening the spacing between GGAG and AUG abrogates RelA-mediated inhibition. Interestingly, RelA binding of target mRNAs carrying GGAG in the coding sequence or close to AUG facilitates target gene regulation by sRNA partners that lack GGAG. However, translation inhibition caused by RelA binding of mRNAs carrying GGAG at position -7 relative to the AUG renders sRNA-mRNA basepairing regulation ineffective. Our study indicates that by binding RNAs carrying GGAG the ribosome-associated RelA protein inhibits translation of specific newly synthesized incoming mRNAs or enables basepairing regulation by their respective sRNA partners, thereby introducing a new regulatory concept for the bacterial response.

An mRNA-mRNA Interaction Couples Expression of a Virulence Factor and Its Chaperone in Listeria monocytogenes

Bacterial pathogens often employ RNA regulatory elements located in the 5' untranslated regions (UTRs) to control gene expression. Using a comparative structural analysis, we examine the structure of 5' UTRs at a global scale in the pathogenic bacterium Listeria monocytogenes under different conditions. In addition to discovering an RNA thermoswitch and detecting simultaneous interaction of ribosomes and small RNAs with mRNA, we identify structural changes in the 5' UTR of an mRNA encoding the post-translocation chaperone PrsA2 during infection conditions. We demonstrate that the 5' UTR of the prsA2 mRNA base pairs with the 3' UTR of the full-length hly mRNA encoding listeriolysin O, thus preventing RNase J1-mediated degradation of the prsA2 transcript. Mutants lacking the hly-prsA2 interaction exhibit reduced virulence properties. This work highlights an additional level of RNA regulation, where the mRNA encoding a chaperone is stabilized by the mRNA encoding its substrate.

The World of Stable Ribonucleoproteins and Its Mapping With Grad-Seq and Related Approaches

Macromolecular complexes of proteins and RNAs are essential building blocks of cells. These stable supramolecular particles can be viewed as minimal biochemical units whose structural organization, i.e., the way the RNA and the protein interact with each other, is directly linked to their biological function. Whether those are dynamic regulatory ribonucleoproteins (RNPs) or integrated molecular machines involved in gene expression, the comprehensive knowledge of these units is critical to our understanding of key molecular mechanisms and cell physiology phenomena. Such is the goal of diverse complexomic approaches and in particular of the recently developed gradient profiling by sequencing (Grad-seq). By separating cellular protein and RNA complexes on a density gradient and quantifying their distributions genome-wide by mass spectrometry and deep sequencing, Grad-seq charts global landscapes of native macromolecular assemblies. In this review, we propose a function-based ontology of stable RNPs and discuss how Grad-seq and related approaches transformed our perspective of bacterial and eukaryotic ribonucleoproteins by guiding the discovery of new RNA-binding proteins and unusual classes of noncoding RNAs. We highlight some methodological aspects and developments that permit to further boost the power of this technique and to look for exciting new biology in understudied and challenging biological models.

A deep learning approach to programmable RNA switches

Engineered RNA elements are programmable tools capable of detecting small molecules, proteins, and nucleic acids. Predicting the behavior of these synthetic biology components remains a challenge, a situation that could be addressed through enhanced pattern recognition from deep learning. Here, we investigate Deep Neural Networks (DNN) to predict toehold switch function as a canonical riboswitch model in synthetic biology. To facilitate DNN training, we synthesize and characterize in vivo a dataset of 91,534 toehold switches spanning 23 viral genomes and 906 human transcription factors. DNNs trained on nucleotide sequences outperform (R2 = 0.43-0.70) previous state-of-the-art thermodynamic and kinetic models (R2 = 0.04-0.15) and allow for human-understandable attention-visualizations (VIS4Map) to identify success and failure modes. This work shows that deep learning approaches can be used for functionality predictions and insight generation in RNA synthetic biology.

Autogenous regulation in vivo of the rpmE gene encoding ribosomal protein L31 (bL31), a key component of the protein-protein intersubunit bridge B1b

Bacterial ribosomal proteins (r-proteins) encoded by nonessential genes often carry out very important tasks in translation. In particular, this is the case of a small basic bacteria-specific r-protein L31 (bL31). Recent studies revealed a crucial role of bL31 in formation of the protein-protein intersubunit bridge B1b and hence its contribution to ribosome dynamics. Our goal was to study in vivo regulation of the rpmE operon encoding bL31. We used a previously developed approach based on chromosomally integrated fusions with the lacZ reporter. E. coli rpmE is transcribed from two promoter regions, and translation of both mRNA transcripts was shown to be feedback regulated by bL31, indicating that the autogenous operator is located within the shorter transcript. The bL31-mediated control of rpmE is gene-specific, as no regulation was found for rpmE-unrelated reporters. Thus, bL31, as many other r-proteins, possesses dual activity in living cells, acting both as an integral ribosome component and an autogenous repressor. Phylogenetic studies revealed the presence of a highly conserved stem-loop structure in the rpmE 5'UTR, a presumable translational operator targeted by bL31, which was further confirmed by site-directed mutagenesis. This stable operator stem-loop separates an AU-rich translational enhancer from a Shine-Dalgarno element, which is a rare case of a noncontiguous translation initiation region. Sequence/structure computational approaches classify bL31 as an RNA-binding protein, consistent with its repressor function discovered here. Mutational analysis of bL31 showed that its unstructured amino-terminal part enriched in lysine is necessary for the repressor activity.

MetaRibo-Seq measures translation in microbiomes

No method exists to measure large-scale translation of genes in uncultured organisms in microbiomes. To overcome this limitation, we develop MetaRibo-Seq, a method for simultaneous ribosome profiling of tens to hundreds of organisms in microbiome samples. MetaRibo-Seq was benchmarked against gold-standard Ribo-Seq in a mock microbial community and applied to five different human fecal samples. Unlike RNA-Seq, Ribo-Seq signal of a predicted gene suggests it encodes a translated protein. We demonstrate two applications of this technique: First, MetaRibo-Seq identifies small genes, whose identification until now has been challenging. For example, MetaRibo-Seq identifies 2,091 translated, previously unannotated small protein families from five fecal samples, more than doubling the number of small proteins predicted to exist in this niche. Second, the combined application of RNA-Seq and MetaRibo-Seq identifies differences in the translation of transcripts. In summary, MetaRibo-Seq enables comprehensive translational profiling in microbiomes and identifies previously unannotated small proteins.

Defining the functions of individual organisms or communities within microbiomes is a challenging task. Here, the authors develop MetaRibo-Seq, a method for simultaneous high-throughput ribosome profiling of organisms in uncultured microbiome samples.

Translational regulation of environmental adaptation in bacteria

Bacteria must rapidly respond to both intracellular and environmental changes to survive. One critical mechanism to rapidly detect and adapt to changes in environmental conditions is control of gene expression at the level of protein synthesis. At each of the three major steps of translation-initiation, elongation, and termination-cells use stimuli to tune translation rate and cellular protein concentrations. For example, changes in nutrient concentrations in the cell can lead to translational responses involving mechanisms such as dynamic folding of riboswitches during translation initiation or the synthesis of alarmones, which drastically alter cell physiology. Moreover, the cell can fine-tune the levels of specific protein products using programmed ribosome pausing or inducing frameshifting. Recent studies have improved understanding and revealed greater complexity regarding long-standing paradigms describing key regulatory steps of translation such as start-site selection and the coupling of transcription and translation. In this review, we describe how bacteria regulate their gene expression at the three translational steps and discuss how translation is used to detect and respond to changes in the cellular environment. Finally, we appraise the costs and benefits of regulation at the translational level in bacteria.

Strategies for enhancing gene expression in Escherichia coli

Regulation of gene expression is fundamental for cellular function. Upon manipulation of the mechanism of gene expression in Escherichia coli, various bioproducts have been developed that are valuable industrially and medically in the last four decades. To efficiently produce bioproducts, numerous molecular tools are used for enhancing expression at the transcriptional and translational levels. Our recent discovery identified a new approach that enhances the gene expression in E. coli using the gene sequence of the eukaryote, Dictyostelium discoideum. In this review, we highlight the current molecular strategies used for high-level gene expression techniques commonly utilized in basic and applied microbiology.

A genetic selection reveals functional metastable structures embedded in a toxin-encoding mRNA

Post-transcriptional regulation plays important roles to fine-tune gene expression in bacteria. In particular, regulation of type I toxin-antitoxin (TA) systems is achieved through sophisticated mechanisms involving toxin mRNA folding. Here, we set up a genetic approach to decipher the molecular underpinnings behind the regulation of a type I TA in Helicobacter pylori. We used the lethality induced by chromosomal inactivation of the antitoxin to select mutations that suppress toxicity. We found that single point mutations are sufficient to allow cell survival. Mutations located either in the 5’ untranslated region or within the open reading frame of the toxin hamper its translation by stabilizing stem-loop structures that sequester the Shine-Dalgarno sequence. We propose that these short hairpins correspond to metastable structures that are transiently formed during transcription to avoid premature toxin expression. This work uncovers the co-transcriptional inhibition of translation as an additional layer of TA regulation in bacteria.

Unstructured 5'-tails act through ribosome standby to override inhibitory structure at ribosome binding sites

Initiation is the rate-limiting step in translation. It is well-known that stable structure at a ribosome binding site (RBS) impedes initiation. The ribosome standby model of de Smit and van Duin, based on studies of the MS2 phage coat cistron, proposed how high translation rates can be reconciled with stable, inhibitory structures at an RBS. Here, we revisited the coat protein system and assessed the translation efficiency from its sequestered RBS by introducing standby mutations. Further experiments with gfp reporter constructs assessed the effects of 5′-tails—as standby sites—with respect to length and sequence contributions. In particular, combining in vivo and in vitro assays, we can show that tails of CA-dinucleotide repeats—and to a lesser extent, AU-repeats—dramatically increase translation rates. Tails of increasing length reach maximal rate-enhancing effects at 16–18 nucleotides. These standby tails are single-stranded and do not exert their effect by structure changes in the neighboring RBS stem–loop. In vitro translation and toeprinting assays furthermore demonstrate that standby effects are exerted at the level of translation initiation. Finally, as expected, destabilizing mutations within the coat RBS indicate an interplay with the effects of standby tails.

The ribosomal protein S1-dependent standby site in tisB mRNA consists of a single-stranded region and a 5' structure element

Ribosome standby is a mechanism that allows translation initiation at ribosome-binding sites that display stable, inhibitory structures. It involves initiator-tRNA-independent 30S subunit binding to single-stranded RNA regions, and the subsequent relocation to the sequestered ribosome-binding sites (RBS). Direct evidence for 30S preloading had previously been elusive. We report here on a detailed characterization of the standby site in tisB mRNA. 30S subunits bind to a single-stranded region and a 5′-stem-loop structure, as shown by fluorescence anisotropy experiments and footprint mapping by cross-linking–immunoprecipitation experiments. Ribosomal protein S1, on its own and in the context of the 30S ribosome, binds to the standby site. This is required for standby-dependent translation, likely reflecting S1-dependent directional unfolding over more than ≈100 nt to reach the sequestered RBS.

In bacteria, stable RNA structures that sequester ribosome-binding sites (RBS) impair translation initiation, and thus protein output. In some cases, ribosome standby can overcome inhibition by structure: 30S subunits bind sequence-nonspecifically to a single-stranded region and, on breathing of the inhibitory structure, relocate to the RBS for initiation. Standby can occur over long distances, as in the active, +42 tisB mRNA, encoding a toxin. This mRNA is translationally silenced by an antitoxin sRNA, IstR-1, that base pairs to the standby site. In tisB and other cases, a direct interaction between 30S subunits and a standby site has remained elusive. Based on fluorescence anisotropy experiments, ribosome toeprinting results, in vitro translation assays, and cross-linking–immunoprecipitation (CLIP) in vitro, carried out on standby-proficient and standby-deficient tisB mRNAs, we provide a thorough characterization of the tisB standby site. 30S subunits and ribosomal protein S1 alone display high-affinity binding to standby-competent fluorescein-labeled +42 mRNA, but not to mRNAs that lack functional standby sites. Ribosomal protein S1 is essential for standby, as 30∆S1 subunits do not support standby-dependent toeprints and TisB translation in vitro. S1 alone- and 30S-CLIP followed by RNA-seq mapping shows that the functional tisB standby site consists of the expected single-stranded region, but surprisingly, also a 5′-end stem-loop structure. Removal of the latter by 5′-truncations, or disruption of the stem, abolishes 30S binding and standby activity. Based on the CLIP-read mapping, the long-distance standby effect in +42 tisB mRNA (∼100 nt) is tentatively explained by S1-dependent directional unfolding toward the downstream RBS.

Introduction to Genome Biology and Diversity

Organisms display astonishing levels of cell and molecular diversity, including genome size, shape, and architecture. In this chapter, we review how the genome can be viewed as both a structural and an informational unit of biological diversity and explicitly define our intended meaning of genetic information. A brief overview of the characteristic features of bacterial, archaeal, and eukaryotic cell types and viruses sets the stage for a review of the differences in organization, size, and packaging strategies of their genomes. We include a detailed review of genetic elements found outside the primary chromosomal structures, as these provide insights into how genomes are sometimes viewed as incomplete informational entities. Lastly, we reassess the definition of the genome in light of recent advancements in our understanding of the diversity of genomic structures and the mechanisms by which genetic information is expressed within the cell. Collectively, these topics comprise a good introduction to genome biology for the newcomer to the field and provide a valuable reference for those developing new statistical or computation methods in genomics. This review also prepares the reader for anticipated transformations in thinking as the field of genome biology progresses.

Expansion of the SOS regulon of Vibrio cholerae through extensive transcriptome analysis and experimental validation

BACKGROUND

The SOS response is an almost ubiquitous response of cells to genotoxic stresses. The full complement of genes in the SOS regulon for Vibrio species has only been addressed through bioinformatic analyses predicting LexA binding box consensus and in vitro validation. Here, we perform whole transcriptome sequencing from Vibrio cholerae treated with mitomycin C as an SOS inducer to characterize the SOS regulon and other pathways affected by this treatment.

RESULTS

Comprehensive transcriptional profiling allowed us to define the full landscape of promoters and transcripts active in V. cholerae. We performed extensive transcription start site (TSS) mapping as well as detection/quantification of the coding and non-coding RNA (ncRNA) repertoire in strain N16961. To improve TSS detection, we developed a new technique to treat RNA extracted from cells grown in various conditions. This allowed for identification of 3078 TSSs with an average 5'UTR of 116 nucleotides, and peak distribution between 16 and 64 nucleotides; as well as 629 ncRNAs. Mitomycin C treatment induced transcription of 737 genes and 28 ncRNAs at least 2 fold, while it repressed 231 genes and 17 ncRNAs. Data analysis revealed that in addition to the core genes known to integrate the SOS regulon, several metabolic pathways were induced. This study allowed for expansion of the Vibrio SOS regulon, as twelve genes (ubiEJB, tatABC, smpA, cep, VC0091, VC1190, VC1369-1370) were found to be co-induced with their adjacent canonical SOS regulon gene(s), through transcriptional read-through. Characterization of UV and mitomycin C susceptibility for mutants of these newly identified SOS regulon genes and other highly induced genes and ncRNAs confirmed their role in DNA damage rescue and protection.

CONCLUSIONS

We show that genotoxic stress induces a pervasive transcriptional response, affecting almost 20% of the V. cholerae genes. We also demonstrate that the SOS regulon is larger than previously known, and its syntenic organization is conserved among Vibrio species. Furthermore, this specific co-localization is found in other γ-proteobacteria for genes recN-smpA and rmuC-tatABC, suggesting SOS regulon conservation in this phylum. Finally, we comment on the limitations of widespread NGS approaches for identification of all RNA species in bacteria.

Transcription and translation of the sigG gene is tuned for proper execution of the switch from early to late gene expression in the developing Bacillus subtilis spore

A cascade of alternative sigma factors directs developmental gene expression during spore formation by the bacterium Bacillus subtilis. As the spore develops, a tightly regulated switch occurs in which the early-acting sigma factor σ^F is replaced by the late-acting sigma factor σ^G. The gene encoding σ^G (sigG) is transcribed by σ^F and by σ^G itself in an autoregulatory loop; yet σ^G activity is not detected until σ^F-dependent gene expression is complete. This separation in σ^F and σ^G activities has been suggested to be due at least in part to a poorly understood intercellular checkpoint pathway that delays sigG expression by σ^F. Here we report the results of a careful examination of sigG expression during sporulation. Unexpectedly, our findings argue against the existence of a regulatory mechanism to delay sigG transcription by σ^F and instead support a model in which sigG is transcribed by σ^F with normal timing, but at levels that are very low. This low-level expression of sigG is the consequence of several intrinsic features of the sigG regulatory and coding sequence—promoter spacing, secondary structure potential of the mRNA, and start codon identity—that dampen its transcription and translation. Especially notable is the presence of a conserved hairpin in the 5’ leader sequence of the sigG mRNA that occludes the ribosome-binding site, reducing translation by up to 4-fold. Finally, we demonstrate that misexpression of sigG from regulatory and coding sequences lacking these features triggers premature σ^G activity in the forespore during sporulation, as well as inappropriate σ^G activity during vegetative growth. Altogether, these data indicate that transcription and translation of the sigG gene is tuned to prevent vegetative expression of σ^G and to ensure the precise timing of the switch from σ^F to σ^G in the developing spore.

Global changes in gene expression occur during normal cellular growth and development, as well as during cancer cell transformation and bacterial pathogenesis. In this study we have investigated the molecular mechanisms that drive the switch from early to late developmental gene expression during spore formation by the model bacterium Bacillus subtilis. At early times, gene expression in the developing spore is directed by the transcription factor σ^F; at later times σ^F is replaced by σ^G. An important, yet poorly understood aspect of this σ^F-to-σ^G transition is how σ^G activation is delayed until the early, σ^F-directed phase of gene expression is complete. Here we have carefully examined expression of the gene encoding σ^G, sigG, and found that its transcription and translation are ordinarily dampened by several features of its regulatory and coding sequences. Moreover, we have found that this “tuning” of sigG expression is required for proper timing of the switch to σ^G. These results reframe our understanding of how sigG is regulated during B. subtilis sporulation and, more broadly, advance our understanding of how global changes in gene expression can be precisely executed at the molecular/genetic level.

Quick change: post-transcriptional regulation in Pseudomonas

Pseudomonas species have evolved dynamic and intricate regulatory networks to fine-tune gene expression, with complex regulation occurring at every stage in the processing of genetic information. This approach enables Pseudomonas to generate precise individual responses to the environment in order to improve their fitness and resource economy. The weak correlations we observe between RNA and protein abundance highlight the significant regulatory contribution of a series of intersecting post-transcriptional pathways, influencing mRNA stability, translational activity and ribosome function, to Pseudomonas environmental responses. This review examines our current understanding of three major post-transcriptional regulatory systems in Pseudomonas spp.; Gac/Rsm, Hfq and RimK, and presents an overview of new research frontiers, emerging genome-wide methodologies, and their potential for the study of global regulatory responses in Pseudomonas.

Transcriptome sequencing of the human pathogen Corynebacterium diphtheriae NCTC 13129 provides detailed insights into its transcriptional landscape and into DtxR-mediated transcriptional regulation

Background

The human pathogen Corynebacterium diphtheriae is the causative agent of diphtheria. In the 1990s a large diphtheria outbreak in Eastern Europe was caused by the strain C. diphtheriae NCTC 13129. Although the genome was sequenced more than a decade ago, not much is known about its transcriptome. Our aim was to use transcriptome sequencing (RNA-Seq) to close this knowledge gap and gain insights into the transcriptional landscape of a C. diphtheriae tox⁺ strain.

Results

We applied two different RNA-Seq techniques, one to retrieve 5′-ends of primary transcripts and the other to characterize the whole transcriptional landscape in order to gain insights into various features of the C. diphtheriae NCTC 13129 transcriptome. By examining the data we identified 1656 transcription start sites (TSS), of which 1202 were assigned to genes and 454 to putative novel transcripts. By using the TSS data promoter regions recognized by the housekeeping sigma factor σ^A and its motifs were analyzed in detail, revealing a well conserved −10 but an only weakly conserved −35 motif, respectively. Furthermore, with the TSS data 5’-UTR lengths were explored. The observed 5’-UTRs range from zero length (leaderless transcripts), which make up 20% of all genes, up to over 450 nt long leaders, which may harbor regulatory functions. The C. diphtheriae transcriptome consists of 471 operons which are further divided into 167 sub-operon structures. In a differential expression analysis approach, we discovered that genetic disruption of the iron-sensing transcription regulator DtxR, which controls expression of diphtheria toxin (DT), causes a strong influence on general gene expression. Nearly 15% of the genome is differentially transcribed, indicating that DtxR might have other regulatory functions in addition to regulation of iron metabolism and DT. Furthermore, our findings shed light on the transcriptional landscape of the DT encoding gene tox and present evidence for two tox antisense RNAs, which point to a new way of transcriptional regulation of toxin production.

Conclusions

This study presents extensive insights into the transcriptome of C. diphtheriae and provides a basis for future studies regarding gene characterization, transcriptional regulatory networks, and regulation of the tox gene in particular.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4481-8) contains supplementary material, which is available to authorized users.

When eukaryotes and prokaryotes look alike: the case of regulatory RNAs

The discovery that all living entities express many RNAs beyond mRNAs, tRNAs and rRNAs has been a surprise in the past two decades. In fact, regulatory RNAs (regRNAs) are plentiful, and we report stunning parallels between their mechanisms and functions in prokaryotes and eukaryotes. For instance, prokaryotic CRISPR (clustered regularly interspaced short palindromic repeats) defense systems are functional analogs to eukaryotic RNA interference processes that preserve the cell against foreign nucleic acid elements. Regulatory RNAs shape the genome in many ways: by controlling mobile element transposition in both domains, via regulation of plasmid counts in prokaryotes, or by directing epigenetic modifications of DNA and associated proteins in eukaryotes. RegRNAs control gene expression extensively at transcriptional and post-transcriptional levels, with crucial roles in fine-tuning cell environmental responses, including intercellular interactions. Although the lengths, structures and outcomes of the regRNAs in all life kingdoms are disparate, they act through similar patterns: by guiding effectors to target molecules or by sequestering macromolecules to hamper their functions. In addition, their biogenesis processes have a lot in common. This unifying vision of regRNAs in all living cells from bacteria to humans points to the possibility of fruitful exchanges between fundamental and applied research in both domains.

Detection of 224 candidate structured RNAs by comparative analysis of specific subsets of intergenic regions

The discovery of structured non-coding RNAs (ncRNAs) in bacteria can reveal new facets of biology and biochemistry. Comparative genomics analyses executed by powerful computer algorithms have successfully been used to uncover many novel bacterial ncRNA classes in recent years. However, this general search strategy favors the discovery of more common ncRNA classes, whereas progressively rarer classes are correspondingly more difficult to identify. In the current study, we confront this problem by devising several methods to select subsets of intergenic regions that can concentrate these rare RNA classes, thereby increasing the probability that comparative sequence analysis approaches will reveal their existence. By implementing these methods, we discovered 224 novel ncRNA classes, which include ROOL RNA, an RNA class averaging 581 nt and present in multiple phyla, several highly conserved and widespread ncRNA classes with properties that suggest sophisticated biochemical functions and a multitude of putative cis-regulatory RNA classes involved in a variety of biological processes. We expect that further research on these newly found RNA classes will reveal additional aspects of novel biology, and allow for greater insights into the biochemistry performed by ncRNAs.

Search for 5'-leader regulatory RNA structures based on gene annotation aided by the RiboGap database

The discovery of noncoding RNAs (ncRNAs) and their importance for gene regulation led us to develop bioinformatics tools to pursue the discovery of novel ncRNAs. Finding ncRNAs de novo is challenging, first due to the difficulty of retrieving large numbers of sequences for given gene activities, and second due to exponential demands on calculation needed for comparative genomics on a large scale. Recently, several tools for the prediction of conserved RNA secondary structure were developed, but many of them are not designed to uncover new ncRNAs, or are too slow for conducting analyses on a large scale. Here we present various approaches using the database RiboGap as a primary tool for finding known ncRNAs and for uncovering simple sequence motifs with regulatory roles. This database also can be used to easily extract intergenic sequences of eubacteria and archaea to find conserved RNA structures upstream of given genes. We also show how to extend analysis further to choose the best candidate ncRNAs for experimental validation.