Emergence of New sRNAs in Enteric Bacteria is Associated with Low Expression and Rapid Evolution

A Peroxide-Responding sRNA Evolved from a Peroxidase mRNA

Small RNAs (sRNAs) are important gene regulators in bacteria, but it is unclear how new sRNAs originate and become part of regulatory networks that coordinate bacterial response to environmental stimuli. Using a covariance modeling-based approach, we analyzed the presence of hundreds of sRNAs in more than a thousand genomes across Enterobacterales, a bacterial order with a confluence of factors that allows robust genome-scale sRNA analyses: several well-studied organisms with fairly conserved genome structures, an established phylogeny, and substantial nucleotide diversity within a narrow evolutionary space. We discovered that a majority of sRNAs arose recently, and uncovered protein-coding genes as a potential source from which new sRNAs arise. A detailed investigation of the emergence of OxyS, a peroxide-responding sRNA, revealed that it evolved from a fragment of a peroxidase messenger RNA. Importantly, although it replaced the ancestral peroxidase, OxyS continues to be part of the ancestral peroxide-response regulon, indicating that an sRNA that arises from a protein-coding gene would inherently be part of the parental protein's regulatory network. This new insight provides a fresh framework for understanding sRNA origin and regulatory integration in bacteria.

Small RNAs that target G-rich sequences are generated by diverse biogenesis pathways in Epsilonproteobacteria

Bacterial small RNAs (sRNAs) are widespread post-transcriptional regulators controlling bacterial stress responses and virulence. Nevertheless, little is known about how they arise and evolve. Homologues can be difficult to identify beyond the strain level using sequence-based approaches, and similar functionalities can arise by convergent evolution. Here, we found that the virulence-associated CJnc190 sRNA of the foodborne pathogen Campylobacter jejuni resembles the RepG sRNA from the gastric pathogen Helicobacter pylori. However, while both sRNAs bind G-rich sites in their target mRNAs using a C/U-rich loop, they largely differ in their biogenesis. RepG is transcribed from a stand-alone gene and does not require processing, whereas CJnc190 is transcribed from two promoters as precursors that are processed by RNase III and also has a cis-encoded antagonist, CJnc180. By comparing CJnc190 homologues in diverse Campylobacter species, we show that RNase III-dependent processing of CJnc190 appears to be a conserved feature even outside of C. jejuni. We also demonstrate the CJnc180 antisense partner is expressed in C. coli, yet here might be derived from the 3'UTR of the upstream flagella-related gene. Our analysis of G-tract targeting sRNAs in Epsilonproteobacteria demonstrates that similar sRNAs can have markedly different biogenesis pathways.

Burning the Candle at Both Ends: Have Exoribonucleases Driven Divergence of Regulatory RNA Mechanisms in Bacteria?

Regulatory RNAs have emerged as ubiquitous gene regulators in all bacterial species studied to date. The combination of sequence-specific RNA interactions and malleable RNA structure has allowed regulatory RNA to adopt different mechanisms of gene regulation in a diversity of genetic backgrounds. In the model GammaproteobacteriaEscherichia coli and Salmonella, the regulatory RNA chaperone Hfq appears to play a global role in gene regulation, directly controlling ∼20 to 25% of the entire transcriptome. While the model FirmicutesBacillus subtilis and Staphylococcus aureus encode a Hfq homologue, its role has been significantly depreciated. These bacteria also have marked differences in RNA turnover. E. coli and Salmonella degrade RNA through internal endonucleolytic and 3′→5′ exonucleolytic cleavage that appears to allow transient accumulation of mRNA 3′ UTR cleavage fragments that contain stabilizing 3′ structures. In contrast, B. subtilis and S. aureus are able to exonucleolytically attack internally cleaved RNA from both the 5′ and 3′ ends, efficiently degrading mRNA 3′ UTR fragments. Here, we propose that the lack of 5′→3′ exoribonuclease activity in Gammaproteobacteria has allowed the accumulation of mRNA 3′ UTR ends as the “default” setting. This in turn may have provided a larger pool of unconstrained RNA sequences that has fueled the expansion of Hfq function and small RNA (sRNA) regulation in E. coli and Salmonella. Conversely, the exoribonuclease RNase J may be a significant barrier to the evolution of 3′ UTR sRNAs in B. subtilis and S. aureus that has limited the pool of RNA ligands available to Hfq and other sRNA chaperones, depreciating their function in these model Firmicutes.

