A genome-wide survey of sRNAs in the symbiotic nitrogen-fixing alpha-proteobacterium Sinorhizobium meliloti

A tripartite bacterial-fungal-plant symbiosis in the mycorrhiza-shaped microbiome drives plant growth and mycorrhization

BACKGROUND

Plant microbiomes play crucial roles in nutrient cycling and plant growth, and are shaped by a complex interplay between plants, microbes, and the environment. The role of bacteria as mediators of the 400-million-year-old partnership between the majority of land plants and, arbuscular mycorrhizal (AM) fungi is still poorly understood. Here, we test whether AM hyphae-associated bacteria influence the success of the AM symbiosis.

RESULTS

Using partitioned microcosms containing field soil, we discovered that AM hyphae and roots selectively assemble their own microbiome from the surrounding soil. In two independent experiments, we identified several bacterial genera, including Devosia, that are consistently enriched on AM hyphae. Subsequently, we isolated 144 pure bacterial isolates from a mycorrhiza-rich sample of extraradical hyphae and isolated Devosia sp. ZB163 as root and hyphal colonizer. We show that this AM-associated bacterium synergistically acts with mycorrhiza on the plant root to strongly promote plant growth, nitrogen uptake, and mycorrhization.

CONCLUSIONS

Our results highlight that AM fungi do not function in  isolation and that the plant-mycorrhiza symbiont can recruit beneficial bacteria that support the symbiosis. Video Abstract.

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.

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.

Reprograming of sRNA target specificity by the leader peptide peTrpL in response to antibiotic exposure

Trans-acting regulatory RNAs have the capacity to base pair with more mRNAs than generally detected under defined conditions, raising the possibility that sRNA target specificities vary depending on the specific metabolic or environmental conditions. In Sinorhizobium meliloti, the sRNA rnTrpL is derived from a tryptophan (Trp) transcription attenuator located upstream of the Trp biosynthesis gene trpE(G). The sRNA rnTrpL contains a small ORF, trpL, encoding the 14-aa leader peptide peTrpL. If Trp is available, efficient trpL translation causes transcription termination and liberation of rnTrpL, which subsequently acts to downregulate the trpDC operon, while peTrpL is known to have a Trp-independent role in posttranscriptional regulation of antibiotic resistance mechanisms. Here, we show that tetracycline (Tc) causes rnTrpL accumulation independently of Trp availability. In the presence of Tc, rnTrpL and peTrpL act collectively to destabilize rplUrpmA mRNA encoding ribosomal proteins L21 and L27. The three molecules, rnTrpL, peTrpL, and rplUrpmA mRNA, form an antibiotic-dependent ribonucleoprotein complex (ARNP). In vitro reconstitution of this ARNP in the presence of competing trpD and rplU transcripts revealed that peTrpL and Tc cause a shift of rnTrpL specificity towards rplU, suggesting that sRNA target prioritization may be readjusted in response to changing environmental conditions.

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.

Genome-wide identification of Azospirillum brasilense Sp245 small RNAs responsive to nitrogen starvation and likely involvement in plant-microbe interactions

BACKGROUND

Small RNAs (sRNAs) are non-coding RNAs known to regulate various biological functions such as stress adaptation, metabolism, virulence as well as pathogenicity across a wide range of bacteria, mainly by controlling mRNA stabilization or regulating translation. Identification and functional characterization of sRNAs has been carried out in various plant growth-promoting bacteria and they have been shown to help the cells cope up with environmental stress. No study has been carried out to uncover these regulatory molecules in the diazotrophic alpha-proteobacterium Azospirillum brasilense Sp245 to date.

RESULTS

Expression-based sRNA identification (RNA-seq) revealed the first list of ~ 468 sRNA candidate genes in A. brasilense Sp245 that were differentially expressed in nitrogen starvation versus non-starved conditions. In parallel, in silico tools also identified 2 of the above as candidate sRNAs. Altogether, putative candidates were stringently curated from RNA-seq data based on known sRNA parameters (size, location, secondary structure, and abundance). In total, ~ 59 significantly expressed sRNAs were identified in this study of which 53 are potentially novel sRNAs as they have no Rfam and BSRD homologs. Sixteen sRNAs were randomly selected and validated for differential expression, which largely was found to be in congruence with the RNA-seq data.

CONCLUSIONS

Differential expression of 468 A. brasilense sRNAs was indicated by RNA-seq data, a subset of which was confirmed by expression analysis. Four of the significantly expressed sRNAs were not observed in nitrogen starvation while 16 sRNAs were found to be exclusively expressed in nitrogen depletion. Putative candidate sRNAs identified have potential mRNA targets primarily involved in stress (abiotic and biotic) adaptability; regulation of bacterial cellular, biological and molecular pathways such as nitrogen fixation, polyhydroxybutyrate synthesis, chemotaxis, biofilm formation and transcriptional regulation. In addition to directly influencing bacteria, some of these sRNAs also have targets influencing plant-microbe interactions through adhesion of bacteria to plant roots directly, suppressing host response, inducing plant defence and signalling.

