Differential evolution in 3'UTRs leads to specific gene expression in Staphylococcus

RNAIII is linked with the pentose phosphate pathway through the activation of RpiRc in Staphylococcus aureus

Staphylococcus aureus RNAIII is a dual-function regulatory RNA that controls the expression of multiple virulence genes and especially the transition from adhesion to the production of exotoxins. However, its contribution to S. aureus central metabolism remains unclear. Using MS2-affinity purification coupled with RNA sequencing, we uncovered more than 50 novel RNAIII-mRNA interactions. Among them, we demonstrate that RNAIII is a major activator of the rpiRc gene, encoding a regulator of the pentose phosphate pathway (PPP). RNAIII binds the 5' UTR of rpiRc mRNA to favor ribosome loading, leading to an increased expression of RpiRc and, subsequently, of two PPP enzymes. Finally, we show that RNAIII and RpiRc are involved in S. aureus fitness in media supplemented with various carbohydrate sources related to PPP and glycolysis. Collectively, our data depict an unprecedented phenotype associated with the RNAIII regulon, especially the direct implication of RNAIII in central metabolic activity modulation. These findings show that the contribution of RNAIII in Staphylococcus aureus adaptation goes far beyond what was initially reported.

IMPORTANCE

Staphylococcus aureus is a major human pathogen involved in acute and chronic infections. Highly recalcitrant to antibiotic treatment, persistent infections are mostly associated with the loss of RNAIII expression, a master RNA regulator responsible for the switch from colonization to infection. Here, we used the MS2 affinity purification coupled with RNA sequencing approach to identify novel mRNA targets of RNAIII and uncover novel functions. We demonstrate that RNAIII is an activator of the expression of genes involved in the pentose phosphate pathway and is implicated in the adjustment of bacterial fitness as a function of carbohydrate sources. Taken together, our results demonstrate an unprecedented role of RNAIII that goes beyond the knowledge gained so far and contributes to a better understanding of the role of RNAIII in bacterial adaptation expression and the coordination of a complex regulatory network.

Leishmania Ribosomal Protein (RP) paralogous genes compensate each other's expression maintaining protein native levels

In the protozoan parasite Leishmania, most genes encoding for ribosomal proteins (RPs) are present as two or more copies in the genome. However, their untranslated regions (UTRs) are predominantly divergent and might be associated with a distinct regulation of the expression of paralogous genes. Herein, we investigated the expression profiles of two RPs (S16 and L13a) encoded by duplicated genes in Leishmania major. The genes encoding for the S16 protein possess identical coding sequences (CDSs) and divergent UTRs, whereas the CDSs of L13a diverge by two amino acids and by their UTRs. Using CRISPR/Cas9 genome editing, we generated knockout (Δ) and endogenously tagged transfectants for each paralog of L13a and S16 genes. Combining tagged and Δ cell lines we found evidence of differential expression of both RPS16 and RPL13a isoforms throughout parasite development, with one isoform consistently more abundant than its respective copy. In addition, compensatory expression was observed for each paralog upon deletion of the corresponding isoform, suggesting functional conservation between these proteins. This differential expression pattern relates to post-translational processes, given compensation occurs at the level of the protein, with no alterations detected at transcript level. Ribosomal profiles for RPL13a indicate a standard behavior for these paralogues suggestive of interaction with heavy RNA-protein complexes, as already reported for other RPs in trypanosomatids. We identified paralog-specific bound to their 3'UTRs which may be influential in regulating paralog expression. In support, we identified conserved cis-elements within the 3'UTRs of RPS16 and RPL13a; cis-elements exclusive to the UTR of the more abundant paralog or to the less abundant ones were identified.

The 3' UTR of vigR is required for virulence in Staphylococcus aureus and has expanded through STAR sequence repeat insertions

