Operon mRNAs are organized into ORF-centric structures that predict translation efficiency

The diversity of Shine-Dalgarno sequences sheds light on the evolution of translation initiation

Shine-Dalgarno (SD) sequences, the core element of prokaryotic ribosome-binding sites, facilitate mRNA translation by base-pair interaction with the anti-SD (aSD) sequence of 16S rRNA. In contrast to this paradigm, an inspection of thousands of prokaryotic species unravels tremendous SD sequence diversity both within and between genomes, whereas aSD sequences remain largely static. The pattern has led many to suggest unidentified mechanisms for translation initiation. Here we review known translation-initiation pathways in prokaryotes. Moreover, we seek to understand the cause and consequence of SD diversity through surveying recent advances in biochemistry, genomics, and high-throughput genetics. These findings collectively show: (1) SD:aSD base pairing is beneficial but nonessential to translation initiation. (2) The 5' untranslated region of mRNA evolves dynamically and correlates with organismal phylogeny and ecological niches. (3) Ribosomes have evolved distinct usage of translation-initiation pathways in different species. We propose a model portraying the SD diversity shaped by optimization of gene expression, adaptation to environments and growth demands, and the species-specific prerequisite of ribosomes to initiate translation. The model highlights the coevolution of ribosomes and mRNA features, leading to functional customization of the translation apparatus in each organism.

How Many Messenger RNAs Can Be Translated by the START Mechanism?

Translation initiation is a key step in the protein synthesis stage of the gene expression pathway of all living cells. In this important process, ribosomes have to accurately find the AUG start codon in order to ensure the integrity of the proteome. “Structure Assisted RNA Translation”, or “START”, has been proposed to use stable secondary structures located in the coding sequence to augment start site selection by steric hindrance of the progression of pre-initiation complex on messenger RNA. This implies that such structures have to be located downstream and at on optimal distance from the AUG start codon (i.e., downstream nucleotide +16). In order to assess the importance of the START mechanism in the overall mRNA translation process, we developed a bioinformatic tool to screen coding sequences for such stable structures in a 50 nucleotide-long window spanning the nucleotides from +16 to +65. We screened eight bacterial genomes and six eukaryotic genomes. We found stable structures in 0.6–2.5% of eukaryotic coding sequences. Among these, approximately half of them were structures predicted to form G-quadruplex structures. In humans, we selected 747 structures. In bacteria, the coding sequences from Gram-positive bacteria contained 2.6–4.2% stable structures, whereas the structures were less abundant in Gram-negative bacteria (0.2–2.7%). In contrast to eukaryotes, putative G-quadruplex structures are very rare in the coding sequence of bacteria. Altogether, our study reveals that the START mechanism seems to be an ancient strategy to facilitate the start codon recognition that is used in different kingdoms of life.

ΦX174 Attenuation by Whole-Genome Codon Deoptimization

Natural selection acting on synonymous mutations in protein-coding genes influences genome composition and evolution. In viruses, introducing synonymous mutations in genes encoding structural proteins can drastically reduce viral growth, providing a means to generate potent, live-attenuated vaccine candidates. However, an improved understanding of what compositional features are under selection and how combinations of synonymous mutations affect viral growth is needed to predictably attenuate viruses and make them resistant to reversion. We systematically recoded all nonoverlapping genes of the bacteriophage ΦX174 with codons rarely used in its Escherichia coli host. The fitness of recombinant viruses decreases as additional deoptimizing mutations are made to the genome, although not always linearly, and not consistently across genes. Combining deoptimizing mutations may reduce viral fitness more or less than expected from the effect size of the constituent mutations and we point out difficulties in untangling correlated compositional features. We test our model by optimizing the same genes and find that the relationship between codon usage and fitness does not hold for optimization, suggesting that wild-type ΦX174 is at a fitness optimum. This work highlights the need to better understand how selection acts on patterns of synonymous codon usage across the genome and provides a convenient system to investigate the genetic determinants of virulence.

The Ongoing Quest to Crack the Genetic Code for Protein Production

Understanding the genetic design principles that determine protein production remains a major challenge. Although the key principles of gene expression were discovered 50 years ago, additional factors are still being uncovered. Both protein-coding and non-coding sequences harbor elements that collectively influence the efficiency of protein production by modulating transcription, mRNA decay, and translation. The influences of many contributing elements are intertwined, which complicates a full understanding of the individual factors. In natural genes, a functional balance between these factors has been obtained in the course of evolution, whereas for genetic-engineering projects, our incomplete understanding still limits optimal design of synthetic genes. However, notable advances have recently been made, supported by high-throughput analysis of synthetic gene libraries as well as by state-of-the-art biomolecular techniques. We discuss here how these advances further strengthen understanding of the gene expression process and how they can be harnessed to optimize protein production.

Structure-seq2 probing of RNA structure upon amino acid starvation reveals both known and novel RNA switches in Bacillus subtilis

RNA structure influences numerous processes in all organisms. In bacteria, these processes include transcription termination and attenuation, small RNA and protein binding, translation initiation, and mRNA stability, and can be regulated via metabolite availability and other stresses. Here we use Structure-seq2 to probe the in vivo RNA structurome of Bacillus subtilis grown in the presence and absence of amino acids. Our results reveal that amino acid starvation results in lower overall dimethyl sulfate (DMS) reactivity of the transcriptome, indicating enhanced protection owing to protein binding or RNA structure. Starvation-induced changes in DMS reactivity correlated inversely with transcript abundance changes. This correlation was particularly pronounced in genes associated with the stringent response and CodY regulons, which are involved in adaptation to nutritional stress, suggesting that RNA structure contributes to transcript abundance change in regulons involved in amino acid metabolism. Structure-seq2 accurately reported on four known amino acid-responsive riboswitches: T-box, SAM, glycine, and lysine riboswitches. Additionally, we discovered a transcription attenuation mechanism that reduces yfmG expression when amino acids are added to the growth medium. We also found that translation of a leader peptide (YfmH) encoded just upstream of yfmG regulates yfmG expression. Our results are consistent with a model in which a slow rate of yfmH translation caused by limitation of the amino acids encoded in YfmH prevents transcription termination in the yfmG leader region by favoring formation of an overlapping antiterminator structure. This novel RNA switch offers a way to simultaneously monitor the levels of multiple amino acids.

A possible universal role for mRNA secondary structure in bacterial translation revealed using a synthetic operon

