Secondary Structure across the Bacterial Transcriptome Reveals Versatile Roles in mRNA Regulation and Function

m6A modification is incorporated into bacterial mRNA without specific functional benefit

N 6-Methyladenosine (m6A), the most abundant modification in eukaryotic messenger RNAs (mRNAs), has also been found at a low level in bacterial mRNAs. However, enzyme(s) that introduce m6A modification on mRNAs in bacteria remain elusive. In this work, we combine deep-sequencing approaches that identify m6A sites with in vitro biochemical studies to identify putative m6A methyltransferases that would modify Escherichia coli mRNAs. We tested four uncharacterized candidates predicted to encode proteins with putative methyltransferase domains, whose deletion decreased the m6A level. However, in vitro analysis with the purified putative methyltransferases revealed that none of them installs m6A on mRNA. Exposure to heat and oxidative stress also changed the m6A level; however, we found no clear correlation between the m6A change and the specific stress. Considering two deep-sequencing approaches with different resolution, we found that m6A methylation on bacterial mRNAs is very low and appears randomly introduced. These results suggest that, in contrast to eukaryotes, the m6A modification in bacterial mRNA lacks a direct enzymatic recognition mechanism and has no clear biological function.

Increasing recombinant protein production in E. coli via FACS-based selection of N-terminal coding DNA libraries

Here, we present a previously undescribed approach to modify N-terminal sequences of recombinant proteins to increase their production yield in Escherichia coli. Prior research has demonstrated that the nucleotides immediately following the start codon can significantly influence protein expression. However, the impact of these sequences is construct-specific and is not universally applicable to all proteins. Most of the previous research has been limited to selecting from a few rationally designed sequences. In contrast, we used a directed evolution-based methodology, screening large numbers of diversified sequences derived from DNA libraries coding for the N-termini of investigated proteins. To facilitate the identification of cells with increased expression of the target construct, we cloned a GFP gene at the C-terminus of the expressed genes and used fluorescent activated cell sorting (FACS) to separate cells based on their fluorescence. By following this systematic workflow, we successfully elevated the yield of soluble recombinant proteins of multiple constructs up to over 30-fold.

The human mitochondrial mRNA structurome reveals mechanisms of gene expression

The human mitochondrial genome encodes crucial oxidative phosphorylation system proteins, pivotal for aerobic energy transduction. They are translated from nine monocistronic and two bicistronic transcripts whose native structures remain unexplored, posing a gap in understanding mitochondrial gene expression. In this work, we devised the mitochondrial dimethyl sulfate mutational profiling with sequencing (mitoDMS-MaPseq) method and applied detection of RNA folding ensembles using expectation-maximization (DREEM) clustering to unravel the native mitochondrial messenger RNA (mt-mRNA) structurome in wild-type (WT) and leucine-rich pentatricopeptide repeat-containing protein (LRPPRC)-deficient cells. Our findings elucidate LRPPRC's role as a holdase contributing to maintaining mt-mRNA folding and efficient translation. mt-mRNA structural insights in WT mitochondria, coupled with metabolic labeling, unveil potential mRNA-programmed translational pausing and a distinct programmed ribosomal frameshifting mechanism. Our data define a critical layer of mitochondrial gene expression regulation. These mt-mRNA folding maps provide a reference for studying mt-mRNA structures in diverse physiological and pathological contexts.

RNase E searches for cleavage sites in RNA by linear diffusion: direct evidence from single-molecule FRET

The ability of obstacles in cellular transcripts to protect downstream but not upstream sites en masse from attack by RNase E has prompted the hypothesis that this mRNA-degrading endonuclease may scan 5'-monophosphorylated RNA linearly for cleavage sites, starting at the 5' end. However, despite its proposed regulatory importance, the migration of RNase E on RNA has never been directly observed. We have now used single-molecule FRET to monitor the dynamics of this homotetrameric enzyme on RNA. Our findings reveal that RNase E slides along unpaired regions of RNA without consuming a molecular source of energy such as ATP and that its forward progress can be impeded when it encounters a large structural discontinuity. This movement, which is bidirectional, occurs in discrete steps of variable length and requires an RNA ligand much longer than needed to occupy a single RNase E subunit. These results indicate that RNase E scans for cleavage sites by one-dimensional diffusion and suggest a possible molecular mechanism.

A new look at the fluorescent protein-based approach for identifying optimal coding sequence for recombinant protein expression in E. coli

Due to the degeneracy of the genetic code, most amino acids are encoded by several codons. The choice among synonymous codons at the N-terminus of genes has a profound effect on protein expression in Escherichia coli. This is often explained by the different contributions of synonymous codons to mRNA secondary structure formation. Strong secondary structures at the 5'-end of mRNA interfere with ribosome binding and affect the process of translation initiation. In silico optimization of the gene 5'-end can significantly increase the level of protein expression; however, this method is not always effective due to the uncertainty of the exact mechanism by which synonymous substitutions affect expression; thus, it may produce nonoptimal variants as well as miss some of the best producers. In this paper, an alternative approach is proposed based on screening a partially randomized library of expression constructs comprising hundreds of selected synonymous variants. The effect of such substitutions was evaluated using the gene of interest fused to the reporter gene of the fluorescent protein with subsequent screening for the most promising candidates according to the reporter's signal intensity. The power of the approach is demonstrated by a significant increase in the prokaryotic expression of three proteins: canine cystatin C, human BCL2-associated athanogene 3 and human cardiac troponin I. This simple approach was suggested which may provide an efficient, easy, and inexpensive optimization method for poorly expressed proteins in bacteria.

Direct and indirect control of Rho-dependent transcription termination by the Escherichia coli lysC riboswitch

Bacterial riboswitches are molecular structures that play a crucial role in controlling gene expression to maintain cellular balance. The Escherichia coli lysC riboswitch has been previously shown to regulate gene expression through translation initiation and mRNA decay. Recent research suggests that lysC gene expression is also influenced by Rho-dependent transcription termination. Through a series of in silico, in vitro, and in vivo experiments, we provide experimental evidence that the lysC riboswitch directly and indirectly modulates Rho transcription termination. Our study demonstrates that Rho-dependent transcription termination plays a significant role in the cotranscriptional regulation of lysC expression. Together with previous studies, our work suggests that lysC expression is governed by a lysine-sensing riboswitch that regulates translation initiation, transcription termination, and mRNA degradation. Notably, both Rho and RNase E target the same region of the RNA molecule, implying that RNase E may degrade Rho-terminated transcripts, providing a means to selectively eliminate these incomplete messenger RNAs. Overall, this study sheds light on the complex regulatory mechanisms used by bacterial riboswitches, emphasizing the role of transcription termination in the control of gene expression and mRNA stability.

RNA-seq reveals multifaceted gene expression response to Fab production in Escherichia coli fed-batch processes with particular focus on ribosome stalling

BACKGROUND

Escherichia coli is a cost-effective expression system for production of antibody fragments like Fabs. Various yield improvement strategies have been applied, however, Fabs remain challenging to produce. This study aimed to characterize the gene expression response of commonly used E. coli strains BL21(DE3) and HMS174(DE3) to periplasmic Fab expression using RNA sequencing (RNA-seq). Two Fabs, Fabx and FTN2, fused to a post-translational translocation signal sequence, were produced in carbon-limited fed-batch cultivations.

RESULTS