Whole-genome Identification of Transcriptional Start Sites by Differential RNA-seq in Bacteria

Gene transcription in bacteria often starts some nucleotides upstream of the start codon. Identifying the specific Transcriptional Start Site (TSS) is essential for genetic manipulation, as in many cases upstream of the start codon there are sequence elements that are involved in gene expression regulation. Taken into account the classical gene structure, we are able to identify two kinds of transcriptional start site: primary and secondary. A primary transcriptional start site is located some nucleotides upstream of the translational start site, while a secondary transcriptional start site is located within the gene encoding sequence. Here, we present a step by step protocol for genome-wide transcriptional start sites determination by differential RNA-sequencing (dRNA-seq) using the enteric pathogen Shigella flexneri serotype 5a strain M90T as model. However, this method can be employed in any other bacterial species of choice. In the first steps, total RNA is purified from bacterial cultures using the hot phenol method. Ribosomal RNA (rRNA) is specifically depleted via hybridization probes using a commercial kit. A 5'-monophosphate-dependent exonuclease (TEX)-treated RNA library enriched in primary transcripts is then prepared for comparison with a library that has not undergone TEX-treatment, followed by ligation of an RNA linker adaptor of known sequence allowing the determination of TSS with single nucleotide precision. Finally, the RNA is processed for Illumina sequencing library preparation and sequenced as purchased service. TSS are identified by in-house bioinformatic analysis. Our protocol is cost-effective as it minimizes the use of commercial kits and employs freely available software.

Complete genome sequence and annotation of the laboratory reference strain Shigella flexneri serotype 5a M90T and genome-wide transcriptional start site determination

Background

Shigella is a Gram-negative facultative intracellular bacterium that causes bacillary dysentery in humans. Shigella invades cells of the colonic mucosa owing to its virulence plasmid-encoded Type 3 Secretion System (T3SS), and multiplies in the target cell cytosol. Although the laboratory reference strain S. flexneri serotype 5a M90T has been extensively used to understand the molecular mechanisms of pathogenesis, its complete genome sequence is not available, thereby greatly limiting studies employing high-throughput sequencing and systems biology approaches.

Results

