Genome rearrangements can make and break small RNA genes

Promoter recruitment drives the emergence of proto-genes in a long-term evolution experiment with Escherichia coli

The phenomenon of de novo gene birth-the emergence of genes from non-genic sequences-has received considerable attention due to the widespread occurrence of genes that are unique to particular species or genomes. Most instances of de novo gene birth have been recognized through comparative analyses of genome sequences in eukaryotes, despite the abundance of novel, lineage-specific genes in bacteria and the relative ease with which bacteria can be studied in an experimental context. Here, we explore the genetic record of the Escherichia coli long-term evolution experiment (LTEE) for changes indicative of "proto-genic" phases of new gene birth in which non-genic sequences evolve stable transcription and/or translation. Over the time span of the LTEE, non-genic regions are frequently transcribed, translated and differentially expressed, with levels of transcription across low-expressed regions increasing in later generations of the experiment. Proto-genes formed downstream of new mutations result either from insertion element activity or chromosomal translocations that fused preexisting regulatory sequences to regions that were not expressed in the LTEE ancestor. Additionally, we identified instances of proto-gene emergence in which a previously unexpressed sequence was transcribed after formation of an upstream promoter, although such cases were rare compared to those caused by recruitment of preexisting promoters. Tracing the origin of the causative mutations, we discovered that most occurred early in the history of the LTEE, often within the first 20,000 generations, and became fixed soon after emergence. Our findings show that proto-genes emerge frequently within evolving populations, can persist stably, and can serve as potential substrates for new gene formation.

Promoter capture drives the emergence of proto-genes in Escherichia coli

The phenomenon of de novo gene birth-the emergence of genes from non-genic sequences-has received considerable attention due to the widespread occurrence of genes that are unique to particular species or genomes. Most instances of de novo gene birth have been recognized through comparative analyses of genome sequences in eukaryotes, despite the abundance of novel, lineage-specific genes in bacteria and the relative ease with which bacteria can be studied in an experimental context. Here, we explore the genetic record of the Escherichia coli Long-Term Evolution Experiment (LTEE) for changes indicative of "proto-genic" phases of new gene birth in which non-genic sequences evolve stable transcription and/or translation. Over the time-span of the LTEE, non-genic regions are frequently transcribed, translated and differentially expressed, thereby serving as raw material for new gene emergence. Most proto-genes result either from insertion element activity or chromosomal translocations that fused pre-existing regulatory sequences to regions that were not expressed in the LTEE ancestor. Additionally, we identified instances of proto-gene emergence in which a previously unexpressed sequence was transcribed after formation of an upstream promoter. Tracing the origin of the causative mutations, we discovered that most occurred early in the history of the LTEE, often within the first 20,000 generations, and became fixed soon after emergence. Our findings show that proto-genes emerge frequently within evolving populations, persist stably, and can serve as potential substrates for new gene formation.

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.

A small RNA is functional in Escherichia fergusonii despite containing a large insertion

Bacterial small RNAs (sRNAs) are important regulators of gene expression; however, the impact of natural mutations on sRNA functions has not been studied extensively. Here we show that the sRNA MgrR contains a unique 53 bp insertion in Escherichia fergusonii, a close relative of Escherichia coli and Salmonella enterica. The insertion is a repetitive extragenic palindromic (REP) sequence that could block transcription, but full-length MgrR is produced in E. fergusonii, showing that the insertion has not affected sRNA production. Additionally, despite containing the large insertion, the sRNA appears to be functional because deletion of mgrR made E. fergusonii more susceptible to H₂O₂. The molecular details of MgrR's roles in H₂O₂defence are yet to be defined, but our results suggest that having an alternative function allowed the sRNA to be retained in E. fergusonii despite it sustaining a large, potentially disruptive mutation.

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.

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.

The Ingenuity of Bacterial Genomes