Natural selection shapes maintenance of orthologous sRNAs in divergent host-restricted bacterial genomes

Historically it has been difficult to study the evolution of bacterial small RNAs (sRNAs) across distantly related species. For example, identifying homologs of sRNAs is often difficult in genomes that have undergone multiple structural rearrangements. Also, some types of regulatory sRNAs evolve at rapid rates. The high degree of genomic synteny among divergent host-restricted bacterial lineages, including intracellular symbionts, is conducive to sRNA maintenance and homolog identification. In turn, symbiont genomes can provide us with novel insights into sRNA evolution. Here, we examine the sRNA expression profile of the obligate symbiont of psyllids, Carsonella ruddii, which has one of the smallest cellular genomes described. Using RNA-seq, we identified 36 and 32 antisense sRNAs (asRNAs) expressed by Carsonella from the psyllids Bactericera cockerelli (Carsonella-BC) and Diaphorina citri (Carsonella-DC), respectively. The majority of these asRNAs were associated with genes that are involved in essential amino acid biosynthetic pathways. Eleven of the asRNAs were conserved in both Carsonella lineages and the majority were maintained by selection. Notably, five of the corresponding coding sequences are also the targets of conserved asRNAs in a distantly related insect symbiont, Buchnera. We detected differential expression of two asRNAs for genes involved in arginine and leucine biosynthesis occurring between two distinct Carsonella-BC life stages. Using asRNAs identified in Carsonella, Buchnera, and Profftella which are all endosymbionts, and Escherichia coli, we determined that regions upstream of these asRNAs encode unique conserved patterns of AT/GC richness, GC skew, and sequence motifs which may be involved in asRNA regulation.

Modulation of Bacterial Fitness and Virulence Through Antisense RNAs

Regulatory RNAs contribute to gene expression control in bacteria. Antisense RNAs (asRNA) are a class of regulatory RNAs that are transcribed from opposite strands of their target genes. Typically, these untranslated transcripts bind to cognate mRNAs and rapidly regulate gene expression at the post-transcriptional level. In this article, we review asRNAs that modulate bacterial fitness and increase virulence. We chose examples that underscore the variety observed in nature including, plasmid- and chromosome-encoded asRNAs, a riboswitch-regulated asRNA, and asRNAs that require other RNAs or RNA-binding proteins for stability and activity. We explore how asRNAs improve bacterial fitness and virulence by modulating plasmid acquisition and maintenance, regulating transposon mobility, increasing resistance against bacteriophages, controlling flagellar production, and regulating nutrient acquisition. We conclude with a brief discussion on how this knowledge is helping to inform current efforts to develop new therapeutics.

Divergently Transcribed ncRNAs in Escherichia coli: Refinement of the Transcription Starts Assumes Functional Diversification

Non-coding regulatory RNAs (ncRNAs) comprise specialized group of essential genetically encoded biological molecules involved in the wide variety of cellular metabolic processes. The progressive increase in the number of newly identified ncRNAs and the defining of their genome location indicate their predominant nesting in intergenic regions and expression under the control of their own regulatory elements. At the same time, the regulation of ncRNA's transcription cannot be considered in isolation from the processes occurring in the immediate genetic environment. A number of experimental data indicate the notable impact of positional regulation of gene expression mediated by dynamic temporal DNA rearrangements accompanying transcription events in the vicinity of neighboring genes. This issue can be perceived as particularly significant for divergently transcribed ncRNAs being actually subjected to double regulatory pressure. Based on available results of RNAseq experiments for Escherichia coli, we screened out divergent ncRNAs and the adjacent genes for the exact positions of transcription start sites (TSSs) and relative efficiency of RNA production. This analysis revealed extension or shortening of some previously annotated ncRNAs resulting in modified secondary structure, confirmed stable expression of four ncRNAs annotated earlier as putative, and approved the possibility of expression of divergently transcribed ncRNAs containing repetitive extragenic palindromic (REP) elements. The biogenesis of secreted ncRNAs from divergently transcribed ffs, chiX, ralA, and ryhB is discussed taking into account positions of TSSs. Refinement of TSSs for the neighboring genes renders some ncRNAs as true antisense overlapping with 5'UTR of divergently transcribed mRNAs.