Production of Fabx impeded cell growth substantially stronger than FTN2 and yields of both Fabs differed considerably. The most noticeable, common changes in Fab-producing cells suggested by our RNA-seq data concern the cell envelope. The Cpx and Psp stress responses, both connected to inner membrane integrity, were activated, presumably by recombinant protein aggregation and impairment of the Sec translocon. The data additionally suggest changes in lipopolysaccharide synthesis, adjustment of membrane permeability, and peptidoglycan maturation and remodeling. Moreover, all Fab-producing strains showed depletion of Mg2+, indicated by activation of the PhoQP two-component signal transduction system during the early stage and sulfur and phosphate starvation during the later stage of the process. Furthermore, our data revealed ribosome stalling, caused by the Fabx amino acid sequence, as a contributor to low Fabx yields. Increased Fabx yields were obtained by a site-specific amino acid exchange replacing the stalling sequence. Contrary to expectations, cell growth was not impacted by presence or removal of the stalling sequence. Considering ribosome rescue is a conserved mechanism, the substantial differences observed in gene expression between BL21(DE3) and HMS174(DE3) in response to ribosome stalling on the recombinant mRNA were surprising.

CONCLUSIONS

Through characterization of the gene expression response to Fab production under industrially relevant cultivation conditions, we identified potential cell engineering targets. Thereby, we hope to enable rational approaches to improve cell fitness and Fab yields. Furthermore, we highlight ribosome stalling caused by the amino acid sequence of the recombinant protein as a possible challenge during recombinant protein production.

Genome-Scale Analysis of the Structure and Function of RNA Pathways and Networks in Pseudomonas aeruginosa

In recent years, several genome-wide approaches based on RNA sequencing (RNA-seq) have been developed. These methods allow a comprehensive and dynamic view of the structure and function of the multi-layered RNA pathways and networks. Many of these approaches, including the promising one of single-cell transcriptome analysis, have been successfully applied to Pseudomonas aeruginosa. However, we are only at the beginning because only a few surrounding conditions have been considered. Here, we aim to illustrate the different types of approaches based on RNA-seq that will lead us in the future to a better understanding of the dynamics of RNA biology in P. aeruginosa.

Decoding stoichiometric protein synthesis in E. coli through translation rate parameters

E. coli is one of the most widely used organisms for understanding the principles of cellular and molecular genetics. However, we are yet to understand the origin of several experimental observations related to the regulation of gene expression in E. coli. One of the prominent examples in this context is the proportional synthesis in multiprotein complexes where all of their obligate subunits are produced in proportion to their stoichiometry. In this work, by combining the next-generation sequencing data with the stochastic simulations of protein synthesis, we explain the origin of proportional protein synthesis in multicomponent complexes. We find that the estimated initiation rates for the translation of all subunits in those complexes are proportional to their stoichiometry. This constraint on protein synthesis kinetics enforces proportional protein synthesis without requiring any feedback mechanism. We also find that the translation initiation rates in E. coli are influenced by the coding sequence length and the enrichment of A and C nucleotides near the start codon. Thus, this study rationalizes the role of conserved and nonrandom features of genes in regulating the translation kinetics and unravels a key principle of the regulation of protein synthesis.

The human mitochondrial mRNA structurome reveals mechanisms of gene expression

The mammalian mitochondrial genome encodes thirteen oxidative phosphorylation system proteins, crucial in aerobic energy transduction. These proteins are translated from 9 monocistronic and 2 bicistronic transcripts, whose native structures remain unexplored, leaving fundamental molecular determinants of mitochondrial gene expression unknown. To address this gap, we developed a mitoDMS-MaPseq approach and used DREEM clustering to resolve the native human mitochondrial mt-mRNA structurome. We gained insights into mt-mRNA biology and translation regulatory mechanisms, including a unique programmed ribosomal frameshifting for the ATP8/ATP6 transcript. Furthermore, absence of the mt-mRNA maintenance factor LRPPRC led to a mitochondrial transcriptome structured differently, with specific mRNA regions exhibiting increased or decreased structuredness. This highlights the role of LRPPRC in maintaining mRNA folding to promote mt-mRNA stabilization and efficient translation. In conclusion, our mt-mRNA folding maps reveal novel mitochondrial gene expression mechanisms, serving as a detailed reference and tool for studying them in different physiological and pathological contexts.

Sequence specificity defines the effectiveness of PPMOs targeting Pseudomonas aeruginosa

Development of new therapeutics against antibiotic resistant pathogenic bacteria is recognized as a priority across the globe. We have reported using peptide-conjugated phosphorodiamidate morpholino oligomers (PPMOs) as species-specific antibiotics. The oligo sequences, 11 bases are designed to be complementary to specific essential genes near the Shine-Dalgarno site and inhibit translation. Here, we analyzed target specificity and the impact of genetic mutations on lead PPMOs targeting the rpsJ or acpP gene of Pseudomonas aeruginosa. Mutants in P. aeruginosa PAO1 were generated with four, two, or one base-pair mutations within the 11-base target sequence of the rpsJ gene. All mutants exhibited increased MICs compared to wild-type PAO1 when treated with the RpsJ PPMO, and the increase in the MICs was proportional to the number of base-pair mutations. Among single base-pair mutants, mutations in the middle of the sequence were more impactful than mutations in 5' or 3' end of the sequence. The increased MICs shown by the rpsJ mutants could be reversed by PPMOs designed to target the mutated rpsJ sequence. BALB/c mice infected intratracheally with mutants demonstrated increased lung burden when treated with RpsJ PPMO compared to wild-type PAO1-infected mice treated with RpsJ PPMO. Treating mice with a PPMOs designed to specifically target the mutant sequence was more effective against these mutant strains. These experiments confirm target specificity of two lead P. aeruginosa PPMOs and illustrate one potential mechanism of resistance that could emerge from an antisense approach.

Transcriptome-wide probing reveals RNA thermometers that regulate translation of glycerol permease genes in Bacillus subtilis

RNA structure regulates bacterial gene expression by several distinct mechanisms via environmental and cellular stimuli, one of which is temperature. While some genome-wide studies have focused on heat shock treatments and the subsequent transcriptomic changes, soil bacteria are less likely to experience such rapid and extreme temperature changes. Though RNA thermometers (RNATs) have been found in 5' untranslated leader regions (5' UTRs) of heat shock and virulence-associated genes, this RNA-controlled mechanism could regulate other genes as well. Using Structure-seq2 and the chemical probe dimethyl sulfate (DMS) at four growth temperatures ranging from 23°C to 42°C, we captured a dynamic response of the Bacillus subtilis transcriptome to temperature. Our transcriptome-wide results show RNA structural changes across all four temperatures and reveal nonmonotonic reactivity trends with increasing temperature. Then, focusing on subregions likely to contain regulatory RNAs, we examined 5' UTRs to identify large, local reactivity changes. This approach led to the discovery of RNATs that control the expression of glpF (glycerol permease) and glpT (glycerol-3-phosphate permease); expression of both genes increased with increased temperature. Results with mutant RNATs indicate that both genes are controlled at the translational level. Increased import of glycerols at high temperatures could provide thermoprotection to proteins.

Effect of mRNA/tRNA mutations on translation speed: Implications for human diseases

Recent discoveries establish tRNAs as central regulators of mRNA translation dynamics, and therefore cotranslational folding and function of the encoded protein. The tRNA pool, whose composition and abundance change in a cell- and tissue-dependent manner, is the main factor which determines mRNA translation velocity. In this review, we discuss a group of pathogenic mutations, in the coding sequences of either protein-coding genes or in tRNA genes, that alter mRNA translation dynamics. We also summarize advances in tRNA biology that have uncovered how variations in tRNA levels on account of genetic mutations affect protein folding and function, and thereby contribute to phenotypic diversity in clinical manifestations.

