Identification of differentially expressed small non-coding RNAs in the legume endosymbiont Sinorhizobium meliloti by comparative genomics

An Integrated Affinity Chromatography-Based Approach to Unravel the sRNA Interactome in Nitrogen-Fixing Rhizobia

The activity mechanism and function of bacterial base-pairing small non-coding RNA regulators (sRNAs) are largely shaped by their main interacting cellular partners, i.e., proteins and mRNAs. We describe here an MS2 affinity chromatography-based procedure adapted to unravel the sRNA interactome in nitrogen-fixing legume endosymbiotic bacteria. The method consists of tagging of the bait sRNA at its 5'-end with the MS2 aptamer followed by pulse overexpression and immobilization of the chimeric transcript from cell lysates by an MS2-MBP fusion protein conjugated to an amylose resin. The sRNA-binding proteins and target mRNAs are further profiled by mass spectrometry and RNAseq, respectively.

A double-negative feedback loop between NtrBC and a small RNA rewires nitrogen metabolism in legume symbionts

The nitrogen (N) status transduced via the NtrBC two-component system is a major signaling cue in the root nodule endosymbiosis of diazotrophic rhizobia with legumes. NtrBC is upregulated in the N-limiting rhizosphere environment at the onset of nodulation but silenced in nodules to favor the assimilation of the fixed N into plant biomass. We reported that the trans-acting sRNA NfeR1 (Nodule Formation Efficiency RNA) broadly influences the symbiotic performance of the α-rhizobium Sinorhizobium meliloti. Here, we show that NfeR1 is indeed an N-responsive sRNA that fine-tunes NtrBC output during the symbiotic transition. Biochemical and genetic approaches unveiled that NtrC and the LysR-type symbiotic regulator LsrB bind at distinct nearby sites in the NfeR1 promoter, acting antagonistically as repressor and activator of transcription, respectively. This complex transcriptional control specifies peak NfeR1 steady-state levels in N-starved and endosymbiotic bacteria. Furthermore, NfeR1 base pairs the translation initiation region of the histidine kinase coding mRNA ntrB, causing a decrease in both NtrB and NtrC abundance as assessed by double-plasmid genetic assays. In the context of endogenous regulation, NfeR1-mediated ntrBC silencing most likely amends the effective strength of the known operon autorepression exerted by NtrC. Accordingly, a lack of NfeR1 shifts the wild-type NtrBC output, restraining the fitness of free-living rhizobia under N stress and plant growth upon nodulation. The mixed NtrBC-NfeR1 double-negative feedback loop is thus an unprecedented adaptive network motif that helps α-rhizobia adjust N metabolism to the demands of an efficient symbiosis with legume plants. IMPORTANCE Root nodule endosymbioses between diazotrophic rhizobia and legumes provide the largest input of combined N to the biosphere, thus representing an alternative to harmful chemical fertilizers for sustainable crop production. Rhizobia have evolved intricate strategies to coordinate N assimilation for their own benefit with N2 fixation to sustain plant growth. The rhizobial N status is transduced by the NtrBC two-component system, the seemingly ubiquitous form of N signal transduction in Proteobacteria. Here, we show that the regulatory sRNA NfeR1 (nodule formation efficiency RNA) of the alfalfa symbiont Sinorhizobium meliloti is transcribed from a complex promoter repressed by NtrC in a N-dependent manner and feedback silences ntrBC by complementary base-pairing. These findings unveil a more prominent role of NtrC as a transcriptional repressor than hitherto anticipated and a novel RNA-based mechanism for NtrBC regulation. The NtrBC-NfeR1 double-negative feedback loop accurately rewires symbiotic S. meliloti N metabolism and is likely conserved in α-rhizobia.

Identification of sigma factor 54-regulated small non-coding RNAs by employing genome-wide and transcriptome-based methods in rhizobium strains

Rhizobium-legume symbiosis is considered as the major contributor of biological nitrogen fixation. Bacterial small non-coding RNAs are crucial regulators in several cellular adaptation processes that occur due to the changes in metabolism, physiology, or the external environment. Identifying and analysing the conditional specific/sigma factor-54 regulated sRNAs provides a better understanding of sRNA regulation/mechanism in symbiotic association. In the present study, we have identified sigma factor 54-regulated sRNAs from the genome of six rhizobium strains and from the RNA-seq data of free-living and symbiotic conditions of Bradyrhizobium diazoefficiens USDA 110 to identify the novel putative sRNAs that are over expressed during the regulation of nitrogen fixation. A total of 1351 sRNAs were predicted from the genome of six rhizobium strains and 1375 sRNAs were predicted from the transcriptome data of B. diazoefficiens USDA 110. Analysis of target mRNA for these novel sRNAs was inferred to target several nodulation and nitrogen fixation genes including nodC, nodJ, nodY, nodJ, nodM, nodW, nodZ, nifD, nifN, nifQ, fixK, fixL, fdx, nolB, and several cytochrome proteins. In addition, sRNAs of B. diazoefficiens USDA 110 which targeted the regulatory genes of nitrogen fixation were confirmed by wet-lab experiments with semi-quantitative reverse transcription polymerase chain reaction. Predicted target mRNAs were functionally classified based on the COG analysis and GO annotations. The genome-wide and transcriptome-based integrated methods have led to the identification of several sRNAs involved in the nodulation and symbiosis. Further validation of the functional role of these sRNAs can help in exploring the role of sRNAs in nitrogen metabolism during free-living and symbiotic association with legumes.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03394-x.

Do small RNAs unlock the below ground microbiome-plant interaction mystery?

Over the past few decades, regulatory RNAs, such as small RNAs (sRNAs), have received increasing attention in the context of host-microbe interactions due to their diverse roles in controlling various biological processes in eukaryotes. In addition, studies have identified an increasing number of sRNAs with novel functions across a wide range of bacteria. What is not well understood is why cells regulate gene expression through post-transcriptional mechanisms rather than at the initiation of transcription. The finding of a multitude of sRNAs and their identified associated targets has allowed further investigation into the role of sRNAs in mediating gene regulation. These foundational data allow for further development of hypotheses concerning how a precise control of gene activity is accomplished through the combination of transcriptional and post-transcriptional regulation. Recently, sRNAs have been reported to participate in interkingdom communication and signalling where sRNAs originating from one kingdom are able to target or control gene expression in another kingdom. For example, small RNAs of fungal pathogens that silence plant genes and vice-versa plant sRNAs that mediate bacterial gene expression. However, there is currently a lack of evidence regarding sRNA-based inter-kingdom signalling across more than two interacting organisms. A habitat that provides an excellent opportunity to investigate interconnectivity is the plant rhizosphere, a multifaceted ecosystem where plants and associated soil microbes are known to interact. In this paper, we discuss how the interconnectivity of bacteria, fungi, and plants within the rhizosphere may be mediated by bacterial sRNAs with a particular focus on disease suppressive and non-suppressive soils. We discuss the potential roles sRNAs may play in the below-ground world and identify potential areas of future research, particularly in reference to the regulation of plant immunity genes by bacterial and fungal communities in disease-suppressive and non-disease-suppressive soils.

Synthetase of the methyl donor S-adenosylmethionine from nitrogen-fixing α-rhizobia can bind functionally diverse RNA species