Infections caused by methicillin-resistant Staphylococcus aureus (MRSA) are alarmingly common, and treatment is confined to last-line antibiotics. Vancomycin is the treatment of choice for MRSA bacteremia, and treatment failure is often associated with vancomycin-intermediate S. aureus isolates. The regulatory 3' UTR of the vigR mRNA contributes to vancomycin tolerance and upregulates the autolysin IsaA. Using MS2-affinity purification coupled with RNA sequencing, we find that the vigR 3' UTR also regulates dapE, a succinyl-diaminopimelate desuccinylase required for lysine and peptidoglycan synthesis, suggesting a broader role in controlling cell wall metabolism and vancomycin tolerance. Deletion of the 3' UTR increased virulence, while the isaA mutant is completely attenuated in a wax moth larvae model. Sequence and structural analyses of vigR indicated that the 3' UTR has expanded through the acquisition of Staphylococcus aureus repeat insertions that contribute sequence for the isaA interaction seed and may functionalize the 3' UTR.

RlmQ: a newly discovered rRNA modification enzyme bridging RNA modification and virulence traits in Staphylococcus aureus

rRNA modifications play crucial roles in fine-tuning the delicate balance between translation speed and accuracy, yet the underlying mechanisms remain elusive. Comparative analyses of the rRNA modifications in taxonomically distant bacteria could help define their general, as well as species-specific, roles. In this study, we identified a new methyltransferase, RlmQ, in Staphylococcus aureus responsible for the Gram-positive specific m7G2601, which is not modified in Escherichia coli (G2574). We also demonstrate the absence of methylation on C1989, equivalent to E. coli C1962, which is methylated at position 5 by the Gram-negative specific RlmI methyltransferase, a paralog of RlmQ. Both modifications (S. aureus m7G2601 and E. coli m5C1962) are situated within the same tRNA accommodation corridor, hinting at a potential shared function in translation. Inactivation of S. aureus rlmQ causes the loss of methylation at G2601 and significantly impacts growth, cytotoxicity, and biofilm formation. These findings unravel the intricate connections between rRNA modifications, translation, and virulence in pathogenic Gram-positive bacteria.

TRS: a method for determining transcript termini from RNAtag-seq sequencing data

In bacteria, determination of the 3' termini of transcripts plays an essential role in regulation of gene expression, affecting the functionality and stability of the transcript. Several experimental approaches were developed to identify the 3' termini of transcripts, however, these were applied only to a limited number of bacteria and growth conditions. Here we present a straightforward approach to identify 3' termini from widely available RNA-seq data without the need for additional experiments. Our approach relies on the observation that the RNAtag-seq sequencing protocol results in overabundance of reads mapped to transcript 3' termini. We present TRS (Termini by Read Starts), a computational pipeline exploiting this property to identify 3' termini in RNAtag-seq data, and show that the identified 3' termini are highly reliable. Since RNAtag-seq data are widely available for many bacteria and growth conditions, our approach paves the way for studying bacterial transcription termination in an unprecedented scope.

Using a whole genome co-expression network to inform the functional characterisation of predicted genomic elements from Mycobacterium tuberculosis transcriptomic data

A whole genome co-expression network was created using Mycobacterium tuberculosis transcriptomic data from publicly available RNA-sequencing experiments covering a wide variety of experimental conditions. The network includes expressed regions with no formal annotation, including putative short RNAs and untranslated regions of expressed transcripts, along with the protein-coding genes. These unannotated expressed transcripts were among the best-connected members of the module sub-networks, making up more than half of the 'hub' elements in modules that include protein-coding genes known to be part of regulatory systems involved in stress response and host adaptation. This data set provides a valuable resource for investigating the role of non-coding RNA, and conserved hypothetical proteins, in transcriptomic remodelling. Based on their connections to genes with known functional groupings and correlations with replicated host conditions, predicted expressed transcripts can be screened as suitable candidates for further experimental validation.

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.

Translatome and transcriptome co-profiling reveals a role of TPRXs in human zygotic genome activation