Led-Seq: ligation-enhanced double-end sequence-based structure analysis of RNA

Structural analysis of RNA is an important and versatile tool to investigate the function of this type of molecules in the cell as well as in vitro. Several robust and reliable procedures are available, relying on chemical modification inducing RT stops or nucleotide misincorporations during reverse transcription. Others are based on cleavage reactions and RT stop signals. However, these methods address only one side of the RT stop or misincorporation position. Here, we describe Led-Seq, a new approach based on lead-induced cleavage of unpaired RNA positions, where both resulting cleavage products are investigated. The RNA fragments carrying 2', 3'-cyclic phosphate or 5'-OH ends are selectively ligated to oligonucleotide adapters by specific RNA ligases. In a deep sequencing analysis, the cleavage sites are identified as ligation positions, avoiding possible false positive signals based on premature RT stops. With a benchmark set of transcripts in Escherichia coli, we show that Led-Seq is an improved and reliable approach based on metal ion-induced phosphodiester hydrolysis to investigate RNA structures in vivo.

Role of the 5' end phosphorylation state for small RNA stability and target RNA regulation in bacteria

In enteric bacteria, several small RNAs (sRNAs) including MicC employ endoribonuclease RNase E to stimulate target RNA decay. A current model proposes that interaction of the sRNA 5' monophosphate (5'P) with the N-terminal sensing pocket of RNase E allosterically activates cleavage of the base-paired target in the active site. In vivo evidence supporting this model is lacking. Here, we engineered a genetic tool allowing us to generate 5' monophosphorylated sRNAs of choice in a controllable manner in the cell. Four sRNAs were tested and none performed better in target destabilization when 5' monophosphorylated. MicC retains full activity even when RNase E is defective in 5'P sensing, whereas regulation is lost upon removal of its scaffolding domain. Interestingly, sRNAs MicC and RyhB that originate with a 5' triphosphate group are dramatically destabilized when 5' monophosphorylated, but stable when in 5' triphosphorylated form. In contrast, the processing-derived sRNAs CpxQ and SroC, which carry 5'P groups naturally, are highly stable. Thus, the 5' phosphorylation state determines stability of naturally triphosphorylated sRNAs, but plays no major role for target RNA destabilization in vivo. In contrast, the RNase E C-terminal half is crucial for MicC-mediated ompD decay, suggesting that interaction with Hfq is mandatory.

RNase Y Autoregulates Its Synthesis in Bacillus subtilis

The instability of messenger RNA is crucial to the control of gene expression. In Bacillus subtilis, RNase Y is the major decay-initiating endoribonuclease. Here, we show how this key enzyme regulates its own synthesis by modulating the longevity of its mRNA. Autoregulation is achieved through cleavages in two regions of the rny (RNase Y) transcript: (i) within the first ~100 nucleotides of the open reading frame, immediately inactivating the mRNA for further rounds of translation; (ii) cleavages in the rny 5' UTR, primarily within the 5'-terminal 50 nucleotides, creating entry sites for the 5' exonuclease J1 whose progression is blocked around position -15 of the rny mRNA, potentially by initiating ribosomes. This links the functional inactivation of the transcript by RNase J1 to translation efficiency, depending on the ribosome occupancy at the translation initiation site. By these mechanisms, RNase Y can initiate degradation of its own mRNA when the enzyme is not occupied with degradation of other RNAs and thus prevent its overexpression beyond the needs of RNA metabolism.

Escherichia coli CspA stimulates translation in the cold of its own mRNA by promoting ribosome progression

Escherichia coli CspA is an RNA binding protein that accumulates during cold-shock and stimulates translation of several mRNAs-including its own. Translation in the cold of cspA mRNA involves a cis-acting thermosensor element, which enhances ribosome binding, and the trans-acting action of CspA. Using reconstituted translation systems and probing experiments we show that, at low temperature, CspA specifically promotes the translation of the cspA mRNA folded in the conformation less accessible to the ribosome, which is formed at 37°C but is retained upon cold shock. CspA interacts with its mRNA without inducing large structural rearrangements, but allowing the progression of the ribosomes during the transition from translation initiation to translation elongation. A similar structure-dependent mechanism may be responsible for the CspA-dependent translation stimulation observed with other probed mRNAs, for which the transition to the elongation phase is progressively facilitated during cold acclimation with the accumulation of CspA.

Graded impact of obstacle size on scanning by RNase E

In countless bacterial species, the lifetimes of most mRNAs are controlled by the regulatory endonuclease RNase E, which preferentially degrades RNAs bearing a 5' monophosphate and locates cleavage sites within them by scanning linearly from the 5' terminus along single-stranded regions. Consequently, its rate of cleavage at distal sites is governed by any obstacles that it may encounter along the way, such as bound proteins or ribosomes or base pairing that is coaxial with the path traversed by this enzyme. Here, we report that the protection afforded by such obstacles is dependent on the size and persistence of the structural discontinuities they create, whereas the molecular composition of obstacles to scanning is of comparatively little consequence. Over a broad range of sizes, incrementally larger discontinuities are incrementally more protective, with corresponding effects on mRNA stability. The graded impact of such obstacles suggests possible explanations for why their effect on scanning is not an all-or-none phenomenon dependent simply on whether the size of the resulting discontinuity exceeds the step length of RNase E.

The Influence of the Nucleotide Composition of Genes and Gene Regulatory Elements on the Efficiency of Protein Expression in Escherichia coli

Recombinant proteins expressed in Escherichia coli are widely used in biochemical research and industrial processes. At the same time, achieving higher protein expression levels and correct protein folding still remains the key problem, since optimization of nutrient media, growth conditions, and methods for induction of protein synthesis do not always lead to the desired result. Often, low protein expression is determined by the sequences of the expressed genes and their regulatory regions. The genetic code is degenerated; 18 out of 20 amino acids are encoded by more than one codon. Choosing between synonymous codons in the coding sequence can significantly affect the level of protein expression and protein folding due to the influence of the gene nucleotide composition on the probability of formation of secondary mRNA structures that affect the ribosome binding at the translation initiation phase, as well as the ribosome movement along the mRNA during elongation, which, in turn, influences the mRNA degradation and the folding of the nascent protein. The nucleotide composition of the mRNA untranslated regions, in particular the promoter and Shine-Dalgarno sequences, also affects the efficiency of mRNA transcription, translation, and degradation. In this review, we describe the genetic principles that determine the efficiency of protein production in Escherichia coli.

Principles of gene regulation quantitatively connect DNA to RNA and proteins in bacteria

Protein concentrations are set by a complex interplay between gene-specific regulatory processes and systemic factors, including cell volume and shared gene expression machineries. Elucidating this interplay is crucial for discerning and designing gene regulatory systems. We quantitatively characterized gene-specific and systemic factors that affect transcription and translation genome-wide for Escherichia coli across many conditions. The results revealed two design principles that make regulation of gene expression insulated from concentrations of shared machineries: RNA polymerase activity is fine-tuned to match translational output, and translational characteristics are similar across most messenger RNAs (mRNAs). Consequently, in bacteria, protein concentration is set primarily at the promoter level. A simple mathematical formula relates promoter activities and protein concentrations across growth conditions, enabling quantitative inference of gene regulation from omics data.

The prokaryotic activity of the IGR IRESs is mediated by ribosomal protein S1