Are Antisense Proteins in Prokaryotes Functional?

Many prokaryotic RNAs are transcribed from loci outside of annotated protein coding genes. Across bacterial species hundreds of short open reading frames antisense to annotated genes show evidence of both transcription and translation, for instance in ribosome profiling data. Determining the functional fraction of these protein products awaits further research, including insights from studies of molecular interactions and detailed evolutionary analysis. There are multiple lines of evidence, however, that many of these newly discovered proteins are of use to the organism. Condition-specific phenotypes have been characterized for a few. These proteins should be added to genome annotations, and the methods for predicting them standardized. Evolutionary analysis of these typically young sequences also may provide important insights into gene evolution. This research should be prioritized for its exciting potential to uncover large numbers of novel proteins with extremely diverse potential practical uses, including applications in synthetic biology and responding to pathogens.

Bacteriophages as sources of small non-coding RNA molecules

Bacteriophages play an essential role in the transferring of genes that contribute to the bacterial virulence and whose products are dangerous to human health. Interestingly, phages carrying virulence genes are mostly temperate and in contrast to lytic phages undergo both lysogenic and lytic cycles. Importantly, expression of the majority of phage genes and subsequent production of phage encoded proteins is suppressed during lysogeny. The expression of the majority of phage genes is tightly linked to lytic development. Among others, small non-coding RNAs (sRNAs) of phage origin are involved in the regulation of phage gene expression and thus play an important role in both phage and host development. In the case of bacteria, sRNAs affect processes such as virulence, colonization ability, motility and cell growth or death. In turn, in the case of phages, they play essential roles during the early stage of infection, maintaining the state of lysogeny and silencing the expression of late structural genes, thereby regulating the transition between phage life cycles. Interestingly, sRNAs have been identified in both lytic and temperate phages and they have been discussed in this work according to this classification. Particular attention was paid to viral sRNAs resembling eukaryotic microRNAs.

A CsrA-Binding, trans-Acting sRNA of Coxiella burnetii Is Necessary for Optimal Intracellular Growth and Vacuole Formation during Early Infection of Host Cells

Coxiella burnetii is an obligate intracellular gammaproteobacterium and zoonotic agent of Q fever. We previously identified 15 small noncoding RNAs (sRNAs) of C. burnetii One of them, CbsR12 (Coxiella burnetii small RNA 12), is highly transcribed during axenic growth and becomes more prominent during infection of cultured mammalian cells. Secondary structure predictions of CbsR12 revealed four putative CsrA-binding sites in stem loops with consensus AGGA/ANGGA motifs. We subsequently determined that CbsR12 binds to recombinant C. burnetii CsrA-2, but not CsrA-1, proteins in vitro Moreover, through a combination of in vitro and cell culture assays, we identified several in trans mRNA targets of CbsR12. Of these, we determined that CbsR12 binds and upregulates translation of carA transcripts coding for carbamoyl phosphate synthetase A, an enzyme that catalyzes the first step of pyrimidine biosynthesis. In addition, CbsR12 binds and downregulates translation of metK transcripts coding for S-adenosylmethionine synthetase, a component of the methionine cycle. Furthermore, we found that CbsR12 binds to and downregulates the quantity of cvpD transcripts, coding for a type IVB effector protein, in mammalian cell culture. Finally, we found that CbsR12 is necessary for expansion of Coxiella-containing vacuoles and affects growth rates in a dose-dependent manner in the early phase of infecting THP-1 cells. This is the first characterization of a trans-acting sRNA of C. burnetii and the first example of a bacterial sRNA that regulates both CarA and MetK synthesis. CbsR12 is one of only a few identified trans-acting sRNAs that interacts with CsrA.IMPORTANCE Regulation of metabolism and virulence in C. burnetii is not well understood. Here, we show that C. burnetii small RNA 12 (CbsR12) is highly transcribed in the metabolically active large-cell variant compared to the nonreplicative small-cell variant. We show that CbsR12 directly regulates several genes involved in metabolism, along with a type IV effector gene, in trans In addition, we demonstrate that CbsR12 binds to CsrA-2 in vitro and induces autoaggregation and biofilm formation when transcribed ectopically in Escherichia coli, consistent with other CsrA-sequestering sRNAs. These results implicate CbsR12 in the indirect regulation of a number of genes via CsrA-mediated regulatory activities. The results also support CbsR12 as a crucial regulatory component early on in a mammalian cell infection.