The genomes of bacteria contain fewer genes and substantially less noncoding DNA than those of eukaryotes, and as a result, they have much less raw material to invent new traits. Yet, bacteria are vastly more taxonomically diverse, numerically abundant, and globally successful in colonizing new habitats compared to eukaryotes. Although bacterial genomes are generally considered to be optimized for efficient growth and rapid adaptation, nonadaptive processes have played a major role in shaping the size, contents, and compact organization of bacterial genomes and have allowed the establishment of deleterious traits that serve as the raw materials for genetic innovation.

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

The evolution of gene expression regulation has contributed to species differentiation. The 3' untranslated regions (3'UTRs) of mRNAs include regulatory elements that modulate gene expression; however, our knowledge of their implications in the divergence of bacterial species is currently limited. In this study, we performed genome-wide comparative analyses of mRNAs encoding orthologous proteins from the genus Staphylococcus and found that mRNA conservation was lost mostly downstream of the coding sequence (CDS), indicating the presence of high sequence diversity in the 3'UTRs of orthologous genes. Transcriptomic mapping of different staphylococcal species confirmed that 3'UTRs were also variable in length. We constructed chimeric mRNAs carrying the 3'UTR of orthologous genes and demonstrated that 3'UTR sequence variations affect protein production. This suggested that species-specific functional 3'UTRs might be specifically selected during evolution. 3'UTR variations may occur through different processes, including gene rearrangements, local nucleotide changes, and the transposition of insertion sequences. By extending the conservation analyses to specific 3'UTRs, as well as the entire set of Escherichia coli and Bacillus subtilis mRNAs, we showed that 3'UTR variability is widespread in bacteria. In summary, our work unveils an evolutionary bias within 3'UTRs that results in species-specific non-coding sequences that may contribute to bacterial diversity.

Prevalence of small base-pairing RNAs derived from diverse genomic loci

Small RNAs (sRNAs) that act by base-pairing have been shown to play important roles in fine-tuning the levels and translation of their target transcripts across a variety of model and pathogenic organisms. Work from many different groups in a wide range of bacterial species has provided evidence for the importance and complexity of sRNA regulatory networks, which allow bacteria to quickly respond to changes in their environment. However, despite the expansive literature, much remains to be learned about all aspects of sRNA-mediated regulation, particularly in bacteria beyond the well-characterized Escherichia coli and Salmonella enterica species. Here we discuss what is known, and what remains to be learned, about the identification of regulatory base-pairing RNAs produced from diverse genomic loci including how their expression is regulated. This article is part of a Special Issue entitled: RNA and gene control in bacteria edited by Dr. M. Guillier and F. Repoila.

RNA-Sequencing Analyses of Small Bacterial RNAs and their Emergence as Virulence Factors in Host-Pathogen Interactions

Proteins have long been considered to be the most prominent factors regulating so-called invasive genes involved in host-pathogen interactions. The possible role of small non-coding RNAs (sRNAs), either intracellular, secreted or packaged in outer membrane vesicles (OMVs), remained unclear until recently. The advent of high-throughput RNA-sequencing (RNA-seq) techniques has accelerated sRNA discovery. RNA-seq radically changed the paradigm on bacterial virulence and pathogenicity to the point that sRNAs are emerging as an important, distinct class of virulence factors in both gram-positive and gram-negative bacteria. The potential of OMVs, as protectors and carriers of these functional, gene regulatory sRNAs between cells, has also provided an additional layer of complexity to the dynamic host-pathogen relationship. Using a non-exhaustive approach and through examples, this review aims to discuss the involvement of sRNAs, either free or loaded in OMVs, in the mechanisms of virulence and pathogenicity during bacterial infection. We provide a brief overview of sRNA origin and importance, and describe the classical and more recent methods of identification that have enabled their discovery, with an emphasis on the theoretical lower limit of RNA sizes considered for RNA sequencing and bioinformatics analyses.