Internal ribosome entry sites (IRESs) are RNA elements capable of initiating translation on an internal portion of a messenger RNA. The intergenic region (IGR) IRES of the Dicistroviridae virus family folds into a triple pseudoknot tertiary structure, allowing it to recruit the ribosome and initiate translation in a structure dependent manner. This IRES has also been reported to drive translation in Escherichia coli and to date is the only described translation initiation signal that functions across domains of life. Here we show that unlike in the eukaryotic context the tertiary structure of the IGR IRES is not required for prokaryotic ribosome recruitment. In E. coli IGR IRES translation efficiency is dependent on ribosomal protein S1 in conjunction with an AG-rich Shine-Dalgarno-like element, supporting a model where the translational activity of the IGR IRESs is due to S1-mediated canonical prokaryotic translation.

Ribosome profiling enhances understanding of mycobacterial translation

A recent addition to the -omics toolkit, ribosome profiling, enables researchers to gain insight into the process and regulation of translation by mapping fragments of mRNA protected from nuclease digestion by ribosome binding. In this review, we discuss how ribosome profiling applied to mycobacteria has led to discoveries about translational regulation. Using case studies, we show that the traditional view of “canonical” translation mechanisms needs expanding to encompass features of mycobacterial translation that are more widespread than previously recognized. We also discuss the limitations of the method and potential future developments that could yield further insight into the fundamental biology of this important human pathogen.

Precise tuning of bacterial translation initiation by non-equilibrium 5'-UTR unfolding observed in single mRNAs

Noncoding, structured 5′-untranslated regions (5′-UTRs) of bacterial messenger RNAs (mRNAs) can control translation efficiency by forming structures that either recruit or repel the ribosome. Here we exploit a 5′-UTR embedded preQ₁-sensing, pseudoknotted translational riboswitch to probe how binding of a small ligand controls recruitment of the bacterial ribosome to the partially overlapping Shine-Dalgarno (SD) sequence. Combining single-molecule fluorescence microscopy with mutational analyses, we find that the stability of 30S ribosomal subunit binding is inversely correlated with the free energy needed to unfold the 5′-UTR during mRNA accommodation into the mRNA binding cleft. Ligand binding to the riboswitch stabilizes the structure to both antagonize 30S recruitment and accelerate 30S dissociation. Proximity of the 5′-UTR and stability of the SD:anti-SD interaction both play important roles in modulating the initial 30S-mRNA interaction. Finally, depletion of small ribosomal subunit protein S1, known to help resolve structured 5′-UTRs, further increases the energetic penalty for mRNA accommodation. The resulting model of rapid standby site exploration followed by gated non-equilibrium unfolding of the 5′-UTR during accommodation provides a mechanistic understanding of how translation efficiency is governed by riboswitches and other dynamic structure motifs embedded upstream of the translation initiation site of bacterial mRNAs.

Binding of 30S Ribosome Induces Single-stranded Conformation Within and Downstream of the Expression Platform in a Translational Riboswitch

Translational riboswitches are bacterial gene regulatory elements found in the 5'-untranslated region of mRNAs. They operate through a conformational refolding reaction that is triggered by a concentration change of a modulating small molecular ligand. The translation initiation region (TIR) is either released from or incorporated into base pairing interactions through the conformational switch. Hence, initiation of translation is regulated by the accessibility of the Shine-Dalgarno sequence and start codon. Interaction with the 30S ribosome is indispensable for the structural switch between functional OFF and ON states. However, on a molecular level it is still not fully resolved how the ribosome is accommodated near or at the translation initiation region in the context of translational riboswitches. The standby model of translation initiation postulates a binding site where the mRNA enters the ribosome and where it resides until the initiation site becomes unstructured and accessible. We here investigated the adenine-sensing riboswitch from Vibrio vulnificus. By application of a ¹⁹F labelling strategy for NMR spectroscopy that utilizes ligation techniques to synthesize differentially ¹⁹F labelled riboswitch molecules we show that nucleotides directly downstream of the riboswitch domain are first involved in productive interaction with the 30S ribosomal subunit. Upon the concerted action of ligand and the ribosomal protein rS1 the TIR becomes available and subsequently the 30S ribosome can slide towards the TIR. It will be interesting to see whether this is a general feature in translational riboswitches or if riboswitches exist where this region is structured and represent yet another layer of regulation.

Quantifying shifts in natural selection on codon usage between protein regions: a population genetics approach

Background

Codon usage bias (CUB), the non-uniform usage of synonymous codons, occurs across all domains of life. Adaptive CUB is hypothesized to result from various selective pressures, including selection for efficient ribosome elongation, accurate translation, mRNA secondary structure, and/or protein folding. Given the critical link between protein folding and protein function, numerous studies have analyzed the relationship between codon usage and protein structure. The results from these studies have often been contradictory, likely reflecting the differing methods used for measuring codon usage and the failure to appropriately control for confounding factors, such as differences in amino acid usage between protein structures and changes in the frequency of different structures with gene expression.

Results

Here we take an explicit population genetics approach to quantify codon-specific shifts in natural selection related to protein structure in S. cerevisiae and E. coli. Unlike other metrics of codon usage, our approach explicitly separates the effects of natural selection, scaled by gene expression, and mutation bias while naturally accounting for a region’s amino acid usage. Bayesian model comparisons suggest selection on codon usage varies only slightly between helix, sheet, and coil secondary structures and, similarly, between structured and intrinsically-disordered regions. Similarly, in contrast to prevous findings, we find selection on codon usage only varies slightly at the termini of helices in E. coli. Using simulated data, we show this previous work indicating “non-optimal” codons are enriched at the beginning of helices in S. cerevisiae was due to failure to control for various confounding factors (e.g. amino acid biases, gene expression, etc.), and rather than selection to modulate cotranslational folding.

Conclusions

Our results reveal a weak relationship between codon usage and protein structure, indicating that differences in selection on codon usage between structures are slight. In addition to the magnitude of differences in selection between protein structures being slight, the observed shifts appear to be idiosyncratic and largely codon-specific rather than systematic reversals in the nature of selection. Overall, our work demonstrates the statistical power and benefits of studying selective shifts on codon usage or other genomic features from an explicitly evolutionary approach. Limitations of this approach and future potential research avenues are discussed.

Supplementary Information

The online version contains supplementary material available at (10.1186/s12864-022-08635-0).

Identification of Promoter Region Markers Associated With Altered Expression of Resistance-Nodulation-Division Antibiotic Efflux Pumps in Acinetobacter baumannii

Genetic alterations leading to the constitutive upregulation of specific efflux pumps contribute to antibacterial resistance in multidrug resistant bacteria. The identification of such resistance markers remains one of the most challenging tasks of genome-level resistance predictors. In this study, 487 non-redundant genetic events were identified in upstream zones of three operons coding for resistance-nodulation-division (RND) efflux pumps of 4,130 Acinetobacter baumannii isolates. These events included insertion sequences, small indels, and single nucleotide polymorphisms. In some cases, alterations explicitly modified the expression motifs described for these operons, such as the promoter boxes, operators, and Shine-Dalgarno sequences. In addition, changes in DNA curvature and mRNA secondary structures, which are structural elements that regulate expression, were also calculated. According to their influence on RND upregulation, the catalog of upstream modifications were associated with "experimentally verified," "presumed," and "probably irrelevant" degrees of certainty. For experimental verification, DNA of upstream sequences independently carrying selected markers, three for each RND operon, were fused to a luciferase reporter plasmid system. Five out of the nine selected markers tested showed significant increases in expression with respect to the wild-type sequence control. In particular, a 25-fold expression increase was observed with the ISAba1 insertion sequence upstream the adeABC pump. Next, overexpression of each of the three multi-specific RND pumps was linked to their respective antibacterial substrates by a deep A. baumannii literature screen. Consequently, a data flow framework was then developed to link genomic upregulatory RND determinants to potential antibiotic resistance. Assignment of potential increases in minimal inhibitory concentrations at the "experimentally verified" level was permitted for 42 isolates to 7-8 unrelated antibacterial agents including tigecycline, which is overlooked by conventional resistome predictors. Thus, our protocol may represent a time-saving filter step prior to laborious confirmation experiments for efflux-driven resistance. Altogether, a computational-experimental pipeline containing all components required for identifying the upstream regulatory resistome is proposed. This schema may provide the foundational stone for the elaboration of tools approaching antibiotic efflux that complement routine resistome predictors for preventing antimicrobial therapy failure against difficult-to-threat bacteria.