EsR240, a non-coding sRNA, is required for the resistance of Edwardsiella tarda to stresses in macrophages and for virulence

Bacterial small non-coding RNAs (sRNAs) are gene expression modulators that respond to environmental changes and pathogenic conditions. In this study, 13 novel sRNAs were identified in the intracellular pathogen, Edwardsiella tarda (E. tarda) ET13 strain, based on RNA sequencing and bioinformatic analyses. Eight of the 13 putative sRNAs from the ET13 strain were transcribed (as indicated by RT-PCR) following exposure to different stresses. The transcription levels of three sRNAs (EsR128, EsR139 and EsR240) were all highly induced under these stress conditions. Northern blot hybridization was employed to verify that EsR240 was expressed in the ET13 strain under both logarithmic and stationary growth phases, and that it formed a single copy transcript in the chromosomes of the ET13 strain. The precise start and end points of EsR240 were determined using 5'and 3' RACE. The conservation of EsR240 was in agreement with the characteristics of sRNA, as indicated by a BLAST analysis. Furthermore, the survival rates of EsR240 mutant were lower than the rates of the wild type ET13 under stress conditions. When the infection time was extended 4 or 6 h, the CFUs of the wild type bacteria increased more significantly within macrophages compared to the mutant. When the intra-peritoneal (i.p.) route of infection was used in mice, the bacterial loads of the tissues in the mice infected with the wild type bacteria were significantly higher than in the mice infected with the mutants. The virulence of the EsR240 mutant was 6.79-fold lower than the wild type bacterium based on the LD₅₀. In addition, the IntaRNA program was used to predict the target genes of EsR240. Out of the top 10 predicted target genes, 9 genes were regulated by EsR240. These target genes may encode FtsH protease modulator YccA, Na⁺ and H⁺ antiporters, FtsX-like permease family protein, glycoside hydrolases or various other proteins. Therefore, EsR240 may positively regulate its target genes in E. tarda to maintain intracellular survival within host macrophages and to increase its virulence.

Transcriptional noise and exaptation as sources for bacterial sRNAs

Understanding how new genes originate and integrate into cellular networks is key to understanding evolution. Bacteria present unique opportunities for both the natural history and experimental study of gene origins, due to their large effective population sizes, rapid generation times, and ease of genetic manipulation. Bacterial small non-coding RNAs (sRNAs), in particular, many of which operate through a simple antisense regulatory logic, may serve as tractable models for exploring processes of gene origin and adaptation. Understanding how and on what timescales these regulatory molecules arise has important implications for understanding the evolution of bacterial regulatory networks, in particular, for the design of comparative studies of sRNA function. Here, we introduce relevant concepts from evolutionary biology and review recent work that has begun to shed light on the timescales and processes through which non-functional transcriptional noise is co-opted to provide regulatory functions. We explore possible scenarios for sRNA origin, focusing on the co-option, or exaptation, of existing genomic structures which may provide protected spaces for sRNA evolution.

Cross-Regulation between Bacteria and Phages at a Posttranscriptional Level