Translational regulation plays a critical role during the oocyte-to-embryo transition (OET) and zygotic genome activation (ZGA). Here, we integrated ultra-low-input Ribo-seq with mRNA-seq to co-profile the translatome and transcriptome in human oocytes and early embryos. Comparison with mouse counterparts identified widespread differentially translated genes functioning in epigenetic reprogramming, transposon defense, and small RNA biogenesis, in part driven by species-specific regulatory elements in 3' untranslated regions. Moreover, PRD-like homeobox transcription factors, including TPRXL, TPRX1, and TPRX2, are highly translated around ZGA. TPRX1/2/L knockdown leads to defective ZGA and preimplantation development. Ectopically expressed TPRXs bind and activate key ZGA genes in human embryonic stem cells. These data reveal the conservation and divergence of translation landscapes during OET and identify critical regulators of human ZGA.

Staphylococcus aureus ftnA 3'-Untranslated Region Modulates Ferritin Production Facilitating Growth Under Iron Starvation Conditions

Iron acquisition and modulation of its intracellular concentration are critical for the development of all living organisms. So far, several proteins have been described to be involved in iron homeostasis. Among them, ferritins act as the major iron storage proteins, sequestering internalized iron and modulating its concentration inside bacterial cells. We previously described that the deletion of the 3’-untranslated region (3’UTR) of the ftnA gene, which codes for ferritin in Staphylococcus aureus, increased the ftnA mRNA and ferritin levels. Here, we show that the ferritin levels are affected by RNase III and PNPase, which target the ftnA 3’UTR. Rifampicin mRNA stability experiments revealed that the half-life of the ftnA mRNA is affected by both RNase III and the ftnA 3’UTR. A transcriptional fusion of the ftnA 3’UTR to the gfp reporter gene decreased green fluorescent protein (GFP) expression, indicating that the ftnA 3’UTR could work as an independent module. Additionally, a chromosomal deletion of the ftnA 3’UTR impaired S. aureus growth under conditions of iron starvation. Overall, this work highlights the biological relevance of the ftnA 3’UTR for iron homeostasis in S. aureus.

Bioinformatics Analysis of the Interaction of miRNAs and piRNAs with Human mRNA Genes Having di- and Trinucleotide Repeats

The variability of nucleotide repeats is considered one of the causes of diseases, but their biological function is not understood. In recent years, the interaction of miRNAs and piRNAs with the mRNAs of genes responsible for developing neurodegenerative and oncological diseases and diabetes have been actively studied. We explored candidate genes with nucleotide repeats to predict associations with miRNAs and piRNAs. The parameters of miRNAs and piRNA binding sites with mRNAs of human genes having nucleotide repeats were determined using the MirTarget program. This program defines the start of the initiation of miRNA and piRNA binding to mRNAs, the localization of miRNA and piRNA binding sites in the 5′-untranslated region (5′UTR), coding sequence (CDS) and 3′-untranslated region (3′UTR); the free energy of binding; and the schemes of nucleotide interactions of miRNAs and piRNAs with mRNAs. The characteristics of miRNAs and piRNA binding sites with mRNAs of 73 human genes were determined. The 5′UTR, 3′UTR and CDS of the mRNAs of genes are involved in the development of neurodegenerative, oncological and diabetes diseases with GU, AC dinucleotide and CCG, CAG, GCC, CGG, CGC trinucleotide repeats. The associations of miRNAs, piRNAs and candidate target genes could be recommended for developing methods for diagnosing diseases, including neurodegenerative diseases, oncological diseases and diabetes.

Synthetic 3'-UTR valves for optimal metabolic flux control in Escherichia coli

As the design of genetic circuitry for synthetic biology becomes more sophisticated, diverse regulatory bioparts are required. Despite their importance, well-characterized 3′-untranslated region (3′-UTR) bioparts are limited. Thus, transcript 3′-ends require further investigation to understand the underlying regulatory role and applications of the 3′-UTR. Here, we revisited the use of Term-Seq in the Escherichia coli strain K-12 MG1655 to enhance our understanding of 3′-UTR regulatory functions and to provide a diverse collection of tunable 3′-UTR bioparts with a wide termination strength range. Comprehensive analysis of 1,629 transcript 3′-end positions revealed multiple 3′-termini classes generated through transcription termination and RNA processing. The examination of individual Rho-independent terminators revealed a reduction in downstream gene expression over a wide range, which led to the design of novel synthetic metabolic valves that control metabolic fluxes in branched pathways. These synthetic metabolic valves determine the optimal balance of heterologous pathways for maximum target biochemical productivity. The regulatory strategy using 3′-UTR bioparts is advantageous over promoter- or 5′-UTR-based transcriptional control as it modulates gene expression at transcription levels without trans-acting element requirements (e.g. transcription factors). Our results provide a foundational platform for 3′-UTR engineering in synthetic biology applications.