Probing RNA structures and functions by solvent accessibility: an overview from experimental and computational perspectives

Characterizing RNA structures and functions have mostly been focused on 2D, secondary and 3D, tertiary structures. Recent advances in experimental and computational techniques for probing or predicting RNA solvent accessibility make this 1D representation of tertiary structures an increasingly attractive feature to explore. Here, we provide a survey of these recent developments, which indicate the emergence of solvent accessibility as a simple 1D property, adding to secondary and tertiary structures for investigating complex structure-function relations of RNAs.

Analysis of mRNA Decay Intermediates in Bacillus subtilis 3' Exoribonuclease and RNA Helicase Mutant Strains

The Bacillus subtilis genome encodes four 3' exoribonucleases: polynucleotide phosphorylase (PNPase), RNase R, RNase PH, and YhaM. Previous work showed that PNPase, encoded by the pnpA gene, is the major 3' exonuclease involved in mRNA turnover; in a pnpA deletion strain, numerous mRNA decay intermediates accumulate. Whether B. subtilis mRNA decay occurs in the context of a degradosome complex is controversial. In this study, global mapping of mRNA decay intermediate 3' ends within coding sequences was performed in strains that were either deleted for or had an inactivating point mutation in the pnpA gene. The patterns of 3'-end accumulation in these strains were highly similar, which may have implications for the role of a degradosome in mRNA decay. A comparison with mapped 3' ends in a strain lacking CshA, the major RNA helicase, indicated that many mRNAs require both PNPase and CshA for efficient decay. Transcriptome sequencing (RNA-seq) analysis of strains lacking RNase R suggested that this enzyme did not play a major role in mRNA turnover in the wild-type strain. Strains were constructed that contained only one of the four known 3' exoribonucleases. When RNase R was the only 3' exonuclease present, it was able to degrade a model mRNA efficiently, showing processive decay even through a strong stem-loop structure that inhibits PNPase processivity. Strains containing only RNase PH or only YhaM were also insensitive to this RNA secondary structure, suggesting the existence of another, as-yet-unidentified, 3' exoribonuclease. IMPORTANCE The ability to rapidly change bacterial gene expression programs in response to environmental conditions is highly dependent on the efficient turnover of mRNA. While much is known about the regulation of gene expression at the transcriptional and translational levels, much less is known about the intermediate step of mRNA decay. Here, we mapped the 3' ends of mRNA decay intermediates in strains that were missing the major 3' exoribonuclease PNPase or the RNA helicase CshA. We also assessed the roles of three other B. subtilis 3' exonucleases in the mRNA decay process. The data confirm the primary role of PNPase in mRNA turnover and suggest the involvement of one or more unknown RNases.

Pervasive translation in Mycobacterium tuberculosis

Most bacterial ORFs are identified by automated prediction algorithms. However, these algorithms often fail to identify ORFs lacking canonical features such as a length of >50 codons or the presence of an upstream Shine-Dalgarno sequence. Here, we use ribosome profiling approaches to identify actively translated ORFs in Mycobacterium tuberculosis. Most of the ORFs we identify have not been previously described, indicating that the M. tuberculosis transcriptome is pervasively translated. The newly described ORFs are predominantly short, with many encoding proteins of ≤50 amino acids. Codon usage of the newly discovered ORFs suggests that most have not been subject to purifying selection, and hence are unlikely to contribute to cell fitness. Nevertheless, we identify 90 new ORFs (median length of 52 codons) that bear the hallmarks of purifying selection. Thus, our data suggest that pervasive translation of short ORFs in Mycobacterium tuberculosis serves as a rich source for the evolution of new functional proteins.

Specific length and structure rather than high thermodynamic stability enable regulatory mRNA stem-loops to pause translation

Translating ribosomes unwind mRNA secondary structures by three basepairs each elongation cycle. Despite the ribosome helicase, certain mRNA stem-loops stimulate programmed ribosomal frameshift by inhibiting translation elongation. Here, using mutagenesis, biochemical and single-molecule experiments, we examine whether high stability of three basepairs, which are unwound by the translating ribosome, is critical for inducing ribosome pauses. We find that encountering frameshift-inducing mRNA stem-loops from the E. coli dnaX mRNA and the gag-pol transcript of Human Immunodeficiency Virus (HIV) hinders A-site tRNA binding and slows down ribosome translocation by 15-20 folds. By contrast, unwinding of first three basepairs adjacent to the mRNA entry channel slows down the translating ribosome by only 2-3 folds. Rather than high thermodynamic stability, specific length and structure enable regulatory mRNA stem-loops to stall translation by forming inhibitory interactions with the ribosome. Our data provide the basis for rationalizing transcriptome-wide studies of translation and searching for novel regulatory mRNA stem-loops.

DsrA Modulates Central Carbon Metabolism and Redox Balance by Directly Repressing pflB Expression in Salmonella Typhimurium

Bacterial small RNAs (sRNAs) function as vital regulators in response to various environmental stresses by base pairing with target mRNAs. The sRNA DsrA, an important posttranscriptional regulator, has been reported to play a crucial role in defense against oxidative stress in Salmonella enterica serovar Typhimurium, but its regulatory mechanism remains unclear. The transcriptome sequencing (RNA-seq) results in this study showed that the genes involved in glycolysis, pyruvate metabolism, the tricarboxylic acid (TCA) cycle, and NADH-dependent respiration exhibited significantly different expression patterns between S. Typhimurium wild type (WT) and the dsrA deletion mutant (ΔdsrA strain) before and after H₂O₂ treatment. This indicated the importance of DsrA in regulating central carbon metabolism (CCM) and NAD(H) homeostasis of S. Typhimurium. To reveal the direct target of DsrA action, fusion proteins of six candidate genes (acnA, srlE, tdcB, nuoH, katG, and pflB) with green fluorescent protein (GFP) were constructed, and the fluorescence analysis showed that the expression of pflB encoding pyruvate-formate lyase was repressed by DsrA. Furthermore, site-directed mutagenesis and RNase E-dependent experiments showed that the direct base pairing of DsrA with pflB mRNA could recruit RNase E to degrade pflB mRNA and reduce the stability of pflB mRNA. In addition, the NAD⁺/NADH ratio in WT-ppflB-pdsrA was significantly lower than that in WT-ppflB, suggesting that the repression of pflB by DsrA could contribute greatly to the redox balance in S. Typhimurium. Taken together, a novel target of DsrA was identified, and its regulatory role was clarified, which demonstrated that DsrA could modulate CCM and redox balance by directly repressing pflB expression in S. Typhimurium.