The study of bacteriophages (phages) and prophages has provided key insights into almost every cellular process as well as led to the discovery of unexpected new mechanisms and the development of valuable tools. This is exemplified for RNA-based regulation. For instance, the characterization and exploitation of the antiphage CRISPR (clustered regularly interspaced short palindromic repeat) systems is revolutionizing molecular biology. Phage-encoded proteins such as the RNA-binding MS2 protein, which is broadly used to isolate tagged RNAs, also have been developed as valuable tools. Hfq, the RNA chaperone protein central to the function of many base-pairing small RNAs (sRNAs), was first characterized as a bacterial host factor required for Qβ phage replication. The ongoing studies of RNAs are continuing to reveal regulatory connections between infecting phages, prophages, and bacteria and to provide novel insights. There are bacterial and prophage sRNAs that regulate prophage genes, which impact bacterial virulence as well as bacterial cell killing. Conversely, phage- and prophage-encoded sRNAs modulate the expression of bacterial genes modifying metabolism. An interesting subcategory of the prophage-encoded sRNAs are sponge RNAs that inhibit the activities of bacterial-encoded sRNAs. Phages also affect posttranscriptional regulation in bacteria through proteins that inhibit or alter the activities of key bacterial proteins involved in posttranscriptional regulation. However, what is most exciting about phage and prophage research, given the millions of phage-encoded genes that have not yet been characterized, is the vast potential for discovering new RNA regulators and novel mechanisms and for gaining insight into the evolution of regulatory RNAs.

Origin, Evolution, and Loss of Bacterial Small RNAs

Despite the central role of bacterial noncoding small RNAs (sRNAs) in posttranscriptional regulation, little is understood about their evolution. Here we compile what has been studied to date and trace a life cycle of sRNAs-from their mechanisms of emergence, through processes of change and frequent neofunctionalization, to their loss from bacterial lineages. Because they possess relatively unrestrictive structural requirements, we find that sRNA origins are varied, and include de novo emergence as well as formation from preexisting genetic elements via duplication events and horizontal gene transfer. The need for only partial complementarity to their mRNA targets facilitates apparent rapid change, which also contributes to significant challenges in tracing sRNAs across broad evolutionary distances. We document that recently emerged sRNAs in particular evolve quickly, mirroring dynamics observed in microRNAs, their functional analogs in eukaryotes. Mutations in mRNA-binding regions, transcriptional regulator or sigma factor binding sites, and protein-binding regions are all likely sources of shifting regulatory roles of sRNAs. Finally, using examples from the few evolutionary studies available, we examine cases of sRNA loss and describe how these may be the result of adaptive in addition to neutral processes. We highlight the need for more-comprehensive analyses of sRNA evolutionary patterns as a means to improve novel sRNA detection, enhance genome annotation, and deepen our understanding of regulatory networks in bacteria.

rRNA Mimicry in RNA Regulation of Gene Expression

The rRNA is the largest and most abundant RNA in bacterial and archaeal cells. It is also one of the best-characterized RNAs in terms of its structural motifs and sequence variation. Production of ribosome components including >50 ribosomal proteins (r-proteins) consumes significant cellular resources. Thus, RNA cis-regulatory structures that interact with r-proteins to repress further r-protein synthesis play an important role in maintaining appropriate stoichiometry between r-proteins and rRNA. Classically, such mRNA structures were thought to directly mimic the rRNA. However, more than 30 years of research has demonstrated that a variety of different recognition and regulatory paradigms are present. This review will demonstrate how structural mimicry between the rRNA and mRNA cis-regulatory structures may take many different forms. The collection of mRNA structures that interact with r-proteins to regulate r-protein operons are best characterized in Escherichia coli, but are increasingly found within species from nearly all phyla of bacteria and several archaea. Furthermore, they represent a unique opportunity to assess the plasticity of RNA structure in the context of RNA-protein interactions. The binding determinants imposed by r-proteins to allow regulation can be fulfilled in many ways. Some r-protein-interacting mRNAs are immediately obvious as rRNA mimics from primary sequence similarity, others are identifiable only after secondary or tertiary structure determination, and some show no obvious similarity. In addition, across different bacterial species a host of different mechanisms of action have been characterized, showing that there is no simple one-size-fits-all solution.