Regulation of Heterogenous LexA Expression in Staphylococcus aureus by an Antisense RNA Originating from Transcriptional Read-Through upon Natural Mispairings in the sbrB Intrinsic Terminator

Bacterial genomes are pervasively transcribed, generating a wide variety of antisense RNAs (asRNAs). Many of them originate from transcriptional read-through events (TREs) during the transcription termination process. Previous transcriptome analyses revealed that the lexA gene from Staphylococcus aureus, which encodes the main SOS response regulator, is affected by the presence of an asRNA. Here, we show that the lexA antisense RNA (lexA-asRNA) is generated by a TRE on the intrinsic terminator (TT_sbrB) of the sbrB gene, which is located downstream of lexA, in the opposite strand. Transcriptional read-through occurs by a natural mutation that destabilizes the TT_sbrB structure and modifies the efficiency of the intrinsic terminator. Restoring the mispairing mutation in the hairpin of TT_sbrB prevented lexA-asRNA transcription. The level of lexA-asRNA directly correlated with cellular stress since the expressions of sbrB and lexA-asRNA depend on the stress transcription factor SigB. Comparative analyses revealed strain-specific nucleotide polymorphisms within TT_sbrB, suggesting that this TT could be prone to accumulating natural mutations. A genome-wide analysis of TREs suggested that mispairings in TT hairpins might provide wider transcriptional connections with downstream genes and, ultimately, transcriptomic variability among S. aureus strains.

The 3'UTR-derived sRNA RsaG coordinates redox homeostasis and metabolism adaptation in response to glucose-6-phosphate uptake in Staphylococcus aureus

Staphylococcus aureus RsaG is a 3'-untranslated region (3'UTR) derived sRNA from the conserved uhpT gene encoding a glucose-6-phosphate (G6P) transporter expressed in response to extracellular G6P. The transcript uhpT-RsaG undergoes degradation from 5'- to 3'-end by the action of the exoribonucleases J1/J2, which are blocked by a stable hairpin structure at the 5'-end of RsaG, leading to its accumulation. RsaG together with uhpT is induced when bacteria are internalized into host cells or in the presence of mucus-secreting cells. Using MS2-affinity purification coupled with RNA sequencing, several RNAs were identified as targets including mRNAs encoding the transcriptional factors Rex, CcpA, SarA, and the sRNA RsaI. Our data suggested that RsaG contributes to the control of redox homeostasis and adjusts metabolism to changing environmental conditions. RsaG uses different molecular mechanisms to stabilize, degrade, or repress the translation of its mRNA targets. Although RsaG is conserved only in closely related species, the uhpT 3'UTR of the ape pathogen S. simiae harbors an sRNA, whose sequence is highly different, and which does not respond to G6P levels. Our results hypothesized that the 3'UTRs from UhpT transporter encoding mRNAs could have rapidly evolved to enable adaptation to host niches.

Direct RNA Sequencing Unfolds the Complex Transcriptome of Vibrio parahaemolyticus