Function of bacterial small non-coding RNAs (sRNAs) and overall RNA metabolism is largely shaped by a vast diversity of RNA-protein interactions. However, in non-model bacteria with defined non-coding transcriptomes the sRNA interactome remains almost unexplored. We used affinity chromatography to capture proteins associated in vivo with MS2-tagged trans-sRNAs that regulate nutrient uptake (AbcR2 and NfeR1) and cell cycle (EcpR1) mRNAs by antisense-based translational inhibition in the nitrogen-fixing α-rhizobia Sinorhizobium meliloti. The three proteomes were rather distinct, with that of EcpR1 particularly enriched in cell cycle-related enzymes, whilst sharing several transcription/translation-related proteins recurrently identified associated with sRNAs. Strikingly, MetK, the synthetase of the major methyl donor S-adenosylmethionine, was reliably recovered as a binding partner of the three sRNAs, which reciprocally co-immunoprecipitated with a FLAG-tagged MetK variant. Induced (over)expression of the trans-sRNAs and MetK depletion did not influence canonical riboregulatory traits, `for example, protein titration or sRNA stability, respectively. An in vitro filter assay confirmed binding of AbcR2, NfeR1 and EcpR1 to MetK and further revealed interaction of the protein with other non-coding and coding transcripts but not with the 5S rRNA. These findings uncover a broad specificity for RNA binding as an unprecedented feature of this housekeeping prokaryotic enzyme.

Evolutionary conservation of RNA sequence and structure

An RNA structure prediction from a single‐sequence RNA folding program is not evidence for an RNA whose structure is important for function. Random sequences have plausible and complex predicted structures not easily distinguishable from those of structural RNAs. How to tell when an RNA has a conserved structure is a question that requires looking at the evolutionary signature left by the conserved RNA. This question is important not just for long noncoding RNAs which usually lack an identified function, but also for RNA binding protein motifs which can be single stranded RNAs or structures. Here we review recent advances using sequence and structural analysis to determine when RNA structure is conserved or not. Although covariation measures assess structural RNA conservation, one must distinguish covariation due to RNA structure from covariation due to independent phylogenetic substitutions. We review a statistical test to measure false positives expected under the null hypothesis of phylogenetic covariation alone (specificity). We also review a complementary test that measures power, that is, expected covariation derived from sequence variation alone (sensitivity). Power in the absence of covariation signals the absence of a conserved RNA structure. We analyze artifacts that falsely identify conserved RNA structure such as the misuse of programs that do not assess significance, the use of inappropriate statistics confounded by signals other than covariation, or misalignments that induce spurious covariation. Among artifacts that obscure the signal of a conserved RNA structure, we discuss the inclusion of pseudogenes in alignments which increase power but destroy covariation.

This article is categorized under:

RNA Structure and Dynamics > RNA Structure, Dynamics and Chemistry

RNA Evolution and Genomics > Computational Analyses of RNA

RNA Evolution and Genomics > RNA and Ribonucleoprotein Evolution

Pervasive RNA Regulation of Metabolism Enhances the Root Colonization Ability of Nitrogen-Fixing Symbiotic α-Rhizobia

The rhizosphere and rhizoplane are nutrient-rich but selective environments for the root microbiome. Here, we deciphered a posttranscriptional network regulated by the homologous trans-small RNAs (sRNAs) AbcR1 and AbcR2, which rewire the metabolism of the nitrogen-fixing α-rhizobium Sinorhizobium meliloti during preinfection stages of symbiosis with its legume host alfalfa. The LysR-type regulator LsrB, which transduces the cell redox state, is indispensable for AbcR1 expression in actively dividing bacteria, whereas the stress-induced transcription of AbcR2 depends on the alternative σ factor RpoH1. MS2 affinity purification coupled with RNA sequencing unveiled exceptionally large and overlapping AbcR1/2 mRNA interactomes, jointly representing ⁓6% of the S. meliloti protein-coding genes. Most mRNAs encode transport/metabolic proteins whose translation is silenced by base pairing to two distinct anti-Shine Dalgarno motifs that function independently in both sRNAs. A metabolic model-aided analysis of the targetomes predicted changes in AbcR1/2 expression driven by shifts in carbon/nitrogen sources, which were confirmed experimentally. Low AbcR1/2 levels in some defined media anticipated overexpression growth phenotypes linked to the silencing of specific mRNAs. As a proof of principle, we confirmed AbcR1/2-mediated downregulation of the l-amino acid AapQ permease. AbcR1/2 interactomes are well represented in rhizosphere-related S. meliloti transcriptomic signatures. Remarkably, a lack of AbcR1 specifically compromised the ability of S. meliloti to colonize the root rhizoplane. The AbcR1 regulon likely ranks the utilization of available substrates to optimize metabolism, thus conferring on S. meliloti an advantage for efficient rhizosphere/rhizoplane colonization. AbcR1 regulation is predicted to be conserved in related α-rhizobia, which opens unprecedented possibilities for engineering highly competitive biofertilizers. IMPORTANCE Nitrogen-fixing root nodule symbioses between rhizobia and legume plants provide more than half of the combined nitrogen incorporated annually into terrestrial ecosystems, rendering plant growth independent of environmentally unfriendly chemical fertilizers. The success of symbiosis depends primarily on the capacity of rhizobia to establish competitive populations in soil and rhizosphere environments. Here, we provide insights into the regulation and architecture of an extensive RNA posttranscriptional network that fine-tunes the metabolism of the alfalfa symbiont S. meliloti, thereby enhancing the ability of this beneficial bacterium to colonize nutrient-rich but extremely selective niches, such as the rhizosphere of its host plant. This pervasive RNA regulation of metabolism is a major adaptive mechanism, predicted to operate in diverse rhizobial species. Because RNA regulation relies on modifiable base-pairing interactions, our findings open unexplored avenues for engineering the legumes rhizobiome within sustainable agricultural practices.

Similarities and differences between 6S RNAs from Bradyrhizobium japonicum and Sinorhizobium meliloti

6S RNA, a conserved and abundant small non-coding RNA found in most bacteria, regulates gene expression by inhibiting RNA polymerase (RNAP) holoenzyme. 6S RNAs from α-proteobacteria have been studied poorly so far. Here, we present a first in-depth analysis of 6S RNAs from two α-proteobacteria species, Bradyrhizobium japonicum and Sinorhizobium meliloti. Although both belong to the order Rhizobiales and are typical nitrogen-fixing symbionts of legumes, their 6S RNA expression profiles were found to differ: B. japonicum 6S RNA accumulated in the stationary phase, thus being reminiscent of Escherichia coli 6S RNA, whereas S. meliloti 6S RNA level peaked at the transition to the stationary phase, similarly to Rhodobacter sphaeroides 6S RNA. We demonstrated in vitro that both RNAs have hallmarks of 6S RNAs: they bind to the σ⁷⁰-type RNAP holoenzyme and serve as templates for de novo transcription of so-called product RNAs (pRNAs) ranging in length from ∼13 to 24 nucleotides, with further evidence of the synthesis of even longer pRNAs. Likewise, stably bound pRNAs were found to rearrange the 6S RNA structure to induce its dissociation from RNAP. Compared with B. japonicum 6S RNA, considerable conformational heterogeneity was observed for S. meliloti 6S RNA and its complexes with pRNAs, even though the two 6S RNAs share ∼75% sequence identity. Overall, our findings suggest that the two rhizobial 6S RNAs have diverged with respect to their regulatory impact on gene expression throughout the bacterial life cycle.

Riboregulation in Nitrogen-Fixing Endosymbiotic Bacteria

Small non-coding RNAs (sRNAs) are ubiquitous components of bacterial adaptive regulatory networks underlying stress responses and chronic intracellular infection of eukaryotic hosts. Thus, sRNA-mediated regulation of gene expression is expected to play a major role in the establishment of mutualistic root nodule endosymbiosis between nitrogen-fixing rhizobia and legume plants. However, knowledge about this level of genetic regulation in this group of plant-interacting bacteria is still rather scarce. Here, we review insights into the rhizobial non-coding transcriptome and sRNA-mediated post-transcriptional regulation of symbiotic relevant traits such as nutrient uptake, cell cycle, quorum sensing, or nodule development. We provide details about the transcriptional control and protein-assisted activity mechanisms of the functionally characterized sRNAs involved in these processes. Finally, we discuss the forthcoming research on riboregulation in legume symbionts.

Navigation through the twists and turns of RNA sequencing technologies: Application to bacterial regulatory RNAs

Discovered in the 1980s, small regulatory RNAs (sRNAs) are now considered key actors in virtually all aspects of bacterial physiology and virulence. Together with transcriptional and translational regulatory proteins, they integrate and often are hubs of complex regulatory networks, responsible for bacterial response/adaptation to various perceived stimuli. The recent development of powerful RNA sequencing technologies has facilitated the identification and characterization of sRNAs (length, structure and expression conditions) and their RNA targets in several bacteria. Nevertheless, it could be very difficult for non-experts to understand the advantages and drawbacks related to each offered option and, consequently, to make an informed choice. Therefore, the main goal of this review is to provide a guide to navigate through the twists and turns of high-throughput RNA sequencing technologies, with a specific focus on those applied to the study of sRNAs. This article is part of a Special Issue entitled: RNA and gene control in bacteria edited by Dr. M. Guillier and F. Repoila.

Targeting Hsp20 Using the Novel Small Non-coding RNA DnrH Regulates Heat Tolerance in Deinococcus radiodurans

Small non-coding RNAs (ncRNAs) are a class of regulatory molecules, which remain understudied in bacteria. In the extremophilic bacterium Deinococcus radiodurans, although hundreds of ncRNAs have been identified, few have been characterized in detail. In this study, we report the identification and characterization of a novel heat-inducible ncRNA named DnrH. Heat tolerance analysis showed that deleting DnrH significantly inhibited viability in response to high temperature conditions. Comparative phenotypic and qRT-PCR analyses of a DnrH mutant (ΔDnrH) and wild-type (WT) D. radiodurans suggested that DnrH is potentially involved in regulating the expression of the heat shock-related gene Hsp20. Microscale thermophoresis and genetic complementation showed that a 28-nucleotide (nt) sequence in the stem-loop structure of DnrH (143–170 nt) pairs with its counterpart in the coding region of Hsp20 mRNA (91–117 nt) via a 22 nt region. In vivo, mutation of the 22-nt region in the D. radiodurans genome led to a reduction in heat tolerance similar to that observed in the DnrH-mutant. Our results show that DnrH positively influences heat tolerance by increasing the transcription of Hsp20 mRNA, demonstrating, for the first time, a ncRNA that directly controls the expression of a heat stress-resistance gene. This work provides new insight into the heat stress response mechanism of D. radiodurans as well as other extremophiles that express similar Hsp20 proteins.

Polyamine biosynthesis and biological roles in rhizobia

Polyamines are ubiquitous molecules containing two or more amino groups that fulfill varied and often essential physiological and regulatory roles in all organisms. In the symbiotic nitrogen-fixing bacteria known as rhizobia, putrescine and homospermidine are invariably produced while spermidine and norspermidine synthesis appears to be restricted to the alfalfa microsymbiont Sinorhizobium meliloti. Studies with rhizobial mutants deficient in the synthesis of one or more polyamines have shown that these compounds are important for growth, stress resistance, motility, exopolysaccharide production and biofilm formation. In this review, we describe these studies and examine how polyamines are synthesized and regulated in rhizobia.

Identification of Small RNA-Protein Partners in Plant Symbiotic Bacteria

The identification of the protein partners of bacterial small noncoding RNAs (sRNAs) is essential to understand the mechanistic principles and functions of riboregulation in prokaryotic cells. Here, we describe an optimized affinity chromatography protocol that enables purification of in vivo formed sRNA-protein complexes in Sinorhizobium meliloti, a genetically tractable nitrogen-fixing plant symbiotic bacterium. The procedure requires the tagging of the desired sRNA with the MS2 aptamer, which is affinity-captured by the MS2-MBP protein conjugated to an amylose resin. As proof of principle, we show recovery of the RNA chaperone Hfq associated to the strictly Hfq-dependent AbcR2 trans-sRNA. This method can be applied for the investigation of sRNA-protein interactions on a broad range of genetically tractable α-proteobacteria.

Primary Characterization of Small RNAs in Symbiotic Nitrogen-Fixing Bacteria

High-throughput transcriptome profiling (RNAseq) has uncovered large and heterogeneous populations of small noncoding RNA species (sRNAs) with potential regulatory roles in bacteria. A large fraction of sRNAs are differentially regulated and rely on protein-assisted antisense interactions to trans-encoded target mRNAs to fine-tune posttranscriptional reprogramming of gene expression in response to external cues. However, annotation and function of sRNAs are still largely overlooked in nonmodel bacteria with complex lifestyles. Here, we describe experimental protocols successfully applied for the accurate annotation, expression profiling and target mRNA identification of trans-acting sRNAs in the nitrogen-fixing α-rhizobium Sinorhizobium meliloti. The protocols presented here can be similarly applied for the characterization of trans-sRNAs in genetically tractable α-proteobacteria of agronomical or clinical relevance interacting with eukaryotic hosts.

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.

Use of RNA Immunoprecipitation Method for Determining Sinorhizobium meliloti RNA-Hfq Protein Associations In Vivo

Background

Soil bacterium Sinorhizobium meliloti (S. meliloti) forms an endosymbiotic partnership with Medicago truncatula (M. truncatula) roots which results in root nodules. The bacteria live within root nodules where they function to fix atmospheric N₂ and supply the host plant with reduced nitrogen. The bacterial RNA-binding protein Hfq (Hfq) is an important regulator for the effectiveness of the nitrogen fixation. RNA immunoprecipitation (RIP) method is a powerful method for detecting the association of Hfq protein with specific RNA in cultured bacteria, yet a RIP method for bacteria living in root nodules remains to be described.

Results

A modified S. meliloti gene encoding a His-tagged Hfq protein (Hfq^His) was placed under the regulation of the native Hfq gene promoter (P_hfqsm). The trans produced Hfq^His protein was accumulated at its nature levels during all stages of the symbiosis, allowing RNAs that associated with the given protein to be immunoprecipitated with the anti-His antibody against the protein from root nodule lysates. RNAs that associated with the protein were selectively enriched in the immunoprecipitated sample. The RNAs were recovered by a simple method using heat and subsequently analyzed by RT-PCR. The nature of PCR products was determined by DNA sequencing. Hfq association with specific RNAs can be analyzed at different conditions (e. g. young or older root nodules) and/or in wild-type versus mutant strains.

Conclusions

This article describes the RIP method for determining Sinorhizobium meliloti RNA-Hfq associations in vivo. It is also applicable to other rhizobia living in planta, although some tissue-specific modification related to sample disruption and homogenization may be needed.

Challenges of ligand identification for the second wave of orphan riboswitch candidates

Orphan riboswitch candidates are noncoding RNA motifs whose representatives are believed to function as genetic regulatory elements, but whose target ligands have yet to be identified. The study of certain orphans, particularly classes that have resisted experimental validation for many years, has led to the discovery of important biological pathways and processes once their ligands were identified. Previously, we highlighted details for four of the most common and intriguing orphan riboswitch candidates. This facilitated the validation of riboswitches for the signaling molecules c-di-AMP, ZTP, and ppGpp, the metal ion Mn2+, and the metabolites guanidine and PRPP. Such studies also yield useful linkages between the ligands sensed by the riboswitches and numerous biochemical pathways. In the current report, we describe the known characteristics of 30 distinct classes of orphan riboswitch candidates - some of which have remained unsolved for over a decade. We also discuss the prospects for uncovering novel biological insights via focused studies on these RNAs. Lastly, we make recommendations for experimental objectives along the path to finding ligands for these mysterious RNAs.

An account of evolutionary specialization: the AbcR small RNAs in the Rhizobiales

The AbcR small RNAs (sRNAs) are a fascinating example of two highly conserved sRNAs that differ tremendously at the functional level among organisms. From their transcriptional activation to their regulatory capabilities, the AbcR sRNAs exhibit varying characteristics in three well-studied bacteria belonging to the Rhizobiales order: the plant symbiont Sinorhizobium meliloti, the plant pathogen Agrobacterium tumefaciens, and the animal pathogen Brucella abortus. This review outlines the similarities and differences of the AbcR sRNAs between each of these organisms, and discusses reasons as to why this group of sRNAs has diverged in their genetic organization and regulatory functions across species. In the end, this review will shed light on how regulatory systems, although seemingly conserved among bacteria, can vary based on the environmental niche and lifestyle of an organism.

RNA silencing in plant symbiotic bacteria: Insights from a protein-centric view

Extensive work in model enterobacteria has evidenced that the RNA chaperone Hfq and several endoribonucleases, such as RNase E or RNase III, serve pivotal roles in small RNA-mediated post-transcriptional silencing of gene expression. Characterization of these protein hubs commonly provide global functional and mechanistic insights into complex sRNA regulatory networks. The legume endosymbiont Sinorhizobium meliloti is a non-classical model bacterium with a very complex lifestyle in which riboregulation is expected to play important adaptive functions. Here, we discuss current knowledge about RNA silencing in S. meliloti from the perspective of the activity of Hfq and a recently discovered endoribonuclease (YbeY) exhibiting unprecedented catalytic versatility for the cleavage of single- and double-stranded RNA molecules.

A Genome-Wide Prediction and Identification of Intergenic Small RNAs by Comparative Analysis in Mesorhizobium huakuii 7653R

In bacteria, small non-coding RNAs (sRNAs) are critical regulators of cellular adaptation to changes in metabolism, physiology, or the external environment. In the last decade, more than 2000 of sRNA families have been reported in the Rfam database and have been shown to exert various regulatory functions in bacterial transcription and translation. However, little is known about sRNAs and their functions in Mesorhizobium. Here, we predicted putative sRNAs in the intergenic regions (IGRs) of M. huakuii 7653R by genome-wide comparisons with four related Mesorhizobial strains. The expression and transcribed regions of candidate sRNAs were analyzed using a set of high-throughput RNA deep sequencing data. In all, 39 candidate sRNAs were found, with 5 located in the symbiotic megaplasmids and 34 in the chromosome of M. huakuii 7653R. Of these, 24 were annotated as functional sRNAs in the Rfam database and 15 were recognized as putative novel sRNAs. The expression of nine selected sRNAs was confirmed by Northern blotting, and most of the nine selected sRNAs were highly expressed in 28 dpi nodules and under symbiosis-mimicking conditions. For those putative novel sRNAs, functional categorizations of their target genes were performed by analyzing the enriched GO terms. In addition, MH_s15 was shown to be an abundant and conserved sRNA.

A 6-Nucleotide Regulatory Motif within the AbcR Small RNAs of Brucella abortus Mediates Host-Pathogen Interactions

In Brucella abortus, two small RNAs (sRNAs), AbcR1 and AbcR2, are responsible for regulating transcripts encoding ABC-type transport systems. AbcR1 and AbcR2 are required for Brucella virulence, as a double chromosomal deletion of both sRNAs results in attenuation in mice. Although these sRNAs are responsible for targeting transcripts for degradation, the mechanism utilized by the AbcR sRNAs to regulate mRNA in Brucella has not been described. Here, two motifs (M1 and M2) were identified in AbcR1 and AbcR2, and complementary motif sequences were defined in AbcR-regulated transcripts. Site-directed mutagenesis of M1 or M2 or of both M1 and M2 in the sRNAs revealed transcripts to be targeted by one or both motifs. Electrophoretic mobility shift assays revealed direct, concentration-dependent binding of both AbcR sRNAs to a target mRNA sequence. These experiments genetically and biochemically characterized two indispensable motifs within the AbcR sRNAs that bind to and regulate transcripts. Additionally, cellular and animal models of infection demonstrated that only M2 in the AbcR sRNAs is required for Brucella virulence. Furthermore, one of the M2-regulated targets, BAB2_0612, was found to be critical for the virulence of B. abortus in a mouse model of infection. Although these sRNAs are highly conserved among Alphaproteobacteria, the present report displays how gene regulation mediated by the AbcR sRNAs has diverged to meet the intricate regulatory requirements of each particular organism and its unique biological niche.

Small RNAs (sRNAs) are important components of bacterial regulation, allowing organisms to quickly adapt to changes in their environments. The AbcR sRNAs are highly conserved throughout the Alphaproteobacteria and negatively regulate myriad transcripts, many encoding ABC-type transport systems. In Brucella abortus, AbcR1 and AbcR2 are functionally redundant, as only a double abcR1 abcR2 (abcR1/2) deletion results in attenuation in vitro and in vivo. In the present study, we confirmed that the AbcR sRNAs have redundant regulatory functions and defined two six-nucleotide motifs, M1 and M2, that the AbcR sRNAs utilize to control gene expression. Importantly, only M2 was linked to B. abortus virulence. Further investigation of M2-regulated targets identified BAB2_0612 as critical for colonization of B. abortus in mice, highlighting the significance of AbcR M2-regulated transcripts for Brucella infection. Overall, our findings define the molecular mechanism of the virulence-associated AbcR system in the pathogenic bacterium B. abortus.

A conserved α-proteobacterial small RNA contributes to osmoadaptation and symbiotic efficiency of rhizobia on legume roots

Small non-coding RNAs (sRNAs) are expected to have pivotal roles in the adaptive responses underlying symbiosis of nitrogen-fixing rhizobia with legumes. Here, we provide primary insights into the function and activity mechanism of the Sinorhizobium meliloti trans-sRNA NfeR1 (Nodule Formation Efficiency RNA). Northern blot probing and transcription tracking with fluorescent promoter-reporter fusions unveiled high nfeR1 expression in response to salt stress and throughout the symbiotic interaction. The strength and differential regulation of nfeR1 transcription are conferred by a motif, which is conserved in nfeR1 promoter regions in α-proteobacteria. NfeR1 loss-of-function compromised osmoadaptation of free-living bacteria, whilst causing misregulation of salt-responsive genes related to stress adaptation, osmolytes catabolism and membrane trafficking. Nodulation tests revealed that lack of NfeR1 affected competitiveness, infectivity, nodule development and symbiotic efficiency of S. meliloti on alfalfa roots. Comparative computer predictions and a genetic reporter assay evidenced a redundant role of three identical NfeR1 unpaired anti Shine-Dalgarno motifs for targeting and downregulation of translation of multiple mRNAs from transporter genes. Our data provide genetic evidence of the hyperosmotic conditions of the endosymbiotic compartments. NfeR1-mediated gene regulation in response to this cue could contribute to coordinate nutrient uptake with the metabolic reprogramming concomitant to symbiotic transitions.

Regulation of Polyhydroxybutyrate Accumulation in Sinorhizobium meliloti by the Trans-Encoded Small RNA MmgR

Riboregulation has a major role in the fine-tuning of multiple bacterial processes. Among the RNA players, trans-encoded untranslated small RNAs (sRNAs) regulate complex metabolic networks by tuning expression from multiple target genes in response to numerous signals. In Sinorhizobium meliloti, over 400 sRNAs are expressed under different stimuli. The sRNA MmgR (standing for Makes more granules Regulator) has been of particular interest to us since its sequence and structure are highly conserved among the alphaproteobacteria and its expression is regulated by the amount and quality of the bacterium's available nitrogen source. In this work, we explored the biological role of MmgR in S. meliloti 2011 by characterizing the effect of a deletion of the internal conserved core of mmgR (mmgR^Δ33-51). This mutation resulted in larger amounts of polyhydroxybutyrate (PHB) distributed into more intracellular granules than are found in the wild-type strain. This phenotype was expressed upon cessation of balanced growth owing to nitrogen depletion in the presence of surplus carbon (i.e., at a carbon/nitrogen molar ratio greater than 10). The normal PHB accumulation was complemented with a wild-type mmgR copy but not with unrelated sRNA genes. Furthermore, the expression of mmgR limited PHB accumulation in the wild type, regardless of the magnitude of the C surplus. Quantitative proteomic profiling and quantitative reverse transcription-PCR (qRT-PCR) revealed that the absence of MmgR results in a posttranscriptional overexpression of both PHB phasin proteins (PhaP1 and PhaP2). Together, our results indicate that the widely conserved alphaproteobacterial MmgR sRNA fine-tunes the regulation of PHB storage in S. melilotiIMPORTANCE High-throughput RNA sequencing has recently uncovered an overwhelming number of trans-encoded small RNAs (sRNAs) in diverse prokaryotes. In the nitrogen-fixing alphaproteobacterial symbiont of alfalfa root nodules Sinorhizobium meliloti, only four out of hundreds of identified sRNA genes have been functionally characterized. Thus, uncovering the biological role of sRNAs currently represents a major issue and one that is particularly challenging because of the usually subtle quantitative regulation contributed by most characterized sRNAs. Here, we have characterized the function of the broadly conserved alphaproteobacterial sRNA gene mmgR in S. meliloti Our results strongly suggest that mmgR encodes a negative regulator of the accumulation of polyhydroxybutyrate, the major carbon and reducing power storage polymer in S. meliloti cells growing under conditions of C/N overbalance.

A statistical test for conserved RNA structure shows lack of evidence for structure in lncRNAs

Many functional RNAs have an evolutionarily conserved secondary structure. Conservation of RNA base pairing induces pairwise covariations in sequence alignments. We developed a computational method, R-scape (RNA Structural Covariation Above Phylogenetic Expectation), that quantitatively tests whether covariation analysis supports the presence of a conserved RNA secondary structure. R-scape analysis finds no statistically significant support for proposed secondary structures of the long noncoding RNAs HOTAIR, SRA, and Xist.

Regulation of fixLJ by Hfq Controls Symbiotically Important Genes in Sinorhizobium meliloti

The RNA-binding chaperone Hfq plays critical roles in the establishment and functionality of the symbiosis between Sinorhizobium meliloti and its legume hosts. A mutation in hfq reduces symbiotic efficiency resulting in a Fix^- phenotype, characterized by the inability of the bacterium to fix nitrogen. At least in part, this is due to the ability of Hfq to regulate the fixLJ operon, which encodes a sensor kinase-response regulator pair that controls expression of the nitrogenase genes. The ability of Hfq to bind fixLJ in vitro and in planta was demonstrated with gel shift and coimmunoprecipitation experiments. Two (ARN)₂ motifs in the fixLJ message were the likely sites through which Hfq exerted its posttranscriptional control. Consistent with the regulatory effects of Hfq, downstream genes controlled by FixLJ (such as nifK, noeB) were also subject to Hfq regulation in planta.

Thermodynamic matchers for the construction of the cuckoo RNA family

RNA family models describe classes of functionally related, non-coding RNAs based on sequence and structure conservation. The most important method for modeling RNA families is the use of covariance models, which are stochastic models that serve in the discovery of yet unknown, homologous RNAs. However, the performance of covariance models in finding remote homologs is poor for RNA families with high sequence conservation, while for families with high structure but low sequence conservation, these models are difficult to built in the first place. A complementary approach to RNA family modeling involves the use of thermodynamic matchers. Thermodynamic matchers are RNA folding programs, based on the established thermodynamic model, but tailored to a specific structural motif. As thermodynamic matchers focus on structure and folding energy, they unfold their potential in discovering homologs, when high structure conservation is paired with low sequence conservation. In contrast to covariance models, construction of thermodynamic matchers does not require an input alignment, but requires human design decisions and experimentation, and hence, model construction is more laborious. Here we report a case study on an RNA family that was constructed by means of thermodynamic matchers. It starts from a set of known but structurally different members of the same RNA family. The consensus secondary structure of this family consists of 2 to 4 adjacent hairpins. Each hairpin loop carries the same motif, CCUCCUCCC, while the stems show high variability in their nucleotide content. The present study describes (1) a novel approach for the integration of the structurally varying family into a single RNA family model by means of the thermodynamic matcher methodology, and (2) provides the results of homology searches that were conducted with this model in a wide spectrum of bacterial species.

The Spot 42 RNA: A regulatory small RNA with roles in the central metabolism

The Spot 42 RNA is a 109 nucleotide long (in Escherichia coli) noncoding small regulatory RNA (sRNA) encoded by the spf (spot fourty-two) gene. spf is found in gamma-proteobacteria and the majority of experimental work on Spot 42 RNA has been performed using E. coli, and recently Aliivibrio salmonicida. In the cell Spot 42 RNA plays essential roles as a regulator in carbohydrate metabolism and uptake, and its expression is activated by glucose, and inhibited by the cAMP-CRP complex. Here we summarize the current knowledge on Spot 42, and present the natural distribution of spf, show family-specific secondary structural features of Spot 42, and link highly conserved structural regions to mRNA target binding.

Transcriptomic Analysis of Sinorhizobium meliloti and Medicago truncatula Symbiosis Using Nitrogen Fixation-Deficient Nodules

The bacterium Sinorhizobium meliloti interacts symbiotically with legume plant hosts such as Medicago truncatula to form nitrogen-fixing root nodules. During symbiosis, plant and bacterial cells differentiate in a coordinated manner, resulting in specialized plant cells that contain nitrogen-fixing bacteroids. Both plant and bacterial genes are required at each developmental stage of symbiosis. We analyzed gene expression in nodules formed by wild-type bacteria on six plant mutants with defects in nitrogen fixation. We observed differential expression of 482 S. meliloti genes with functions in cell envelope homeostasis, cell division, stress response, energy metabolism, and nitrogen fixation. We simultaneously analyzed gene expression in M. truncatula and observed differential regulation of host processes that may trigger bacteroid differentiation and control bacterial infection. Our analyses of developmentally arrested plant mutants indicate that plants use distinct means to control bacterial infection during early and late symbiotic stages.

Small RNA Deep-Sequencing Analyses Reveal a New Regulator of Virulence in Agrobacterium fabrum C58

Novel ways of regulating Ti plasmid functions were investigated by studying small RNAs (sRNAs) that are known to act as posttranscriptional regulators in plant pathogenic bacteria. sRNA-seq analyses of Agrobacterium fabrum C58 allowed us to identify 1,108 small transcripts expressed in several growth conditions that could be sRNAs. A quarter of them were confirmed by bioinformatics or by biological experiments. Antisense RNAs represent 24% of the candidates and they are over-represented on the pTi (with 62% of pTi sRNAs), suggesting differences in the regulatory mechanisms between the essential and accessory replicons. Moreover, a large number of these pTi antisense RNAs are transcribed opposite to those genes involved in virulence. Others are 5'- and 3'-untranslated region RNAs and trans-encoded RNAs. We have validated, by rapid amplification of cDNA ends polymerase chain reaction, the transcription of 14 trans-encoded RNAs, among which RNA1111 is expressed from the pTiC58. Its deletion decreased the aggressiveness of A. fabrum C58 on tomatoes, tobaccos, and kalanchoe, suggesting that this sRNA activates virulence. The identification of its putative target mRNAs (6b gene, virC2, virD3, and traA) suggests that this sRNA may coordinate two of the major pTi functions, the infection of plants and its dissemination among bacteria.

A stress-induced small RNA modulates alpha-rhizobial cell cycle progression

Mechanisms adjusting replication initiation and cell cycle progression in response to environmental conditions are crucial for microbial survival. Functional characterization of the trans-encoded small non-coding RNA (trans-sRNA) EcpR1 in the plant-symbiotic alpha-proteobacterium Sinorhizobium meliloti revealed a role of this class of riboregulators in modulation of cell cycle regulation. EcpR1 is broadly conserved in at least five families of the Rhizobiales and is predicted to form a stable structure with two defined stem-loop domains. In S. meliloti, this trans-sRNA is encoded downstream of the divK-pleD operon. ecpR1 belongs to the stringent response regulon, and its expression was induced by various stress factors and in stationary phase. Induced EcpR1 overproduction led to cell elongation and increased DNA content, while deletion of ecpR1 resulted in reduced competitiveness. Computationally predicted EcpR1 targets were enriched with cell cycle-related mRNAs. Post-transcriptional repression of the cell cycle key regulatory genes gcrA and dnaA mediated by mRNA base-pairing with the strongly conserved loop 1 of EcpR1 was experimentally confirmed by two-plasmid differential gene expression assays and compensatory changes in sRNA and mRNA. Evidence is presented for EcpR1 promoting RNase E-dependent degradation of the dnaA mRNA. We propose that EcpR1 contributes to modulation of cell cycle regulation under detrimental conditions.

Microorganisms frequently encounter adverse conditions unfavorable for cell proliferation. They have evolved diverse mechanisms, including transcriptional control and targeted protein degradation, to adjust cell cycle progression in response to environmental cues. Non-coding RNAs are widespread regulators of various cellular processes in all domains of life. In prokaryotes, trans-encoded small non-coding RNAs (trans-sRNAs) contribute to a rapid cellular response to changing environments, but so far have not been directly related to cell cycle regulation. Here, we report the first example of a trans-sRNA (EcpR1) with two experimentally confirmed targets in the core of cell cycle regulation and demonstrate that in the plant-symbiotic alpha-proteobacterium Sinorhizobium meliloti the regulatory mechanism involves base-pairing of this sRNA with the dnaA and gcrA mRNAs. Most trans-sRNAs are restricted to closely related species, but the stress-induced EcpR1 is broadly conserved in the order of Rhizobiales suggesting an evolutionary advantage conferred by ecpR1. It broadens the functional diversity of prokaryotic sRNAs and adds a new regulatory level to the mechanisms that contribute to interlinking stress responses with the cell cycle machinery.

Small regulatory RNAs and the fine-tuning of plant–bacteria interactions

Small regulatory RNAs (sRNAs) play a key role in many physiological and adaptive responses in bacteria. Faced with rapidly changing environments, it is more advantageous for bacteria to use sRNA-mediated responses than regulation by protein transcriptional factors, as sRNAs act at the post-transcriptional level and require less energy and time for their synthesis and turnover. The use of RNA deep sequencing has provided hundreds of sRNA candidates in different bacterial species that interact with plants. Here, we review the most recent results for the involvement of bacterial sRNAs in beneficial as well as deleterious plant–bacteria interactions. We describe the current view for the role of sRNAs, which are suggested to improve competition for both niches and resources in plant-interacting bacteria. These sRNAs also help plant-associated bacteria individually adapt to the rapidly changing conditions to which they are exposed, during different stages of this interaction.

Sibling rivalry: related bacterial small RNAs and their redundant and non-redundant roles

Small RNA molecules (sRNAs) are now recognized as key regulators controlling bacterial gene expression, as sRNAs provide a quick and efficient means of positively or negatively altering the expression of specific genes. To date, numerous sRNAs have been identified and characterized in a myriad of bacterial species, but more recently, a theme in bacterial sRNAs has emerged: the presence of more than one highly related sRNAs produced by a given bacterium, here termed sibling sRNAs. Sibling sRNAs are those that are highly similar at the nucleotide level, and while it might be expected that sibling sRNAs exert identical regulatory functions on the expression of target genes based on their high degree of relatedness, emerging evidence is demonstrating that this is not always the case. Indeed, there are several examples of bacterial sibling sRNAs with non-redundant regulatory functions, but there are also instances of apparent regulatory redundancy between sibling sRNAs. This review provides a comprehensive overview of the current knowledge of bacterial sibling sRNAs, and also discusses important questions about the significance and evolutionary implications of this emerging class of regulators.

RNA-Seq analysis of the multipartite genome of Rhizobium etli CE3 shows different replicon contributions under heat and saline shock

Background

Regulation of transcription is essential for any organism and Rhizobium etli (a multi-replicon, nitrogen-fixing symbiotic bacterium) is no exception. This bacterium is commonly found in the rhizosphere (free-living) or inside of root-nodules of the common bean (Phaseolus vulgaris) in a symbiotic relationship. Abiotic stresses, such as high soil temperatures and salinity, compromise the genetic stability of R. etli and therefore its symbiotic interaction with P. vulgaris. However, it is still unclear which genes are up- or down-regulated to cope with these stress conditions. The aim of this study was to identify the genes and non-coding RNAs (ncRNAs) that are differentially expressed under heat and saline shock, as well as the promoter regions of the up-regulated loci.

Results

Analysing the heat and saline shock responses of R. etli CE3 through RNA-Seq, we identified 756 and 392 differentially expressed genes, respectively, and 106 were up-regulated under both conditions. Notably, the set of genes over-expressed under either condition was preferentially encoded on plasmids, although this observation was more significant for the heat shock response. In contrast, during either saline shock or heat shock, the down-regulated genes were principally chromosomally encoded. Our functional analysis shows that genes encoding chaperone proteins were up-regulated during the heat shock response, whereas genes involved in the metabolism of compatible solutes were up-regulated following saline shock. Furthermore, we identified thirteen and nine ncRNAs that were differentially expressed under heat and saline shock, respectively, as well as eleven ncRNAs that had not been previously identified. Finally, using an in silico analysis, we studied the promoter motifs in all of the non-coding regions associated with the genes and ncRNAs up-regulated under both conditions.

Conclusions

Our data suggest that the replicon contribution is different for different stress responses and that the heat shock response is more complex than the saline shock response. In general, this work exemplifies how strategies that not only consider differentially regulated genes but also regulatory elements of the stress response provide a more comprehensive view of bacterial gene regulation.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-770) contains supplementary material, which is available to authorized users.

Riboregulation in plant-associated α-proteobacteria

The symbiotic α-rhizobia Sinorhizobium meliloti, Bradyrhizobium japonicum, Rhizobium etli and the related plant pathogen Agrobacterium tumefaciens are important model organisms for studying plant-microbe interactions. These metabolically versatile soil bacteria are characterized by complex lifestyles and large genomes. Here we summarize the recent knowledge on their small non-coding RNAs (sRNAs) including conservation, function, and interaction of the sRNAs with the RNA chaperone Hfq. In each of these organisms, an inventory of hundreds of cis- and trans-encoded sRNAs with regulatory potential was uncovered by high-throughput approaches and used for the construction of 39 sRNA family models. Genome-wide analyses of hfq mutants and co-immunoprecipitation with tagged Hfq revealed a major impact of the RNA chaperone on the physiology of plant-associated α-proteobacteria including symbiosis and virulence. Highly conserved members of the SmelC411 family are the AbcR sRNAs, which predominantly regulate ABC transport systems. AbcR1 of A. tumefaciens controls the uptake of the plant-generated signaling molecule GABA and is a central regulator of nutrient uptake systems. It has similar functions in S. meliloti and the human pathogen Brucella abortus. As RNA degradation is an important process in RNA-based gene regulation, a short overview on ribonucleases in plant-associated α-proteobacteria concludes this review.

A multiobjective method for robust identification of bacterial small non-coding RNAs

MOTIVATION

Small non-coding RNAs (sRNAs) have major roles in the post-transcriptional regulation in prokaryotes. The experimental validation of a relatively small number of sRNAs in few species requires developing computational algorithms capable of robustly encoding the available knowledge and using this knowledge to predict sRNAs within and across species.

RESULTS

We present a novel methodology designed to identify bacterial sRNAs by incorporating the knowledge encoded by different sRNA prediction methods and optimally aggregating them as potential predictors. Because some of these methods emphasize specificity, whereas others emphasize sensitivity while detecting sRNAs, their optimal aggregation constitutes trade-off solutions between these two contradictory objectives that enhance their individual merits. Many non-redundant optimal aggregations uncovered by using multiobjective optimization techniques are then combined into a multiclassifier, which ensures robustness during detection and prediction even in genomes with distinct nucleotide composition. By training with sRNAs in Salmonella enterica Typhimurium, we were able to successfully predict sRNAs in Sinorhizobium meliloti, as well as in multiple and poorly annotated species. The proposed methodology, like a meta-analysis approach, may begin to lay a possible foundation for developing robust predictive methods across a wide spectrum of genomic variability.

AVAILABILITY AND IMPLEMENTATION

Scripts created for the experimentation are available at http://m4m.ugr.es/SupInfo/sRNAOS/sRNAOSscripts.zip.

CONTACT

delval@decsai.ugr.es

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Two separate modules of the conserved regulatory RNA AbcR1 address multiple target mRNAs in and outside of the translation initiation region

The small RNA AbcR1 regulates the expression of ABC transporters in the plant pathogen Agrobacterium tumefaciens, the plant symbiont Sinorhizobium meliloti, and the human pathogen Brucella abortus. A combination of proteomic and bioinformatic approaches suggested dozens of AbcR1 targets in A. tumefaciens. Several of these newly discovered targets are involved in the uptake of amino acids, their derivatives, and sugars. Among the latter is the periplasmic sugar-binding protein ChvE, a component of the virulence signal transduction system. We examined 16 targets and their interaction with AbcR1 in close detail. In addition to the previously described mRNA interaction site of AbcR1 (M1), the CopraRNA program predicted a second functional module (M2) as target-binding site. Both M1 and M2 contain single-stranded anti-SD motifs. Using mutated AbcR1 variants, we systematically tested by band shift experiments, which sRNA region is responsible for mRNA binding and gene regulation. On the target site, we find that AbcR1 interacts with some mRNAs in the translation initiation region and with others far into their coding sequence. Our data show that AbcR1 is a versatile master regulator of nutrient uptake systems in A. tumefaciens and related bacteria.

Bioinformatics of prokaryotic RNAs

The genome of most prokaryotes gives rise to surprisingly complex transcriptomes, comprising not only protein-coding mRNAs, often organized as operons, but also harbors dozens or even hundreds of highly structured small regulatory RNAs and unexpectedly large levels of anti-sense transcripts. Comprehensive surveys of prokaryotic transcriptomes and the need to characterize also their non-coding components is heavily dependent on computational methods and workflows, many of which have been developed or at least adapted specifically for the use with bacterial and archaeal data. This review provides an overview on the state-of-the-art of RNA bioinformatics focusing on applications to prokaryotes.

Genome-wide profiling of Hfq-binding RNAs uncovers extensive post-transcriptional rewiring of major stress response and symbiotic regulons in Sinorhizobium meliloti

The RNA chaperone Hfq is a global post-transcriptional regulator in bacteria. Here, we used RNAseq to analyze RNA populations from the legume symbiont Sinorhizobium meliloti that were co-immunoprecipitated (CoIP-RNA) with a FLAG-tagged Hfq in five growth/stress conditions. Hfq-bound transcripts (1315) were largely identified in stressed bacteria and derived from small RNAs (sRNAs), both trans-encoded (6.4%) and antisense (asRNAs; 6.3%), and mRNAs (86%). Pull-down with Hfq recovered a small proportion of annotated S. meliloti sRNAs (14% of trans-sRNAs and 2% of asRNAs) suggesting a discrete impact of this protein in sRNA pathways. Nonetheless, Hfq selectively stabilized CoIP-enriched sRNAs, anticipating that these interactions are functionally significant. Transcription of 26 Hfq-bound sRNAs was predicted to occur from promoters recognized by the major stress σ factors σ(E2) or σ(H1/2). Recovery rates of sRNAs in each of the CoIP-RNA libraries suggest a large impact of Hfq-assisted riboregulation in S. meliloti osmoadaptation. Hfq directly targeted 18% of the predicted S. meliloti mRNAs, which encode functionally diverse proteins involved in transport and metabolism, σ(E2)-dependent stress responses, quorum sensing, flagella biosynthesis, ribosome, and membrane assembly or symbiotic nitrogen fixation. Canonical targeting of the 5' regions of two of the ABC transporter mRNAs by the homologous Hfq-binding AbcR1 and AbcR2 sRNAs leading to inhibition of protein synthesis was confirmed in vivo. We therefore provide a comprehensive resource for the systems-level deciphering of hitherto unexplored S. meliloti stress and symbiotic post-transcriptional regulons and the identification of Hfq-dependent sRNA-mRNA regulatory pairs.

Small RNA sX13: a multifaceted regulator of virulence in the plant pathogen Xanthomonas

Small noncoding RNAs (sRNAs) are ubiquitous posttranscriptional regulators of gene expression. Using the model plant-pathogenic bacterium Xanthomonas campestris pv. vesicatoria (Xcv), we investigated the highly expressed and conserved sRNA sX13 in detail. Deletion of sX13 impinged on Xcv virulence and the expression of genes encoding components and substrates of the Hrp type III secretion (T3S) system. qRT-PCR analyses revealed that sX13 promotes mRNA accumulation of HrpX, a key regulator of the T3S system, whereas the mRNA level of the master regulator HrpG was unaffected. Complementation studies suggest that sX13 acts upstream of HrpG. Microarray analyses identified 63 sX13-regulated genes, which are involved in signal transduction, motility, transcriptional and posttranscriptional regulation and virulence. Structure analyses of in vitro transcribed sX13 revealed a structure with three stable stems and three apical C-rich loops. A computational search for putative regulatory motifs revealed that sX13-repressed mRNAs predominantly harbor G-rich motifs in proximity of translation start sites. Mutation of sX13 loops differentially affected Xcv virulence and the mRNA abundance of putative targets. Using a GFP-based reporter system, we demonstrated that sX13-mediated repression of protein synthesis requires both the C-rich motifs in sX13 and G-rich motifs in potential target mRNAs. Although the RNA-binding protein Hfq was dispensable for sX13 activity, the hfq mRNA and Hfq::GFP abundance were negatively regulated by sX13. In addition, we found that G-rich motifs in sX13-repressed mRNAs can serve as translational enhancers and are located at the ribosome-binding site in 5% of all protein-coding Xcv genes. Our study revealed that sX13 represents a novel class of virulence regulators and provides insights into sRNA-mediated modulation of adaptive processes in the plant pathogen Xanthomonas.

Since the discovery of the first regulatory RNA in 1981, hundreds of small RNAs (sRNAs) have been identified in bacteria. Although sRNA-mediated control of virulence was demonstrated for numerous animal- and human-pathogenic bacteria, sRNAs and their functions in plant-pathogenic bacteria have been enigmatic. We discovered that the sRNA sX13 is a novel virulence regulator of Xanthomonas campestris pv. vesicatoria (Xcv), which causes bacterial spot disease on pepper and tomato. sX13 contributes to the Xcv-plant interaction by promoting the synthesis of an essential pathogenicity factor of Xcv, i. e., the type III secretion system. Thus, in addition to transcriptional regulation, sRNA-mediated posttranscriptional regulation contributes to virulence of plant-pathogenic xanthomonads. To repress target mRNAs carrying G-rich motifs, sX13 employs C-rich loops. Hence, sX13 exhibits striking structural similarity to sRNAs in distantly related human pathogens, e. g., Staphylococcus aureus and Helicobacter pylori, suggesting that structure-driven target regulation via C-rich motifs represents a conserved feature of sRNA-mediated posttranscriptional regulation. Furthermore, sX13 is the first sRNA shown to control the mRNA level of hfq, which encodes a conserved RNA-binding protein required for sRNA activity and virulence in many enteric bacteria.

Independent activity of the homologous small regulatory RNAs AbcR1 and AbcR2 in the legume symbiont Sinorhizobium meliloti

The legume symbiont Sinorhizobium meliloti expresses a plethora of small noncoding RNAs (sRNAs) whose function is mostly unknown. Here, we have functionally characterized two tandemly encoded S. meliloti Rm1021 sRNAs that are similar in sequence and structure. Homologous sRNAs (designated AbcR1 and AbcR2) have been shown to regulate several ABC transporters in the related α-proteobacteria Agrobacterium tumefaciens and Brucella abortus. In Rm1021, AbcR1 and AbcR2 exhibit divergent unlinked regulation and are stabilized by the RNA chaperone Hfq. AbcR1 is transcribed in actively dividing bacteria, either in culture, rhizosphere or within the invasion zone of mature alfalfa nodules. Conversely, AbcR2 expression is induced upon entry into stationary phase and under abiotic stress. Only deletion of AbcR1 resulted into a discrete growth delay in rich medium, but both are dispensable for symbiosis. Periplasmic proteome profiling revealed down-regulation of the branched-chain amino acid binding protein LivK by AbcR1, but not by AbcR2. A double-plasmid reporter assay confirmed the predicted specific targeting of the 5′-untranslated region of the livK mRNA by AbcR1 in vivo. Our findings provide evidences of independent regulatory functions of these sRNAs, probably to fine-tune nutrient uptake in free-living and undifferentiated symbiotic rhizobia.

Computational identification of functional RNA homologs in metagenomic data

A key step toward understanding a metagenomics data set is the identification of functional sequence elements within it, such as protein coding genes and structural RNAs. Relative to protein coding genes, structural RNAs are more difficult to identify because of their reduced alphabet size, lack of open reading frames, and short length. Infernal is a software package that implements “covariance models” (CMs) for RNA homology search, which harness both sequence and structural conservation when searching for RNA homologs. Thanks to the added statistical signal inherent in the secondary structure conservation of many RNA families, Infernal is more powerful than sequence-only based methods such as BLAST and profile HMMs. Together with the Rfam database of CMs, Infernal is a useful tool for identifying RNAs in metagenomics data sets.

Next-generation annotation of prokaryotic genomes with EuGene-P: application to Sinorhizobium meliloti 2011

The availability of next-generation sequences of transcripts from prokaryotic organisms offers the opportunity to design a new generation of automated genome annotation tools not yet available for prokaryotes. In this work, we designed EuGene-P, the first integrative prokaryotic gene finder tool which combines a variety of high-throughput data, including oriented RNA-Seq data, directly into the prediction process. This enables the automated prediction of coding sequences (CDSs), untranslated regions, transcription start sites (TSSs) and non-coding RNA (ncRNA, sense and antisense) genes. EuGene-P was used to comprehensively and accurately annotate the genome of the nitrogen-fixing bacterium Sinorhizobium meliloti strain 2011, leading to the prediction of 6308 CDSs as well as 1876 ncRNAs. Among them, 1280 appeared as antisense to a CDS, which supports recent findings that antisense transcription activity is widespread in bacteria. Moreover, 4077 TSSs upstream of protein-coding or non-coding genes were precisely mapped providing valuable data for the study of promoter regions. By looking for RpoE2-binding sites upstream of annotated TSSs, we were able to extend the S. meliloti RpoE2 regulon by ∼3-fold. Altogether, these observations demonstrate the power of EuGene-P to produce a reliable and high-resolution automatic annotation of prokaryotic genomes.

Complete Genome Sequence of the Alfalfa Symbiont Sinorhizobium/Ensifer meliloti Strain GR4

We present the complete nucleotide sequence of the multipartite genome of Sinorhizobium/Ensifer meliloti GR4, a predominant rhizobial strain in an agricultural field site. The genome (total size, 7.14 Mb) consists of five replicons: one chromosome, two expected symbiotic megaplasmids (pRmeGR4c and pRmeGR4d), and two accessory plasmids (pRmeGR4a and pRmeGR4b).

Insights into the noncoding RNome of nitrogen-fixing endosymbiotic α-proteobacteria

Symbiotic chronic infection of legumes by rhizobia involves transition of invading bacteria from a free-living environment in soil to an intracellular state as differentiated nitrogen-fixing bacteroids within the nodules elicited in the host plant. The adaptive flexibility demanded by this complex lifestyle is likely facilitated by the large set of regulatory proteins encoded by rhizobial genomes. However, proteins are not the only relevant players in the regulation of gene expression in bacteria. Large-scale high-throughput analysis of prokaryotic genomes is evidencing the expression of an unexpected plethora of small untranslated transcripts (sRNAs) with housekeeping or regulatory roles. sRNAs mostly act in response to environmental cues as post-transcriptional regulators of gene expression through protein-assisted base-pairing interactions with target mRNAs. Riboregulation contributes to fine-tune a wide range of bacterial processes which, in intracellular animal pathogens, largely compromise virulence traits. Here, we summarize the incipient knowledge about the noncoding RNome structure of nitrogen-fixing endosymbiotic bacteria as inferred from genome-wide searches for sRNA genes in the alfalfa partner Sinorhizobium meliloti and further comparative genomics analysis. The biology of relevant S. meliloti RNA chaperones (e.g., Hfq) is also reviewed as a first global indicator of the impact of riboregulation in the establishment of the symbiotic interaction.

Quantitative proteomic analysis of the Hfq-regulon in Sinorhizobium meliloti 2011

Riboregulation stands for RNA-based control of gene expression. In bacteria, small non-coding RNAs (sRNAs) are a major class of riboregulatory elements, most of which act at the post-transcriptional level by base-pairing target mRNA genes. The RNA chaperone Hfq facilitates antisense interactions between target mRNAs and regulatory sRNAs, thus influencing mRNA stability and/or translation rate. In the α-proteobacterium Sinorhizobium meliloti strain 2011, the identification and detection of multiple sRNAs genes and the broadly pleitropic phenotype associated to the absence of a functional Hfq protein both support the existence of riboregulatory circuits controlling gene expression to ensure the fitness of this bacterium in both free living and symbiotic conditions. In order to identify target mRNAs subject to Hfq-dependent riboregulation, we have compared the proteome of an hfq mutant and the wild type S. meliloti by quantitative proteomics following protein labelling with ¹⁵N. Among 2139 univocally identified proteins, a total of 195 proteins showed a differential abundance between the Hfq mutant and the wild type strain; 65 proteins accumulated ≥2-fold whereas 130 were downregulated (≤0.5-fold) in the absence of Hfq. This profound proteomic impact implies a major role for Hfq on regulation of diverse physiological processes in S. meliloti, from transport of small molecules to homeostasis of iron and nitrogen. Changes in the cellular levels of proteins involved in transport of nucleotides, peptides and amino acids, and in iron homeostasis, were confirmed with phenotypic assays. These results represent the first quantitative proteomic analysis in S. meliloti. The comparative analysis of the hfq mutant proteome allowed identification of novel strongly Hfq-regulated genes in S. meliloti.

Computational analysis of noncoding RNAs

Noncoding RNAs have emerged as important key players in the cell. Understanding their surprisingly diverse range of functions is challenging for experimental and computational biology. Here, we review computational methods to analyze noncoding RNAs. The topics covered include basic and advanced techniques to predict RNA structures, annotation of noncoding RNAs in genomic data, mining RNA-seq data for novel transcripts and prediction of transcript structures, computational aspects of microRNAs, and database resources.

Dual RpoH sigma factors and transcriptional plasticity in a symbiotic bacterium

Sinorhizobium meliloti can live as a soil saprophyte and can engage in a nitrogen-fixing symbiosis with plant roots. To succeed in such diverse environments, the bacteria must continually adjust gene expression. Transcriptional plasticity in eubacteria is often mediated by alternative sigma (σ) factors interacting with core RNA polymerase. The S. meliloti genome encodes 14 of these alternative σ factors, including two putative RpoH ("heat shock") σ factors. We used custom Affymetrix symbiosis chips to characterize the global transcriptional response of S. meliloti rpoH1, rpoH2, and rpoH1 rpoH2 mutants during heat shock and stationary-phase growth. Under these conditions, expression of over 300 genes is dependent on rpoH1 and rpoH2. We mapped transcript start sites of 69 rpoH-dependent genes using 5' RACE (5' rapid amplification of cDNA ends), which allowed us to determine putative RpoH1-dependent, RpoH2-dependent, and dual-promoter (RpoH1- and RpoH2-dependent) consensus sequences that were each used to search the genome for other potential direct targets of RpoH. The inferred S. meliloti RpoH promoter consensus sequences share features of Escherichia coli RpoH promoters but lack extended -10 motifs.

Identification of two small regulatory RNAs linked to virulence in Brucella abortus 2308

Hfq is an RNA-binding protein that functions in post-transcriptional gene regulation by mediating interactions between mRNAs and small regulatory RNAs (sRNAs). Two proteins encoded by BAB1_1794 and BAB2_0612 are highly over-produced in a Brucella abortus hfq mutant compared with the parental strain, and recently, expression of orthologues of these proteins in Agrobacterium tumefaciens was shown to be regulated by two sRNAs, called AbcR1 and AbcR2. Orthologous sRNAs (likewise designated AbcR1 and AbcR2) have been identified in B. abortus 2308. In Brucella, abcR1 and abcR2 single mutants are not defective in their ability to survive in cultured murine macrophages, but an abcR1 abcR2 double mutant exhibits significant attenuation in macrophages. Additionally, the abcR1 abcR2 double mutant displays significant attenuation in a mouse model of chronic Brucella infection. Quantitative proteomics and microarray analyses revealed that the AbcR sRNAs predominantly regulate genes predicted to be involved in amino acid and polyamine transport and metabolism, and Northern blot analyses indicate that the AbcR sRNAs accelerate the degradation of the target mRNAs. In an Escherichia coli two-plasmid reporter system, overexpression of either AbcR1 or AbcR2 was sufficient for regulation of target mRNAs, indicating that the AbcR sRNAs from B. abortus 2308 perform redundant regulatory functions.