It's a small, small world: unravelling the role and evolution of small RNAs in organelle and endosymbiont genomes

Organelles and host-restricted bacterial symbionts are characterized by having highly reduced genomes that lack many key regulatory genes and elements. Thus, it has been hypothesized that the eukaryotic nuclear genome is primarily responsible for regulating these symbioses. However, with the discovery of organelle- and symbiont-expressed small RNAs (sRNAs) there is emerging evidence that these sRNAs may play a role in gene regulation as well. Here, we compare the diversity of organelle and bacterial symbiont sRNAs recently identified using genome-enabled '-omic' technologies and discuss their potential role in gene regulation. We also discuss how the genome architecture of small genomes may influence the evolution of these sRNAs and their potential function. Additionally, these new studies suggest that some sRNAs are conserved within organelle and symbiont taxa and respond to changes in the environment and/or their hosts. In summary, these results suggest that organelle and symbiont sRNAs may play a role in gene regulation in addition to nuclear-encoded host mechanisms.

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.

Functional annotation and distribution overview of RNA families in 27 Streptococcus agalactiae genomes

Background

Streptococcus agalactiae, also known as Group B Streptococcus (GBS), is a Gram-positive bacterium that colonizes the gastrointestinal and genitourinary tract of humans. This bacterium has also been isolated from various animals, such as fish and cattle. Non-coding RNAs (ncRNAs) can act as regulators of gene expression in bacteria, such as Streptococcus pneumoniae and Streptococcus pyogenes. However, little is known about the genomic distribution of ncRNAs and RNA families in S. agalactiae.

Results

Comparative genome analysis of 27 S. agalactiae strains showed more than 5 thousand genomic regions identified and classified as Core, Exclusive, and Shared genome sequences. We identified 27 to 89 RNA families per genome distributed over these regions, from these, 25 were in Core regions while Shared and Exclusive regions showed variations amongst strains. We propose that the amount and type of ncRNA present in each genome can provide a pattern to contribute in the identification of the clonal types.

Conclusions

The identification of RNA families provides an insight over ncRNAs, sRNAs and ribozymes function, that can be further explored as targets for antibiotic development or studied in gene regulation of cellular processes. RNA families could be considered as markers to determine infection capabilities of different strains. Lastly, pan-genome analysis of GBS including the full range of functional transcripts provides a broader approach in the understanding of this pathogen.

Electronic supplementary material

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

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

Non-coding small RNAs (sRNAs) are critical to post-transcriptional gene regulation in bacteria. However, unlike for protein-coding genes, the evolutionary forces that shape sRNAs are not understood. We investigated sRNAs in enteric bacteria and discovered that recently emerged sRNAs evolve at significantly faster rates than older sRNAs. Concomitantly, younger sRNAs are expressed at significantly lower levels than older sRNAs. This process could potentially facilitate the integration of newly emerged sRNAs into bacterial regulatory networks. Furthermore, it has previously been difficult to trace the evolutionary histories of sRNAs because rapid evolution obscures their original sources. We overcame this challenge by identifying a recently evolved sRNA in Escherichia coli, which allowed us to determine that novel sRNAs could emerge from vestigial bacteriophage genes, the first known source for sRNA origination.

Competing for Iron: Duplication and Amplification of the isd Locus in Staphylococcus lugdunensis HKU09-01 Provides a Competitive Advantage to Overcome Nutritional Limitation