In bacteria, translation re-initiation is crucial for synthesizing proteins encoded by genes that are organized into operons. The mechanisms regulating translation re-initiation remain, however, poorly understood. We now describe the ribosome termination structure (RTS), a conserved and stable mRNA secondary structure localized immediately downstream of stop codons, and provide experimental evidence for its role in governing re-initiation efficiency in a synthetic Escherichia coli operon. We further report that RTSs are abundant, being associated with 18%-65% of genes in 128 analyzed bacterial genomes representing all phyla, and are selectively depleted when translation re-initiation is advantageous yet selectively enriched so as to insulate translation when re-initiation is deleterious. Our results support a potentially universal role for the RTS in controlling translation termination-insulation and re-initiation across bacteria.

Ribosome recycling is not critical for translational coupling in Escherichia coli

We used ribosome profiling to characterize the biological role of ribosome recycling factor (RRF) in Escherichia coli. As expected, RRF depletion leads to enrichment of post-termination 70S complexes in 3′-UTRs. We also observe that elongating ribosomes are unable to complete translation because they are blocked by non-recycled ribosomes at stop codons. Previous studies have suggested a role for recycling in translational coupling within operons; if a ribosome remains bound to an mRNA after termination, it may re-initiate downstream. We found, however, that RRF depletion did not significantly affect coupling efficiency in reporter assays or in ribosome density genome-wide. These findings argue that re-initiation is not a major mechanism of translational coupling in E. coli. Finally, RRF depletion has dramatic effects on the activity of ribosome rescue factors tmRNA and ArfA. Our results provide a global view of the effects of the loss of ribosome recycling on protein synthesis in E. coli.

On the Way to Understanding the Interplay between the RNA Structure and Functions in Cells: A Genome-Wide Perspective

RNAs adopt specific structures in order to perform their biological activities. The structure of RNA is an important layer of gene expression regulation, and can impact a plethora of cellular processes, starting with transcription, RNA processing, and translation, and ending with RNA turnover. The development of high-throughput technologies has enabled a deeper insight into the sophisticated interplay between the structure of the cellular transcriptome and the living cells environment. In this review, we present the current view on the RNA structure in vivo resulting from the most recent transcriptome-wide studies in different organisms, including mammalians, yeast, plants, and bacteria. We focus on the relationship between the mRNA structure and translation, mRNA stability and degradation, protein binding, and RNA posttranscriptional modifications.

Intact RNA structurome reveals mRNA structure-mediated regulation of miRNA cleavage in vivo

MicroRNA (miRNA)-mediated cleavage is involved in numerous essential cellular pathways. miRNAs recognize target RNAs via sequence complementarity. In addition to complementarity, in vitro and in silico studies have suggested that RNA structure may influence the accessibility of mRNAs to miRNA-induced silencing complexes (miRISCs), thereby affecting RNA silencing. However, the regulatory mechanism of mRNA structure in miRNA cleavage remains elusive. We investigated the role of in vivo RNA secondary structure in miRNA cleavage by developing the new CAP-STRUCTURE-seq method to capture the intact mRNA structurome in Arabidopsis thaliana. This approach revealed that miRNA target sites were not structurally accessible for miRISC binding prior to cleavage in vivo. Instead, we found that the unfolding of the target site structure plays a key role in miRISC activity in vivo. We found that the single-strandedness of the two nucleotides immediately downstream of the target site, named Target Adjacent nucleotide Motif, can promote miRNA cleavage but not miRNA binding, thus decoupling target site binding from cleavage. Our findings demonstrate that mRNA structure in vivo can modulate miRNA cleavage, providing evidence of mRNA structure-dependent regulation of biological processes.

Lead-seq: transcriptome-wide structure probing in vivo using lead(II) ions