Conventional bacterial genome annotation provides information about coding sequences but ignores untranslated regions and operons. However, untranslated regions contain important regulatory elements as well as targets for many regulatory factors, such as small RNAs. Operon maps are also essential for functional gene analysis. In the last decade, considerable progress has been made in the study of bacterial transcriptomes through transcriptome sequencing (RNA-seq). Given the compact nature of bacterial genomes, many challenges still cannot be resolved through short reads generated using classical RNA-seq because of fragmentation and loss of the full-length information. Direct RNA sequencing is a technology that sequences the native RNA directly without information loss or bias. Here, we employed direct RNA sequencing to annotate the Vibrio parahaemolyticus transcriptome with its full features, including transcription start sites (TSSs), transcription termination sites, and operon maps. A total of 4,103 TSSs were identified. In comparison to short-read sequencing, full-length information provided a deeper view of TSS classification, showing that most internal and antisense TSSs were actually a result of gene overlap. Sequencing the transcriptome of V. parahaemolyticus grown with bile allowed us to study the landscape of pathogenicity island Vp-PAI. Some genes in this region were reannotated, providing more accurate annotation to increase precision in their characterization. Quantitative detection of operons in V. parahaemolyticus showed high complexity in some operons, shedding light on a greater extent of regulation within the same operon. Our study using direct RNA sequencing provides a quantitative and high-resolution landscape of the V. parahaemolyticus transcriptome.

IMPORTANCEVibrio parahaemolyticus is a halophilic bacterium found in the marine environment. Outbreaks of gastroenteritis resulting from seafood poisoning by these pathogens have risen over the past 2 decades. Upon ingestion by humans—often through the consumption of raw or undercooked seafood—V. parahaemolyticus senses the host environment and expresses numerous genes, the products of which synergize to synthesize and secrete toxins that can cause acute gastroenteritis. To understand the regulation of such adaptive response, mRNA transcripts must be mapped accurately. However, due to the limitations of common sequencing methods, not all features of bacterial transcriptomes are always reported. We applied direct RNA sequencing to analyze the V. parahaemolyticus transcriptome. Mapping the full features of the transcriptome is anticipated to enhance our understanding of gene regulation in this bacterium and provides a data set for future work. Additionally, this study reveals a deeper view of a complicated transcriptome landscape, demonstrating the importance of applying such methods to other bacterial models.

Siblings or doppelgängers? Deciphering the evolution of structured cis-regulatory RNAs beyond homology

Structured cis-regulatory RNAs have evolved across all domains of life, highlighting the utility and plasticity of RNA as a regulatory molecule. Homologous RNA sequences and structures often have similar functions, but homology may also be deceiving. The challenges that derive from trying to assign function to structure and vice versa are not trivial. Bacterial riboswitches, viral and eukaryotic IRESes, CITEs, and 3′ UTR elements employ an array of mechanisms to exert their effects. Bioinformatic searches coupled with biochemical and functional validation have elucidated some shared and many unique ways cis-regulators are employed in mRNA transcripts. As cis-regulatory RNAs are resolved in greater detail, it is increasingly apparent that shared homology can mask the full spectrum of mRNA cis-regulator functional diversity. Furthermore, similar functions may be obscured by lack of obvious sequence similarity. Thus looking beyond homology is crucial for furthering our understanding of RNA-based regulation.

Bacterial 3'UTRs: A Useful Resource in Post-transcriptional Regulation

Bacterial messenger RNAs (mRNAs) are composed of 5′ and 3′ untranslated regions (UTRs) that flank the coding sequences (CDSs). In eukaryotes, 3′UTRs play key roles in post-transcriptional regulatory mechanisms. Shortening or deregulation of these regions is associated with diseases such as cancer and metabolic disorders. Comparatively, little is known about the functions of 3′UTRs in bacteria. Over the past few years, 3′UTRs have emerged as important players in the regulation of relevant bacterial processes such as virulence, iron metabolism, and biofilm formation. This MiniReview is an update for the different 3′UTR-mediated mechanisms that regulate gene expression in bacteria. Some of these include 3′UTRs that interact with the 5′UTR of the same transcript to modulate translation, 3′UTRs that are targeted by specific ribonucleases, RNA-binding proteins and small RNAs (sRNAs), and 3′UTRs that act as reservoirs of trans-acting sRNAs, among others. In addition, recent findings regarding a differential evolution of bacterial 3′UTRs and its impact in the species-specific expression of orthologous genes are also discussed.

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.