IMPORTANCE Small RNA DsrA plays an important role in defending against oxidative stress in bacteria. In this study, we identified a novel target (pflB, encoding pyruvate-formate lyase) of DsrA and demonstrated its potential regulatory mechanism in S. Typhimurium by transcriptome analysis. In silico prediction revealed a direct base pairing between DsrA and pflB mRNA, which was confirmed in site-directed mutagenesis experiments. The interaction of DsrA-pflB mRNA could greatly contribute to the regulation of central carbon metabolism and intracellular redox balance in S. Typhimurium. These findings provided a better understanding of the critical roles of small RNA in central metabolism and stress responses in foodborne pathogens.

mRNAs sequestered in stress granules recover nearly completely for translation

Stress granules (SGs) are membrane-less condensates composed of RNA and protein that assemble in response to stress stimuli and disassemble when stress is lifted. Both assembly and disassembly are tightly controlled processes, yet, it remains elusive whether mRNAs in SGs completely recover for translation following stress relief. Using RNA-seq of translating fractions in human cell line, we found that higher fraction of the m⁶A-modified mRNAs recovered for translation compared to unmodified mRNAs, i.e. 95% vs 84%, respectively. Considering structural mRNA analysis, we found that the m⁶A modification enhances structuring at nucleotides in its close vicinity. Our results suggest that SG-sequestered mRNAs disassemble nearly completely from SGs and the m⁶A modification may display some advantage to the mRNAs in their recovery for translation likely by m⁶A-driven structural stabilization.

Combinations of slow-translating codon clusters can increase mRNA half-life in Saccharomyces cerevisiae