The dynamic conformation of RNA molecules within living cells is key to their function. Recent advances in probing the RNA structurome in vivo, including the use of SHAPE (Selective 2'-Hydroxyl Acylation analyzed by Primer Extension) or kethoxal reagents or DMS (dimethyl sulfate), provided unprecedented insights into the architecture of RNA molecules in the living cell. Here, we report the establishment of lead probing in a global RNA structuromics approach. In order to elucidate the transcriptome-wide RNA landscape in the enteric pathogen Yersinia pseudotuberculosis, we combined lead(II) acetate-mediated cleavage of single-stranded RNA regions with high-throughput sequencing. This new approach, termed 'Lead-seq', provides structural information independent of base identity. We show that the method recapitulates secondary structures of tRNAs, RNase P RNA, tmRNA, 16S rRNA and the rpsT 5'-untranslated region, and that it reveals global structural features of mRNAs. The application of Lead-seq to Y. pseudotuberculosis cells grown at two different temperatures unveiled the first temperature-responsive in vivo RNA structurome of a bacterial pathogen. The translation of candidate genes derived from this approach was confirmed to be temperature regulated. Overall, this study establishes Lead-seq as complementary approach to interrogate intracellular RNA structures on a global scale.

Recommendations for bacterial ribosome profiling experiments based on bioinformatic evaluation of published data

Ribosome profiling (RIBO-Seq) has improved our understanding of bacterial translation, including finding many unannotated genes. However, protocols for RIBO-Seq and corresponding data analysis are not yet standardized. Here, we analyzed 48 RIBO-Seq samples from nine studies of Escherichia coli K12 grown in lysogeny broth medium and particularly focused on the size-selection step. We show that for conventional expression analysis, a size range between 22 and 30 nucleotides is sufficient to obtain protein-coding fragments, which has the advantage of removing many unwanted rRNA and tRNA reads. More specific analyses may require longer reads and a corresponding improvement in rRNA/tRNA depletion. There is no consensus about the appropriate sequencing depth for RIBO-Seq experiments in prokaryotes, and studies vary significantly in total read number. Our analysis suggests that 20 million reads that are not mapping to rRNA/tRNA are required for global detection of translated annotated genes. We also highlight the influence of drug-induced ribosome stalling, which causes bias at translation start sites. The resulting accumulation of reads at the start site may be especially useful for detecting weakly expressed genes. As different methods suit different questions, it may not be possible to produce a "one-size-fits-all" ribosome profiling data set. Therefore, experiments should be carefully designed in light of the scientific questions of interest. We propose some basic characteristics that should be reported with any new RIBO-Seq data sets. Careful attention to the factors discussed should improve prokaryotic gene detection and the comparability of ribosome profiling data sets.

RSVdb: a comprehensive database of transcriptome RNA structure

RNA fulfills a crucial regulatory role in cells by folding into a complex RNA structure. To date, a chemical compound, dimethyl sulfate (DMS), has been developed to probe the RNA structure at the transcriptome level effectively. We proposed a database, RSVdb (https://taolab.nwafu.edu.cn/rsvdb/), for the browsing and visualization of transcriptome RNA structures. RSVdb, including 626 225 RNAs with validated DMS reactivity from 178 samples in eight species, supports four main functions: information retrieval, research overview, structure prediction and resource download. Users can search for species, studies, transcripts and genes of interest; browse the quality control of sequencing data and statistical charts of RNA structure information; preview and perform online prediction of RNA structures in silico and under DMS restraint of different experimental treatments and download RNA structure data for species and studies. Together, RSVdb provides a reference for RNA structure and will support future research on the function of RNA structure at the transcriptome level.

Impact of C-terminal amino acid composition on protein expression in bacteria

The C-terminal sequence of a protein is involved in processes such as efficiency of translation termination and protein degradation. However, the general relationship between features of this C-terminal sequence and levels of protein expression remains unknown. Here, we identified C-terminal amino acid biases that are ubiquitous across the bacterial taxonomy (1,582 genomes). We showed that the frequency is higher for positively charged amino acids (lysine, arginine), while hydrophobic amino acids and threonine are lower. We then studied the impact of C-terminal composition on protein levels in a library of Mycoplasma pneumoniae mutants, covering all possible combinations of the two last codons. We found that charged and polar residues, in particular lysine, led to higher expression, while hydrophobic and aromatic residues led to lower expression, with a difference in protein levels up to fourfold. We further showed that modulation of protein degradation rate could be one of the main mechanisms driving these differences. Our results demonstrate that the identity of the last amino acids has a strong influence on protein expression levels.

RNA Polymerase's Relationship with the Ribosome: Not So Physical, Most of the Time

In bacteria, the rates of transcription elongation and translation elongation are coordinated, changing together in response to growth conditions. It has been proposed that this is due to physical coupling of RNA polymerase and the lead ribosome on nascent mRNA, an interaction important for preventing premature transcription termination by Rho factor. Recent studies challenge this view and provide evidence that coordination is indirect, mediated in Escherichia coli by the alarmone (p)ppGpp. Here, we discuss these new findings and how they shape our understanding of the functional relationship between RNA polymerase and the ribosome as well as the basis of transcriptional polarity.

Translational initiation in E. coli occurs at the correct sites genome-wide in the absence of mRNA-rRNA base-pairing

Shine-Dalgarno (SD) motifs are thought to play an important role in translational initiation in bacteria. Paradoxically, ribosome profiling studies in E. coli show no correlation between the strength of an mRNA’s SD motif and how efficiently it is translated. Performing profiling on ribosomes with altered anti-Shine-Dalgarno sequences, we reveal a genome-wide correlation between SD strength and ribosome occupancy that was previously masked by other contributing factors. Using the antibiotic retapamulin to trap initiation complexes at start codons, we find that the mutant ribosomes select start sites correctly, arguing that start sites are hard-wired for initiation through the action of other mRNA features. We show that A-rich sequences upstream of start codons promote initiation. Taken together, our genome-wide study reveals that SD motifs are not necessary for ribosomes to determine where initiation occurs, though they do affect how efficiently initiation occurs.

Synonymous codon substitutions perturb cotranslational protein folding in vivo and impair cell fitness

In the cell, proteins are synthesized from N to C terminus and begin to fold during translation. Cotranslational folding mechanisms are therefore linked to elongation rate, which varies as a function of synonymous codon usage. However, synonymous codon substitutions can affect many distinct cellular processes, which has complicated attempts to deconvolve the extent to which synonymous codon usage can promote or frustrate proper protein folding in vivo. Although previous studies have shown that some synonymous changes can lead to different final structures, other substitutions will likely be more subtle, perturbing predominantly the protein folding pathway without radically altering the final structure. Here we show that synonymous codon substitutions encoding a single essential enzyme lead to dramatically slower cell growth. These mutations do not prevent active enzyme formation; instead, they predominantly alter the protein folding mechanism, leading to enhanced degradation in vivo. These results support a model in which synonymous codon substitutions can impair cell fitness by significantly perturbing cotranslational protein folding mechanisms, despite the chaperoning provided by the cellular protein homeostasis network.

How the initiating ribosome copes with ppGpp to translate mRNAs

During host colonization, bacteria use the alarmones (p)ppGpp to reshape their proteome by acting pleiotropically on DNA, RNA, and protein synthesis. Here, we elucidate how the initiating ribosome senses the cellular pool of guanosine nucleotides and regulates the progression towards protein synthesis. Our results show that the affinity of guanosine triphosphate (GTP) and the inhibitory concentration of ppGpp for the 30S-bound initiation factor IF2 vary depending on the programmed mRNA. The TufA mRNA enhanced GTP affinity for 30S complexes, resulting in improved ppGpp tolerance and allowing efficient protein synthesis. Conversely, the InfA mRNA allowed ppGpp to compete with GTP for IF2, thus stalling 30S complexes. Structural modeling and biochemical analysis of the TufA mRNA unveiled a structured enhancer of translation initiation (SETI) composed of two consecutive hairpins proximal to the translation initiation region (TIR) that largely account for ppGpp tolerance under physiological concentrations of guanosine nucleotides. Furthermore, our results show that the mechanism enhancing ppGpp tolerance is not restricted to the TufA mRNA, as similar ppGpp tolerance was found for the SETI-containing Rnr mRNA. Finally, we show that IF2 can use pppGpp to promote the formation of 30S initiation complexes (ICs), albeit requiring higher factor concentration and resulting in slower transitions to translation elongation. Altogether, our data unveil a novel regulatory mechanism at the onset of protein synthesis that tolerates physiological concentrations of ppGpp and that bacteria can exploit to modulate their proteome as a function of the nutritional shift happening during stringent response and infection.

A short translational ramp determines the efficiency of protein synthesis

Translation initiation is a major rate-limiting step for protein synthesis. However, recent studies strongly suggest that the efficiency of protein synthesis is additionally regulated by multiple factors that impact the elongation phase. To assess the influence of early elongation on protein synthesis, we employed a library of more than 250,000 reporters combined with in vitro and in vivo protein expression assays. Here we report that the identity of the amino acids encoded by codons 3 to 5 impact protein yield. This effect is independent of tRNA abundance, translation initiation efficiency, or overall mRNA structure. Single-molecule measurements of translation kinetics revealed pausing of the ribosome and aborted protein synthesis on codons 4 and 5 of distinct amino acid and nucleotide compositions. Finally, introduction of preferred sequence motifs only at specific codon positions improves protein synthesis efficiency for recombinant proteins. Collectively, our data underscore the critical role of early elongation events in translational control of gene expression.

Several factors contribute to the efficiency of protein expression. Here the authors show that the identity of amino acids encoded by codons at position 3–5 significantly impact translation efficiency and protein expression levels.

From reporters to endogenous genes: the impact of the first five codons on translation efficiency in Escherichia coli

Translation initiation is a critical step in the regulation of protein synthesis, and it is subjected to different control mechanisms, such as 5' UTR secondary structure and initiation codon context, that can influence the rates at which initiation and consequentially translation occur. For some genes, translation elongation also affects the rate of protein synthesis. With a GFP library containing nearly all possible combinations of nucleotides from the 3rd to the 5th codon positions in the protein coding region of the mRNA, it was previously demonstrated that some nucleotide combinations increased GFP expression up to four orders of magnitude. While it is clear that the codon region from positions 3 to 5 can influence protein expression levels of artificial constructs, its impact on endogenous proteins is still unknown. Through bioinformatics analysis, we identified the nucleotide combinations of the GFP library in Escherichia coli genes and examined the correlation between the expected levels of translation according to the GFP data with the experimental measures of protein expression. We observed that E. coli genes were enriched with the nucleotide compositions that enhanced protein expression in the GFP library, but surprisingly, it seemed to affect the translation efficiency only marginally. Nevertheless, our data indicate that different enterobacteria present similar nucleotide composition enrichment as E. coli, suggesting an evolutionary pressure towards the conservation of short translational enhancer sequences.

On the role of mRNA secondary structure in bacterial translation

Messenger RNA (mRNA) is no longer considered as a mere informational molecule whose sole function is to convey the genetic information specified by DNA to the ribosome. Beyond this primary function, mRNA also contains additional instructions that influence the way and the extent to which this message is translated by the ribosome into protein(s). Indeed, owing to its intrinsic propensity to quickly and dynamically fold and form higher order structures, mRNA exhibits a second layer of structural information specified by the sequence itself. Besides influencing transcription and mRNA stability, this additional information also affects translation, and more precisely the frequency of translation initiation, the choice of open reading frame by recoding, the elongation speed, and the folding of the nascent protein. Many studies in bacteria have shown that mRNA secondary structure participates to the rapid adaptation of these versatile organisms to changing environmental conditions by efficiently tuning translation in response to diverse signals, such as the presence of ligands, regulatory proteins, or small RNAs. This article is categorized under: Regulatory RNAs/RNAi/Riboswitches > Regulatory RNAs RNA Structure and Dynamics > Influence of RNA Structure in Biological Systems Translation > Translation Regulation.

Global analysis of protein synthesis in Flavobacterium johnsoniae reveals the use of Kozak-like sequences in diverse bacteria

In all cells, initiation of translation is tuned by intrinsic features of the mRNA. Here, we analyze translation in Flavobacterium johnsoniae, a representative of the Bacteroidetes. Members of this phylum naturally lack Shine-Dalgarno (SD) sequences in their mRNA, and yet their ribosomes retain the conserved anti-SD sequence. Translation initiation is tuned by mRNA secondary structure and by the identities of several key nucleotides upstream of the start codon. Positive determinants include adenine at position -3, reminiscent of the Kozak sequence of Eukarya. Comparative analysis of Escherichia coli reveals use of the same Kozak-like sequence to enhance initiation, suggesting an ancient and widespread mechanism. Elimination of contacts between A-3 and the conserved β-hairpin of ribosomal protein uS7 fails to diminish the contribution of A-3 to initiation, suggesting an indirect mode of recognition. Also, we find that, in the Bacteroidetes, the trinucleotide AUG is underrepresented in the vicinity of the start codon, which presumably helps compensate for the absence of SD sequences in these organisms.

Translational coupling via termination-reinitiation in archaea and bacteria

The genomes of many prokaryotes contain substantial fractions of gene pairs with overlapping stop and start codons (ATGA or TGATG). A potential benefit of overlapping gene pairs is translational coupling. In 720 genomes of archaea and bacteria representing all major phyla, we identify substantial, albeit highly variable, fractions of co-directed overlapping gene pairs. Various patterns are observed for the utilization of the SD motif for de novo initiation at upstream genes versus reinitiation at overlapping gene pairs. We experimentally test the predicted coupling in 9 gene pairs from the archaeon Haloferax volcanii and 5 gene pairs from the bacterium Escherichia coli. In 13 of 14 cases, translation of both genes is strictly coupled. Mutational analysis of SD motifs located upstream of the downstream genes indicate that the contribution of the SD to translational coupling widely varies from gene to gene. The nearly universal, abundant occurrence of overlapping gene pairs suggests that tight translational coupling is widespread in archaea and bacteria.

Archaea and bacteria often have gene pairs with overlapping stop and start codons, suggesting translational coupling. Here, Huber et al. analyse overlapping gene pairs from 720 genomes, and validate translational coupling via termination-reinitiation for 14 gene pairs in Haloferax volcanii and Escherichia coli.

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.

High-throughput determination of RNA structures

RNA performs and regulates a diverse range of cellular processes, with new functional roles being uncovered at a rapid pace. Interest is growing in how these functions are linked to RNA structures that form in the complex cellular environment. A growing suite of technologies that use advances in RNA structural probes, high-throughput sequencing and new computational approaches to interrogate RNA structure at unprecedented throughput are beginning to provide insights into RNA structures at new spatial, temporal and cellular scales.

Multifaceted Stoichiometry Control of Bacterial Operons Revealed by Deep Proteome Quantification

More than half of the protein-coding genes in bacteria are organized in polycistronic operons composed of two or more genes. It remains under debate whether the operon organization maintains the stoichiometric expression of the genes within an operon. In this study, we performed a label-free data-independent acquisition hyper reaction monitoring mass-spectrometry (HRM-MS) experiment to quantify the Escherichia coli proteome in exponential phase and quantified 93.6% of the cytosolic proteins, covering 67.9% and 56.0% of the translating polycistronic operons in BW25113 and MG1655 strains, respectively. We found that the translational regulation contributes largely to the proteome complexity: the shorter operons tend to be more tightly controlled for stoichiometry than longer operons; the operons which mainly code for complexes is more tightly controlled for stoichiometry than the operons which mainly code for metabolic pathways. The gene interval (distance between adjacent genes in one operon) may serve as a regulatory factor for stoichiometry. The catalytic efficiency might be a driving force for differential expression of enzymes encoded in one operon. These results illustrated the multifaceted nature of the operon regulation: the operon unified transcriptional level and gene-specific translational level. This multi-level regulation benefits the host by optimizing the efficiency of the productivity of metabolic pathways and maintenance of different types of protein complexes.

Programming Structured DNA Assemblies to Probe Biophysical Processes

Structural DNA nanotechnology is beginning to emerge as a widely accessible research tool to mechanistically study diverse biophysical processes. Enabled by scaffolded DNA origami in which a long single strand of DNA is weaved throughout an entire target nucleic acid assembly to ensure its proper folding, assemblies of nearly any geometric shape can now be programmed in a fully automatic manner to interface with biology on the 1-100-nm scale. Here, we review the major design and synthesis principles that have enabled the fabrication of a specific subclass of scaffolded DNA origami objects called wireframe assemblies. These objects offer unprecedented control over the nanoscale organization of biomolecules, including biomolecular copy numbers, presentation on convex or concave geometries, and internal versus external functionalization, in addition to stability in physiological buffer. To highlight the power and versatility of this synthetic structural biology approach to probing molecular and cellular biophysics, we feature its application to three leading areas of investigation: light harvesting and nanoscale energy transport, RNA structural biology, and immune receptor signaling, with an outlook toward unique mechanistic insight that may be gained in these areas in the coming decade.

A systematically-revised ribosome profiling method for bacteria reveals pauses at single-codon resolution

In eukaryotes, ribosome profiling provides insight into the mechanism of protein synthesis at the codon level. In bacteria, however, the method has been more problematic and no consensus has emerged for how to best prepare profiling samples. Here, we identify the sources of these problems and describe new solutions for arresting translation and harvesting cells in order to overcome them. These improvements remove confounding artifacts and improve the resolution to allow analyses of ribosome behavior at the codon level. With a clearer view of the translational landscape in vivo, we observe that filtering cultures leads to translational pauses at serine and glycine codons through the reduction of tRNA aminoacylation levels. This observation illustrates how bacterial ribosome profiling studies can yield insight into the mechanism of protein synthesis at the codon level and how these mechanisms are regulated in response to changes in the physiology of the cell.

Protocol for Ribosome Profiling in Bacteria

Ribosome profiling provides information on the position of ribosomes on mRNA on a genomic scale. Although this information is often used to detect changes in gene expression under different conditions, it also has great potential for yielding insight into the mechanism and regulation of protein synthesis itself. First developed in yeast, ribosome profiling involves the isolation and sequencing of ribosome-protected mRNA fragments generated by nuclease treatment. Since the application of ribosome profiling in bacteria has been problematic, we report here a systematically optimized protocol for E. coli that we have used with success for other bacteria as well. Cells are harvested by flash-freezing cultures directly in liquid nitrogen. After lysis, translation is arrested by high magnesium concentration without the use of antibiotics. These improvements eliminate artifacts induced by harvesting cells by centrifugation or filtration and by use of chloramphenicol to arrest translation. These improvements are especially appropriate for studies where the exact position of the ribosome is critical, and not merely the number of ribosomes per message, such as studies aimed at monitoring differences in local elongation rates.

Genome-wide probing RNA structure with the modified DMS-MaPseq in Arabidopsis

Transcripts have intrinsic propensity to form stable secondary structure that is fundamental to regulate RNA transcription, splicing, translation, RNA localization and turnover. Numerous methods that integrate chemical reactions with next-generation sequencing (NGS) have been applied to study in vivo RNA structure, providing new insights into RNA biology. Dimethyl sulfate (DMS) probing coupled with mutational profiling through NGS (DMS-MaPseq) is a newly developed method for revealing genome-wide or target-specific RNA structure. Herein, we present our experimental protocol of a modified DMS-MaPseq method for plant materials. The DMS treatment condition was optimized, and library preparation procedures were simplified. We also provided custom scripts for bioinformatic analysis of genome-wide DMS-MaPseq data. Bioinformatic results showed that our method could generate high-quality and reproducible data. Further, we assessed sequencing depth and coverage for genome-wide RNA structure profiling in Arabidopsis, and provided two examples of in vivo structure of mobile RNAs. We hope that our modified DMS-MaPseq method will serve as a powerful tool for analyzing in vivo RNA structurome in plants.

Revealing lncRNA Structures and Interactions by Sequencing-Based Approaches

Long noncoding RNAs (lncRNAs) have emerged as significant players in almost every level of gene function and regulation. Thus, characterizing the structures and interactions of lncRNAs is essential for understanding their mechanistic roles in cells. Through a combination of (bio)chemical approaches and automated capillary and high-throughput sequencing (HTS), the complexity and diversity of RNA structures and interactions has been revealed in the transcriptomes of multiple species. These methods have uncovered important biological insights into the mechanistic and functional roles of lncRNA in gene expression and RNA metabolism, as well as in development and disease. In this review, we summarize the latest sequencing strategies to reveal RNA structure, RNA-RNA, RNA-DNA, and RNA-protein interactions, and highlight the recent applications of these approaches to map functional lncRNAs. We discuss the advantages and limitations of these strategies, and provide recommendations to further advance methodologies capable of mapping RNA structure and interactions in order to discover new biology of lncRNAs and decipher their molecular mechanisms and implication in diseases.

Ribosome Profiling of Synechocystis Reveals Altered Ribosome Allocation at Carbon Starvation

Ribosome profiling accesses the translational step of gene expression via deep sequencing of ribosome-protected mRNA footprints. Pairing of ribosome profiling and transcriptomics data provides a translational efficiency for each gene. Here, the translatome and transcriptome of the model cyanobacterium Synechocystis were compared under carbon-replete and carbon starvation conditions. The latter may be experienced when cyanobacteria are cultivated in poorly mixed bioreactors or engineered to be product-secreting cell factories. A small fraction of genes (<200), including stress response genes, showed changes in translational efficiency during carbon starvation, indicating condition-dependent translation-level regulation. We observed ribosome occupancy in untranslated regions, possibly due to an alternative translation initiation mechanism in Synechocystis. The higher proportion of ribosomes residing in untranslated regions during carbon starvation may be a mechanism to quickly inactivate superfluous ribosomes. This work provides the first ribosome profiling data for cyanobacteria and reveals new regulation strategies for coping with nutrient limitation.

Technique Development for Probing RNA Structure In Vivo and Genome-Wide

How organisms perceive and respond to their surroundings is one of the great questions in biology. It is clear that RNA plays key roles in sensing. Cellular and environmental cues that RNA responds to include temperature, ions, metabolites, and biopolymers. Recent advances in next-generation sequencing and in vivo chemical probing have provided unprecedented insights into RNA folding in vivo and genome-wide. Patterns of chemical reactivity have implicated control of gene expression by RNA and aided prediction of RNA structure. Central to these advances has been development of molecular biological and chemical techniques. Key advances are improvements in the quality, cost, and throughput of library preparation; availability of a wider array of chemicals for probing RNA structure in vivo; and robustness and user friendliness of data analysis. Insights from probing transcriptomes and future directions are provided.

RNA, Action through Interactions

As transcription of the human genome is quite pervasive, it is possible that many novel functions of the noncoding genome have yet to be identified. Often the noncoding genome's functions are carried out by their RNA transcripts, which may rely on their structures and/or extensive interactions with other molecules. Recent technology developments are transforming the fields of RNA biology from studying one RNA at a time to transcriptome-wide mapping of structures and interactions. Here, we highlight the recent advances in transcriptome-wide RNA interaction analysis. These technologies revealed surprising versatility of RNA to participate in diverse molecular systems. For example, tens of thousands of RNA-RNA interactions have been revealed in cultured cells as well as in mouse brain, including interactions between transposon-produced transcripts and mRNAs. In addition, most transcription start sites in the human genome are associated with noncoding RNA transcribed from other genomic loci. These recent discoveries expanded our understanding of RNAs' roles in chromatin organization, gene regulation, and intracellular signaling.

Prokaryotic coding regions have little if any specific depletion of Shine-Dalgarno motifs

The Shine-Dalgarno motif occurs in front of prokaryotic start codons, and is complementary to the 3’ end of the 16S ribosomal RNA. Hybridization between the Shine-Dalgarno sequence and the anti-Shine-Dalgarno region of the16S rRNA (CCUCCU) directs the ribosome to the start AUG of the mRNA for translation. Shine-Dalgarno-like motifs (AGGAGG in E. coli) are depleted from open reading frames of most prokaryotes. This may be because hybridization of the 16S rRNA at Shine-Dalgarnos inside genes would slow translation or induce internal initiation. However, we analyzed 128 species from diverse phyla where the 16S rRNA gene(s) lack the anti-Shine-Dalgarno sequence, and so the 16S rRNA is incapable of interacting with Shine-Dalgarno-like sequences. Despite this lack of an anti-Shine-Dalgarno, half of these species still displayed depletion of Shine-Dalgarno-like sequences when analyzed by previous methods. Depletion of the same G-rich sequences was seen by these methods even in eukaryotes, which do not use the Shine-Dalgarno mechanism. We suggest previous methods are partly detecting a non-specific depletion of G-rich sequences. Alternative informatics approaches show that most prokaryotes have only slight, if any, specific depletion of Shine-Dalgarno-like sequences from open reading frames. Together with recent evidence that ribosomes do not pause at ORF-internal Shine-Dalgarno motifs, these results suggest the presence of ORF-internal Shine-Dalgarno-like motifs may be inconsequential, perhaps because internal regions of prokaryotic mRNAs may be structurally “shielded” from translation initiation.

Re-annotation of 12,495 prokaryotic 16S rRNA 3' ends and analysis of Shine-Dalgarno and anti-Shine-Dalgarno sequences

We examined 20,648 prokaryotic unique taxids with respect to the annotation of the 3' end of the 16S rRNA, which contains the anti-Shine-Dalgarno sequence. We used the sequence of highly conserved helix 45 of the 16S rRNA as a guide. By this criterion, 8,153 annotated 3' ends correctly included the anti-Shine-Dalgarno sequence, but 12,495 were foreshortened or otherwise mis-annotated, missing part or all of the anti-Shine-Dalgarno sequence, which immediately follows helix 45. We re-annotated, giving a total of 20,648 16S rRNA 3' ends. The vast majority indeed contained a consensus anti-Shine-Dalgarno sequence, embedded in a highly conserved 13 base "tail". However, 128 exceptional organisms had either a variant anti-Shine-Dalgarno, or no recognizable anti-Shine-Dalgarno, in their 16S rRNA(s). For organisms both with and without an anti-Shine-Dalgarno, we identified the Shine-Dalgarno motifs actually enriched in front of each organism's open reading frames. This showed to what extent the Shine-Dalgarno motifs correlated with anti-Shine Dalgarno motifs. In general, organisms whose rRNAs lacked a perfect anti-Shine-Dalgarno motif also lacked a recognizable Shine-Dalgarno. For organisms whose 16S rRNAs contained a perfect anti-Shine-Dalgarno motif, a variety of results were obtained. We found one genus, Alteromonas, where several taxids apparently maintain two different types of 16S rRNA genes, with different, but conserved, antiSDs. The fact that some organisms do not seem to have or use Shine-Dalgarno motifs supports the idea that prokaryotes have other robust mechanisms for recognizing start codons for translation.

Analyses of mRNA structure dynamics identify embryonic gene regulatory programs

RNA folding plays a crucial role in RNA function. However, knowledge of the global structure of the transcriptome is limited to cellular systems at steady state, thus hindering the understanding of RNA structure dynamics during biological transitions and how it influences gene function. Here, we characterized mRNA structure dynamics during zebrafish development. We observed that on a global level, translation guides structure rather than structure guiding translation. We detected a decrease in structure in translated regions and identified the ribosome as a major remodeler of RNA structure in vivo. In contrast, we found that 3' untranslated regions (UTRs) form highly folded structures in vivo, which can affect gene expression by modulating microRNA activity. Furthermore, dynamic 3'-UTR structures contain RNA-decay elements, such as the regulatory elements in nanog and ccna1, two genes encoding key maternal factors orchestrating the maternal-to-zygotic transition. These results reveal a central role of RNA structure dynamics in gene regulatory programs.

Accessibility of the Shine-Dalgarno Sequence Dictates N-Terminal Codon Bias in E. coli

Despite considerable efforts, no physical mechanism has been shown to explain N-terminal codon bias in prokaryotic genomes. Using a systematic study of synonymous substitutions in two endogenous E. coli genes, we show that interactions between the coding region and the upstream Shine-Dalgarno (SD) sequence modulate the efficiency of translation initiation, affecting both intracellular mRNA and protein levels due to the inherent coupling of transcription and translation in E. coli. We further demonstrate that far-downstream mutations can also modulate mRNA levels by occluding the SD sequence through the formation of non-equilibrium secondary structures. By contrast, a non-endogenous RNA polymerase that decouples transcription and translation largely alleviates the effects of synonymous substitutions on mRNA levels. Finally, a complementary statistical analysis of the E. coli genome specifically implicates avoidance of intra-molecular base pairing with the SD sequence. Our results provide general physical insights into the coding-level features that optimize protein expression in prokaryotes.

Evolutionary Convergence of Pathway-Specific Enzyme Expression Stoichiometry

Coexpression of proteins in response to pathway-inducing signals is the founding paradigm of gene regulation. However, it remains unexplored whether the relative abundance of co-regulated proteins requires precise tuning. Here, we present large-scale analyses of protein stoichiometry and corresponding regulatory strategies for 21 pathways and 67-224 operons in divergent bacteria separated by 0.6-2 billion years. Using end-enriched RNA-sequencing (Rend-seq) with single-nucleotide resolution, we found that many bacterial gene clusters encoding conserved pathways have undergone massive divergence in transcript abundance and architectures via remodeling of internal promoters and terminators. Remarkably, these evolutionary changes are compensated post-transcriptionally to maintain preferred stoichiometry of protein synthesis rates. Even more strikingly, in eukaryotic budding yeast, functionally analogous proteins that arose independently from bacterial counterparts also evolved to convergent in-pathway expression. The broad requirement for exact protein stoichiometries despite regulatory divergence provides an unexpected principle for building biological pathways both in nature and for synthetic activities.

Extensive reshaping of bacterial operons by programmed mRNA decay

Bacterial operons synchronize the expression of multiple genes by placing them under the control of a shared promoter. It was previously shown that polycistronic transcripts can undergo differential RNA decay, leaving some genes within the polycistron more stable than others, but the extent of regulation by differential mRNA decay or its evolutionary conservation remains unknown. Here, we find that a substantial fraction of E. coli genes display non-uniform mRNA stoichiometries despite being coded from the same operon. We further show that these altered operon stoichiometries are shaped post-transcriptionally by differential mRNA decay, which is regulated by RNA structures that protect specific regions in the transcript from degradation. These protective RNA structures are generally coded within the protein-coding regions of the regulated genes and are frequently evolutionarily conserved. Furthermore, we provide evidence that differences in ribosome densities across polycistronic transcript segments, together with the conserved structural RNA elements, play a major role in the differential decay process. Our results highlight a major role for differential mRNA decay in shaping bacterial transcriptomes.

Bacteria utilize operonic transcription to synchronize the expression of multiple consecutive genes. However, this strategy lacks the ability to fine-tune the expression of specific operon members, which is often biologically important. In this report, we integrate multiple transcriptome-wide RNA-sequencing methods to show that bacteria commonly employ differential mRNA decay rates for genes residing within the same operon, generating differential transcript abundances for equally transcribed operon members, at steady state. By comparing the transcriptomes of different bacteria, we show that differential decay not only regulates the expression levels of hundreds of genes but also often evolutionarily conserved, providing support for its biological importance. By mapping the RNA termini positions at steady-state, we show that stabilized operon segments are protected from different RNases through a combination of protective RNA structures, which surprisingly, are often encoded within protein-coding regions and are evolutionarily conserved. In addition, we provide evidence that differential ribosome densities over the regulated operons guide the initial events in the differential decay mechanism. Our results highlight differential mRNA decay as a major shaping force of bacterial transcriptomes and gene regulatory programs.

Insights into the Fundamental Physiology of the Uncultured Fe-Oxidizing Bacterium Leptothrix ochracea

Leptothrix ochracea is known for producing large volumes of iron oxyhydroxide sheaths that alter wetland biogeochemistry. For over a century, these delicate structures have fascinated microbiologists and geoscientists. Because L. ochracea still resists long-term in vitro culture, the debate regarding its metabolic classification dates back to 1885. We developed a novel culturing technique for L. ochracea using in situ natural waters and coupled this with single-cell genomics and nanoscale secondary-ion mass spectrophotometry (nanoSIMS) to probe L. ochracea's physiology. In microslide cultures L. ochracea doubled every 5.7 h and had an absolute growth requirement for ferrous iron, the genomic capacity for iron oxidation, and a branched electron transport chain with cytochromes putatively involved in lithotrophic iron oxidation. Additionally, its genome encoded several electron transport chain proteins, including a molybdopterin alternative complex III (ACIII), a cytochrome bd oxidase reductase, and several terminal oxidase genes. L. ochracea contained two key autotrophic proteins in the Calvin-Benson-Bassham cycle, a form II ribulose bisphosphate carboxylase, and a phosphoribulose kinase. L. ochracea also assimilated bicarbonate, although calculations suggest that bicarbonate assimilation is a small fraction of its total carbon assimilation. Finally, L. ochracea's fundamental physiology is a hybrid of those of the chemolithotrophic Gallionella-type iron-oxidizing bacteria and the sheathed, heterotrophic filamentous metal-oxidizing bacteria of the Leptothrix-Sphaerotilus genera. This allows L. ochracea to inhabit a unique niche within the neutrophilic iron seeps.IMPORTANCE Leptothrix ochracea was one of three groups of organisms that Sergei Winogradsky used in the 1880s to develop his hypothesis on chemolithotrophy. L. ochracea continues to resist cultivation and appears to have an absolute requirement for organic-rich waters, suggesting that its true physiology remains unknown. Further, L. ochracea is an ecological engineer; a few L. ochracea cells can generate prodigious volumes of iron oxyhydroxides, changing the ecosystem's geochemistry and ecology. Therefore, to determine L. ochracea's basic physiology, we employed new single-cell techniques to demonstrate that L. ochracea oxidizes iron to generate energy and, despite having predicted genes for autotrophic growth, assimilates a fraction of the total CO₂ that autotrophs do. Although not a true chemolithoautotroph, L. ochracea's physiological strategy allows it to be flexible and to extensively colonize iron-rich wetlands.

A Stress Response that Monitors and Regulates mRNA Structure Is Central to Cold Shock Adaptation

Temperature influences the structural and functional properties of cellular components, necessitating stress responses to restore homeostasis following temperature shift. Whereas the circuitry controlling the heat shock response is well understood, that controlling the E. coli cold shock adaptation program is not. We found that during the growth arrest phase (acclimation) that follows shift to low temperature, protein synthesis increases, and open reading frame (ORF)-wide mRNA secondary structure decreases. To identify the regulatory system controlling this process, we screened for players required for increased translation. We identified a two-member mRNA surveillance system that enables recovery of translation during acclimation: RNase R assures appropriate mRNA degradation and the Csps dynamically adjust mRNA secondary structure to globally modulate protein expression level. An autoregulatory switch in which Csps tune their own expression to cellular demand enables dynamic control of global translation. The universality of Csps in bacteria suggests broad utilization of this control mechanism.

Rice In Vivo RNA Structurome Reveals RNA Secondary Structure Conservation and Divergence in Plants

RNA secondary structure plays a critical role in gene regulation. Rice (Oryza sativa) is one of the most important food crops in the world. However, RNA structure in rice has scarcely been studied. Here, we have successfully generated in vivo Structure-seq libraries in rice. We found that the structural flexibility of mRNAs might associate with the dynamics of biological function. Higher N⁶-methyladenosine (m⁶A) modification tends to have less RNA structure in 3' UTR, whereas GC content does not significantly affect in vivo mRNA structure to maintain efficient biological processes such as translation. Comparative analysis of RNA structurome between rice and Arabidopsis revealed that higher GC content does not lead to stronger structure and less RNA structural flexibility. Moreover, we found a weak correlation between sequence and structure conservation of the orthologs between rice and Arabidopsis. The conservation and divergence of both sequence and in vivo RNA structure corresponds to diverse and specific biological processes. Our results indicate that RNA secondary structure might offer a separate layer of selection to the sequence between monocot and dicot. Therefore, our study implies that RNA structure evolves differently in various biological processes to maintain robustness in development and adaptational flexibility during angiosperm evolution.

Pervasive Regulatory Functions of mRNA Structure Revealed by High-Resolution SHAPE Probing

mRNAs can fold into complex structures that regulate gene expression. Resolving such structures de novo has remained challenging and has limited our understanding of the prevalence and functions of mRNA structure. We use SHAPE-MaP experiments in living E. coli cells to derive quantitative, nucleotide-resolution structure models for 194 endogenous transcripts encompassing approximately 400 genes. Individual mRNAs have exceptionally diverse architectures, and most contain well-defined structures. Active translation destabilizes mRNA structure in cells. Nevertheless, mRNA structure remains similar between in-cell and cell-free environments, indicating broad potential for structure-mediated gene regulation. We find that the translation efficiency of endogenous genes is regulated by unfolding kinetics of structures overlapping the ribosome binding site. We discover conserved structured elements in 35% of UTRs, several of which we validate as novel protein binding motifs. RNA structure regulates every gene studied here in a meaningful way, implying that most functional structures remain to be discovered.

Functional 5' UTR mRNA structures in eukaryotic translation regulation and how to find them

RNA molecules can fold into intricate shapes that can provide an additional layer of control of gene expression beyond that of their sequence. In this Review, we discuss the current mechanistic understanding of structures in 5' untranslated regions (UTRs) of eukaryotic mRNAs and the emerging methodologies used to explore them. These structures may regulate cap-dependent translation initiation through helicase-mediated remodelling of RNA structures and higher-order RNA interactions, as well as cap-independent translation initiation through internal ribosome entry sites (IRESs), mRNA modifications and other specialized translation pathways. We discuss known 5' UTR RNA structures and how new structure probing technologies coupled with prospective validation, particularly compensatory mutagenesis, are likely to identify classes of structured RNA elements that shape post-transcriptional control of gene expression and the development of multicellular organisms.

Pausing of Chloroplast Ribosomes Is Induced by Multiple Features and Is Linked to the Assembly of Photosynthetic Complexes

Many mRNAs contain pause sites that briefly interrupt the progress of translation. Specific features that induce ribosome pausing have been described; however, their individual contributions to pause-site formation, and the overall biological significance of ribosome pausing, remain largely unclear. We have taken advantage of the compact genome of chloroplasts to carry out a plastid genome-wide survey of pause sites, as a basis for studying the impact of pausing on posttranslational processes. Based on ribosomal profiling of Arabidopsis (Arabidopsis thaliana) chloroplast mRNAs, we demonstrate that a combination of factors-mRNA secondary structure, internal Shine-Dalgarno sequences, and positively charged amino acids in the nascent peptide chain-explains 95% of the major pause sites on plastid mRNAs, whereas codon usage has little impact. The distribution of the pause sites is nonrandom and conforms to distinct patterns in the vicinity of sequences coding for transmembrane domains, which depend on their orientation within the membrane as well as being next to sequences coding for cofactor binding sites. We found strong indications that the mechanisms causing ribosomal pausing and at least some of the ribosomes pause sites are conserved between distantly related plant species. In addition, the positions of features that cause pausing are well conserved in photoautotrophic plants, but less so in their nonphotosynthetic, parasitic relatives, implying that the synthesis and assembly of photosynthetic multiprotein complexes requires localized ribosome pausing.

Bacterial RNA Biology on a Genome Scale

Bacteria are an exceedingly diverse group of organisms whose molecular exploration is experiencing a renaissance. While the classical view of bacterial gene expression was relatively simple, the emerging view is more complex, encompassing extensive post-transcriptional control involving riboswitches, RNA thermometers, and regulatory small RNAs (sRNAs) associated with the RNA-binding proteins CsrA, Hfq, and ProQ, as well as CRISPR/Cas systems that are programmed by RNAs. Moreover, increasing interest in members of the human microbiota and environmental microbial communities has highlighted the importance of understudied bacterial species with largely unknown transcriptome structures and RNA-based control mechanisms. Collectively, this creates a need for global RNA biology approaches that can rapidly and comprehensively analyze the RNA composition of a bacterium of interest. We review such approaches with a focus on RNA-seq as a versatile tool to investigate the different layers of gene expression in which RNA is made, processed, regulated, modified, translated, and turned over.

Local genic base composition impacts protein production and cellular fitness

The maintenance of a G + C content that is higher than the mutational input to a genome provides support for the view that selection serves to increase G + C contents in bacteria. Recent experimental evidence from Escherichia coli demonstrated that selection for increasing G + C content operates at the level of translation, but the precise mechanism by which this occurs is unknown. To determine the substrate of selection, we asked whether selection on G + C content acts across all sites within a gene or is confined to particular genic regions or nucleotide positions. We systematically altered the G + C contents of the GFP gene and assayed its effects on the fitness of strains harboring each variant. Fitness differences were attributable to the base compositional variation in the terminal portion of the gene, suggesting a connection to the folding of a specific protein feature. Variants containing sequence features that are thought to result in rapid translation, such as low G + C content and high levels of codon adaptation, displayed highly reduced growth rates. Taken together, our results show that purifying selection acting against A and T mutations most likely results from their tendency to increase the rate of translation, which can perturb the dynamics of protein folding.

Model-based genome-wide determination of RNA chain elongation rates in Escherichia coli

Dynamics in the process of transcription are often simplified, yet they play an important role in transcript folding, translation into functional protein and DNA supercoiling. While the modulation of the speed of transcription of individual genes and its role in regulation and proper protein folding has been analyzed in depth, the functional relevance of differences in transcription speeds as well as the factors influencing it have not yet been determined on a genome-wide scale. Here we determined transcription speeds for the majority of E. coli genes based on experimental data. We find large differences in transcription speed between individual genes and a strong influence of both cellular location as well as the relative importance of genes for cellular function on transcription speeds. Investigating factors influencing transcription speeds we observe both codon composition as well as factors associated to DNA topology as most important factors influencing transcription speeds. Moreover, we show that differences in transcription speeds are sufficient to explain the timing of regulatory responses during environmental shifts and highlight the importance of the consideration of transcription speeds in the design of experiments measuring transcriptomic responses to perturbations.