We have sequenced, assembled, annotated and manually curated the full genome of S. flexneri 5a M90T. This yielded two complete circular contigs, the chromosome and the virulence plasmid (pWR100). To obtain the genome sequence, we have employed long-read PacBio DNA sequencing followed by polishing with Illumina RNA-seq data. This provides a new hybrid strategy to prepare gapless, highly accurate genome sequences, which also cover AT-rich tracks or repetitive sequences that are transcribed. Furthermore, we have performed genome-wide analysis of transcriptional start sites (TSS) and determined the length of 5′ untranslated regions (5′-UTRs) at typical culture conditions for the inoculum of in vitro infection experiments. We identified 6723 primary TSS (pTSS) and 7328 secondary TSS (sTSS). The S. flexneri 5a M90T annotated genome sequence and the transcriptional start sites are integrated into RegulonDB (http://regulondb.ccg.unam.mx) and RSAT (http://embnet.ccg.unam.mx/rsat/) databases to use their analysis tools in the S. flexneri 5a M90T genome.

Conclusions

We provide the first complete genome for S. flexneri serotype 5a, specifically the laboratory reference strain M90T. Our work opens the possibility of employing S. flexneri M90T in high-quality systems biology studies such as transcriptomic and differential expression analyses or in genome evolution studies. Moreover, the catalogue of TSS that we report here can be used in molecular pathogenesis studies as a resource to know which genes are transcribed before infection of host cells. The genome sequence, together with the analysis of transcriptional start sites, is also a valuable tool for precise genetic manipulation of S. flexneri 5a M90T. Further, we present a new hybrid strategy to prepare gapless, highly accurate genome sequences. Unlike currently used hybrid strategies combining long- and short-read DNA sequencing technologies to maximize accuracy, our workflow using long-read DNA sequencing and short-read RNA sequencing provides the added value of using non-redundant technologies, which yield distinct, exploitable datasets.

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.

Transcription attenuation-derived small RNA rnTrpL regulates tryptophan biosynthesis gene expression in trans

Ribosome-mediated transcription attenuation is a basic posttranscriptional regulation mechanism in bacteria. Liberated attenuator RNAs arising in this process are generally considered nonfunctional. In Sinorhizobium meliloti, the tryptophan (Trp) biosynthesis genes are organized into three operons, trpE(G), ppiD-trpDC-moaC-moeA, and trpFBA-accD-folC, of which only the first one, trpE(G), contains a short ORF (trpL) in the 5′-UTR and is regulated by transcription attenuation. Under conditions of Trp sufficiency, transcription is terminated between trpL and trpE(G), and a small attenuator RNA, rnTrpL, is produced. Here, we show that rnTrpL base-pairs with trpD and destabilizes the polycistronic trpDC mRNA, indicating rnTrpL-mediated downregulation of the trpDC operon in trans. Although all three trp operons are regulated in response to Trp availability, only in the two operons trpE(G) and trpDC the Trp-mediated regulation is controlled by rnTrpL. Together, our data show that the trp attenuator coordinates trpE(G) and trpDC expression posttranscriptionally by two fundamentally different mechanisms: ribosome-mediated transcription attenuation in cis and base-pairing in trans. Also, we present evidence that rnTrpL-mediated regulation of trpDC genes expression in trans is conserved in Agrobacterium and Bradyrhizobium, suggesting that the small attenuator RNAs may have additional conserved functions in the control of bacterial gene expression.

Characterization of novel small RNAs (sRNAs) contributing to the desiccation response of Salmonella enterica serovar Typhimurium

Noncoding RNA (ncRNA) modulation of gene expression has now been ubiquitously observed across all domains of life. An increasingly apparent role of ncRNAs is to coordinate changes in gene expressions in response to environmental stress. Salmonella enterica, a common food-born pathogen, is known for its striking ability to survive, adapt, and thrive in various unfavourable environments which makes it a particularly difficult pathogen to eliminate as well as an interesting model in which to study ncRNA contributions to cellular stress response. Mounting evidence now suggests that small RNAs (sRNAs) represent key regulators of Salmonella stress adaptation. Approximately 50-500 nucleotides in length, sRNAs regulate gene expression through complementary base pairing with molecular targets and have recently been suggested to outnumber protein-coding genes in bacteria. In this work, we employ small RNA transcriptome sequencing to characterize changes in the sRNA profiles of Salmonella in response to desiccation. In all, we identify 102 previously annotated sRNAs significantly differentially expressed during desiccation; and excitingly, 71 novel sRNAs likewise differentially expressed. Small transcript northern blotting and qRT-PCRs confirm the identities and expressions of several of our novel sRNAs, and computational analyses indicate the majority are highly conserved and structurally related to characterized sRNAs. Predicted sRNA targets include several proteins necessary for desiccation survival and this, in part, suggests a role for desiccation-regulated sRNAs in this stress response. Furthermore, we find individual knock-outs of two of the novel sRNAs identified herein, either sRNA1320429 or sRNA3981754, significantly impairs the ability of Salmonella to survive desiccation, confirming their involvements (and suggesting the potential involvements of other sRNAs we identify in this work) in the Salmonella response to desiccation.

Whole genome sequence of moderate halophilic marine bacterium Marinobacter litoralis SW-45: Abundance of non-coding RNAs

A report on the de novo Whole Genome Sequence (WGS) of Marinobacter litoralis SW-45, a moderately salt-tolerant bacterium isolated from the seawater in Malaysia is presented. The strain has a genome size of 3.45 Mb and is capable of producing halophilic lipase, protease and esterase enzymes. Computational prediction of non-coding RNA (ncRNA) genes in M. litoralis SW-45 was performed using standalone software known as the non-coding RNA characterization (nocoRNAc). In addition, a phylogenetic tree showing the evolutionary relationship between the strain and other members of the genus Marinobacter was constructed using 16SrRNA sequence information. A total of 385 ncRNA transcripts, 1124 terminator region, and 2350 Stress Induced Duplex Destabilization sites were predicted. The current WGS shotgun project has provided the relevant genetic information that may be useful for the strain's improvement studies. This manuscript gives the first description of M. litoralis with a complete genome.

Highly conserved nucleotide motifs present in the 5'UTR of the heme-receptor gene shmR are required for HmuP-dependent expression of shmR in Ensifer meliloti

Heme may represent a major iron-source for bacteria. In the symbiotic nitrogen-fixing bacterium Ensifer meliloti 1021, iron acquisition from heme depends on the outer-membrane heme-receptor ShmR. Expression of shmR gene is repressed by iron in a RirA dependent manner while under iron-limitation its expression requires the small protein HmuP. In this work, we identified highly conserved nucleotide motifs present upstream the shmR gene. These motifs are widely distributed among Alpha and Beta Proteobacteria, and correlate with the presence of HmuP coding sequences in bacterial genomes. According to data presented in this work, we named these new motifs as HmuP-responsive elements (HPREs). In the analyzed genomes, the HPREs were always present upstream of genes encoding putative heme-receptors. Moreover, in those Alpha and Beta Proteobacteria where transcriptional start sites for shmR homologs are known, HPREs were located in the 5'UTR region. In this work we show that in E. meliloti 1021, HPREs are involved in HmuP-dependent shmR expression. Moreover, we show that changes in sequence composition of the HPREs correlate with changes in a predicted RNA secondary structure element and affect shmR gene expression.

Multidisciplinary approaches for studying rhizobium–legume symbioses

The rhizobium–legume symbiosis is a major source of fixed nitrogen (ammonia) in the biosphere. The potential for this process to increase agricultural yield while reducing the reliance on nitrogen-based fertilizers has generated interest in understanding and manipulating this process. For decades, rhizobium research has benefited from the use of leading techniques from a very broad set of fields, including population genetics, molecular genetics, genomics, and systems biology. In this review, we summarize many of the research strategies that have been employed in the study of rhizobia and the unique knowledge gained from these diverse tools, with a focus on genome- and systems-level approaches. We then describe ongoing synthetic biology approaches aimed at improving existing symbioses or engineering completely new symbiotic interactions. The review concludes with our perspective of the future directions and challenges of the field, with an emphasis on how the application of a multidisciplinary approach and the development of new methods will be necessary to ensure successful biotechnological manipulation of the symbiosis.

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.

Widespread Antisense Transcription in Prokaryotes

Although bacterial genomes are usually densely protein-coding, genome-wide mapping approaches of transcriptional start sites revealed that a significant fraction of the identified promoters drive the transcription of noncoding RNAs. These can be trans-acting RNAs, mainly originating from intergenic regions and, in many studied examples, possessing regulatory functions. However, a significant fraction of these noncoding RNAs consist of natural antisense transcripts (asRNAs), which overlap other transcriptional units. Naturally occurring asRNAs were first observed to play a role in bacterial plasmid replication and in bacteriophage λ more than 30 years ago. Today's view is that asRNAs abound in all three domains of life. There are several examples of asRNAs in bacteria with clearly defined functions. Nevertheless, many asRNAs appear to result from pervasive initiation of transcription, and some data point toward global functions of such widespread transcriptional activity, explaining why the search for a specific regulatory role is sometimes futile. In this review, we give an overview about the occurrence of antisense transcription in bacteria, highlight particular examples of functionally characterized asRNAs, and discuss recent evidence pointing at global relevance in RNA processing and transcription-coupled DNA repair.

Pervasive Transcription of Mitochondrial, Plastid, and Nucleomorph Genomes across Diverse Plastid-Bearing Species

Organelle genomes exhibit remarkable diversity in content, structure, and size, and in their modes of gene expression, which are governed by both organelle- and nuclear-encoded machinery. Next generation sequencing (NGS) has generated unprecedented amounts of genomic and transcriptomic data, which can be used to investigate organelle genome transcription. However, most of the available eukaryotic RNA-sequencing (RNA-seq) data are used to study nuclear transcription only, even though large numbers of organelle-derived reads can typically be mined from these experiments. Here, we use publicly available RNA-seq data to assess organelle genome transcription in 59 diverse plastid-bearing species. Our RNA mapping analyses unraveled pervasive (full or near-full) transcription of mitochondrial, plastid, and nucleomorph genomes. In all cases, 85% or more of the organelle genome was recovered from the RNA data, including noncoding (intergenic and intronic) regions. These results reinforce the idea that organelles transcribe all or nearly all of their genomic material and are dependent on post-transcriptional processing of polycistronic transcripts. We explore the possibility that transcribed intergenic regions are producing functional noncoding RNAs, and that organelle genome noncoding content might provide raw material for generating regulatory RNAs.

Functional Genomics Approaches to Studying Symbioses between Legumes and Nitrogen-Fixing Rhizobia

Biological nitrogen fixation gives legumes a pronounced growth advantage in nitrogen-deprived soils and is of considerable ecological and economic interest. In exchange for reduced atmospheric nitrogen, typically given to the plant in the form of amides or ureides, the legume provides nitrogen-fixing rhizobia with nutrients and highly specialised root structures called nodules. To elucidate the molecular basis underlying physiological adaptations on a genome-wide scale, functional genomics approaches, such as transcriptomics, proteomics, and metabolomics, have been used. This review presents an overview of the different functional genomics approaches that have been performed on rhizobial symbiosis, with a focus on studies investigating the molecular mechanisms used by the bacterial partner to interact with the legume. While rhizobia belonging to the alpha-proteobacterial group (alpha-rhizobia) have been well studied, few studies to date have investigated this process in beta-proteobacteria (beta-rhizobia).

Regulation Mediated by N-Acyl Homoserine Lactone Quorum Sensing Signals in the Rhizobium-Legume Symbiosis

Soil-dwelling bacteria collectively referred to as rhizobia synthesize and perceive N-acyl-homoserine lactone (AHL) signals to regulate gene expression in a population density-dependent manner. AHL-mediated signaling in these bacteria regulates several functions which are important for the establishment of nitrogen-fixing symbiosis with legume plants. Moreover, rhizobial AHL act as interkingdom signals triggering plant responses that impact the plant-bacteria interaction. Both the regulatory mechanisms that control AHL synthesis in rhizobia and the set of bacterial genes and associated traits under quorum sensing (QS) control vary greatly among the rhizobial species. In this article, we focus on the well-known QS system of the alfalfa symbiont Sinorhizobium(Ensifer)meliloti. Bacterial genes, environmental factors and transcriptional and posttranscriptional regulatory mechanisms that control AHL production in this Rhizobium, as well as the effects of the signaling molecule on bacterial phenotypes and plant responses will be reviewed. Current knowledge of S. meliloti QS will be compared with that of other rhizobia. Finally, participation of the legume host in QS by interfering with rhizobial AHL perception through the production of molecular mimics will also be addressed.

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.

An sRNA and Cold Shock Protein Homolog-Based Feedforward Loop Post-transcriptionally Controls Cell Cycle Master Regulator CtrA

Adjustment of cell cycle progression is crucial for bacterial survival and adaptation under adverse conditions. However, the understanding of modulation of cell cycle control in response to environmental changes is rather incomplete. In α-proteobacteria, the broadly conserved cell cycle master regulator CtrA underlies multiple levels of control, including coupling of cell cycle and cell differentiation. CtrA levels are known to be tightly controlled through diverse transcriptional and post-translational mechanisms. Here, small RNA (sRNA)-mediated post-transcriptional regulation is uncovered as an additional level of CtrA fine-tuning. Computational predictions as well as transcriptome and proteome studies consistently suggested targeting of ctrA and the putative cold shock chaperone cspA5 mRNAs by the trans-encoded sRNA (trans-sRNA) GspR (formerly SmelC775) in several Sinorhizobium species. GspR strongly accumulated in the stationary growth phase, especially in minimal medium (MM) cultures. Lack of the gspR locus confers a fitness disadvantage in competition with the wild type, while its overproduction hampers cell growth, suggesting that this riboregulator interferes with cell cycle progression. An eGFP-based reporter in vivo assay, involving wild-type and mutant sRNA and mRNA pairs, experimentally confirmed GspR-dependent post-transcriptional down-regulation of ctrA and cspA5 expression, which most likely occurs through base-pairing to the respective mRNA. The energetically favored secondary structure of GspR is predicted to comprise three stem-loop domains, with stem-loop 1 and stem-loop 3 targeting ctrA and cspA5 mRNA, respectively. Moreover, this work reports evidence for post-transcriptional control of ctrA by CspA5. Thus, this regulation and GspR-mediated post-transcriptional repression of ctrA and cspA5 expression constitute a coherent feed-forward loop, which may enhance the negative effect of GspR on CtrA levels. This novel regulatory circuit involving the riboregulator GspR, CtrA, and a cold shock chaperone may contribute to fine-tuning of ctrA expression.

Transcriptomic Studies of the Effect of nod Gene-Inducing Molecules in Rhizobia: Different Weapons, One Purpose

Simultaneous quantification of transcripts of the whole bacterial genome allows the analysis of the global transcriptional response under changing conditions. RNA-seq and microarrays are the most used techniques to measure these transcriptomic changes, and both complement each other in transcriptome profiling. In this review, we exhaustively compiled the symbiosis-related transcriptomic reports (microarrays and RNA sequencing) carried out hitherto in rhizobia. This review is specially focused on transcriptomic changes that takes place when five rhizobial species, Bradyrhizobium japonicum (=diazoefficiens) USDA 110, Rhizobium leguminosarum biovar viciae 3841, Rhizobium tropici CIAT 899, Sinorhizobium (=Ensifer) meliloti 1021 and S. fredii HH103, recognize inducing flavonoids, plant-exuded phenolic compounds that activate the biosynthesis and export of Nod factors (NF) in all analysed rhizobia. Interestingly, our global transcriptomic comparison also indicates that each rhizobial species possesses its own arsenal of molecular weapons accompanying the set of NF in order to establish a successful interaction with host legumes.

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.

Pervasive, Genome-Wide Transcription in the Organelle Genomes of Diverse Plastid-Bearing Protists

Organelle genomes are among the most sequenced kinds of chromosome. This is largely because they are small and widely used in molecular studies, but also because next-generation sequencing technologies made sequencing easier, faster, and cheaper. However, studies of organelle RNA have not kept pace with those of DNA, despite huge amounts of freely available eukaryotic RNA-sequencing (RNA-seq) data. Little is known about organelle transcription in nonmodel species, and most of the available eukaryotic RNA-seq data have not been mined for organelle transcripts. Here, we use publicly available RNA-seq experiments to investigate organelle transcription in 30 diverse plastid-bearing protists with varying organelle genomic architectures. Mapping RNA-seq data to organelle genomes revealed pervasive, genome-wide transcription, regardless of the taxonomic grouping, gene organization, or noncoding content. For every species analyzed, transcripts covered ≥85% of the mitochondrial and/or plastid genomes (all of which were ≤105 kb), indicating that most of the organelle DNA—coding and noncoding—is transcriptionally active. These results follow earlier studies of model species showing that organellar transcription is coupled and ubiquitous across the genome, requiring significant downstream processing of polycistronic transcripts. Our findings suggest that noncoding organelle DNA can be transcriptionally active, raising questions about the underlying function of these transcripts and underscoring the utility of publicly available RNA-seq data for recovering complete genome sequences. If pervasive transcription is also found in bigger organelle genomes (>105 kb) and across a broader range of eukaryotes, this could indicate that noncoding organelle RNAs are regulating fundamental processes within eukaryotic cells.

Sinorhizobium meliloti YbeY is an endoribonuclease with unprecedented catalytic features, acting as silencing enzyme in riboregulation

Structural and biochemical features suggest that the almost ubiquitous bacterial YbeY protein may serve catalytic and/or Hfq-like protective functions central to small RNA (sRNA)-mediated regulation and RNA metabolism. We have biochemically and genetically characterized the YbeY ortholog of the legume symbiont Sinorhizobium meliloti (SmYbeY). Co-immunoprecipitation (CoIP) with a FLAG-tagged SmYbeY yielded a poor enrichment in RNA species, compared to Hfq CoIP-RNA uncovered previously by a similar experimental setup. Purified SmYbeY behaved as a monomer that indistinctly cleaved single- and double-stranded RNA substrates, a unique ability among bacterial endoribonucleases. SmYbeY-mediated catalysis was supported by the divalent metal ions Mg²⁺, Mn²⁺ and Ca²⁺, which influenced in a different manner cleavage efficiency and reactivity patterns, with Ca²⁺ specifically blocking activity on double-stranded and some structured RNA molecules. SmYbeY loss-of-function compromised expression of core energy and RNA metabolism genes, whilst promoting accumulation of motility, late symbiotic and transport mRNAs. Some of the latter transcripts are known Hfq-binding sRNA targets and might be SmYbeY substrates. Genetic reporter and in vitro assays confirmed that SmYbeY is required for sRNA-mediated down-regulation of the amino acid ABC transporter prbA mRNA. We have thus discovered a bacterial endoribonuclease with unprecedented catalytic features, acting also as gene silencing enzyme.

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 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.

RNase E and RNase J are needed for S-adenosylmethionine homeostasis in Sinorhizobium meliloti

The ribonucleases (RNases) E and J play major roles in E. coli and Bacillus subtilis, respectively, and co-exist in Sinorhizobium meliloti. We analysed S. meliloti 2011 mutants with mini-Tn5 insertions in the corresponding genes rne and rnj and found many overlapping effects. We observed similar changes in mRNA levels, including lower mRNA levels of the motility and chemotaxis related genes flaA, flgB and cheR and higher levels of ndvA (important for glucan export). The acyl-homoserine lactone (AHL) levels were also higher during exponential growth in both RNase mutants, despite no increase in the expression of the sinI AHL synthase gene. Furthermore, several RNAs from both mutants migrated aberrantly in denaturing gels at 300 V but not under stronger denaturing conditions at 1300 V. The similarities between the two mutants could be explained by increased levels of the key methyl donor S-adenosylmethionine (SAM), since this may result in faster AHL synthesis leading to higher AHL accumulation as well as in uncontrolled methylation of macromolecules including RNA, which may strengthen RNA secondary structures. Indeed, we found that in both mutants the N6-methyladenosine content was increased almost threefold and the SAM level was increased at least sevenfold. Complementation by induced ectopic expression of the respective RNase restored the AHL and SAM levels in each of the mutants. In summary, our data show that both RNase E and RNase J are needed for SAM homeostasis in S. meliloti.

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.

Genome-wide identification and characterization of small RNAs in Rhodobacter capsulatus and identification of small RNAs affected by loss of the response regulator CtrA

Small non-coding RNAs (sRNAs) are involved in the control of numerous cellular processes through various regulatory mechanisms, and in the past decade many studies have identified sRNAs in a multitude of bacterial species using RNA sequencing (RNA-seq). Here, we present the first genome-wide analysis of sRNA sequencing data in Rhodobacter capsulatus, a purple nonsulfur photosynthetic alphaproteobacterium. Using a recently developed bioinformatics approach, sRNA-Detect, we detected 422 putative sRNAs from R. capsulatus RNA-seq data. Based on their sequence similarity to sRNAs in a sRNA collection, consisting of published putative sRNAs from 23 additional bacterial species, and RNA databases, the sequences of 124 putative sRNAs were conserved in at least one other bacterial species; and, 19 putative sRNAs were assigned a predicted function. We bioinformatically characterized all putative sRNAs and applied machine learning approaches to calculate the probability of a nucleotide sequence to be a bona fide sRNA. The resulting quantitative model was able to correctly classify 95.2% of sequences in a validation set. We found that putative cis-targets for antisense and partially overlapping sRNAs were enriched with protein-coding genes involved in primary metabolic processes, photosynthesis, compound binding, and with genes forming part of macromolecular complexes. We performed differential expression analysis to compare the wild type strain to a mutant lacking the response regulator CtrA, an important regulator of gene expression in R. capsulatus, and identified 18 putative sRNAs with differing levels in the two strains. Finally, we validated the existence and expression patterns of four novel sRNAs by Northern blot analysis.

Massive parallel insertion site sequencing of an arrayed Sinorhizobium meliloti signature-tagged mini-Tn 5 transposon mutant library

Transposon mutagenesis in conjunction with identification of genomic transposon insertion sites is a powerful tool for gene function studies. We have implemented a protocol for parallel determination of transposon insertion sites by Illumina sequencing involving a hierarchical barcoding method that allowed for tracking back insertion sites to individual clones of an arrayed signature-tagged transposon mutant library. This protocol was applied to further characterize a signature-tagged mini-Tn 5 mutant library comprising about 12,000 mutants of the symbiotic nitrogen-fixing alphaproteobacterium Sinorhizobium meliloti (Pobigaylo et al., 2006; Appl. Environ. Microbiol. 72, 4329-4337). Previously, insertion sites have been determined for 5000 mutants of this library. Combining an adapter-free, inverse PCR method for sequencing library preparation with next generation sequencing, we identified 4473 novel insertion sites, increasing the total number of transposon mutants with known insertion site to 9562. The number of protein-coding genes that were hit at least once by a transposon increased by 1231 to a total number of 3673 disrupted genes, which represents 59% of the predicted protein-coding genes in S. meliloti.

A putative 3-hydroxyisobutyryl-CoA hydrolase is required for efficient symbiotic nitrogen fixation in Sinorhizobium meliloti and Sinorhizobium fredii NGR234

We report that the smb20752 gene of the alfalfa symbiont Sinorhizobium meliloti is a novel symbiotic gene required for full N₂ -fixation. Deletion of smb20752 resulted in lower nitrogenase activity and smaller nodules without impacting overall nodule morphology. Orthologs of smb20752 were present in all alpha and beta rhizobia, including the ngr_b20860 gene of Sinorhizobium fredii NGR234. A ngr_b20860 mutant formed Fix^- determinate nodules that developed normally to a late stage of the symbiosis on the host plants Macroptilium atropurpureum and Vigna unguiculata. However an early symbiotic defect was evident during symbiosis with Leucaena leucocephala, producing Fix^- indeterminate nodules. The smb20752 and ngr_b20860 genes encode putative 3-hydroxyisobutyryl-CoA (HIB-CoA) hydrolases. HIB-CoA hydrolases are required for l-valine catabolism and appear to prevent the accumulation of toxic metabolic intermediates, particularly methacrylyl-CoA. Evidence presented here and elsewhere (Curson et al., , PLoS ONE 9:e97660) demonstrated that Smb20752 and NGR_b20860 can also prevent metabolic toxicity, are required for l-valine metabolism, and play an undefined role in 3-hydroxybutyrate catabolism. We present evidence that the symbiotic defect of the HIB-CoA hydrolase mutants is independent of the inability to catabolize l-valine and suggest it relates to the toxicity resulting from metabolism of other compounds possibly related to 3-hydroxybutyric acid.

Expression of the Sinorhizobium meliloti small RNA gene mmgR is controlled by the nitrogen source

Small non-coding regulatory RNAs (sRNAs) are key players in post-transcriptional regulation of gene expression. Hundreds of sRNAs have been identified in Sinorhizobium meliloti, but their biological function remains unknown for most of them. In this study, we characterized the expression pattern of the gene encoding the 77-nt sRNA MmgR in S. meliloti strain 2011. A chromosomal transcriptional reporter fusion (PmmgR-gfp) showed that the mmgR promoter is active along different stages of the interaction with alfalfa roots. In pure cultures, PmmgR-gfp activity paralleled the sRNA abundance indicating that mmgR expression is primarily controlled at the level of transcriptional initiation. PmmgR-gfp activity was higher during growth in rhizobial defined medium (RDM) than in TY medium. Furthermore, PmmgR-gfp was induced at 60 min after shifting growing cells from TY to RDM medium, i.e. shorter than the cell doubling time. In defined RDM medium containing NO3 (-), both PmmgR-gfp and MmgR level were repressed by the addition of tryptone or single amino acids, suggesting that mmgR expression depends on the cellular nitrogen (N) status. In silico analysis failed to detect conserved motifs upstream the promoter RNA polymerase binding site, but revealed a strongly conserved motif centered at -28 that may be linked to the observed regulatory pattern by the N source.

Bacterial antisense RNAs are mainly the product of transcriptional noise

cis-Encoded antisense RNAs (asRNAs) are widespread along bacterial transcriptomes. However, the role of most of these RNAs remains unknown, and there is an ongoing discussion as to what extent these transcripts are the result of transcriptional noise. We show, by comparative transcriptomics of 20 bacterial species and one chloroplast, that the number of asRNAs is exponentially dependent on the genomic AT content and that expression of asRNA at low levels exerts little impact in terms of energy consumption. A transcription model simulating mRNA and asRNA production indicates that the asRNA regulatory effect is only observed above certain expression thresholds, substantially higher than physiological transcript levels. These predictions were verified experimentally by overexpressing nine different asRNAs in Mycoplasma pneumoniae. Our results suggest that most of the antisense transcripts found in bacteria are the consequence of transcriptional noise, arising at spurious promoters throughout the genome.

Isolation and characterization of active promoters from Gluconacetobacter diazotrophicus strain PAL5 using a promoter-trapping plasmid

Gluconacetobacter diazotrophicus is a nitrogen-fixing, endophytic bacterium that has the potential to promote plant growth and increase yield. Genetically modified strains might get more benefits to host plants, including through expression of useful proteins, such as Cry toxins from B. thuringiensis, or enzymes involved in phytohormone production, proteins with antagonistic activity for phytopathogens, or that improve nutrient utilization by the plant. For that, expression systems for G. diazotrophicus are needed, which requires active promoters fused to foreign (or innate) genes. This article describes the construction of a G. diazotrophicus PAL5 promoter library using a promoter-less lacZ-bearing vector, and the identification of six active promoters through β-galactosidase activity assays, sequencing and localization in the bacterial genome. The characterized promoters, which are located on distinct regions of the bacterial genome and encoding either sense or antisense transcripts, present variable expression strengths and might be used in the future for expressing useful proteins.

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 stress-related, rhizobial small RNA RcsR1 destabilizes the autoinducer synthase encoding mRNA sinI in Sinorhizobium meliloti

Quorum sensing is a cell density-dependent communication system of bacteria relying on autoinducer molecules. During the analysis of the post-transcriptional regulation of quorum sensing in the nitrogen fixing plant symbiont Sinorhizobium meliloti, we predicted and verified a direct interaction between the 5'-UTR of sinI mRNA encoding the autoinducer synthase and a small RNA (sRNA), which we named RcsR1. In vitro, RcsR1 prevented cleavage in the 5'-UTR of sinI by RNase E and impaired sinI translation. In line with low ribosomal occupancy and transcript destabilization upon binding of RcsR1 to sinI, overproduction of RcsR1 in S. meliloti resulted in lower level and shorter half-life of sinI mRNA, and in decreased autoinducer amount. Although RcsR1 can influence quorum sensing via sinI, its level did not vary at different cell densities, but decreased under salt stress and increased at low temperature. We found that RcsR1 and its stress-related expression pattern, but not the interaction with sinI homologs, are conserved in Sinorhizobium, Rhizobium and Agrobacterium. Consistently, overproduction of RcsR1 in S. meliloti and Agrobacterium tumefaciens inhibited growth at high salinity. We identified conserved targets of RcsR1 and showed that most conserved interactions and the effect on growth under salt stress are mediated by the first stem-loop of RcsR1, while its central part is responsible for the species-specific interaction with sinI. We conclude that RcsR1 is an ancient, stress-related riboregulator in rhizobia and propose that it links stress responses to quorum sensing in S. meliloti.

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.

RNA regulators of host immunity and pathogen adaptive responses in the oral cavity

The recent explosion of RNA-seq studies has resulted in a newfound appreciation for the importance of riboregulatory RNAs in the posttranscriptional control of eukaryotic and prokaryotic genetic networks. The current review will explore the role of trans-riboregulatory RNAs in various adaptive responses of host and pathogen in the oral cavity.

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.

Whole-genome mapping of 5' RNA ends in bacteria by tagged sequencing: a comprehensive view in Enterococcus faecalis

Enterococcus faecalis is the third cause of nosocomial infections. To obtain the first snapshot of transcriptional organizations in this bacterium, we used a modified RNA-seq approach enabling to discriminate primary from processed 5' RNA ends. We also validated our approach by confirming known features in Escherichia coli. We mapped 559 transcription start sites (TSSs) and 352 processing sites (PSSs) in E. faecalis. A blind motif search retrieved canonical features of SigA- and SigN-dependent promoters preceding transcription start sites mapped. We discovered 85 novel putative regulatory RNAs, small- and antisense RNAs, and 72 transcriptional antisense organizations. Presented data constitute a significant insight into bacterial RNA landscapes and a step toward the inference of regulatory processes at transcriptional and post-transcriptional levels in a comprehensive manner.

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.

Differential RNA-seq: the approach behind and the biological insight gained

RNA-sequencing has revolutionized the quantitative and qualitative analysis of transcriptomes in both prokaryotes and eukaryotes. It provides a generic approach for gene expression profiling, annotation of transcript boundaries and operons, as well as identifying novel transcripts including small noncoding RNA molecules and antisense RNAs. We recently developed a differential RNA-seq (dRNA-seq) method which in addition to the above, yields information as to whether a given RNA is a primary or processed transcript. Originally applied to describe the primary transcriptome of the gastric pathogen Helicobacter pylori, dRNA-seq has since provided global maps of transcriptional start sites in diverse species, informed new biology in the CRISPR-Cas9 system, advanced to a tool for comparative transcriptomics, and inspired simultaneous RNA-seq of pathogen and host.

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.