The presence of a single cluster of nonoptimal codons was found to decrease a transcript's half-life through the interaction of the ribosome-associated quality control machinery with stalled ribosomes in Saccharomyces cerevisiae The impact of multiple nonoptimal codon clusters on a transcript's half-life, however, is unknown. Using a kinetic model, we predict that inserting a second nonoptimal cluster near the 5' end can lead to synergistic effects that increase a messenger RNA's (mRNA's) half-life in S. cerevisiae Specifically, the 5' end cluster suppresses the formation of ribosome queues, reducing the interaction of ribosome-associated quality control factors with stalled ribosomes. We experimentally validate this prediction by introducing two nonoptimal clusters into three different genes and find that their mRNA half-life increases up to fourfold. The model also predicts that in the presence of two clusters, the cluster closest to the 5' end is the primary determinant of mRNA half-life. These results suggest the "translational ramp," in which nonoptimal codons are located near the start codon and increase translational efficiency, may have the additional biological benefit of allowing downstream slow-codon clusters to be present without decreasing mRNA half-life. These results indicate that codon usage bias plays a more nuanced role in controlling cellular protein levels than previously thought.

Multifaceted impact of a nucleoside monophosphate kinase on 5'-end-dependent mRNA degradation in bacteria

A key pathway for mRNA degradation in bacterial cells begins with conversion of the initial 5'-terminal triphosphate to a monophosphate, a modification that renders transcripts more vulnerable to attack by ribonucleases whose affinity for monophosphorylated 5' ends potentiates their catalytic efficacy. In Escherichia coli, the only proteins known to be important for controlling degradation via this pathway are the RNA pyrophosphohydrolase RppH, its heteromeric partner DapF, and the 5'-monophosphate-assisted endonucleases RNase E and RNase G. We have now identified the metabolic enzyme cytidylate kinase as another protein that affects rates of 5'-end-dependent mRNA degradation in E. coli. It does so by utilizing two distinct mechanisms to influence the 5'-terminal phosphorylation state of RNA, each dependent on the catalytic activity of cytidylate kinase and not its mere presence in cells. First, this enzyme acts in conjunction with DapF to stimulate the conversion of 5' triphosphates to monophosphates by RppH. In addition, it suppresses the direct synthesis of monophosphorylated transcripts that begin with cytidine by reducing the cellular concentration of cytidine monophosphate, thereby disfavoring the 5'-terminal incorporation of this nucleotide by RNA polymerase during transcription initiation. Together, these findings suggest dual signaling pathways by which nucleotide metabolism can impact mRNA degradation in bacteria.

Impact of pseudouridylation, substrate fold, and degradosome organization on the endonuclease activity of RNase E

The conserved endoribonuclease RNase E dominates the dynamic landscape of RNA metabolism and underpins control mediated by small regulatory RNAs in diverse bacterial species. We explored the enzyme's hydrolytic mechanism, allosteric activation, and interplay with partner proteins in the multicomponent RNA degradosome assembly of Escherichia coli. RNase E cleaves single-stranded RNA with preference to attack the phosphate located at the 5' nucleotide preceding uracil, and we corroborate key interactions that select that base. Unexpectedly, RNase E activity is impeded strongly when the recognized uracil is isomerized to 5-ribosyluracil (pseudouridine), from which we infer the detailed geometry of the hydrolytic attack process. Kinetics analyses support models for recognition of secondary structure in substrates by RNase E and for allosteric autoregulation. The catalytic power of the enzyme is boosted when it is assembled into the multienzyme RNA degradosome, most likely as a consequence of substrate capture and presentation. Our results rationalize the origins of substrate preferences of RNase E and illuminate its catalytic mechanism, supporting the roles of allosteric domain closure and cooperation with other components of the RNA degradosome complex.

Analysis of 11,430 recombinant protein production experiments reveals that protein yield is tunable by synonymous codon changes of translation initiation sites

Recombinant protein production is a key process in generating proteins of interest in the pharmaceutical industry and biomedical research. However, about 50% of recombinant proteins fail to be expressed in a variety of host cells. Here we show that the accessibility of translation initiation sites modelled using the mRNA base-unpairing across the Boltzmann's ensemble significantly outperforms alternative features. This approach accurately predicts the successes or failures of expression experiments, which utilised Escherichia coli cells to express 11,430 recombinant proteins from over 189 diverse species. On this basis, we develop TIsigner that uses simulated annealing to modify up to the first nine codons of mRNAs with synonymous substitutions. We show that accessibility captures the key propensity beyond the target region (initiation sites in this case), as a modest number of synonymous changes is sufficient to tune the recombinant protein expression levels. We build a stochastic simulation model and show that higher accessibility leads to higher protein production and slower cell growth, supporting the idea of protein cost, where cell growth is constrained by protein circuits during overexpression.

Ribosome as a Translocase and Helicase

During protein synthesis, ribosome moves along mRNA to decode one codon after the other. Ribosome translocation is induced by a universally conserved protein, elongation factor G (EF-G) in bacteria and elongation factor 2 (EF-2) in eukaryotes. EF-G-induced translocation results in unwinding of the intramolecular secondary structures of mRNA by three base pairs at a time that renders the translating ribosome a processive helicase. Professor Alexander Sergeevich Spirin has made numerous seminal contributions to understanding the molecular mechanism of translocation. Here, we review Spirin's insights into the ribosomal translocation and recent advances in the field that stemmed from Spirin's pioneering work. We also discuss key remaining challenges in studies of translocase and helicase activities of the ribosome.

Quantitative Control for Stoichiometric Protein Synthesis

Bacterial protein synthesis rates have evolved to maintain preferred stoichiometries at striking precision, from the components of protein complexes to constituents of entire pathways. Setting relative protein production rates to be well within a factor of two requires concerted tuning of transcription, RNA turnover, and translation, allowing many potential regulatory strategies to achieve the preferred output. The last decade has seen a greatly expanded capacity for precise interrogation of each step of the central dogma genome-wide. Here, we summarize how these technologies have shaped the current understanding of diverse bacterial regulatory architectures underpinning stoichiometric protein synthesis. We focus on the emerging expanded view of bacterial operons, which encode diverse primary and secondary mRNA structures for tuning protein stoichiometry. Emphasis is placed on how quantitative tuning is achieved. We discuss the challenges and open questions in the application of quantitative, genome-wide methodologies to the problem of precise protein production. Expected final online publication date for the Annual Review of Microbiology, Volume 75 is October 2021. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Estimation of peptide elongation times from ribosome profiling spectra

Ribosome profiling spectra bear rich information on translation control and dynamics. Yet, due to technical biases in library generation, extracting quantitative measures of discrete translation events has remained elusive. Using maximum likelihood statistics and data set from Escherichia coli we develop a robust method for neutralizing technical biases (e.g. base specific RNase preferences in ribosome-protected mRNA fragments (RPF) generation), which allows for correct estimation of translation times at single codon resolution. Furthermore, we validated the method with available datasets from E. coli treated with antibiotic to inhibit isoleucyl-tRNA synthetase, and two datasets from Saccharomyces cerevisiae treated with two RNases with distinct cleavage signatures. We demonstrate that our approach accounts for RNase cleavage preferences and provides bias-corrected translation times estimates. Our approach provides a solution to the long-standing problem of extracting reliable information about peptide elongation times from highly noisy and technically biased ribosome profiling spectra.

smORFer: a modular algorithm to detect small ORFs in prokaryotes

Emerging evidence places small proteins (≤50 amino acids) more centrally in physiological processes. Yet, their functional identification and the systematic genome annotation of their cognate small open-reading frames (smORFs) remains challenging both experimentally and computationally. Ribosome profiling or Ribo-Seq (that is a deep sequencing of ribosome-protected fragments) enables detecting of actively translated open-reading frames (ORFs) and empirical annotation of coding sequences (CDSs) using the in-register translation pattern that is characteristic for genuinely translating ribosomes. Multiple identifiers of ORFs that use the 3-nt periodicity in Ribo-Seq data sets have been successful in eukaryotic smORF annotation. They have difficulties evaluating prokaryotic genomes due to the unique architecture (e.g. polycistronic messages, overlapping ORFs, leaderless translation, non-canonical initiation etc.). Here, we present a new algorithm, smORFer, which performs with high accuracy in prokaryotic organisms in detecting putative smORFs. The unique feature of smORFer is that it uses an integrated approach and considers structural features of the genetic sequence along with in-frame translation and uses Fourier transform to convert these parameters into a measurable score to faithfully select smORFs. The algorithm is executed in a modular way, and dependent on the data available for a particular organism, different modules can be selected for smORF search.

Towards the characterization of the hidden world of small proteins in Staphylococcus aureus, a proteogenomics approach

Small proteins play essential roles in bacterial physiology and virulence, however, automated algorithms for genome annotation are often not yet able to accurately predict the corresponding genes. The accuracy and reliability of genome annotations, particularly for small open reading frames (sORFs), can be significantly improved by integrating protein evidence from experimental approaches. Here we present a highly optimized and flexible bioinformatics workflow for bacterial proteogenomics covering all steps from (i) generation of protein databases, (ii) database searches and (iii) peptide-to-genome mapping to (iv) visualization of results. We used the workflow to identify high quality peptide spectrum matches (PSMs) for small proteins (≤ 100 aa, SP100) in Staphylococcus aureus Newman. Protein extracts from S. aureus were subjected to different experimental workflows for protein digestion and prefractionation and measured with highly sensitive mass spectrometers. In total, 175 proteins with up to 100 aa (SP100) were identified. Out of these 24 (ranging from 9 to 99 aa) were novel and not contained in the used genome annotation.144 SP100 are highly conserved and were found in at least 50% of the publicly available S. aureus genomes, while 127 are additionally conserved in other staphylococci. Almost half of the identified SP100 were basic, suggesting a role in binding to more acidic molecules such as nucleic acids or phospholipids.

Diversity and prevalence of ANTAR RNAs across actinobacteria

Background

Computational approaches are often used to predict regulatory RNAs in bacteria, but their success is limited to RNAs that are highly conserved across phyla, in sequence and structure. The ANTAR regulatory system consists of a family of RNAs (the ANTAR-target RNAs) that selectively recruit ANTAR proteins. This protein-RNA complex together regulates genes at the level of translation or transcriptional elongation. Despite the widespread distribution of ANTAR proteins in bacteria, their target RNAs haven’t been identified in certain bacterial phyla such as actinobacteria.

Results

Here, by using a computational search model that is tuned to actinobacterial genomes, we comprehensively identify ANTAR-target RNAs in actinobacteria. These RNA motifs lie in select transcripts, often overlapping with the ribosome binding site or start codon, to regulate translation. Transcripts harboring ANTAR-target RNAs majorly encode proteins involved in the transport and metabolism of cellular metabolites like sugars, amino acids and ions; or encode transcription factors that in turn regulate diverse genes.

Conclusion

In this report, we substantially diversify and expand the family of ANTAR RNAs across bacteria. These findings now provide a starting point to investigate the actinobacterial processes that are regulated by ANTAR.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12866-021-02234-x.

In Vivo and In Vitro Genome-Wide Profiling of RNA Secondary Structures Reveals Key Regulatory Features in Plasmodium falciparum

It is widely accepted that the structure of RNA plays important roles in a number of biological processes, such as polyadenylation, splicing, and catalytic functions. Dynamic changes in RNA structure are able to regulate the gene expression programme and can be used as a highly specific and subtle mechanism for governing cellular processes. However, the nature of most RNA secondary structures in Plasmodium falciparum has not been determined. To investigate the genome-wide RNA secondary structural features at single-nucleotide resolution in P. falciparum, we applied a novel high-throughput method utilizing the chemical modification of RNA structures to characterize these structures. Structural data from parasites are in close agreement with the known 18S ribosomal RNA secondary structures of P. falciparum and can help to predict the in vivo RNA secondary structure of a total of 3,396 transcripts in the ring-stage and trophozoite-stage developmental cycles. By parallel analysis of RNA structures in vivo and in vitro during the Plasmodium parasite ring-stage and trophozoite-stage intraerythrocytic developmental cycles, we identified some key regulatory features. Recent studies have established that the RNA structure is a ubiquitous and fundamental regulator of gene expression. Our study indicate that there is a critical connection between RNA secondary structure and mRNA abundance during the complex biological programme of P. falciparum. This work presents a useful framework and important results, which may facilitate further research investigating the interactions between RNA secondary structure and the complex biological programme in P. falciparum. The RNA secondary structure characterized in this study has potential applications and important implications regarding the identification of RNA structural elements, which are important for parasite infection and elucidating host-parasite interactions and parasites in the environment.

Mfd regulates RNA polymerase association with hard-to-transcribe regions in vivo, especially those with structured RNAs

RNA polymerase (RNAP) encounters various roadblocks during transcription. These obstacles can impede RNAP movement and influence transcription, ultimately necessitating the activity of RNAP-associated factors. One such factor is the bacterial protein Mfd, a highly conserved DNA translocase and evolvability factor that interacts with RNAP. Although Mfd is thought to function primarily in the repair of DNA lesions that stall RNAP, increasing evidence suggests that it may also be important for transcription regulation. However, this is yet to be fully characterized. To shed light on Mfd's in vivo functions, we identified the chromosomal regions where it associates. We analyzed Mfd's impact on RNAP association and transcription regulation genome-wide. We found that Mfd represses RNAP association at many chromosomal regions. We found that these regions show increased RNAP pausing, suggesting that they are hard to transcribe. Interestingly, we noticed that the majority of the regions where Mfd regulates transcription contain highly structured regulatory RNAs. The RNAs identified regulate a myriad of biological processes, ranging from metabolism to transfer RNA regulation to toxin-antitoxin (TA) functions. We found that cells lacking Mfd are highly sensitive to toxin overexpression. Finally, we found that Mfd promotes mutagenesis in at least one toxin gene, suggesting that its function in regulating transcription may promote evolution of certain TA systems and other regions containing strong RNA secondary structures. We conclude that Mfd is an RNAP cofactor that is important, and at times critical, for transcription regulation at hard-to-transcribe regions, especially those that express structured regulatory RNAs.

Non-optimal codon usage preferences of coronaviruses determine their promiscuity for infecting multiple hosts

Circulating animal coronaviruses occasionally infect humans. The SARS-CoV-2 is responsible for the current worldwide outbreak of COVID-19 that has resulted in 2 112 844 deaths as of late January 2021. We compared genetic code preferences in 496 viruses, including 34 coronaviruses and 242 corresponding hosts, to uncover patterns that distinguish single- and 'promiscuous' multiple-host-infecting viruses. Based on a codon usage preference score, promiscuous viruses were shown to significantly employ nonoptimal codons, namely codons that involve 'wobble' binding to anticodons, as compared to single-host viruses. The codon adaptation index (CAI) and the effective number of codons (ENC) were calculated for all viruses and hosts. Promiscuous viruses were less adapted hosts vs single-host viruses (P-value = 4.392e-11). All coronaviruses exploit nonoptimal codons to infect multiple hosts. We found that nonoptimal codon preferences at the beginning of viral coding sequences enhance the translational efficiency of viral proteins within the host. Finally, coronaviruses lack endogenous RNA degradation motifs to a significant degree, thereby increasing viral mRNA burden and infection load. To conclude, we found that promiscuously infecting coronaviruses prefer nonoptimal codon usage to remove degradation motifs from their RNAs and to dramatically increase their viral RNA production rates.

Systematic Quantification of Sequence and Structural Determinants Controlling mRNA stability in Bacterial Operons

mRNA degradation is a central process that affects all gene expression levels, and yet, the determinants that control mRNA decay rates remain poorly characterized. Here, we applied a synthetic biology, learn-by-design approach to elucidate the sequence and structural determinants that control mRNA stability in bacterial operons. We designed, constructed, and characterized 82 operons in Escherichia coli, systematically varying RNase binding site characteristics, translation initiation rates, and transcriptional terminator efficiencies in the 5' untranslated region (UTR), intergenic, and 3' UTR regions, followed by measuring their mRNA levels using reverse transcription quantitative polymerase chain reaction (RT-qPCR) assays during exponential growth. We show that introducing long single-stranded RNA into 5' UTRs reduced mRNA levels by up to 9.4-fold and that lowering translation rates reduced mRNA levels by up to 11.8-fold. We also found that RNase binding sites in intergenic regions had much lower effects on mRNA levels. Surprisingly, changing the transcriptional termination efficiency or introducing long single-stranded RNA into 3' UTRs had no effect on upstream mRNA levels. From these measurements, we developed and validated biophysical models of ribosome protection and RNase activity with excellent quantitative agreement. We also formulated design rules to rationally control a mRNA's stability, facilitating the automated design of engineered genetic systems with desired functionalities.

A snapshot of translation in Mycobacterium tuberculosis during exponential growth and nutrient starvation revealed by ribosome profiling

Mycobacterium tuberculosis, which causes tuberculosis, can undergo prolonged periods of non-replicating persistence in the host. The mechanisms underlying this are not fully understood, but translational regulation is thought to play a role. A large proportion of mRNA transcripts expressed in M. tuberculosis lack canonical bacterial translation initiation signals, but little is known about the implications of this for fine-tuning of translation. Here, we perform ribosome profiling to characterize the translational landscape of M. tuberculosis under conditions of exponential growth and nutrient starvation. Our data reveal robust, widespread translation of non-canonical transcripts and point toward different translation initiation mechanisms compared to canonical Shine-Dalgarno transcripts. During nutrient starvation, patterns of ribosome recruitment vary, suggesting that regulation of translation in this pathogen is more complex than originally thought. Our data represent a rich resource for others seeking to understand translational regulation in bacterial pathogens.

Fine-Tuning Multi-Gene Clusters via Well-Characterized Gene Expression Regulatory Elements: Case Study of the Arginine Synthesis Pathway in C. glutamicum

Promoters and ribosome binding sites (RBSs) are routinely applied in gene expression regulation, but their orthogonality and combinatorial effects have not yet been systematically studied in Corynebacterium glutamicum. Here, 17 core promoters and 29 RBSs in C. glutamicum were characterized, which exhibited 470-fold and 430-fold in transcriptional and translational activity, respectively. By comparing the expression of two reporter genes regulated by multiple RBSs, the RBS efficacy showed significant dependence on the gene context, besides the RBSs' strength, reflecting the poor orthogonality of RBSs. Bicistron-modified RBS (referred as bc-RBS) was adapted to C. glutamicum, which improved RBS reliability. By coupling a series of promoters with RBSs/bc-RBSs, a much broader regulation range that spanned 4 orders of magnitude was observed compared with that of a sole element, and the contribution to gene expression of RBS was more than that of promoter. Finally, promoters and RBSs were applied as built-in elements to fine-tune the gene cluster in the arginine synthesis pathway in C. glutamicum. Compared with the original strain, more arginine (1.61-fold) or citrulline (2.35-fold) was accumulated in a 7 L bioreactor by strains with the gene expression regulation system rationally engineered. We demonstrated that, via combination of well-characterized gene elements, and overall consideration for both transcription and translation, the biosynthesis pathway can be effectively balanced, and the yield of a target metabolite can be further improved.

Translational Control by Ribosome Pausing in Bacteria: How a Non-uniform Pace of Translation Affects Protein Production and Folding

Protein homeostasis of bacterial cells is maintained by coordinated processes of protein production, folding, and degradation. Translational efficiency of a given mRNA depends on how often the ribosomes initiate synthesis of a new polypeptide and how quickly they read the coding sequence to produce a full-length protein. The pace of ribosomes along the mRNA is not uniform: periods of rapid synthesis are separated by pauses. Here, we summarize recent evidence on how ribosome pausing affects translational efficiency and protein folding. We discuss the factors that slow down translation elongation and affect the quality of the newly synthesized protein. Ribosome pausing emerges as important factor contributing to the regulatory programs that ensure the quality of the proteome and integrate the cellular and environmental cues into regulatory circuits of the cell.

Characterizing Genetic Parts and Devices Using RNA Sequencing

Synthetic genetic circuits are composed of many parts that must interact and function together to produce a desired pattern of gene expression. A challenge when assembling circuits is that genetic parts often behave differently within a circuit, potentially impacting the desired functionality. Existing debugging methods based on fluorescent reporter proteins allow for only a few internal states to be monitored simultaneously, making diagnosis of the root cause impossible for large systems. Here, we present a tool called the Genetic Analyzer which uses RNA sequencing data to simultaneously characterize all transcriptional parts (e.g., promoters and terminators) and devices (e.g., sensors and logic gates) in complex genetic circuits. This provides a complete picture of the inner workings of a genetic circuit enabling faults to be easily identified and fixed. We construct a complete workflow to coordinate the execution of the various data processing and analysis steps and explain the options available when adapting these for the characterization of new systems.