Towards a Dynamic Interaction Network of Life to unify and expand the evolutionary theory

The classic Darwinian theory and the Synthetic evolutionary theory and their linear models, while invaluable to study the origins and evolution of species, are not primarily designed to model the evolution of organisations, typically that of ecosystems, nor that of processes. How could evolutionary theory better explain the evolution of biological complexity and diversity? Inclusive network-based analyses of dynamic systems could retrace interactions between (related or unrelated) components. This theoretical shift from a Tree of Life to a Dynamic Interaction Network of Life, which is supported by diverse molecular, cellular, microbiological, organismal, ecological and evolutionary studies, would further unify evolutionary biology.

Identification of novel MITEs (miniature inverted-repeat transposable elements) in Coxiella burnetii: implications for protein and small RNA evolution

Background

Coxiella burnetii is a Gram-negative gammaproteobacterium and zoonotic agent of Q fever. C. burnetii’s genome contains an abundance of pseudogenes and numerous selfish genetic elements. MITEs (miniature inverted-repeat transposable elements) are non-autonomous transposons that occur in all domains of life and are thought to be insertion sequences (ISs) that have lost their transposase function. Like most transposable elements (TEs), MITEs are thought to play an active role in evolution by altering gene function and expression through insertion and deletion activities. However, information regarding bacterial MITEs is limited.

Results

We describe two MITE families discovered during research on small non-coding RNAs (sRNAs) of C. burnetii. Two sRNAs, Cbsr3 and Cbsr13, were found to originate from a novel MITE family, termed QMITE1. Another sRNA, CbsR16, was found to originate from a separate and novel MITE family, termed QMITE2. Members of each family occur ~ 50 times within the strains evaluated. QMITE1 is a typical MITE of 300-400 bp with short (2-3 nt) direct repeats (DRs) of variable sequence and is often found overlapping annotated open reading frames (ORFs). Additionally, QMITE1 elements possess sigma-70 promoters and are transcriptionally active at several loci, potentially influencing expression of nearby genes. QMITE2 is smaller (150-190 bps), but has longer (7-11 nt) DRs of variable sequences and is mainly found in the 3′ untranslated region of annotated ORFs and intergenic regions. QMITE2 contains a GTAG repetitive extragenic palindrome (REP) that serves as a target for IS1111 TE insertion. Both QMITE1 and QMITE2 display inter-strain linkage and sequence conservation, suggesting that they are adaptive and existed before divergence of C. burnetii strains.

Conclusions

We have discovered two novel MITE families of C. burnetii. Our finding that MITEs serve as a source for sRNAs is novel. QMITE2 has a unique structure and occurs in large or small versions with unique DRs that display linkage and sequence conservation between strains, allowing for tracking of genomic rearrangements. QMITE1 and QMITE2 copies are hypothesized to influence expression of neighboring genes involved in DNA repair and virulence through transcriptional interference and ribonuclease processing.

Electronic supplementary material

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

Impact of bacterial sRNAs in stress responses

Bacterial life is harsh and involves numerous environmental and internal challenges that are perceived as stresses. Consequently, adequate responses to survive, cope with, and counteract stress conditions have evolved. In the last few decades, a class of small, non-coding RNAs (sRNAs) has been shown to be involved as key players in stress responses. This review will discuss — primarily from an enterobacterial perspective — selected stress response pathways that involve antisense-type sRNAs. These include themes of how bacteria deal with severe envelope stress, threats of DNA damage, problems with poisoning due to toxic sugar intermediates, issues of iron homeostasis, and nutrient limitation/starvation. The examples discussed highlight how stress relief can be achieved, and how sRNAs act mechanistically in regulatory circuits. For some cases, we will propose scenarios that may suggest why contributions from post-transcriptional control by sRNAs, rather than transcriptional control alone, appear to be a beneficial and universally selected feature.