Staphylococcus lugdunensis is a coagulase negative bacterial pathogen that is particularly associated with severe cases of infectious endocarditis. Unique amongst the coagulase-negative staphylococci, S. lugdunensis harbors an iron regulated surface determinant locus (isd). This locus facilitates the acquisition of heme as a source of nutrient iron during infection and allows iron limitation caused by “nutritional immunity” to be overcome. The isd locus is duplicated in S. lugdunensis HKU09-01 and we show here that the duplication is intrinsically unstable and undergoes accordion-like amplification and segregation leading to extensive isd copy number variation. Amplification of the locus increased the level of expression of Isd proteins and improved binding of hemoglobin to the cell surface of S. lugdunensis. Furthermore, Isd overexpression provided an advantage when strains were competing for a limited amount of hemoglobin as the sole source of iron. Gene duplications and amplifications (GDA) are events of fundamental importance for bacterial evolution and are frequently associated with antibiotic resistance in many species. As such, GDAs are regarded as evolutionary adaptions to novel selective pressures in hostile environments pointing towards a special importance of isd for S. lugdunensis. For the first time we show an example of a GDA that involves a virulence factor of a Gram-positive pathogen and link the GDA directly to a competitive advantage when the bacteria were struggling with selective pressures mimicking “nutritional immunity”.

Sometimes changing environmental conditions force bacteria to boost protein expression above the level that can be achieved by transcriptional or translational control. Gene duplication and amplification (GDA) represents a simple and effective means to augment protein expression and is frequently associated with the development of resistance against antibacterial agents. As such GDAs can be seen as recent evolutionary adaptions towards novel selective pressures within the environment. We describe herein the GDA of the isd locus in a clinical isolate of Staphylococcus lugdunensis a bacterial species frequently associated with infectious endocarditis. This amplification made individual lineages significantly more successful when competing for the trace element iron supplied in the form of hemoglobin (hb). Human body fluids are actively depleted of iron to prevent bacterial proliferation and hb represents the most abundant iron source during infection. This deepens our understanding about adaption processes of pathogens and demonstrates how selective pressures drive the evolution of pathogens to become more successful. This study suggests that investigating GDAs in clinical isolates can help to identify chromosomal regions that are of special importance for bacterial fitness under changing environmental conditions. The gene products can therefore represent attractive targets for novel antibiotics supporting the natural immune defenses by targeting the vulnerable spots of the bacterial anti-immunity strategies.

Identification and Characterization of Novel Small RNAs in Rickettsia prowazekii

Emerging evidence implicates a critically important role for bacterial small RNAs (sRNAs) as post-transcriptional regulators of physiology, metabolism, stress/adaptive responses, and virulence, but the roles of sRNAs in pathogenic Rickettsia species remain poorly understood. Here, we report on the identification of both novel and well-known bacterial sRNAs in Rickettsia prowazekii, known to cause epidemic typhus in humans. RNA sequencing of human microvascular endothelial cells (HMECs), the preferred targets during human rickettsioses, infected with R. prowazekii revealed the presence of 35 trans-acting and 23 cis-acting sRNAs, respectively. Of these, expression of two trans-acting (Rp_sR17 and Rp_sR60) and one cis-acting (Rp_sR47) novel sRNAs and four well-characterized bacterial sRNAs (RNaseP_bact_a, α-tmRNA, 4.5S RNA, 6S RNA) was further confirmed by Northern blot or RT-PCR analyses. The transcriptional start sites of five novel rickettsial sRNAs and 6S RNA were next determined using 5' RLM-RACE yielding evidence for their independent biogenesis in R. prowazekii. Finally, computational approaches were employed to determine the secondary structures and potential mRNA targets of novel sRNAs. Together, these results establish the presence and expression of sRNAs in R. prowazekii during host cell infection and suggest potential functional roles for these important post-transcriptional regulators in rickettsial biology and pathogenesis.

How do base-pairing small RNAs evolve?

The increasing numbers of characterized base-pairing small RNAs (sRNAs) and the identification of these regulators in a broad range of bacteria are allowing comparisons between species and explorations of sRNA evolution. In this review, we describe some examples of trans-encoded base-pairing sRNAs that are species-specific and others that are more broadly distributed. We also describe examples of sRNA orthologs where different features are conserved. These examples provide the background for a discussion of mechanisms of sRNA evolution and selective pressures on the sRNAs and their mRNA target(s).