Original Research ARTICLE
A Genome-Wide Prediction and Identification of Intergenic Small RNAs by Comparative Analysis in Mesorhizobium huakuii 7653R
- State Key Laboratory of Agricultural Microbiology, Huazhong Agricultural University, Wuhan, China
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.
Small non-coding RNAs (sRNAs) in bacteria are usually 50–500 nucleotides (nts) in length and fulfill their riboregulations together with Hfq (an RNA chaperone) and ribonucleases (RNases) or RNase-containing complexes to mediate the expression and stability of target mRNAs (Repoila and Darfeuille, 2009; Richards and Vanderpool, 2011; Hoe et al., 2013; Morita et al., 2015). With advances in bio-computational and experimental methods, increasing numbers of sRNAs have been found and identified in bacteria (Sridhar and Gunasekaran, 2013; Becker et al., 2014; Stubben et al., 2014; Tsai et al., 2015).
Most sRNAs studies have been carried out in Gram-negative bacteria and have mainly focused on the survey and characterization of trans-encoded RNAs (del Val et al., 2012; Robledo et al., 2015b). sRNAs contain at least two conserved stem-loop structures (Man et al., 2011; Richards and Vanderpool, 2011; Storz et al., 2011; Gottesman and Storz, 2015). Depending on their specific stem-loop structures, sRNAs can set up partial or complete base-pairing with target mRNAs to modulate mRNA posttranscriptionally or can act by sequestering protein (Backofen and Hess, 2010; Bobrovskyy et al., 2015; Schu et al., 2015; Wagner and Romby, 2015; Wang et al., 2015). Alternative stem-loop structures allow sRNAs to combine with different effect factors, such as coenzymes, purines, amino acids and other compounds, and which then cause responses to changes in cell's physiology or environment (Wilderman et al., 2004; Guillier and Gottesman, 2008; Beisel and Storz, 2010; Storz et al., 2011). Moreover, in most trans-encoded sRNAs, the basically consensus Hfq-binding signatures and anti-Shine-Dalgarno sequences were inferred from covariance models in each sRNA family (Monteiro et al., 2012; Sobrero and Valverde, 2012; Wilms et al., 2012a). sRNAs regulate diverse processes, including carbon metabolism, quorum sensing, cell division, virulence, iron uptake, oxidative stress, heat shock and antibiotic resistance so on (Valverde et al., 2008; Berghoff et al., 2009; Repoila and Darfeuille, 2009; Stubben et al., 2014; Wagner and Romby, 2015; Baumgardt et al., 2016). Recently, A nodule formation efficiency sRNA NfeR1 from S. meliloti was reported to effect on osmoadaptation and symbiotic efficiency of rhizobia (Robledo et al., 2017). Moreover, a trans-sRNA EcpR1 was revealed to be broadly conserved in Rhizobiales and to contribute to the modulation of cell cycle regulation under detrimental conditions (Robledo et al., 2015a). MmgR is another trans-encoded small RNA, and highly conserved among the a-proteobacteria. MmgR regulates the cellular Polyhydroxybutyrate accumulation and is controlled by the cellular nitrogen status in S. meliloti (Ceizel Borella et al., 2016; Lagares et al., 2016, 2017). In cyanobacteria, a nitrogen stress-induced RNA (NsiR4) was shown to be involved in the regulation of glutamine synthetase, a key enzyme required for biological nitrogen assimilation (Klahn et al., 2015).
In recent decades, many bio-computational methods have been developed to predict bacterial sRNAs. The RNA sequence homology, the thermodynamically favorable secondary structure, and the conserved and consensus secondary structures are the initial parameters that are commonly adopted to predict bacterial sRNAs in intergenic genomic regions (IGRs), which can be scanned by using software such as QRNA, RNALfold and RNAz (del Val et al., 2007; Livny, 2007; Livny and Waldor, 2007). Subsequently, more complex combined approaches have been used to predict sRNAs by comparative genomics. More recently, transcriptional signs-based methods that include the prediction of transcription factor binding sites, promoters and terminator signals have changed the focus toward the predicting bacterial genomic transcription units (Chang et al., 2010; Pellin et al., 2012; Sridhar and Gunasekaran, 2013; Su et al., 2014; Tsai et al., 2015).
In rhizobia, Sinorhizobium meliloti is the first strain in which sRNAs in IGRs were screened by comparative genomic sequences from eight related alpha-Proteobacteria using the programs eQRNA and RNAz programs as predictive tools. Eight of the original 32 candidates were confirmed to express small transcripts by using Northern blotting experiments (del Val et al., 2007). Subsequent comprehensive genome-wide screening and identification of sRNAs were carried out in a few species of Rhizobials, including Sinorhizobium meliloti, Rhizobium etli and Bradyrhizobium japonicum (Ulve et al., 2007; Valverde et al., 2008; Voss et al., 2009; Schluter et al., 2010; Vercruysse et al., 2010; Madhugiri et al., 2012; Becker et al., 2014; Lopez-Leal et al., 2014; Hahn et al., 2016). Mesorhizobium huakuii is a Gram-negative bacterium that belongs to the Rhizobials of alpha-Proteobacteria. M. huakuii interacts with its specific host plant, Astragalus sinicus L., and performs nitrogen-fixation by forming indeterminate-type nodules. Recently, the sequencing of the entire genome of M. huakuii 7653R was completed. Its genome was found to be composed of a chromosome (6,364,365 bp), and two megaplasmids, pMhu7653Ra (193,835 bp) and pMhu7653Rb (323,475 bp), with 7,205 protein-coding genes (Wang et al., 2014). Subsequently, the transcriptomes of M. huakuii 7653R in bacteroids and free-living cells were analyzed and compared using RNA-seq and microarrays. However, due to the low sequencing depth, the detected genes were limited to mRNAs (≥200 bp)(Peng et al., 2014). In this study, we predicted the existence of novel sRNA genes in the IGRs of M. huakuii 7653R via genome-wide comparisons with four related mesorhizobial strains, including M. huakuii bv. loti MAFF303099, M. ciceri bv. biserrulae WSM1271, M. australicum WSM2073 and M. opportunistum WSM2075. The transcription units of the predicted sRNAs were further analyzed basing on the high-throughput deep sequencing data of the Mesorhizobium huakuii 7653R global transcriptome, and the expression profiles of the nine selected sRNAs under diverse stress conditions were revealed using Northern blotting.
Prediction of Potential sRNAs in the IGRs of M. huakuii 7653R
The IGRs of M. huakuii 7653R with a length ≥50 nt are 5125, and these were compared with four other Mesorhizobial strains by WU-BLASTN to identify conserved candidate sequences in the IGRs. The resulting homologies of IGRs revealed by these genome-wide comparisons were further analyzed. Raw data from ~1,500 sequences were detected after WU-BLAST analysis. Each of the conserved intergenic regions in these candidate sequences was scanned individually using eQRNA and RNAz, and overlapping sequences were determined. Meanwhile, the promoters and Rho-independent terminators were predicted using Promoter 2.0, RNAMotif and Erpin, and the resulting information was used to further assess the sequences mentioned above (Tables S1, S2). Finally, 40 conserved sequences were identified as potential sRNAs (Table 1).
Table 1. The candidate sRNAs in the IGRs of M. huakuii 7653R and the annotated functional sRNAs deduced by comparing the predicted sRNAs of 7653R with annotated sRNAs in the Rfam database.
To find the annotated known sRNAs in 7653R, we performed Rfam searches using the above potential sRNAs. Of the 40 predicted sRNAs, 24 were annotated in the Rfam database, thus, the other 16 sRNAs are regarded as putative novel sRNAs (Table 1). Several known sRNAs, which are highly conserved in bacteria, were successfully predicted, including 6S RNA (MH_32), 4.5S (MH_s40), RNase P class A (MH_s20), and tmRNA (MH_s19), suggesting that the methods used in this study worked well. A few glycine, cobalamin and TPP riboswitches were also detected. Notably, most of the sRNAs (35 of 40) were located on the chromosome, with only a few (three in pMhu7653Ra and two in pMhu7653Rb) located on megaplasmids, indicating possible differential conservation of the IGR sequences between the chromosome and megaplasmids among the tested bacteria. Predicting the sRNAs that locate on the symbiotic megaplasmids seems to be difficult.
Predicted sRNA Transcription Units
As noted above, the genomic comparison identified 40 potential sRNAs. Further analyses were conducted to describe their transcription units and transcriptional start regions using Illumina high-throughput sequencing data. These analyses were performed on RNA isolated from 28- and 50-dpi (days post-infection) symbiotic nodules and free-living cells of 7653R (see Materials and Methods) that were treated with TAP. The depths of RNA-seq were shown in Tables S3, S4.
A comparison of the analysis between the predicted sRNAs and RNA-seq products showed that, apart from 5 riboswitches, 19 functional annotated sRNAs and 15 of 16 putative novel sRNAs were detected in the sequencing dataset, and only one was not expressed in plant nodules or in free-living cultured cells (Table S5). Of the riboswitches (cis-regulators of mRNA that control the translation of mRNA), MH_s14, was observed to be highly expressed in all sequenced samples. In total, other than the 5 riboswitches, 34 of the 35 predicted sRNAs had corresponding RNA-seq sequencing products, and only MH_s6 was absent in the sequenced samples. These results suggest that MH_s6 expression might be induced by different stress conditions rather than the logarithmic growth state or symbiosis or that the MH_s6 transcript levels were below the threshold of detection using RNA-seq. Notably, except for a few sRNAs (MH_s1, MH_8, MH_s10, and MH_s20), the lengths of the sequenced sRNAs in the RNA-seq data were slightly longer (10 of 34 sRNAs) or shorter (20 of 34 sRNAs) than the predicted lengths. The predicted transcriptional start and stop regions of the sRNAs were also compared with the RNA-seq data, and the results showed that 12 and 5 predicted sRNAs had approximate matches to their sequenced start and stop positions, respectively. Thus, more than 50% of the predicted sRNAs showed transcript length discrepancies with their RNA-seq products. This result has also been observed by others (Gonzalez et al., 2008; Valverde et al., 2008; Otto et al., 2012; Mentz et al., 2013), and the main reason may be ascribed to RNA processing or degradation in the cell (Madhugiri et al., 2012; Saramago et al., 2014).
Annotation and Distribution of Known Conserved sRNAs
There are four sRNAs, including 6S RNA, 4.5S RNA, bacterial RNase P class A, and tmRNA, that are highly conserved and regarded as housekeeping sRNAs in bacteria. In this study, these four housekeeping sRNAs were detected by both bioinformatic analysis and RNA-seq, and their expression profiles and expression levels are respectively shown in Figure 1 and Table 2. In agreement with the earlier reports, RNase P (MH_S20), an omnipresent endoribonuclease, was highly expressed under both free-living and symbiotic conditions and displayed the highest expression levels compared with other housekeeping sRNAs (Regalia et al., 2002; Vercruysse et al., 2010). The 4.5S RNA (MH_S40), which forms the signal recognition particle (SRP) with the Ffh protein, was also relatively highly expressed, especially in symbiotic bacteroids (~1.5-fold higher than in free-living cells). The 6S sRNA (MH_S32) has been reported to be more abundant during the stationary phase of cell growth. Its function is to block δ70-dependent transcription in alpha-Proteobacteria. In our sequencing results, 6S was strongly expressed in exponential-phase cells (OD600 = 0.8), and abundant transcripts were also detected in mature (28 dpi) and senescent (50 dpi) nodules. The expression of tmRNA (MH_S19), which rescues stalled ribosomes and tags incomplete polypeptides for degradation, was lower in both free-living and symbiotic conditions, and down-regulated in mature and senescent nodules compared with exponential-phase cells (>2-fold).
Figure 1. Expression profiles of four known conserved small RNAs. RNAseq raw read coverage traces of each sRNA are shown and were derived from the RNAseq libraries of free-living cells (red), 28-dpi (green) nodules, and 50-dpi nodules (blue). The coordinates for the x- and y-axes represent the genome alignment and coverage, respectively. (A) Expression profile of MH_s19 (tm RNA); (B) Expression profile of MH_s20 (RNase P RNA); (C) Expression profile of MH_s32 (6S RNA); (D) Expression profile of MH_s40 (4.5S RNA). In these four known conserved sRNAs, RNase P showed the highest expression levels under both free-living and symbiotic conditions.
Table 2. The differentiation and fold-changes of known conserved sRNA expression levels under different growth conditions based on the M. huakuii 7653R RNA-seq data.
ctRNA (counter-transcribed RNA) is a highly conserved small RNA located in repB–repC intergenic regions. It regulates the replication and incompatibility of repABC-type plasmids (MacLellan et al., 2005; Vercruysse et al., 2010; Rivera-Urbalejo et al., 2015). M. huakuii 7653R possesses two symbiotic megaplasmids required for symbiosis and nitrogen fixation. In our results, the levels of ctRNA (MH_s1 in pMhu7653Ra and MH_s4 in pMhu7653Rb) were very high in free-living cells, approximately 3- to 4-fold higher in 28 dpi and 50 dpi nodules, respectively (Figure 2). Moreover, a cis-encoded highly abundant sequencing product was observed in the repB–repC intergenic regions of pMhu7653Ra (located at 2601–2727) and pMhu7653Rb (located at 2344–2474), respectively, so we speculated that this cis-encoded product could be formed by ctRNA terminating the transcription of repC. The results further showed that ctRNA is trans-transcribed and highly expressed in rhizobia (Yip et al., 2015).
Figure 2. Expression profiles of ctRNA and its opposite strands in pMhu7653Ra (MH_s1) and pMhu7653Rb (MH_s4). The levels of ctRNA were very high in free-living cells in both pMhu7653Ra and pMhu7653Rb. (A) Expression profile of MH_s1 (+); (B) Expression profile of MH_s1 (−); (C) Expression profile of MH_s4 (+); (D) Expression profile of MH_s4 (−). A cis-encoded, highly-abundant sequencing product was observed in the opposite strands of ctRNA, which was located at 2,601–2,727 in pMhu7653Ra and 2,344–2,474 in pMhu7653Rb.
Experimental Verification of Nine Selected sRNAs by Northern Blotting
To confirm the existence of the predicted novel sRNAs and to understand their basic functions, 9 of the 15 potential novel sRNAs were selected randomly and detected by conducting Northern blotting on the total RNA extracted from cultured cells grown under non-stressed and stressed conditions (Figure 3), and 5S RNA was used as a positive control (Figure S1). As shown, the sizes of the nine sRNAs that were tested roughly consistent with those of the RNA sequencing dataset. Two candidate sRNAs (MH_s10 and MH_s22) exhibited two bands, and MH_s25 showed a complex banding pattern. Six of the nine tested sRNAs (MH_3, MH_s10, MH_s15, MH_s22, MH_s25, and MH_s36) displayed strong signals in response to a variety of external stimuli. MH_s7, MH_s11 and MH_s39 were expressed under only rare growth conditions. In addition, after analyzing the sequencing products again, another highly expressed sequence (MH_s25) was detected that located between chromosomes 4531972 and 4532239. Because this transcript product lacked typical conserved structures and 4532144–4532330 is close to the predicted region (4532142–4532343), we regard MH_s25 (4532144–4532330) as a potential sRNAs. Apart from 4532144–4532330 and 4532142–4532343, the further information on other Northern blotting signals of MH_s25 could not be detected in the RNA-seq data. The main reason may be that the bio-informational analysis and compilation were performed to deal with an enormous number of non-coding RNAs in the RNA-seq data, which led to some shorter or less abundant transcripts usually be merged into the longer ones in a same transcription unit, so the sequence information for complex banding pattern appearing in Northern blots was difficult to obtain.
Figure 3. Expression profiles (A–C), northern blots (D–F), and transcription directions (G–I) of nine selected small RNAs. The expression profiles of nine selected sRNAs from the RNAseq raw reads of free-living cells, 28-dpi nodules and 50-dpi nodules are shown on the top of each column. The details are described in Figure 1. Northern blot results are shown in the middle of each column. Lanes 1–9 represent different tested conditions, which include symbiosis (lane 1, 28 dpi nodules), micro-oxygen stress (lane 2, >99% N2), oxidative stress (lane 3, 2 mM H2O2), salt stress (lane 4, 4 M NaCl), acid stress (lane 5, pH 5), alkali stress (lane 6, pH 9), cold (lane 7, 20°C), heat stress (lane 8, 37°C) and the unstressed condition (lane 9, log phase cells) (see Section Materials and Methods). Below the read coverage traces, small RNA gene regions and transcription directions are indicated by arrows of different lengths and colors. Gray arrows represent flanking genes and white arrows represent the sRNA genes.
The expression and abundance of sRNAs were affected by various external stimuli, and these condition-specific expression patterns help to infer the functions of different sRNAs (Vercruysse et al., 2010). In our study, most of the tested sRNAs (MH_3, MH_s10, MH_s15, MH_s22, MH_s25, and MH_s36) were highly abundant in symbiotic, micro-oxygen, oxidative and salt-stress conditions, but few sRNAs accumulated under heat, cold, acid, or alkali stress conditions. MH_s15 and MH_s25 were widely expressed sRNAs and accumulated heavily under most of the tested induction conditions, suggesting that they function in multiple cellar processes. The expression patterns of MH_s7, MH_s11 and MH_s39 were unique, and these sRNAs accumulated under only a few of the tested conditions. Nearly all the detected sRNAs (8 of 9, MH_s11 being the exception) were expressed under symbiotic conditions or symbiosis-mimicking (micro-oxygen and H2O2) stimuli, indicating that a large number of sRNAs are required during rhizobial symbiotic nitrogen fixation. However, more studies are required to characterize the functions and mechanisms of action of these sRNAs.
Conservation and Structural Characterization of Nine Identified sRNAs
To further examine the characteristics of the detected sRNAs, a BLASTN search (NCBI) was performed to analyze the conservation of all the detected sRNAs in M. huakuii 7653R. The results showed that most of the detected sRNAs are highly conserved and only distributed in Mesorhizobium spp. and Bradyrhizobium japonicum, except for MH_s11, which has homologs in Bartonella spp. The detected sRNAs contain highly conserved secondary structures and more than two loops (Figure 4). MH_s15 has highly conserved structures, with a single, clear band shown on the Northern blots, and abundant polyA/U in the loop and tail of its secondary structure, which is the putative Hfq-binding site. In addition, the transcription levels of the genes upstream and downstream of MH_s15 were upregulated in 28 and 50 dpi nodules, respectively, compared with free-living cells (from 1.5- to 7-fold), according to the RNA-seq results. Therefore, the target genes of MH_s15 were analyzed by using the IntaRNA software and the potential conserved interactive regions were shown in Table S6, but further experiments are needed to characterize the interactions between MH_s15 and its target genes and determine the influence of MH_s15 on the symbiosis between M. huakuii 7653R and host.
Figure 4. The secondary structure prediction of nine small RNAs determined by using RNA-fold. The secondary structures are colored according to the base-pairing probabilities. The unpaired regions are colored according to the probabilities of being unpaired. The lowest minimum free energy (MFE) of each sRNA, analyzed using RNA-fold, is shown in the diagram. (A) Secondary structure of MH_s3; (B) Secondary structure of MH_s7; (C) Secondary structure of MH_s10; (D) Secondary structure of MH_s11; (E) Secondary structure of MH_s15; (F) Secondary structure of MH_s22; (G) Secondary structure of MH_s25; (H) Secondary structure of MH_s36; (I) Secondary structure of MH_s39.
Target Genes Prediction of Candidate sRNAs and Functional Analysis
In order to better understand the potential function of the 16 new sRNA candidates identified in this work, the target genes of 16 putative novel sRNAs were predicted by using the webserver IntaRNA. These candidate target genes were subjective to GO annotation and enrichment analysis. The WEGO output of the target genes was shown in Figures S2, S3. The enriched GO terms and functional categorization of candidate small RNA target genes were summarized and demonstrated in a pie chart (Figure 5, Table S8). As it showed in Figure 5A, classified in terms of cellular component, a majority of the predicted target genes is enriched in the function of membrane systems (40%), cell part (34%) and intracellular part (15%). When classified as molecular functions, a total of 38% of the enriched GO terms are related to catalytic activity, and 21% of the target genes are involved in the binding of nucleic acid, nucleoside and nucleotide (Figure 5B). Regarding the enriched GO term of biological process, as shown in Figure 5C, the target genes were distributed widely in various biological functions, from cellular process, to biological regulation, to primary metabolic process, to nitrogen compound metabolic process, and to response to stimulus and stress. Importantly to point out, 48% of the target genes are enriched in all kinds of metabolic processes, such as primary metabolism, second metabolism, nitrogen compound metabolism, alcohol metabolism etc. Remarkably, 9% of the target genes are involved in the nitrogen compound metabolic process, indicating the important roles of the sRNAs playing in the nitrogen fixing and utilization. Moreover, 16 and 17% of the target genes are respectively involved in biological regulation and cellular process, and 2% of the target genes are connected with response to stimulus and stress (Figure 5C).
Figure 5. GO analyses and functional classifications of predicted target genes of 16 putative novel sRNA in Mesorhizobium huakuii 7653R. The target genes of sRNA candidates were predicted by using the program IntaRNA, the GO term annotations were finished by the tool of InterPro, and the program WEGO was used to classify and plot the GO annotations. The enriched GO terms and functional categorization of candidate sRNA target genes were summarized. It showed that they are involved in cellular components (A), molecular functions (B), and biological processes (C). The target genes contained in each functional category are indicated as percentages (in brackets) of the total number of genes with GO annotations.
Genome-wide screens for sRNAs have mainly used bioinformatic prediction, RNA-seq, and tiling array analyses. In general, the number of non-coding transcripts is large, reaching hundreds or thousands of species in an RNA-seq or tiling array experiment (Vercruysse et al., 2010; Otto et al., 2012; Wilms et al., 2012b; Ignatov et al., 2013; Jimenez-Zurdo and Robledo, 2015). Hence, effective and accurate sRNA analysis technology is necessary to extract potential sRNA sequences from the complicated data. In this study, because most trans-encoded sRNAs are not followed by a Rho-independent terminator (Toffano-Nioche et al., 2012; Mentz et al., 2013), we focused on the conserved intergenic regions and stable secondary structures to screen for putative novel sRNAs in the IGRs of the M. huakuii 7653R genome. Based on the results of a comparative analysis of five related strains, we extracted 40 potential sRNAs and used RNA deep sequencing data to further identify their characteristics. Using this strategy, 24 annotated sRNAs with well-characterized functions in bacteria and 15 putative novel sRNAs were detected.
Except for five cis-regulatory riboswitches, 34 of the 35 predicted sRNAs in M. huakuii 7653R were consistent with the RNA-seq data. thus, the technology used for predicting sRNAs based on structural conservation is a reliable means of revealing sRNAs in bacteria. However, the method has some flaws. Apart from the generation of a small number of false positives (Sridhar and Gunasekaran, 2013), there were obvious size differences between the predicted sRNAs and the RNA-seq products which was also observed in previous articles (Tsai et al., 2015). This result may be related to the processing or degradation of sRNAs and the corresponding RNA-seq products (Arraiano et al., 2010; Vercruysse et al., 2010; Saramago et al., 2014). More importantly, the lengths of the sequenced non-coding RNAs showed some differences between free-living cells and plant nodules, where the non-coding RNAs in free-living cells were usually slightly longer than those in nodules. This result indicates that the sizes of sRNAs are influenced by the growth conditions. Moreover, the conservation of M. huakuii 7653R IGR sequences might result in some false judgments about sRNA transcription units. M. huakuii 7653R is a special case in rhizobia and forms nodules only in A. sinicus L. Most of the widespread sRNAs in rhizobia, e.g., Smr45C, Smr15C1/C2(AbcR1/2) and Smr7C so on, have not been discovered in 7653R. Hence, secondary structure conservation was adopted as a main predictive parameter to screen for the putative sRNAs in this study.
In our experiments, the diverse stresses were designed to stimulate the transcription of the sRNAs. Reactive oxygen (H2O2, ) is found to have a key function in a plant's defense system against pathogens. In rhizobia, H2O2 has not been detected in free-living cells or bacteroids, but H2O2 accumulation was observed in the infection stages of the symbiotic interaction and in ultrathin sections of mature 6-week-old nodules. For this reason, H2O2 played an important signaling role in early symbiotic and senescence processes, and rhizobia have an efficient antioxidant defense system to detoxify H2O2 (Herouart et al., 2002; Pauly et al., 2006; Puppo et al., 2013; Janczarek et al., 2015). During symbiotic nitrogen fixation, anoxic stress or a low oxygen concentration is required for nitrogenase activity and the expression of nitrogen-fixing genes in bacteroids. Two oxygen-responsive regulatory systems, FixLJ-FixK and RegSR-NifA, control numerous symbiotic genes in rhizobia. When the oxygen concentration was lowered from 21 to 0.5%, the number of genes controlled by the FixLJ-FixK2 regulatory cassette was shown to increase progressively. Moreover, NifA-dependent genes were expressed only when the oxygen concentrations was below 2% in the gas phase. Therefore, micro-oxygen and H2O2 stimuli were applied to the free-living cells to mimicked symbiotic conditions (Sciotti et al., 2003; Dixon and Kahn, 2004; Lindemann et al., 2007; Mesa et al., 2008). By Northern blotting, we further confirmed the transcription and sizes of the nine predicted sRNAs, and found that the expression of the majority of the detected sRNAs could be induced by symbiotic or symbiosis-mimicking conditions. The expression of sRNAs is usually driven by specific environmental conditions, and in turn, the function of sRNAs is relevant to the stimulus or stress (Caswell et al., 2012; Robledo et al., 2015a). Hence, the present results indicate that a large number of sRNAs are necessary to regulate the expression of genes involved in the differentiation and maturation of bacteroids during rhizobial symbiotic nitrogen fixation. Further, we predicted the target genes of 16 functional-unknown sRNAs including those 9 validated sRNAs, and sorted out the target genes through GO terms analysis by utilizing WEGO, by which a lot of useful information to figure out the potential functions of those novel sRNA candidates is provided. It was shown that a majority of the targets are involved in various metabolic or cellular biological processes, indicating that these sRNAs play fundamental roles in bacterial primary metabolism, enzyme catalytic activity, transport, replication and transcription process etc. The significant roles of the sRNAs playing in rhizobia infection, bacteroid differentiation, nitrogen metabolism and response to stresses are revealed as well in the enriched GO terms. These observations are not strange, since these biological processes actually reflect the life cycle of rhizobia in soil and host plant. For rhizobia, to carry out symbiotic nitrogen fixation in nodule, a series of signaling transduction and mutual molecular interactions between rhizobia and host plant are required. During rhizobia symbiosis and nitrogen fixation process, the size and morphology of rhizobia change significantly, which involves many cellular parts and cellular processes. Also, many transcription factors play essential roles in regulating nodule formation and rhizobia nitrogen fixation. Rhizobia have to survive within the infected host plant cells, during which it has to response to the various stress environments. Besides, symbiotic nitrogen fixation involves many transportation events regarding carbon and energy supplying and nitrogen compound assimilation.
So far, little is known about the role of sRNAs and riboregulation in the control of symbiotic plant-microbe interactions, and the only studies that have been perfo rmed have focused on S. meliloti (del Val et al., 2012; Becker et al., 2014; Jimenez-Zurdo and Robledo, 2015). The homologous AbcR1 and AbcR2 (also known as Sm15/16 or Sm15C1/C2) were identified as trans-acting sRNA to fine-tune nutrient uptake (Torres-Quesada et al., 2013). EcpR1 were considered to control bacterial cell cycle-related genes. The overexpression of EcpR1 leads to cell elongation, and deletion of EcpR1 reduces bacterial competitiveness, but does not influence on symbiosis in S. meliloti (Robledo et al., 2015a). RcsR1 regulates quorum sensing through the autoinducer synthase gene, sinI (Baumgardt et al., 2016), and MmgR controls the cellular Polyhydroxybutyrate (PHB) accumulation, but all these reported sRNAs are dispensable for symbiosis (Ceizel Borella et al., 2016; Lagares et al., 2016, 2017). Recently, a novel sRNA NfeR1 (nodule formation efficiency sRNA) was characterized in S. meliloti. NfeR1 shows high expression in salt stress and throughout the symbiotic interaction. NfeR1 is the firstly reported sRNA to contribute to infectivity, nodule formation and development, and symbiotic nitrogen-fixtion efficiency in rhizobia (Robledo et al., 2017). In this work, MH_s15 was shown to be an abundant sRNA by using Northern blot. In addition, it is worth noting that the transcription levels of the genes upstream and downstream of MH_s15 was upregulated in nodules compared with free-living cells, and an overexpression of MH_s15 resulted in a significantly defected symbiotic phenotype (data not showed). It suggested that the function of MH_s15 would be associated with symbiosis, but more experiments are needed to elucidate its function and mechanism.
In light of the strict host specificity of M. huakuii 7653R for A. sinicus L., and especially most of the sRNA identification and validation studies in endosymbiotic nitrogen- fixing bacteria have been done in species belonging to other genus but not in Mesorhizobium, we believe that our findings in this work provide novel insights into the sRNA regulation of symbiotic plant-microbe interactions.
Materials and Methods
Prediction of sRNAs in the IGRs of M. huakuii 7653R
The comparative analysis of IGRs was carried out among M. austrqalicum WSM2073, M. ciceri WSM1271, M. Loti MAFF303099, M. opportunistum WSM2075 and M. huakuii 7653R. Their genomic sequences and protein annotations were downloaded from the NCBI bacterial databases (ftp://ftp.ncbi.nih.gov/genomes/Bacteria). The predictions of sRNAs were performed as previously described (del Val et al., 2007). IGRs with a length ≥50 nt were scanned and homologous regions were clustered using WU-BLASTN (http://blast.wustl.edu) in the tested genomes with E-values ≤ 1e−5. The alignments were interrogated using eQRNA version 2.0.4 (ftp://ftp.genetics.wustl.edu/pub/eddy/software/qrna.tar.Z) as described previously (del Val et al., 2007). RNAz version 0.1.1 (http://www.tbi.univie.ac.at/~wash/RNAz) was used to screen the conserved intergenic regions in these candidate genomic sequences, and overlapping sequences between eQRNA and RNAz were selected. Subsequently, pairwise alignments in both eQRNA and RNAz were calculated using RNAfold (ViennaRNA-2.1.1) to determine the stable secondary structures, with the stem-length options set to a minimum of 20 bp. The promoters and Rho-independent terminators were analyzed using Promoter 2.0, RNAMotif and Erpin to further assess the resulting sequences. Finally, the nucleotide sequences of potential sRNAs were further inspected according to the results of RNA-seq deep sequencing. Each predicted sRNA was used to query the Rfam database to capture any annotated RNA regulatory elements.
Bacterial Cultivation and Plant Growth Conditions
Rhizobium strains were cultured in tryptone yeast medium (TY) or AMS medium (FAM2 minimal salts medium supplemented with 0.1% vol/vol aniline) at 28°C, and Escherichia coli strains were grown in Luria-Bertani (LB) at 37°C. When required, antibiotics were added at the following concentrations: ampicillin (Ap) 100 μg/ml; gentamicin (Gm) 20 μg/ml; streptomycin (Str) 20 μg/ml and kanamycin (Kan) 50 μg/ml.
In the Northern blotting analyses, free-living rhizobial cells were cultivated in TY at 28°C, and cells were harvested in the log phase (OD600 = 0.4–0.5). When required, the growing log cultures were subjected to the following stress conditions as described previously (Madhugiri et al., 2012; Klahn et al., 2015): oxidative stress (2 mM H2O2 for 10 min), salt stress (4 M NaCl for 30 min), acid stress (pH5 for 90 min), alkali stress (pH9 for 90 min), cold stress (20°C for 90 min) or heat stress (37°C for 90 min). For micro-oxygen stress, log phase cells growing in TY (OD600 = 0.4–0.5) were shifted to a low O2 environment by continuous sparging of the culture with purified N2 for 24 h (N2 ≥ 99%). Controls were non-stressed exponentially growing cultures. For symbiotic conditions, surface-sterilized seeds of A. sinicus L. cultivar XY202 were cultivated in pots filled with sterile sand and irrigated with nitrogen-free Fahraeus solution, and cultivation was performed in a green house with a 16 h light and 8 h darkness cycle at 22°C and 20°C, respectively. Root nodules were harvested 28 or 50 days after M. huakuii 7653R inoculation.
Extraction of Total RNA and Northern Blotting
Total RNA was extracted from M. huakuii 7653R to detect the selected sRNAs by Northern blotting. Following exposure to stress conditions, cells were quickly centrifuged, resuspended in TRIzol Reagent (Roche, USA), and immediately frozen in liquid nitrogen and stored at –80°C. Plant nodules were resuspended in TRIzol Reagent and frozen in liquid nitrogen, then immediately subjected to grinding using a baked mortar and pestle to crush the nodule tissues. Isolation of total RNA and Northern blotting were conducted according to the manufacturer's protocols. Crude RNA extracts were purified with DNaseI to remove potential genomic DNA, and the resulting RNA was stored at –80°C. Northern hybridizations were performed with end-labeled [γ−32P]ATP DNA probes (Table S7). Hybridization signals were monitored using a radioisotope imaging system (Fujifilm, FLA 5100). The 5S rRNA was used as a control, and hybridization band intensities were relatively quantified with the Biorad Quantity One software.
Extraction and Purification of Total RNA and Deep Sequencing of cDNAs
Total RNA was extracted using an RNeasy Mini Kit (Qiagen) following the manufacturer's instructions to carry out the RNA deep sequencing. The RNA integrity number (RIN) of the total RNA samples was checked using an Agilent Bioanalyzer 2100 (Agilent). Qualified total RNA was further purified using an RNeasy micro kit (QIAGEN) and RNase-Free DNase Set (QIAGEN) according to the manufacturer's instructions.
Total RNA, including small-sized RNAs from 28 dpi (days post-infection) nodules, 50 dpi nodules and free-living cells at OD600 = 0.8, was isolated and separated into small (< 100 nt) and long (> 100 nt) fractions using a miRNeasy Mini Kit (Qiagen) or a mirVana miRNA Isolation Kit (Ambion), respectively, according to the manufacturer's instructions. Total purified RNA was then treated with tobacco acid pyrophosphatase (TAP) (2.5 U/5 μg RNA in 50 μl for 2 h) to maintain the start regions of the transcription products. The quantity and quality of the processed RNA samples were assessed using of a NanoDrop ND-1000 spectrophotometer and Agilent 2100 Bioanalyzer. Qualified RNA was ligated to 5′- and 3′-RNA adapters by PCR, then the free adapters were removed and the remaining cDNA was purified. A Qubit® 2.0 Fluorometer and Agilent2100 instrument were used to monitor the abundance and size of cDNA, respectively. Solexa sequencing was performed using an Illumina HiSeq 2500. After that, all sequencing data were mapped by using the genomic sequence of M. huakuii 7653R to obtain the effective reads.
Target Prediction and the Gene Ontology (GO) Analysis of Candidate sRNAs
In order to figure out the potential function of the 16 new sRNA candidates identified in this work, we predicted the target genes for each of these sRNAs by using the webserver IntaRNA (http://rna.informatik.uni-freiburg.de/IntaRNA/Input.jsp). All of the encoded protein sequences of predicted target genes for each candidate sRNA were submitted to InterPro (http://www.ebi.ac.uk/interpro/sequence-search) for protein sequence analysis & classification, and the GO annotations (GO_ID) were extracted from the output. Further, the WEGO (Web Gene Ontology Annotation Plot) webserver (http://wego.genomics.org.cn/cgi-bin/wego/index.pl) was used to display the GO annotation results in histograms. The enriched GO terms of candidate small RNA target genes were summarized and demonstrated in a pie chart.
Other Tools and Software
Searches for promoters in the M. huakuii 7653R genome were performed using Promoter 2.0 (http://www.cbs.dtu.dk/services/Promoter/) and BPROM 184.108.40.2061 (http://www.molquest.com/). Predictions of Rho-independent terminators were implemented with RNAMotif 3.0.7 (http://casegroup.rutgers.edu/casegr-sh-2.5.html), TransTermHP 2.09 (transterm.cbcb.umd.edu/), FindTerm 220.127.116.111 (http://www.molquest.com/) and Erpin 5.4 (http://tagc.univ-mrs.fr/erpin/). Searches for known functional RNA sequences among the predicted sRNAs were conducted within the Rfam database using Rfam, version 12.0 (http://rfam.xfam.org/search) (Nawrocki et al., 2015). RNA secondary structures were analyzed using RNAfold 2.1.9 (http://rna.tbi.univie.ac.at/cgi-bin/RNAWebSuite/RNAfold.cgi). The target genes for the identified sRNAs were predicted using IntaRNA soft (http://rna.informatik.uni-freiburg.de/).
XF and LY conceived and designed the experiments. ZW, ZJ and XF prepared RNA samples and conducted the analysis of the RNA-seq data. ZZ and HB predicted and identified the sRNAs in the M. huakuii 7653R genome. WX and MB performed the functional analysis of the target genes. The article was written and revised by XF, ZW, and LY.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
This work was supported by the funds from the National Key Research and Development Program of China (grant no. 2016YFD0100700), the National Natural Science Foundation of China (grant no. 31670243; 31371549; 31570844), the Fundamental Research Funds for the Central Universities (2016PY025 and 2011QC066), and the Open Project of the State Key Laboratory of Agricultural Microbiology (no. AMLKF201502).
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fmicb.2017.01730/full#supplementary-material
Arraiano, C. M., Andrade, J. M., Domingues, S., Guinote, I. B., Malecki, M., Matos, R. G., et al. (2010). The critical role of RNA processing and degradation in the control of gene expression. FEMS Microbiol. Rev. 34, 883–923. doi: 10.1111/j.1574-6976.2010.00242.x
Baumgardt, K., Smidova, K., Rahn, H., Lochnit, G., Robledo, M., and Evguenieva-Hackenberg, E. (2016). The stress-related, rhizobial small RNA RcsR1 destabilizes the autoinducer synthase encoding mRNA sinI in Sinorhizobium meliloti. RNA Biol. 13, 486–499. doi: 10.1080/15476286.2015.1110673
Becker, A., Overloper, A., Schluter, J. P., Reinkensmeier, J., Robledo, M., Giegerich, R., et al. (2014). Riboregulation in plant-associated alpha-proteobacteria. RNA Biol. 11, 550–562. doi: 10.4161/rna.29625
Berghoff, B. A., Glaeser, J., Sharma, C. M., Vogel, J., and Klug, G. (2009). Photooxidative stress-induced and abundant small RNAs in Rhodobacter sphaeroides. Mol. Microbiol. 74, 1497–1512. doi: 10.1111/j.1365-2958.2009.06949.x
Bobrovskyy, M., Vanderpool, C. K., and Richards, G. R. (2015). Small RNAs regulate primary and secondary metabolism in gram-negative bacteria. Microbiol. Spectr. 3, 1–25. doi: 10.1128/microbiolspec.MBP-0009-2014
Caswell, C. C., Gaines, J. M., Ciborowski, P., Smith, D., Borchers, C. H., Roux, C. M., et al. (2012). Identification of two small regulatory RNAs linked to virulence in Brucella abortus 2308. Mol. Microbiol. 85, 345–360. doi: 10.1111/j.1365-2958.2012.08117.x
Ceizel Borella, G., Lagares, A. Jr., and Valverde, C. (2016). Expression of the Sinorhizobium meliloti small RNA gene mmgR is controlled by the nitrogen source. FEMS Microbiol. Lett. 363, 1–9. doi: 10.1093/femsle/fnw069
Chang, T. H., Wub, L. C., Lin, J. H., Huang, H. D., Liu, B. J., Cheng, K. F., et al. (2010). Prediction of small non-coding RNA in bacterial genomes using support vector machines. Expert Syst. Appl. 37, 5549–5557. doi: 10.1016/j.eswa.2010.02.058
del Val, C., Rivas, E., Torres-Quesada, O., Toro, N., and Jimenez-Zurdo, J. I. (2007). Identification of differentially expressed small non-coding RNAs in the legume endosymbiont Sinorhizobium meliloti by comparative genomics. Mol. Microbiol. 66, 1080–1091. doi: 10.1111/j.1365-2958.2007.05978.x
del Val, C., Romero-Zaliz, R., Torres-Quesada, O., Peregrina, A., Toro, N., and Jimenez-Zurdo, J. I. (2012). A survey of sRNA families in alpha-proteobacteria. RNA Biol. 9, 119–129. doi: 10.4161/rna.18643
Gonzalez, N., Heeb, S., Valverde, C., Kay, E., Reimmann, C., Junier, T., et al. (2008). Genome-wide search reveals a novel GacA-regulated small RNA in Pseudomonas species. BMC Genomics 9:167. doi: 10.1186/1471-2164-9-167
Guillier, M., and Gottesman, S. (2008). The 5′ end of two redundant sRNAs is involved in the regulation of multiple targets, including their own regulator. Nucleic Acids Res. 36, 6781–6794. doi: 10.1093/nar/gkn742
Hahn, J., Tsoy, O. V., Thalmann, S., Cuklina, J., Gelfand, M. S., and Evguenieva-Hackenberg, E. (2016). Small open reading frames, non-coding RNAs and repetitive elements in Bradyrhizobium japonicum USDA 110. PLoS ONE 11:e0165429. doi: 10.1371/journal.pone.0165429
Herouart, D., Baudouin, E., Frendo, P., Harrison, J., Santos, R., Jamet, A., et al. (2002). Reactive oxygen species, nitric oxide and glutathione: a key role in the establishment of the legume-Rhizobium symbiosis? Plant Physiol. Biochem. 40, 619–624. doi: 10.1016/S0981-9428(02)01415-8
Ignatov, D., Malakho, S., Majorov, K., Skvortsov, T., Apt, A., and Azhikina, T. (2013). RNA-Seq analysis of Mycobacterium avium non-coding transcriptome. PLoS ONE 8:e74209. doi: 10.1371/journal.pone.0074209
Janczarek, M., Rachwal, K., Marzec, A., Grzadziel, J., and Palusinska-Szysz, M. (2015). Signal molecules and cell-surface components involved in early stages of the legume-rhizobium interactions. Appl. Soil Ecol. 85, 94–113. doi: 10.1016/j.apsoil.2014.08.010
Klahn, S., Schaal, C., Georg, J., Baumgartner, D., Knippen, G., Hagemann, M., et al. (2015). The sRNA NsiR4 is involved in nitrogen assimilation control in cyanobacteria by targeting glutamine synthetase inactivating factor IF7. Proc. Natl. Acad. Sci. U.S.A. 112, E6243–E6252. doi: 10.1073/pnas.1508412112
Lagares, A. Jr., Borella, G. C., Linne, U., Becker, A., and Valverde, C. (2017). Regulation of polyhydroxybutyrate accumulation in Sinorhizobium meliloti by the trans-encoded small RNA MmgR. J. Bacteriol. 199:e00776-16. doi: 10.1128/JB.00776-16
Lagares, A. Jr., Roux, I., and Valverde, C. (2016). Phylogenetic distribution and evolutionary pattern of an alpha-proteobacterial small RNA gene that controls polyhydroxybutyrate accumulation in Sinorhizobium meliloti. Mol. Phylogenet. Evol. 99, 182–193. doi: 10.1016/j.ympev.2016.03.026
Lindemann, A., Moser, A., Pessi, G., Hauser, F., Friberg, M., Hennecke, H., et al. (2007). New target genes controlled by the Bradyrhizobium japonicum two-component regulatory system RegSR. J. Bacteriol. 189, 8928–8943. doi: 10.1128/JB.01088-07
Livny, J. (2007). Efficient annotation of bacterial genomes for small, noncoding RNAs using the integrative computational tool sRNAPredict2. Methods Mol. Biol. 395, 475–488. doi: 10.1007/978-1-59745-514-5_30
Lopez-Leal, G., Tabche, M. L., Castillo-Ramirez, S., Mendoza-Vargas, A., Ramirez-Romero, M. A., and Davila, G. (2014). RNA-Seq analysis of the multipartite genome of Rhizobium etli CE3 shows different replicon contributions under heat and saline shock. BMC Genomics 15:770. doi: 10.1186/1471-2164-15-770
MacLellan, S. R., Smallbone, L. A., Sibley, C. D., and Finan, T. M. (2005). The expression of a novel antisense gene mediates incompatibility within the large repABC family of alpha-proteobacterial plasmids. Mol. Microbiol. 55, 611–623. doi: 10.1111/j.1365-2958.2004.04412.x
Madhugiri, R., Pessi, G., Voss, B., Hahn, J., Sharma, C. M., Reinhardt, R., et al. (2012). Small RNAs of the Bradyrhizobium/Rhodopseudomonas lineage and their analysis. RNA Biol. 9, 47–58. doi: 10.4161/rna.9.1.18008
Man, S., Cheng, R., Miao, C., Gong, Q., Gu, Y., Lu, X., et al. (2011). Artificial trans-encoded small non-coding RNAs specifically silence the selected gene expression in bacteria. Nucleic Acids Res. 39:e50. doi: 10.1093/nar/gkr034
Mentz, A., Neshat, A., Pfeifer-Sancar, K., Pühler, A., Rückert, C., and Kalinowski, A. J. (2013). Comprehensive discovery and characterization of small RNAs in Corynebacterium glutamicum ATCC 13032. BMC Genomics 14:714. doi: 10.1186/1471-2164-14-714
Mesa, S., Hauser, F., Friberg, M., Malaguti, E., Fischer, H. M., and Hennecke, H. (2008). Comprehensive assessment of the regulons controlled by the FixLJ-FixK2-FixK1 cascade in Bradyrhizobium japonicum. J. Bacteriol. 190, 6568–6579. doi: 10.1128/JB.00748-08
Monteiro, C., Papenfort, K., Hentrich, K., Ahmad, I., Le Guyon, S., Reimann, R., et al. (2012). Hfq and Hfq-dependent small RNAs are major contributors to multicellular development in Salmonella enterica serovar Typhimurium. RNA Biol. 9, 489–502. doi: 10.4161/rna.19682
Morita, T., Ueda, M., Kubo, K., and Aiba, H. (2015). Insights into transcription termination of Hfq-binding sRNAs of Escherichia coli and characterization of readthrough products. RNA 21, 1490–1501. doi: 10.1261/rna.051870.115
Nawrocki, E. P., Burge, S. W., Bateman, A., Daub, J., Eberhardt, R. Y., Eddy, S. R., et al. (2015). Rfam 12.0: updates to the RNA families database. Nucleic Acids Res. 43, 130–137. doi: 10.1093/nar/gku1063
Pauly, N., Pucciariello, C., Mandon, K., Innocenti, G., Jamet, A., Baudouin, E., et al. (2006). Reactive oxygen and nitrogen species and glutathione: key players in the legume-Rhizobium symbiosis. J. Exp. Bot. 57, 1769–1776. doi: 10.1093/jxb/erj184
Pellin, D., Miotto, P., Ambrosi, A., Cirillo, D. M., and Di Serio, C. (2012). A genome-wide identification analysis of small regulatory RNAs in Mycobacterium tuberculosis by RNA-Seq and conservation analysis. PLoS ONE 7:e32723. doi: 10.1371/journal.pone.0032723
Peng, J., Hao, B., Liu, L., Wang, S., Ma, B., Yang, Y., et al. (2014). RNA-Seq and microarrays analyses reveal global differential transcriptomes of Mesorhizobium huakuii 7653R between bacteroids and free-living cells. PLoS ONE 9:e93626. doi: 10.1371/journal.pone.0093626
Puppo, A., Pauly, N., Boscari, A., Mandon, K., and Brouquisse, R. (2013). Hydrogen peroxide and nitric oxide: key regulators of the Legume-Rhizobium and mycorrhizal symbioses. Antioxid. Redox Signal. 18, 2202–2219. doi: 10.1089/ars.2012.5136
Rivera-Urbalejo, A., Perez-Oseguera, A., Carreon-Rodriguez, O. E., and Cevallos, M. A. (2015). Mutations in an antisense RNA, involved in the replication control of a repABC plasmid, that disrupt plasmid incompatibility and mediate plasmid speciation. Plasmid 78, 48–58. doi: 10.1016/j.plasmid.2015.01.004
Robledo, M., Peregrina, A., Millan, V., Garcia-Tomsig, N. I., Torres-Quesada, O., Mateos, P. F., et al. (2017). A conserved alpha-proteobacterial small RNA contributes to osmoadaptation and symbiotic efficiency of rhizobia on legume roots. Environ. Microbiol. 19, 2661–2680. doi: 10.1111/1462-2920.13757
Saramago, M., Barria, C., Dos Santos, R. F., Silva, I. J., Pobre, V., Domingues, S., et al. (2014). The role of RNases in the regulation of small RNAs. Curr. Opin. Microbiol. 18, 105–115. doi: 10.1016/j.mib.2014.02.009
Schluter, J. P., Reinkensmeier, J., Daschkey, S., Evguenieva-Hackenberg, E., Janssen, S., Janicke, S., et al. (2010). A genome-wide survey of sRNAs in the symbiotic nitrogen-fixing alpha-proteobacterium Sinorhizobium meliloti. BMC Genomics 11:245. doi: 10.1186/1471-2164-11-245
Sciotti, M. A., Chanfon, A., Hennecke, H., and Fischer, H. M. (2003). Disparate oxygen responsiveness of two regulatory cascades that control expression of symbiotic genes in Bradyrhizobium japonicum. J. Bacteriol. 185, 5639–5642. doi: 10.1128/JB.185.18.5639-5642.2003
Stubben, C. J., Micheva-Viteva, S. N., Shou, Y., Buddenborg, S. K., Dunbar, J. M., and Hong-Geller, E. (2014). Differential expression of small RNAs from Burkholderia thailandensis in response to varying environmental and stress conditions. BMC Genomics 15:385. doi: 10.1186/1471-2164-15-385
Su, G., Christensen, O. F., Janss, L., and Lund, M. S. (2014). Comparison of genomic predictions using genomic relationship matrices built with different weighting factors to account for locus-specific variances. J. Dairy Sci. 97, 6547–6559. doi: 10.3168/jds.2014-8210
Toffano-Nioche, C., Nguyen, A. N., Kuchly, C., Ott, A., Gautheret, D., Bouloc, P., et al. (2012). Transcriptomic profiling of the oyster pathogen Vibrio splendidus opens a window on the evolutionary dynamics of the small RNA repertoire in the Vibrio genus. RNA 18, 2201–2219. doi: 10.1261/rna.033324.112
Torres-Quesada, O., Millan, V., Nisa-Martinez, R., Bardou, F., Crespi, M., Toro, N., et al. (2013). Independent activity of the homologous small regulatory RNAs AbcR1 and AbcR2 in the legume symbiont Sinorhizobium meliloti. PLoS ONE 8:e68147. doi: 10.1371/journal.pone.0068147
Tsai, C. H., Liao, R., Chou, B., Palumbo, M., and Contreras, L. M. (2015). Genome-wide analyses in bacteria show small-rna enrichment for long and conserved intergenic regions. J. Bacteriol. 197, 40–50. doi: 10.1128/JB.02359-14
Ulve, V. M., Sevin, E. W., Cheron, A., and Barloy-Hubler, F. (2007). Identification of chromosomal alpha-proteobacterial small RNAs by comparative genome analysis and detection in Sinorhizobium meliloti strain 1021. BMC Genomics 8:467. doi: 10.1186/1471-2164-8-467
Valverde, C., Livny, J., Schluter, J. P., Reinkensmeier, J., Becker, A., and Parisi, G. (2008). Prediction of Sinorhizobium meliloti sRNA genes and experimental detection in strain 2011. BMC Genomics 9:416. doi: 10.1186/1471-2164-9-416
Vercruysse, M., Fauvart, M., Cloots, L., Engelen, K., Thijs, I. M., Marchal, K., et al. (2010). Genome-wide detection of predicted non-coding RNAs in Rhizobium etli expressed during free-living and host-associated growth using a high-resolution tiling array. BMC Genomics 11:53. doi: 10.1186/1471-2164-11-53
Voss, B., Holscher, M., Baumgarth, B., Kalbfleisch, A., Kaya, C., Hess, W. R., et al. (2009). Expression of small RNAs in Rhizobiales and protection of a small RNA and its degradation products by Hfq in Sinorhizobium meliloti. Biochem. Biophys. Res. Commun. 390, 331–336. doi: 10.1016/j.bbrc.2009.09.125
Wang, L. J., Wang, W. W., Li, F. D., Zhang, J. H., Wu, J. H., Gong, Q. G., et al. (2015). Structural insights into the recognition of the internal A-rich linker from OxyS sRNA by Escherichia coli Hfq. Nucleic Acids Res. 43, 2400–2411. doi: 10.1093/nar/gkv072
Wang, S., Hao, B., Li, J., Gu, H., Peng, J., Xie, F., et al. (2014). Whole-genome sequencing of Mesorhizobium huakuii 7653R provides molecular insights into host specificity and symbiosis island dynamics. BMC Genomics 15:440. doi: 10.1186/1471-2164-15-440
Wilderman, P. J., Sowa, N. A., FitzGerald, D. J., FitzGerald, P. C., Gottesman, S., Ochsner, U. A., et al. (2004). Identification of tandem duplicate regulatory small RNAs in Pseudomonas aeruginosa involved in iron homeostasis. Proc. Natl. Acad. Sci. U.S.A. 101, 9792–9797. doi: 10.1073/pnas.0403423101
Wilms, I., Moller, P., Stock, A. M., Gurski, R., Lai, E. M., and Narberhaus, F. (2012a). Hfq influences multiple transport systems and virulence in the plant pathogen Agrobacterium tumefaciens. J. Bacteriol. 194, 5209–5217. doi: 10.1128/JB.00510-12
Wilms, I., Overloper, A., Nowrousian, M., Sharma, C. M., and Narberhaus, F. (2012b). Deep sequencing uncovers numerous small RNAs on all four replicons of the plant pathogen Agrobacterium tumefaciens. RNA Biol. 9, 446–457. doi: 10.4161/rna.17212
Keywords: small RNAs, comparative analysis, RNA-seq, Northern blotting, Mesorhizobium huakuii
Citation: Fuli X, Wenlong Z, Xiao W, Jing Z, Baohai H, Zhengzheng Z, Bin-Guang M and Youguo L (2017) A Genome-Wide Prediction and Identification of Intergenic Small RNAs by Comparative Analysis in Mesorhizobium huakuii 7653R. Front. Microbiol. 8:1730. doi: 10.3389/fmicb.2017.01730
Received: 03 July 2017; Accepted: 24 August 2017;
Published: 08 September 2017.
Edited by:Essaid Ait Barka, University of Reims Champagne-Ardenne, France
Reviewed by:Ming Wang, University of California, Riverside, United States
Markus Rienth, Institut National de la Recherche Agronomique (INRA), France
Copyright © 2017 Fuli, Wenlong, Xiao, Jing, Baohai, Zhengzheng, Bin-Guang and Youguo. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Li Youguo, email@example.com