- Research article
- Open Access
Transcriptome profiling of a Rhizobium leguminosarum bv. trifolii rosR mutant reveals the role of the transcriptional regulator RosR in motility, synthesis of cell-surface components, and other cellular processes
BMC Genomics volume 16, Article number: 1111 (2015)
Rhizobium leguminosarum bv. trifolii is a soil bacterium capable of establishing a symbiotic relationship with red clover (Trifolium pratense). The presence of surface polysaccharides and other extracellular components as well as motility and competitiveness are essential traits for both adaptation of this bacterium to changing environmental conditions and successful infection of host plant roots. The R. leguminosarum bv. trifolii rosR gene encodes a protein belonging to the family of Ros/MucR transcriptional regulators, which contain a Cys2His2-type zinc-finger motif and are involved in the regulation of exopolysaccharide synthesis in several rhizobial species. Previously, it was established that a mutation in the rosR gene significantly decreased exopolysaccharide synthesis, increased bacterial sensitivity to some stress factors, and negatively affected infection of clover roots.
RNA-Seq analysis performed for the R. leguminosarum bv. trifolii wild-type strain Rt24.2 and its derivative Rt2472 carrying a rosR mutation identified a large number of genes which were differentially expressed in these two backgrounds. A considerable majority of these genes were up-regulated in the mutant (63.22 %), indicating that RosR functions mainly as a repressor. Transcriptome profiling of the rosR mutant revealed a role of this regulator in several cellular processes, including the synthesis of cell-surface components and polysaccharides, motility, and bacterial metabolism. Moreover, it was established that the Rt2472 strain was characterized by a longer generation time and showed an increased aggregation ability, but was impaired in motility as a result of considerably reduced flagellation of its cells.
The comparative transcriptome analysis of R. leguminosarum bv. trifolii wild-type Rt24.2 and the Rt2472 mutant identified a set of genes belonging to the RosR regulon and confirmed the important role of RosR in the regulatory network. The data obtained in this study indicate that this protein affects several cellular processes and plays an important role in bacterial adaptation to environmental conditions.
Rhizobia are a unique group of bacteria which can either exist as free-living organisms in soil or establish a nitrogen-fixing symbiosis with legumes [1, 2]. This type of plant–microbe interaction is essential for the functioning of the biosphere, since it provides nutrients to plants and increases soil fertility and field crops. However, before rhizobia can find a compatible host plant, they often must survive long periods of time as free-living bacteria in the soil. In such periods, rhizobial cells are exposed to the action of environmental factors such as nutrient limitation, pH, salinity, drought, temperature variation, heavy metals, and oxidative stress [3–9]. In order to adapt to these conditions, rhizobia have developed a wide range of strategies and features that allow them to survive in the soil. One of the most important of these adaptations is the composition of the bacterial envelope. The outer surface of these bacteria contains various polysaccharides (PSs), such as lipopolysaccharide (LPS), capsular polysaccharide (CPS) and cellulose fibrils, as well as the neutral polysaccharide (NP, glucomannan) and gel-forming polysaccharide (GPS) recently reported in Rhizobium leguminosarum [10–12]. NP, cellulose fibrils, and LPS are required for attachment to and biofilm formation on soil particles and plant roots, as well as successful infection of the host plant and adaptation to the conditions prevailing inside the nodule, a specific organ formed on legume roots [13–17].
In addition, exopolysaccharide (EPS), secreted in large amounts by rhizobia and weakly associated with their cells, plays an important protective role against desiccation and other stress conditions occurring in the soil. This polymer is also crucial for the attachment of rhizobia to abiotic as well as biotic surfaces, including plant roots, and establishment of an effective symbiosis [11, 12, 18, 19]. Another type of PS synthesized by these bacteria is β-glucan, which is found in their periplasmic space and is required for adaptation to hypo-osmotic conditions [20–22].
Apart from the influence of envelope components, bacterial adaptation to changing soil conditions and interaction with plants are affected by some of the proteins that rhizobia secrete into the environment, including those exported via Type I secretion system (e.g. NodO, glycanases, adhesins, and substrate-binding proteins) [11, 23–25]. Some of these secreted proteins are required for the attachment of rhizobia to host plant roots, which is a very important step in the initiation of symbiosis. An R. leguminosarum protein called rhicadhesin plays a significant role in bacterial attachment to root hairs, particularly at alkaline pH [10, 11]. Other proteins involved in this process include PlyB glycanase (which affects EPS processing), cadherin-like proteins (calcium-binding adherence proteins), and Rhizobium-adhering proteins (Raps) secreted via Type I secretion system encoded in R. leguminosarum by the prsD and prsE genes [11, 23, 24]. Moreover, some proteins secreted by rhizobia into the environment allow the bacteria to survive in the absence of their legume hosts. Among these, several proteins associated with the uptake of various nutrients have been identified (e.g. proteins binding sorbitol, ribose and other sugars, leucine/isoleucine/valine, arginine, and dipeptides) .
Furthermore, motility of rhizobial cells is an essential trait for their survival in the environment as well as competitiveness [26–31]. In the soil, the availability of carbon sources utilizable by bacteria is limited, since most components of organic matter form complexes with other compounds, thus becoming resistant to bacterial degradation . Because of this, the rhizosphere, a narrow soil region rich in various types of substances secreted by plant roots, is a very attractive niche for bacteria [11, 33, 34]. Therefore, the ability to sense nutrients and move towards them provides a competitive advantage for motile and chemotactic rhizobia over non-motile and non-chemotactic strains . Several studies have shown that non-motile or non-chemotactic bacteria are less adaptive to changing environmental conditions and less competitive when it comes to plant root infection [27, 28, 35–38]. Additionally, swarming of rhizobial cells is affected by exudates from legume seeds .
Rhizobium leguminosarum bv. trifolii is a microsymbiont of red clover (Trifolium pratense). Previously, rosR encoding a regulatory protein which positively influenced EPS production was identified and characterized in this bacterium [40, 41]. A mutation in this gene resulted in a significant decrease in EPS synthesis . R. leguminosarum bv. trifolii RosR belongs to the family of Ros/MucR transcriptional regulators, which are involved in the regulation of EPS synthesis in several rhizobial species, including Sinorhizobium meliloti, Rhizobium etli, and Agrobacterium tumefaciens [43–46]. RosR is a small protein (15.7 kDa) characterized by a Cys2His2-type zinc-finger motif, which is responsible for binding to a 22-bp-long sequence called the RosR-box. Previously, it was established that RosR recognized the RosR-box motif located in the rosR upstream region and negatively regulated the transcription of its own gene . A R. leguminosarum bv. trifolii strain with a rosR mutation formed dry, wrinkled colonies on agar plates, which were significantly different from those formed by the wild-type bacteria. Furthermore, this mutant strain showed some changes in the O-chain of LPS. The rosR mutation also resulted in a decreased biofilm formation, a higher sensitivity to osmotic and oxidative stresses, and a significantly impaired symbiosis with clover plants (in comparison to the wild-type strain, the nodules formed were less numerous and were unable to fix nitrogen) [8, 47]. All these data suggest that RosR is an essential protein involved in several cellular processes and, possibly, an important element of the rhizobial regulatory network. Although previous results concerning the rosR mutant indicate that RosR is engaged in adaptation to stress conditions, it is still unknown how many genes are under the control of this transcriptional regulator.
In this study, we performed a comparative transcriptomic analysis of R. leguminosarum bv. trifolii wild-type 24.2 and its derivative mutant strain Rt2472rosR, which provided novel data on RosR-mediated regulation of gene expression in this bacterium. It was found that RosR influenced the expression of a large number of genes, including those related to the synthesis of cell-surface components and polysaccharides, motility, and metabolic pathways. Among these genes, a significant majority were up-regulated in the rosR mutant, suggesting that RosR functions in R. leguminosarum bv. trifolii cells mainly as a negative regulatory protein.
Results and discussion
RNA-Seq analysis of the wild-type strain Rt24.2 and its derivative rosR mutant Rt2472
Our previous studies of R. leguminosarum bv. trifolii rosR had suggested that this gene could play a global regulatory role in the functioning of rhizobial cells. As a further step in the investigation of this problem, in the present study, we used comparative transcriptomic analysis to establish the set of genes belonging to the RosR regulon. We compared wild-type strain Rt24.2 with its derivative, the Rt2472rosR mutant, obtained via random mutagenesis, which has a mini-Tn5 transposon located inside the rosR coding region, between 151 bp and 152 bp of the open reading frame . First, in order to identify the genetic bases of Rt24.2 and to use it as a reference strain to map the sequence reads from the transcriptome data sets, a draft genomic sequence of this strain was obtained (181 contigs with a total length of 7,653,217 bp, GC content = 60.6 %). The sequence of Rt24.2 was compared with genomes of closely related rhizobial strains (including R. leguminosarum bv. viciae 3841, R. leguminosarum bvs. trifolii WSM2304 and WSM1325, and Rhizobium etli CFN42). In total, 7,374 potential coding sequences were identified in the Rt24.2 genome. Among these, 7,149 sequences (96.95 %) were matched to Rhizobium, 141 (1.91 %) to other bacterial species, and only 84 (1.13 %) were strain-specific and were not matched to any rhizobia.
In order to compare the expression of these genes in the rosR mutant to their expression in the wild-type bacterium, three RNA-Seq libraries were prepared for the two strains grown in the same conditions. Qualitative analyses of RNA samples after depletion of 23S and 16S rRNA indicated that they were of high quality and that a great majority of rRNA had been removed. The mRNAs of the Rt24.2 and Rt2472 strains were sequenced using Illumina MiSeq with SBS technology (general features of the individual runs are shown in Additional file 1A). On average, 5,368,924 reads (SD = 1,005,388) for the wild-type strain and 5,365,247 reads (SD = 1,864,858) for the rosR mutant were unambiguously mapped. This indicated that similar amounts of data were mapped for each strain. The results were subsequently analyzed using the Cuffdiff tool from the Cufflinks package to normalize the data and test for differential expression between the strains tested (Additional file 1B). Box plots generated in order to compare the distributions of RNA-Seq data between biological repetitions of an individual strain and between Rt24.2 and Rt2472 showed similar values for these strains (comparable distributions seen as similar box sizes and box plot whiskers) (Additional file 1C).
Identification of genes differentially expressed in the rosR mutant Rt2472 and the wild-type strain Rt24.2
An analysis of the functional composition of the Rt24.2 transcriptome showed that the most numerously represented categories were related to metabolism processes, especially the functional groups (COGs) G, E, P, and C, which are associated with the transport and metabolism of carbohydrates, amino acids, and inorganic ions as well as energy production and conversion, respectively (Fig. 1). Also, COG classes K (transcription), J (translation/ribosome structure and biogenesis), M (cell wall/membrane/envelope biogenesis), and the poorly characterized class R were highly represented in the R. leguminosarum transcriptome.
In order to establish which genes belonged to the RosR regulon in R. leguminosarum, transcriptional profiles for Rt24.2 and Rt2472 were obtained using the genome-wide RNA-Seq approach. The high reproducibility between the analyzed strains (evidenced by the high correlation factors of read counts) and the relative expression values of the RNA-Seq data allowed us to obtain p-values adjusted with Cuffdiff (Cufflinks package). A histogram of FPKM values (Fragments Per Kilobase of transcript per Million fragments mapped) for the genes, with three biological repetitions for each strain analyzed, is presented in Additional file 2A. A point diagram of FPKM values between the samples and a diagram indicating the dependence of p-values on fold changes are shown in Additional file 2B and C, respectively. Based on the fold changes of gene expressions obtained in the rosR mutant and the wild-type strain (log2 of mut/wt values >2), it was established that 1106 genes were transcribed at significantly different levels in these genetic backgrounds. This indicated that RosR influences the expression of many loci and plays an essential role in the rhizobial regulatory network (Additional file 3). Among the genes in question, a significant majority were up-regulated in the rosR mutant (63.22 %), whereas 36.78 % were down-regulated (Fig. 1a), suggesting that RosR functions mainly as a negative regulator of gene transcription in this bacterium. We successfully classified 89.48 % of the genes into particular COGs (Fig. 1b) . Most of these genes belonged to the following groups: transcription (COG K) (9.33 %), general functions (R) (10.13 %), unknown function (S) (8.05 %), transport and metabolism of amino acids (E) (7.74 %) and carbohydrates (G) (7.18 %), signal transduction mechanisms (T) (5.74 %), cell wall/membrane/envelope biogenesis (M) (5.41 %), and cell motility (N) (4.23 %). Moreover, it was found that for a majority of the COGs (K, M, O, G, C, T, V, P, I, H, R, S and that with hypothetical functions) a significantly larger number of the genes were up-regulated than down-regulated in the rosR mutant in relation to the parental strain (Fig. 1b). In contrast, a great majority of the genes belonging to COGs N and E were expressed at significantly lower levels in the rosR mutant versus the wild-type, indicating that RosR positively affected the transcription of these genes.
When the individual genes belonging to the RosR regulon were analyzed, the highest fold changes in mut/wt expression were found for Rt780_45 encoding a periplasmic protein involved in polysaccharide export (log2 mut/wt = 7.89), Rt622_32 and Rt622_33 encoding subunits of a cytochrome/quinol oxidase (−6.94 and −7.72, respectively), Rt780_44 encoding a glycosyl transferase involved in cell wall biogenesis (7.64), Rt620_2 encoding a transcriptional regulator from the Crp/Fnr family (7.05), and Rt622_21 encoding an adenylate cyclase (-6.32) (Fig. 2). A heatmap for 100 genes with the highest fold change values is shown in Fig. 3. The characteristics and putative functions of these genes are specified in Additional file 4.
Transcription factors and signal transduction mechanisms
A functional category that was highly represented in R. leguminosarum bv. trifolii and differentially expressed in the Rt2472 mutant and the parental strain Rt24.2 was transcription (COG K) (Fig. 1b). A significant majority of the genes from this group were up-regulated (92), with only 25 genes being down-regulated in the rosR mutant. These genes encoded transcriptional factors belonging to various families such as Crp/Fnr, LysR, TetR, LysR, AraC, LuxR, MarR, and GntR. Among these proteins, Rt782_109 (P-II), Rt648_50 (PhoB), Rt648_25 (LacI), Rt620_2 (Crp), Rt620_13 (ROK), and Rt646_15 (RaiR) seem to be the most interesting (Fig. 2). The response regulator PhoB (Rt648_50) together with the sensor kinase PhoR (Rt628_39) form a two-component regulatory system used by rhizobial cells to recognize and adapt to phosphate starvation . The catabolic protein Crp, the ROK protein binding Crp and cAMP, the LacI protein and the adenylate cyclase Rt622_21 are most probably engaged in the regulation of carbon metabolism. Rt782_109 (P-II) is highly homologous to GlnK proteins involved in nitrogen regulation in R. etli and S. meliloti , whereas Rt646_15 is identical to the LuxR-type regulator RaiR of R. leguminosarum bv. viciae engaged in quorum sensing . A transcriptome analysis performed previously for an S. meliloti glnK mutant indicated that GlnK, which belongs to P-II proteins, is an important component of the regulatory cascade involved in nitrogen stress adaptation in rhizobial cells . Also, several genes encoding LuxR-type regulators (e.g. Rt615_20 and Rt643_8), which possess a CheY-like domain responsible for binding chemotaxis proteins were found in this group (Fig. 2).
Moreover, many genes from COG T, involved in signal transduction mechanisms, were differentially expressed in the two genetic backgrounds studied, and a majority of them were down-regulated in the rosR mutant (Figs. 1-2, Additional file 4). Among them, there were several genes coding for methyl-accepting chemotaxis proteins (e.g. Rt761_94 and Rt647_8), chemotaxis proteins Rt628_21 (CheY) and Rt628_27 (CheX), and the response regulator Rt628_22 (CheB), which seem to play a significant role in signaling. All these data indicate that RosR affects the expression of many genes related to the response of rhizobial cells to environmental conditions (nitrogen, phosphate and carbon starvation) and quorum sensing.
Cell motility and chemotaxis
Among the functional groups analyzed in this study, COG N is unique because of the fact that a great majority of the genes (90.6 %) associated with cell motility and chemotaxis were down-regulated in the rosR mutant compared to the wild-type strain (Fig. 1b). This indicated that RosR positively affected this trait in R. leguminosarum. Among these genes were those encoding the transcriptional regulator Rt634_21 (Rem), flagellar biosynthesis proteins (e.g. Rt628_11 (FliN) and Rt628_18 (FliF)), the flagellar motor protein Rt634_17 (MotB), as well as other flagellar proteins (Rt634_4 (FlgG), Rt634_8 (FlgH), and Rt634_5 (FlgA)). Moreover, transcription of Rt634_12 (flaA), Rt634_13 (flaB), Rt634_14 (flaC), and Rt634_15 (flaD) genes encoding flagellins A, B, C, and D, respectively, was significantly decreased in the rosR mutant (Fig. 2).
In addition, several genes involved in chemotaxis were down-regulated in the rosR mutant. Apart from the MCP-type chemoreceptor genes mentioned earlier, also Rt628_24, Rt634_19, and Rt628_25 genes, showing a high homology to cheW, cheD, and cheA, respectively, of R. leguminosarum bv. viciae strains 3841 and VF39 were found in this group (Fig. 2) [27–29].
Previously, Tambalo and others [28, 29] had identified two chemotaxis gene clusters (che1 and che2) involved in the motility of R. leguminosarum and found that the che1 cluster was the major genetic region that controlled swimming and chemotaxis (the che2 cluster had a minor effect on chemotaxis). The first cluster contained such chemotaxis genes as mcpC, cheY, cheA, cheW, cheR, cheB, and cheD. Moreover, flaA, flaB, flaC, flaD, motA, motB, motC, fli, and flg genes, as well as the regulatory genes rem, visR, and visN were located in this region. Protein products of flaA, flaB, and flaC proved to be major components of the flagellar filament . As reported earlier, LuxR-type regulators VisN/VisR and the Rem protein are components of the regulatory cascade involved in the expression of many flagellar, motility, and chemotaxis genes in R. leguminosarum. A great majority of these genes (flaA, flaB, flaC, flaD, motA, motB, mcpD, and the che1 operon) are under positive regulation by this cascade [28, 29]. In this study, we found that the expression of many of these motility and chemotaxis genes was also positively regulated by RosR, which suggests that this protein may be an important element of this regulatory cascade.
Carbon and nitrogen transport and metabolism, and energy production
In the rosR mutant, the expression of a large group of genes related to bacterial metabolism was significantly changed as well (Fig. 1b). This was especially visible in the case of the genes classified to COGs E (97 genes), G (90 genes), and C (54 genes). Many of these genes encoded enzymes and components of various transport systems. In group E, changes in expression were particularly evident in the case of genes encoding two transport systems (Fig. 2, Additional file 4): the expression of genes Rt620_74, Rt620_75, Rt620_76, and Rt620_78 coding for components of the ABC-type branched-chain amino acid transport system was strongly up-regulated, whereas transcription of genes Rt619_162, Rt619_ 163, Rt619_164, and Rt619_165 encoding components of the ABC-type dipeptide/oligopeptide transport system was down-regulated in the rosR mutant.
Similarly in the case of group G, expression of many genes for sugar transport systems was significantly decreased in the mutant background. Among them, the most important role seems to be played by Rt794_91 and Rt794_95 (Fig. 2).
In addition, large differences in expression levels between the mutant and the wild-type cells were found for genes encoding different enzymes (e.g. Rt615_5 coding for a putative phosphoenolpyruvate carboxylase, Rt642_41 coding for a NADH dehydrogenase, and Rt622_33 coding for a cytochrome/quinol oxidase) (Fig. 2).
The differences observed in the expression of group G, E, and C genes suggested large disturbances in the cellular metabolism of the rosR mutant. These results are in agreement with the data obtained previously in Biolog tests which indicated that the rosR mutant utilized fewer carbon and nitrogen sources, and some of them utilized less effectively than the parental strain .
Several authors have reported that the ability to utilize numerous carbon and energy sources plays a very important role in both adaptation of rhizobia to soil conditions and their competitiveness in host plant infection. For example, catabolism of homoserine, a substantial component of pea root exudate, was linked with competiveness for nodulation of this host plant, as evidenced for R. leguminosarum bv. viciae 3841 . On the other hand, an inability to catabolize galactose led to an increased ability to compete for nodule occupancy in S. meliloti .
Cell-surface component synthesis and cell envelope biogenesis
Several genes associated with cell envelope biogenesis and synthesis of various polysaccharides and other surface components were differentially expressed in the analyzed strains (Fig. 1b). In this functional group (COG M), high up-regulation in the rosR mutant was established for the following loci: Rt780_45 and Rt780_39 encoding polysaccharide biosynthesis proteins, and Rt780_44 and Rt780_43 encoding glycosyl transferases (Fig. 2). These genes are located in the same genetic region and are probably engaged in LPS biosynthesis. Similarly, expression of some genes possibly involved in LPS or cell-wall biosynthesis, which are clustered in another genomic region called Pss-II, was highly increased in the mutant background. Among these were Rt620_3 and Rt620_4 coding for glycosyl transferases, and Rt620_1 and Rt620_5 coding for polysaccharide biosynthesis proteins, homologous to pssY and pssL2 of the TA1 strain, respectively .
In the wild-type Rt24.2, RosR also decreased the expression of Rt782_16, encoding the cyclic β-1,2-glucan ABC transporter NdvA, and Rt627_76, which shows a high homology to the RL4404 (gelA) gene involved in GPS synthesis in R. leguminosarum bv. viciae 3841 (Fig. 2). β-1,2-glucan is required for adaptation of rhizobia to hypo-osmotic conditions, motility, and efficient symbiosis with host plants, whereas GPS is important for attachment to root hairs and competitive nodule infection [14, 21, 22].
In contrast, several genes potentially involved in peptidoglycan biosynthesis (Rt634_20, Rt768_57 and Rt621_6) and formation of cellulose fibrils (Rt623_74), as well as those encoding other glycosyl transferases of unknown functions (e.g. Rt623_92 and Rt634_11) were down-regulated in the rosR mutant (Fig. 2) .
Previous studies concerning RosR had indicated that this protein was a positive regulator of EPS synthesis in R. leguminosarum . Therefore, we had expected that expression of many of the genes participating in this process, among them a significant majority of those located in the chromosomal cluster named Pss-I , would be positively affected by this regulator. Surprisingly, it was found that only a few of these genes encoding glycosyl transferases (Rt772_4 (pssS), Rt772_10 (pssI), Rt772_12 (pssC), and Rt620_62 (pssA)), were down-regulated in the mutant background, and the level of this regulation was very low (log2 fold change from −1.45 to −1.12) (Fig. 2). This finding suggests that the observed levels of regulation of pss gene transcription by RosR are sufficient for proper synthesis of EPS in R. leguminosarum. However, it cannot be excluded that RosR exerts a positive influence on the expression of other genes linked with EPS synthesis, whose function in this process has not yet been established. Previously, we had found that the expression of pssA, which encodes a glycosyl transferase involved in the first step of EPS synthesis, was positively regulated by RosR, but the observed effect was not high (a ~2-fold difference in the expression between the rosR mutant and the wild-type strain) . A similar observation was made for MucR, a positive regulator of succinoglycan synthesis in S. meliloti which shows a significant similarity to the R. leguminosarum RosR; it was established that this protein had an only slight effect on the expression of the genes involved in this process (namely, it increased the transcription of exoYFQ and exoK and decreased the expression of exoH and exoX genes) .
In summary, the data obtained in this study indicate that RosR regulates the expression of many genes related to the synthesis of various rhizobial polysaccharides. In the wild-type cells, this protein represses the transcription of several genes involved in the biosynthesis of LPS, β-glucan, GPS and/or CPS, while positively regulates the expression of genes which participate in peptidoglycan, cellulose and EPS synthesis. All these components play an essential role in the proper functioning of the bacterial envelope and adaptation of the bacteria to environmental conditions [11, 54].
Furthermore, RosR proved to be a negative regulator of the expression of several genes from COG M, which encode proteins engaged in the processing of extracellular polysaccharides and in bacterial behavior. Among them, the most interesting genes were Rt620_6 encoding a 1,3-1,4-β-glycanase, Rt772_17 encoding component PrsD of the Type I protein secretion system, Rt772_19 coding for an autoaggregation protein, and Rt644_17 coding for a protein required for bacterial attachment to host plants (Fig. 2). The protein product of the Rt772_19 gene exhibited a high homology to RapA1 agglutinin of R. leguminosarum bvs. trifolii and viciae, which belongs to the family of Rhizobium-adhering proteins (Raps) [24, 55].
The comparative RNA-Seq analysis of the rosR mutant and the wild-type strain revealed the role of the RosR protein in several cellular processes. Similar results had been obtained for S. meliloti strain 1021 with a mutation in fadD encoding a long-chain fatty acyl-coenzyme A ligase. In this latter strain, more than a thousand genes were identified as differentially expressed, including those for some metabolic activities, chemotaxis, motility, and iron uptake as well as stress-related genes . Also, an RNA-Seq analysis of the traI and ngrI mutants of S. fredii NGR234 had identified a large set of genes which were differentially expressed in comparison to the parental strain (316 and 466 genes, respectively) . These genes included, among others, those related to quorum sensing, motility and EPS biosynthesis. These data indicate that mutations in some essential genes in rhizobia may affect the expression of a huge number of other genes involved in various cellular processes.
It was established in this study that the mutation in rosR did not affect the genome structure or RNA processing (COGs A and B). However, some genes involved in cell division (COG D) were differentially expressed in the analyzed backgrounds (Fig. 1b). A majority of these genes were down-regulated in the rosR mutant (e.g. Rt792_202 coding for a chromosome segregation protein and Rt626_125 and Rt657_260 genes coding for cell division proteins FtsZ and FtsK, respectively) (Fig. 2).
Properties of rosR mutant cells
In order to verify the data obtained from the RNA-Seq analysis, we characterized several traits of the rosR mutant. To do so, in addition to the wild-type Rt24.2 and the Rt2472 mutant strains, Rt2472 after complementation of the mutation by rosR introduced on the pRC24 plasmid was used. At first, it was established that the Rt2472 mutant grew slower than the Rt24.2 and Rt2472(pRC24) strains on solid as well as in liquid media. Therefore, the doubling time of these strains was determined. In this experiment, it was established that the rosR mutant had a longer generation time (6.0 ± 0.5 h) than the wild-type (4.5 ± 0.5 h) and the Rt2472(pRC24) (5.0 ± 0.5 h) cells. Moreover, Rt2472 growth on agar plates was visible as rough colonies forming clumps, which significantly differed from the mucoid colonies formed by both Rt24.2 and Rt2472(pRC24). Similarly, in the liquid cultures, Rt2472 showed a tendency to aggregate. It was found that the mutant cells aggregated more effectively than the wild-type cells (Rt2472 = 36.69 % ± 2.91, Rt24.2 = 9.74 % ± 0.78, Rt2472(pRC24) = 10.10 % ± 0.85). A more detailed structure of these clumps was established in fluorescent microscopy using gfp-tagged strains (Fig. 4a, Additional file 5). In cultures of Rt24.2 and Rt2472 complementant, single and non-aggregated cells were visible. In contrast, the strain with the rosR mutation formed large aggregates with many grouped cells, which considerably outnumbered the single cells present in the culture. Moreover, large amounts of extracellular matrix surrounding the cells inside the clumps were observed. This phenotype could be explained by the increased expression of genes coding for the agglutinin RapA1, PrsD of the Type I secretion system, as well as glycosyl transferases involved in the synthesis of surface polysaccharides in the mutant cells. Previous studies had reported that RapA1 is a calcium-binding, cell-surface-associated agglutinin, whose function is to adhere to the root hairs . Moreover, a positive correlation between bacterial aggregation and biofilm formation had been observed, in which polysaccharides play a significant role [19, 57].
Apart from the differences described above, the rosR mutant showed significantly slower migration on both semi-solid and solid media in relation to both Rt24.2 and Rt2472(pRC24) (Fig. 4b, Additional file 6). This difference was especially visible on 0.3 % agar plates after 21 days of incubation. In order to establish whether the motility of this mutant was altered by flagella (more precisely, their number and/or structure), cells of the Rt2472, Rt24.2 and Rt2472(pRC24) strains were observed in transmission electron microscopy. The wild-type cells were rod-shaped, with an average length of 1.716 ± 0.38 μm (n = 15) (Rt2472(pRC24) = 1.792 ± 0.41 μm), and were peritrichously flagellated with 2–5 flagella (Fig. 4c). Moreover, these cells were surrounded by a thick and light layer of extracellular material, which was probably composed of EPS (Additional file 7). In contrast, the Rt2472 cells were slightly shorter (1.56 ± 0.15 μm). Furthermore, these cells were surrounded by a thinner but more densely packed extracellular layer than in the wild-type. In addition, rosR mutant cells were frequently non-flagellated and sporadically contained 1-2 flagella which were significantly shorter than those formed by the wild-type cells. This suggested serious disturbances in flagellum formation and/or function in the Rt2472 strain.
The type of cell flagellation is a strain-specific property, as shown by the fact that R. leguminosarum bv. viciae VF39 and S. meliloti 1021 strains were peritrichously flagellated similarly to Rt24.2, whereas R. leguminosarum bv. viciae 3841 had subpolar flagella. Also, the number of flagella might be strain-specific. For example, strain Rl3841 had 1–3 plain subpolar flagella while strain VF39 had 4–7 peritrichous flagella . It had been established that mutations in certain flagellar genes significantly affected rhizobial motility [29, 58]. For instance, a mutation in flaA resulted in non-motile VF39 cells and extremely reduced the motility of Rl3841 cells, whereas mutations in flaB and flaC resulted in shorter flagellar filaments, which reduced swimming and swarming motility in the cells of both of these strains. However, mutations in other than motility-related genes could also affect this bacterial trait. For example, a fabF mutant of Rl3841, deficient in 27-hydroxyoctacosanoate-modified lipopolysaccharide, was impaired in motility, biofilm formation, and desiccation tolerance . Motility of rhizobial cells is also affected by growth conditions such as agar concentration, inoculum size, temperature, and carbon source .
Analysis of transcriptional fusions in the Rt2472 and Rt24.2 strains
Several R. leguminosarum genes with different expression levels in the wild-type vs. the mutant background (rapA1, prsD, pssY, crp1, celA, gelA, ndvA), as well as those, whose expression was not affected by rosR (plyA, rfuA, pssB, exoB, Rl3414, Rl3425, ghy, nodA, ndvB), were chosen as representative genes for the experiment designed to validate the data obtained from transcriptomic analysis. The genes in these two groups showed a wide range of promoter strengths (from low to very high promoter activities). For this analysis, transcriptional fusions containing the promoter regions of these genes upstream of a promoterless lacZ were constructed (Additional file 8). In addition, plasmids with transcriptional fusions of motility-related genes (mcpC, mcpD, flaA, visN and rem), kindly provided by Prof C. Yost from the University of Regina (Canada) and Prof. M. Hynes from the University of Calgary (Canada), were included in the experiment. All these plasmid fusions were introduced into both Rt2472 and Rt24.2, and β-galactosidase or β-glucuronidase activity was measured (Fig. 5). Based on this analysis, it was found that promoters of the ndvA, ndvB, prsD, pssB, celA, gelA, rem, visN, flaA, mcpD, and mcpC genes were strong, since their activities in the wild-type strain were higher than 1500 Miller units. In contrast, a few of the studied genes (pssY, exoB, crp1, rapA1, plyA, rfuA, nodA, Rl3414, Rl3425, and ghy) possessed weak promoters (values below 1000 Miller units). When β-galactosidase activities for the individual fusions were compared between Rt2472 and Rt24.2, significant differences in the levels of transcription were found for crp1, pssY, rapA1, gelA, prsD, ndvA, celA (Fig. 5a), flaA, and rem genes (Fig. 5b). Among them, the highest Rt2472/Rt24.2 ratio was found for the following transcriptional fusions: crp1-lacZ, pssY-lacZ, gelA-lacZ, rapA1-lacZ, and flaA-gusA. On the other hand, similar levels of enzymatic activity in both of the strains tested were established for rfuA, ndvB, plyA, pssB, nodA, Rl3425, Rl3414, and ghy, indicating that RosR did not influence the expression of these genes.
In summary, the results obtained from the analysis of transcriptional fusions confirm that RosR is engaged in the regulation of the expression of rapA1, prsD, flaA, crp1, as well as ndvA, gelA, celA, and pssY involved in the synthesis of various rhizobial PS and other cell-surface components. These data are in conformity with those obtained in the RNA-Seq analysis. Previously, Bittinger and Handelsman, in their study of R. etli CE3, found that a mutation in a gene homologous to R. leguminosarum bv. trifolii rosR also affected the expression of several genes, including those involved in cellular metabolism (e.g. genes encoding coproporphyrinogen III oxidase, trehalose-phosphate synthase, and repressor of sucrose degradation operon) as well as those engaged in the synthesis and modification of polysaccharides (exoY, prsD, genes coding for cellulose synthase, CelR2 regulator of cellulose synthesis, and O-antigen acetylase) . RosR of R. etli, similarly to R. leguminosarum RosR, functions as a negative regulator of gene expression, and this protein proved to be essential for nodulation competitiveness of R. etli on bean plants .
Identification of RosR-box motifs in promoter regions of the genes differentially expressed in the rosR mutant and the wild-type strain
In order to establish which of the large group of genes differentially expressed in the mutant Rt2472 and the wild-type strain Rt24.2 were directly regulated by RosR, additional bioinformatics analyses were performed using a 22-nt-long RosR-box sequence and promoter regions of these genes. Sequence motifs showing a significant similarity to the RosR-box were found in the regulatory regions of many of these genes (195) (Additional file 8). Selected genes from this group and their RosR-box motifs are shown in Table 1. Among them, there were genes related to bacterial chemotaxis and motility (e.g. Rt792_47, Rt792_70, Rt628_20 (cheD), Rt634_26 (flbT), and Rt628_11 (fliN)), transport and metabolism of carbon and nitrogen sources (e.g. Rt620_74, Rt794_99, Rt617_43, Rt617_78, Rt619_162, and Rt794_91), and synthesis of cell-surface components and envelope biogenesis (e.g. Rt620_1 (pssY), Rt620_4, Rt780_45, Rt627_76 (gelA), Rt772_17 (prsD)). Moreover, several genes which encode proteins engaged in signal transduction and transcriptional regulators belonging to various families (Crp/Fnr, TetR, LuxR, LysR, AraC, and GntR) also possessed RosR-box motifs in their promoter regions. In this group, there were, among others, such genes as Rt620_2 (crp), Rt648_25 (lacI), Rt634_21 (rem), Rt782_109 (p-II), Rt646_16 (raiI), Rt648_50 (phoB) and Rt628_39 (phoR), which are involved in carbon, phosphate and nitrogen regulation, as well as motility and quorum sensing in R. leguminosarum [27, 28, 40, 49–51, 61]. These data indicate that RosR possibly controls a large group of these genes directly by binding to the sequence motifs located in their regulatory regions. However, in our study, there were many genes in which we did not identify motifs resembling the RosR-box; the expression of these genes may be indirectly regulated by RosR (probably mediated by the action of some transcriptional regulators encoded by COG K genes). Among the transcriptional regulators of R. leguminosarum mentioned above, the activator Rem, which together with LuxR-type regulators VisN and VisR positively controls the expression of many motility and chemotaxis genes, seems to be an interesting example of this relationship [49–51, 62]. Besides direct regulation of the transcription of several COG N genes, RosR was also found to indirectly affect the expression of some motility genes by modulating rem expression. Thus, RosR could be an important element of a second regulatory system controlling the expression of motility genes that are not under the control of the VisN/VisR-Rem regulatory cascade. The existence of such a system in R. leguminosarum was previously postulated by Tambalo and colleagues .
Another possibility is that RosR regulates some genes of the RosR regulon, for which RosR-box motifs were not found in our study, by acting together with another regulatory protein, as evidenced for the RosR homolog in S. meliloti, MucR . MucR single-handedly regulates the expression of many genes related to the synthesis of succinoglycan and the Nod factor, cell motility, and nitrogen fixation (e.g. exoH, ExoX, exoYFQ, exoK, nodD, and fix genes), and in a complex with WggR (ExpG) is involved in the regulation of galactoglucan synthesis genes [62–65].
In summary, data presented in this study suggest that RosR, similarly to S. meliloti MucR, plays an essential role in the regulatory network in R. leguminosarum, linking regulatory cascades associated with various environmental factors.
In this study, we performed transcriptome profiling of the rosR mutant of R. leguminosarum bv. trifolii strain 24.2 using the next-generation RNA-Seq technology. This analysis allowed us to identify a large set of genes linked to motility, synthesis of cell-surface components, and several metabolic pathways, whose expression was affected by the rosR mutation. It was established that RosR is a global transcriptional regulator, which functions mainly as a repressor and is involved in several cellular processes. On the other hand, this protein is engaged in the positive regulation of many motility-related genes. The data presented in this study expand our knowledge of the role of RosR in the functioning of R. leguminosarum cells and provide insight into the regulatory network of this bacterium. Moreover, our study confirms the observation that the rosR mutation significantly affects the behavior of bacterial cells, as evidenced by their altered morphology, flagellation, motility, and aggregation ability.
Bacterial strains, plasmids, and culture conditions
The bacterial strains and plasmids used in this study are listed in Additional file 9. R. leguminosarum strains were grown in an energy-rich 79CA medium with 1 % glycerol as a carbon source and TY medium at 28 °C on a rotary shaker (160 rpm) as described previously . E. coli strains were grown in Luria-Bertani (LB) medium at 37 °C . When required, antibiotics were used at the following final concentrations: kanamycin, 40 μg ml−1; rifampicin, 40 μg ml−1; ampicillin, 100 μg ml−1; gentamicin, 10 μg ml−1; tetracycline, 10 or 5 μg ml−1; and nalidixic acid, 40 μg ml−1. The insertion of an antibiotic resistance cassette in the Rt2472 genome was stable even in the absence of kanamycin in the medium. Therefore, in order to exclude the influence of antibiotics on bacterial metabolism, they were only used for the growth and storage of the mutant strain on agar plates. Cultures of Rt2472, for a significant majority of the experiments (RNA isolation, motility and aggregation assays, determination of generation time, and flagella visualization) did not contain kanamycin. Tetracycline was used at a concentration of 10 μg ml−1 in two experiments (visualization of aggregation of the Rt24.2, Rt2472 and Rt2472(pBR24) strains carrying a pHC60 vector in fluorescent microscope, and β-galactosidase (β-glucuronidase) assays in the Rt24.2 and Rt2472 strains carrying transcriptional fusions). In the case of Rt2472(pRC24), the dose of tetracycline was reduced to a final concentration of 5 μg ml−1.
Standard molecular techniques were used for DNA isolation, cloning, PCR reactions, ligation, transformation, and sequencing . For PCR amplifications, the REDTaq Ready PCR Reaction Mix (Sigma, USA) was used. PCR products cloned into plasmid vectors were verified by sequencing using the BigDye terminator cycle sequencing kit (Applied Biosystems) and the ABI Prism 310 sequencer. Sequence analyses were done with the FASTA and BLAST programs from the National Center for Biotechnology Information (Bethesda, MD, USA) and the European Bioinformatic Institute (Hinxton, UK), respectively. Promoter prediction in upstream regions of the analyzed genes was done using BDGP Neural Network Promoter Prediction , whereas identification of promoter motifs was performed with the Malign program using S. meliloti CTTGAC-N17-18-CTATAT and Escherichia coli TTGACA-N17-18-TATAAT promoter consensuses [68, 69]. A RosR-box motif (5’-CGGAATCTAGGGGTGGATTTCG-3’), as a query, and the Rt24.2 genomic sequence were used to search for the genes regulated by RosR [68, 70]. Sequence motifs with a high similarity to the RosR-box sequence were identified using position weight matrices (PWM) and subsequently were annotated to the genes. The motifs located within a distance of 300-nt upstream to 150-nt downstream of the translation start site of the genes were selected. All calculations were done using the Biopython motif package and custom python scripts [71, 72].
Total RNA isolation
For total RNA isolation, 25-ml cultures of Rt24.2 and Rt2472 were grown for 24 hours in 79CA medium in a rotary shaker. Then, the cultures were centrifuged for 10 min at 12,000 × g, and the bacterial pellet was suspended in 15 ml Trizol, shaken vigorously for 30 seconds, and incubated for 5 min at room temperature. After this period, 0.2 ml of chloroform was added per each 1 ml of Trizol used. The mixtures were shaken vigorously for 15 seconds, incubated for 8 minutes at room temperature, and subsequently centrifuged for 15 min at 12,000 × g at 4 °C. The RNA present in the water phase (12 ml) was transferred to a new tube and precipitated using 6 ml of isopropanol. Then, the mixtures were incubated at room temperature for 15 min and centrifuged for 15 min at 12,000 × g at 4 °C. The RNA pellets obtained were washed twice with 1 ml 75 % ethanol, mixed using a vortex mixer, centrifuged for 5 min at 7,500 × g at 4 °C, and dried in air for 10 min. Subsequently, the RNA pellets were dissolved in deionized RNase-free water and incubated for 10 min at 55 °C. The concentration of RNA isolated from the Rt24.2 and Rt2472 strains and their quality were determined using a NanoDrop 2000 spectrophotometer (Thermo Scientific, Wilmington, USA). Traces of DNA from the RNA solution were removed using a TURBO DNA-free Kit (Ambion, Life Technologies) according to the manufacturer’s instruction. Possible contamination of RNA by DNA was checked in PCR reactions using 2 primer pairs complementary to pssA (pssAG1f and pssA2r) and pssY (pssY5f and pssY5r) encoding glycosyl transferases in R. leguminosarum bv. trifolii as described earlier (Additional file 10) . Then, rRNA from the total RNAs isolated from the Rt24.2 and Rt2472 strains was removed using a Ribo-Zero Magnetic Kit for Gram-Negative bacteria, containing magnetic beads (Epicentre, Illumina). The rRNA-depleted mRNA was precipitated using ice-cold ethanol at a ratio 3:1 (v/v, ethanol : RNA solution). The sample was gently mixed and incubated at −20 °C for an hour and then centrifuged for 30 minutes at 12,000 × g. The pellet was washed twice using ice-cold 70 % ethanol and centrifuged at 12,000 × g for 5 minutes. The mRNA pellet obtained was dissolved in deionized RNase- and DNase-free water, and the amount of RNA was quantified using a NanoDrop 2000 Spectrophotometer. RNA integrity was assessed using an Agilent Bioanalyzer 2100 and an RNA 6000 Pico Kit (Agilent Technologies).
DNA sequencing and mapping, and analysis of RNA-Seq data
A draft genome of the Rt24.2 strain was sequenced using an Illumina MiSeq run (TruSeq DNA Technology) with paired-end reads. The initial preparation of the reads for analysis, including elimination of adapters and filtration of low quality reads, was performed using Cutadapt software . De novo assembly of the Rt24.2 sequence was performed with the MaSurCA genome assembler using the default values for bacterial genomes . The available switches were set for the assembly of bacterial genomes according to the software manual and the following parameters: ovlMerSize = 30, cgwErrorRate = 0.25, ovlMemory = 4GB, JF_SIZE = 100,000,000. As a result, the genomic sequence of the Rt24.2 strain was obtained in the form of 181 contigs with a total length of 7,653,217 bp (GC content = 60.6 %, average coverage 78.26×). In these sequences, coding regions were identified using the Prodigal program . In total, 7,374 putative coding regions were found, for which amino acid sequences were obtained and, subsequently, matched to the highly homologous proteins available in databases using the Blastp program .
The mRNAs of Rt24.2 and Rt2472 were sequenced in 250-bp cycles using the MiSeq System with SBS technology (Illumina). Three independent biological experiments were performed for each of the two strains. Image analyses and basecalling were performed with the use of MiSeq Control Software (MCS 188.8.131.52) and an RTA component (RTA 1.18.54) (Illumina). The RNA-Seq data for Rt24.2 and Rt2472 were aligned with Bowtie2  (within the Tophat package) using the reference Rt24.2 genome. Mapping parameters were set so that account could be taken of the strand-specific library option (fr-firststrand). The detection of novel splicing forms was disabled (option --no-novel-juncs). The median read number per CDS was above 1,000 (log value >3). For identification of genes with statistically significant differences in the expression levels between the Rt24.2 and Rt2472 strains, the Cufflinks\Cuffdiff software significance threshold value was set to α = 0.05 (using the Benjamini-Hochberg False Discovery Rate (FDR) correction)  (http://compbio.mit.edu/cummeRbund/manual_2_0.html). CDSs with FDR-corrected p-values for different expression between the rosR mutant and the wild-type lower than α = 0.05 were considered significant. As a result of the analysis, a list of genes differentially expressed in the tested strains with fold changes in the mutant strain versus the wild-type strain was obtained; the list also contained normalized expression FPKM (Fragments Per Kilobase of transcript per Million fragments mapped) values along with the Cuffdiff statistic test data.
To establish data quality and compare the expression levels between the Rt24.2 and Rt2472 strains, additional analyses were performed using cummeRbund (strain-specific repetition similarity calculated using a dendrogram with the Jensen-Shannon distance, distribution of FPKM values between the strains and between repetitions for an individual strain tested).
Functional categorization of genes
The Clusters of Orthologous Groups database was used to classify genes differentially expressed in Rt2472 and Rt24.2 into functional categories (COG) .
Construction of plasmids bearing transcriptional fusions with a reporter lacZ gene
In order to construct plasmids containing promoter regions of the analyzed genes cloned upstream of a promoterless lacZ gene, the broad-host-range vector pMP220 was used . Using a genomic DNA of the Rt24.2 strain and 14 primer pairs (Additional file 10), PCR products encompassing the upstream regions of the crp1, plyA, rfuA, pssB, ndvA, pssY, celA, gelA, rapA1, prsD, exoB, Rl3414, Rl3425, and ghy genes were obtained. Subsequently, these amplicons were digested with appropriate restriction enzymes and cloned into the respective sites in the pUC19 (for the plyA, rfuA, pssB, ndvA, pssY, celA, prsD, exoB, Rl3425, and ghy promoters) or in the pMP220 plasmid (for the crp1, gelA, rapA1, and Rl3414 promoters). The inserts of the latter plasmid group were re-cloned into the pUC19 vector. The nucleotide sequences of all these PCR products were verified by sequencing. Then, pUC19-derived fragments were re-cloned into the corresponding sites in pMP220, yielding plasmids which contained promoter fragments of the analyzed genes located upstream of the lacZ gene (Additional file 9). The resulting transcriptional fusions were introduced into E. coli S17-1 [81, 82] and subsequently into the Rt24.2 and Rt2472 strains by bi-parental conjugation. For this purpose, 24-h cultures of E. coli S17-1 derivatives carrying fusions (used as donor strains) and Rt24.2 and Rt2472 (used as recipient strains) were mixed at a 1:10 (v/v) ratio and centrifuged at 6,000 × g for 10 min. The bacterial pellet was suspended in 1 ml of sterilized water, and the mixture obtained was centrifuged for 10 min at 6,000 × g. Then, the pellet was suspended in 0.2 ml of sterilized water, and the bacterial mixture was placed on 79CA agar plates and incubated for 48 h at 28 °C. After this time, the bacteria were collected from the plates to 1 ml of sterilized water, and 0.1-ml aliquots were spread on 79CA agar plates supplemented with rifampicin and tetracycline. The transconjugants were obtained after one-week incubation at 28 °C.
β-galactosidase and β-glucuronidase assays
The Rt24.2 and Rt2472 strains containing the transcriptional fusions were grown for 24 h in 79CA medium supplemented with tetracycline. The assays for β-galactosidase activity were carried out according to the protocol described by Miller  using 2-Nitrophenyl-β-D-galactopyranoside as a substrate (Sigma, USA). In the case of transcriptional fusions containing the reporter gene gusA, β-glucuronidase activity was determined using p-Nitrophenyl-β-D-glucuronide (Sigma, USA) . The values reported are given in Miller units and are averages of at least three independent experiments.
The ability of Rt24.2, Rt2472 and Rt2472(pRC24) cells to move was determined in 79CA medium with 0.3 % and 0.7 % agar. Plates containing 20 ml of these semi-solid or solid media were air-dried for 24 h and then inoculated with 5 μl of 24-h bacterial cultures in liquid 79CA medium of OD600 = 0.4. The plates were incubated at 25 °C for 3 and 21 days, respectively. Afterwards, the diameter of the injection site was measured. The experiment was repeated twice with three biological repetitions per each strain tested.
Determination of generation time
In order to determine the generation time of Rt24.2, Rt2472 and Rt2472(pRC24) cells, 5 ml of liquid 79CA medium was inoculated with the bacteria to an optical density of OD600 = 0.15. Then, serial dilutions of the cultures were done after each 30 min during a 9-h period of growth. To this end, 100-μl aliquots of a particular dilution (from 10−2 to 10−7) were spread on 79CA agar plates. After three days of incubation at 28 °C, colonies were counted for each strain and each dilution tested. The number of bacterial cells present in the culture of an individual strain was calculated using the number of colonies × dilution factor × 10. The experiment was repeated twice with three biological repetitions per analyzed strain.
The autoaggregation assay for the rhizobial strains was performed according to the method of Sorroche et al.  with a slight modification. Briefly, Rt24.2, Rt2472 and Rt2472(pRC24) were cultured in 5 ml of liquid 79CA medium on a rotary shaker for 24 h at 28 °C. The cultures were subsequently left to settle for 24 h without agitation. Afterwards, a 0.3-ml sample was taken from the upper part of the suspension, and its optical density (OD600) was measured (A2). Then, the cultures were vigorously mixed for 30 s, and OD600 was measured (A1). The autoaggregation percentage was calculated as follows: 100 * [1 − (A2/A1)]. The experiment was repeated twice with three biological repetitions per each strain tested.
Transmission electron microscopy
Transmission electron microscopy was performed using a slightly modified procedure described by Tambalo et al. . The Rt24.2, Rt2472 and Rt2472(pRC24) strains were grown on both solid 79CA and TY agar plates. The investigated strains (especially Rt24.2 and Rt2472(pRC24)) produced significantly smaller amounts of extracellular polysaccharides on TY than on 79CA medium. A formvar copper grid was gently placed on the plate containing bacteria and then incubated in 4 % ammonium molybdate for 1.5 min in order to dye the bacteria. The cells of each individual strain were observed using a LEO912AB transmission electron microscope.
Broughton WJ, Zhang F, Perret X, Staehelin C. Signals exchanged between legumes and Rhizobium: agricultural uses and perspectives. Plant Soil. 2003;252:129–37.
Janczarek M, Rachwał K, Marzec A, Grządziel J, Palusińska-Szysz M. Signal molecules and cell-surface components involved in early stages of the legume–rhizobium interactions. Appl Soil Ecol. 2014;85:94–113.
Frey SD, Blum LK. Effect of pH on competition for nodule occupancy by type-i and type-ii strains of Rhizobium leguminosarum bv. phaseoli. Plant Soil. 1994;163:157–64.
Talbi C, Argandoña M, Salvador M, Alché JD, Vargas C, Bedmar EJ, et al. Burkholderia phymatum improves salt tolerance of symbiotic nitrogen fixation in Phaseolus vulgaris. Plant Soil. 2013;367:673–85.
Suzuki N, Rivero RM, Shulaev V, Blumwald E, Mittler R. Abiotic and biotic stress combinations. New Phytol. 2014;203:32–43.
Maynaud G, Brunel B, Mornico D, Durot M, Severac D, Dubois E, et al. Genome-wide transcriptional responses of two metal-tolerant symbiotic Mesorhizobium isolates to zinc and cadmium exposure. BMC Genomics. 2013;14:292.
Barthelemy-Delaux C, Marburger D, Delaux PM, Conley S, Ané JM. Effect of drought on Bradyrhizobium japonicum populations in Midwest soils. Plant Soil. 2014;382:165–73.
Jaszek M, Janczarek M, Kuczyński K, Piersiak T, Grzywnowicz K. The response of the Rhizobium leguminosarum bv. trifolii wild-type and exopolysaccharide-deficient mutants to oxidative stress. Plant Soil. 2014;376:75–94.
López-Leal G, Tabche ML, Castillo-Ramírez S, Mendoza-Vargas A, Ramírez-Romero MA, Dávila G. RNA-Seq analysis of the multipartite genome of Rhizobium etli CE3 shows different replicon contributions under heat and saline shock. BMC Genomics. 2014;15:770.
Laus MC, Logman TJ, Lamers GE, Van Brussel AA, Carlson RW, Kijne JW. A novel polar surface polysaccharide from Rhizobium leguminosarum binds host plant lectin. Mol Microbiol. 2006;59:1704–13.
Downie JA. The roles of extracellular proteins, polysaccharides and signals in the interactions of rhizobia with legume roots. FEMS Microbiol Rev. 2010;34:150–70.
Janczarek M, Rachwał K, Cieśla J, Ginalska G, Bieganowski A. Production of exopolysaccharide by Rhizobium leguminosarum bv. trifolii and its role in bacterial attachment and surface properties. Plant Soil. 2015;388:211–27.
Bogino PC, Oliva Mde L, Sorroche FG, Giordano W. The role of bacterial biofilms and surface components in plant-bacterial associations. Int J Mol Sci. 2013;14:15838–59.
Williams A, Wilkinson A, Krehenbrink M, Russo DM, Zorreguieta M, Downie JA. Glucomannan-mediated attachment of Rhizobium leguminosarum to pea root hairs is required for competitive nodule infection. J Bacteriol. 2008;190:4706–15.
Vanderlinde EM, Yost CK. Mutation of the Sensor Kinase chvG in Rhizobium leguminosarum negatively impacts cellular metabolism, outer membrane stability, and symbiosis. J Bacteriol. 2012;194:768–77.
Russo DM, Abdian PL, Posadas DM, Williams A, Vozza N, Giordano W, et al. Lipopolysaccharide O-chain core region required for cellular cohesion and compaction of in vitro and root biofilms developed by Rhizobium leguminosarum. Appl Environ Microbiol. 2015;81:1013–23.
Sorroche FG, Spesia MB, Zorreguieta A, Giordano W. A positive correlation between bacterial autoaggregation and biofilm formation in native Sinorhizobium meliloti isolates from Argentina. Appl Environ Microbiol. 2012;78:4092–101.
Rinaudi LV, Giordano W. An integrated view of biofilm formation in rhizobia. FEMS Microbiol Lett. 2010;304:1–11.
Balsanelli E, de Baura VA, de Oliveira Pedrosa F, de Souza EM, Monteiro RA. Exopolysaccharide biosynthesis enables mature biofilm formation on abiotic surfaces by Herbaspirillum seropedicae. PLoS ONE 2014; doi:10.1371/journal.pone.0110392.
Choma A, Komaniecka I. Characterization of cyclic β-glucans of Bradyrhizobium by MALDI-TOF mass spectrometry. Carbohydr Res. 2011;346:1945–50.
Gay-Fraret J, Ardissone S, Kambara K, Broughton WJ, Deakin WJ, Le Quéré A. Cyclic-β-glucans of Rhizobium (Sinorhizobium) sp. strain NGR234 are required for hypo-osmotic adaptation, motility, and efficient symbiosis with host plants. FEMS Microbiol Lett. 2012;333:28–36.
Pérez-Mendoza D, Rodríguez-Carvajal MA, Romero-Jiménez L, de Araujo FG, Lloret J, Gallegos MT, et al. Novel mixed-linkage β-glucan activated by c-di-GMP in Sinorhizobium meliloti. Proc Natl Acad Sci U S A. 2015;112:757–65.
Russo DM, Williams A, Edwards A, Posadas DM, Finnie C, Dankert M, et al. Proteins exported via the PrsD-PrsE type I secretion system and the acidic exopolysaccharide are involved in biofilm formation by Rhizobium leguminosarum. J Bacteriol. 2006;188:4474–86.
Krehenbrink M, Downie JA. Identification of protein secretion systems and novel secreted proteins in Rhizobium leguminosarum bv. viciae. BMC Genomics. 2008;9:55.
Janczarek M, Rachwał K. Mutation in the pssA gene involved in exopolysaccharide synthesis leads to several physiological and symbiotic defects in Rhizobium leguminosarum bv. trifolii. Int J Mol Sci. 2013;14:23711–35.
Yost CK, Del Bel KL, Quandt J, Hynes MF. Rhizobium leguminosarum methyl-accepting chemotaxis protein genes are down-regulated in the pea nodule. Arch Microbiol. 2004;182:505–13.
Miller LD, Yost CK, Hynes MF, Alexandre G. The major chemotaxis gene cluster of Rhizobium leguminosarum bv. viciae is essential for competitive nodulation. Mol Microbiol. 2007;63:348–62.
Tambalo DD, Bustard DE, Del Bel KL, Koval SF, Khan MF, Hynes MF. Characterization and functional analysis of seven flagellin genes in Rhizobium leguminosarum. BMC Microbiol. 2010;10:219.
Tambalo DD, Yost CK, Hynes MF. Characterization of swarming motility in Rhizobium leguminosarum bv. viciae. FEMS Microbiol Lett. 2010;307:165–74.
Nogales J, Domínguez-Ferreras A, Amaya-Gómez CV, van Dillewijn P, Cuéllar V, Sanjuán J, et al. Transcriptome profiling of a Sinorhizobium meliloti fadD mutant reveals the role of rhizobactin 1021 biosynthesis and regulation genes in the control of swarming. BMC Genomics. 2010;11:157.
Vicario JC, Dardanelli MS, Giordano W. Swimming and swarming motility properties of peanut-nodulating rhizobia. FEMS Microbiol Lett. 2015;362:1–6.
Morita RY. Bioavailability of energy and its relationship to growth and starvation survival in nature. Can J Microbiol. 1988;34:436–44.
Walker TS, Bais HP, Grotewold E, Vivanco JM. Root exudation and rhizosphere biology. Plant Physiol. 2003;132:44–51.
Thorne SH, Williams HD. Adaptation to nutrient starvation in Rhizobium leguminosarum bv. phaseoli: analysis of survival, stress resistance, and changes in macromolecular synthesis during entry to and exit from stationary phase. J Bacteriol. 1997;179:6894–901.
Yost CK, Rochepeau P, Hynes MF. Rhizobium leguminosarum contains a group of genes that appear to code for methyl-accepting chemotaxis proteins. Microbiology. 1998;144:1945–56.
Liu R, Van Mai T, Schmidt EL. Nodulating competitiveness of a non-motile Tn7 mutant of Bradyrhizobium japonicum in nonsterile soil. Appl Environ Microbiol. 1989;55:1895–900.
Caetano-Anolles G, Crist-Estes DK, Bauer WD. Chemotaxis of Rhizobium meliloti to the plant flavone luteolin requires functional nodulation genes. J Bacteriol. 1988;170:3164–9.
Bauer WD, Caetano-Anollés G. Chemotaxis, induced gene expression and competitiveness in the rhizosphere. Plant Soil. 1990;129:45–52.
Tambalo DD, Vanderlinde EM, Robinson S, Halmillawewa A, Hynes MF, Yost CK. Legume seed extracts and Physcomitrella patens extracts influence swarming behaviour in Rhizobium leguminosarum. Can J Microbiol. 2014;60:15–24.
Janczarek M, Skorupska A. The Rhizobium leguminosarum bv. trifolii RosR: transcriptional regulator involved in exopolysaccharide production. Mol Plant Microbe Interact. 2007;20:867–81.
Janczarek M, Skorupska A. Rhizobium leguminosarum bv. trifolii rosR gene expression is regulated by catabolic repression. FEMS Microbiol Lett. 2009;291:112–9.
Janczarek M, Jaroszuk-Ściseł J, Skorupska A. Multiple copies of rosR and pssA genes enhance exopolysaccharide production, symbiotic competitiveness and clover nodulation in Rhizobium leguminosarum bv. trifolii. Antonie Van Leeuwenhoek. 2009;96:471–86.
D’Souza-Ault MR, Cooley MB, Kado CI. Analysis of the Ros repressor of Agrobacterium virC and virD operons: molecular intercommunication between plasmid and chromosomal genes. J Bacteriol. 1993;175:3486–90.
Keller M, Roxlau A, Weng WM, Schmidt M, Quandt J, Niehaus K, et al. Molecular analysis of the Rhizobium meliloti mucR gene regulating the biosynthesis of the exopolysaccharides succinoglycan and galactoglucan. Mol Plant Microbe Interact. 1995;8:267–77.
Bittinger MA, Milner JL, Saville BJ, Handelsman J. rosR, a determinant of nodulation competitiveness in Rhizobium etli. Mol Plant Microbe Interact. 1997;10:180–6.
Chou AY, Archdeacon J, Kado CI. Agrobacterium transcriptional regulator Ros is a prokaryotic zinc finger protein that regulates the plant oncogene ipt. Proc Natl Acad Sci U S A. 1998;95:5293–8.
Janczarek M, Kutkowska J, Piersiak T, Skorupska A. Rhizobium leguminosarum bv. trifolii rosR is required for interaction with clover, biofilm formation and adaptation to the environment. BMC Microbiol. 2010;10:284.
Tatusov RL, Natale DA, Garkavtsev IV, Tatusova TA, Shankavaram UT, Rao BS, et al. The COG database: new developments in phylogenetic classification of proteins from complete genomes. Nucleic Acids Res. 2001;29:22–8.
Santos-Beneit F. The Pho regulon: a huge regulatory network in bacteria. Front Microbiol. 2015;6:402.
Yurgel SN, Rice J, Kahn ML. Transcriptome analysis of the role of GlnD/GlnBK in nitrogen stress adaptation by Sinorhizobium meliloti Rm1021. PLoS ONE 2013; doi:10.1371/journal.pone.0058028.
Wisniewski-Dyé F, Jones J, Chhabra SR, Downie JA. raiIR genes are part of a quorum-sensing network controlled by CinI and CinR in Rhizobium leguminosarum. J Bacteriol. 2002;184:1597–606.
Vanderlinde EM, Hynes MF, Yost CK. Homoserine catabolism by Rhizobium leguminosarum bv. viciae 3841 requires a plasmid-borne gene cluster that also affects competitiveness for nodulation. Environ Microbiol. 2014;16:205–17.
Geddes BA, Oresnik IJ. Inability to catabolize galactose leads to increased ability to compete for nodule occupancy in Sinorhizobium meliloti. J Bacteriol. 2012;194:5044–53.
Król JE, Mazur A, Marczak M, Skorupska A. Syntenic arrangements of the surface polysaccharide biosynthesis genes in Rhizobium leguminosarum. Genomics. 2007;89:237–47.
Ausmees N, Jacobsson K, Lindberg M. A unipolarly located, cell-surface-associated agglutinin, RapA1, belongs to a family of Rhizobium-adhering proteins (Rap) in Rhizobium leguminosarum bv. trifolii. Microbiology. 2001;147:549–59.
Krysciak D, Grote J, Rodriguez Orbegoso M, Utpatel C, Förstner KU, Li L, et al. RNA sequencing analysis of the broad-host-range strain Sinorhizobium fredii NGR234 identifies a large set of genes linked to quorum sensing-dependent regulation in the background of a traI and ngrI deletion mutant. Appl Environ Microbiol. 2014;80:5655–71.
Gao M, Coggin A, Yagnik K, Teplitski M. Role of specific quorum-sensing signals in the regulation of exopolysaccharide II production within Sinorhizobium meliloti spreading colonies. PLoS ONE 2012; doi:10.1371/journal.pone.0042611.
Bahlawane C, McIntosh M, Krol E, Becker A. Sinorhizobium meliloti regulator MucR couples exopolysaccharide synthesis and motility. Mol Plant Microbe Interact. 2008;21:1498–509.
Vanderlinde EM, Muszyński A, Harrison JJ, Koval SF, Foreman DL, Ceri H, et al. Rhizobium leguminosarum biovar viciae 3841, deficient in 27-hydroxyoctacosanoate-modified lipopolysaccharide, is impaired in desiccation tolerance, biofilm formation and motility. Microbiology. 2009;155:3055–69.
Bittinger MA, Handelsman J. Identification of genes in the RosR regulon of Rhizobium etli. J Bacteriol. 2000;182:1706–13.
Tambalo DD, Del Bel KL, Bustard DE, Greenwood PR, Steedman AE, Hynes MF. Regulation of flagellar, motility and chemotaxis genes in Rhizobium leguminosarum by the VisN/R-Rem cascade. Microbiology. 2010;156:1673–85.
Sourjik V, Muschler P, Scharf B, Schmitt R. VisN and VisR are global regulators of chemotaxis, flagellar, and motility genes in Sinorhizobium (Rhizobium) meliloti. J Bacteriol. 2000;182:782–8.
Bertram-Drogatz PA, Quester I, Becker A, Pühler A. The Sinorhizobium meliloti MucR protein, which is essential for the production of high-molecular-weight succinoglycan exopolysaccharide, binds to short DNA regions upstream of exoH and exoY. Mol Gen Genet. 1998;257:433–41.
Rüberg S, Pühler A, Becker A. Biosynthesis of the exopolysaccharide galactoglucan in Sinorhizobium meliloti is subject to a complex control by the phosphate-dependent regulator PhoB and the proteins ExpG and MucR. Microbiology. 1999;145:603–11.
Mueller K, González JE. Complex regulation of symbiotic functions is coordinated by MucR and quorum sensing in Sinorhizobium meliloti. J Bacteriol. 2011;193:485–96.
Sambrook J, Fritsch EF, Maniatis T. Molecular Cloning: a laboratory manual. New York: Cold Spring Harbour Laboratory; 1989.
Berkeley Drosophila Genome Project (BDGP), Neural Network Promoter Prediction. http://www.fruitfly.org/seq_tools/promoter.html.
A.N. Belozersky Institute of Physico-Chemical Biology, GeneBee, http://www.genebee.msu.su/services/malign.
MacLellan SR, MacLean AM, Finan TM. Promoter prediction in the rhizobia. Microbiology. 2006;152:1751–63.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215:403–10.
Stormo GD. DNA binding sites: representation and discovery. Bioinformatics. 2000;16:16–23.
Cock PJA, Antao T, Chang JT, Chapman BA, Cox CJ, Dalke A, et al. Biopython: freely available Python tools for computational molecular biology and bioinformatics. Bioinformatics. 2009;25:1422–3.
Janczarek M, Kalita M, Skorupska A. New taxonomic markers for identification of Rhizobium leguminosarum and discrimination between closely related species. Arch Microbiol. 2009;191:207–19.
Marcell M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet journal. 2011;1:10–2. ISSN 2226-6089 [S.1.].
Zimin AV, Marçais G, Puiu D, Roberts M, Salzberg SL, Yorke JA. The MaSuRCA genome assembler. Bioinformatics. 2013;29:2669–77.
Hyatt D, Chen G, Locascio P, Land M, Larimer F, Hauser L. Prodigal: prokaryotic gene recognition and translation initiation site identification. BMC Bioinformatics. 2010;11(1):119.
Bowtie - An ultrafast memory-efficient short read aligner. John Hopkins University. http://bowtie-bio.sourceforge.net.
Trapnell C, Roberts A, Goff L, Pertea G, Kim D, Kelley DR, et al. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat Protoc. 2012;7:562–78.
Clusters of Orthologous Groups (COG). http://www.ncbi.nlm.nih.gov/COG/.
Spaink HP, Okker RJH, Wijffelman CA, Pees E, Lugtenberg BJJ. Promoters in the nodulation region of the Rhizobium leguminosarum Sym plasmid pRL1JI. Plant Mol Biol. 1987;9:29–37.
Simon R, Priefer U, Pühler A. A broad host range mobilization system for in vivo genetic engineering: transposon mutagenesis in Gram-negative bacteria. Biotechnology. 1983;1:784–91.
Maj D, Wielbo J, Marek-Kozaczuk M, Skorupska A. Response to flavonoids as a factor influencing competitiveness and symbiotic activity of Rhizobium leguminosarum. Microbiol Res. 2010;165:50–60.
Miller JH. Experiments in molecular genetics. Cold Spring Harbor, NY: Cold Spring Harbor Laboratory Press; 1972.
Cheng HP, Walker GC. Succinoglycan is required for initiation and elongation of the infection threads during nodulation of alfalfa by Rhizobium meliloti. J Bacteriol. 1998;180:5183–91.
This work was supported by the grant of the National Science Centre no 2012/07/B/NZ1/00099. We thank Prof C. Yost from the University of Regina (Canada) and Prof. M. Hynes from the University of Calgary (Canada) for providing transcriptional fusions of motility genes.
The authors declare that they have no competing interests.
KR carried out the experiments concerning the properties of the Rt24.2 and Rt2472 strains, participated in the molecular genetic studies, prepared the figures; EM carried out bioinformatics analyses of genes differentially expressed in Rt2472 versus Rt24.2 and participated in drafting the manuscript; MJ conceived and designed the experiments, participated in the molecular genetic studies, analyzed the data and wrote the manuscript. All the authors have read and approved the final version of the manuscript.
General features of the total sequenced and mapped reads for the R. leguminosarum bv. trifolii wild-type strain Rt24.2 and its derivative, the Rt2472 rosR mutant* (A); (B) Dendrogram showing the similarity of the biological repetitions for Rt24.2 and Rt2472 measured using the Jensen–Shannon distance. The samples for Rt24.2 are marked with wt_0, wt_1, and wt_2, whereas those for Rt2472 are marked with mut_0, mut_1, and mut_3; (C) Box plot of the distribution of FPKM values between three biological repetitions for the Rt24.2 and Rt2472 strains. (PDF 173 kb)
General features of FPKM values obtained for Rt24.2 and Rt2472 genes. (A) Histogram of FPKM values for Rt24.2 and Rt2472 genes. The three biological repetitions for Rt24.2 are marked with wt_0, wt_1, and wt_2, whereas those for Rt2472 are marked with mut_0, mut_1, and mut_3; (B) Point diagram showing FPKM values of the individual genes obtained from Rt24.2 and Rt2472; (C) Point diagram showing log2 fold change values and p-values for the individual genes differentially expressed in Rt2472 and Rt24.2. (PDF 327 kb)
A list of coding sequences in R. leguminosarum bv. trifolii, for which differences in the expression levels between the Rt2472 rosR mutant (mut) and the wild-type strain Rt24.2 (wt) were established (log 2 fold change >2), with classification of the sequences into COG groups and putative functions of the encoded proteins. (XLSX 148 kb)
List of putative CDS showing highest differences in the expression levels between the rosR mutant Rt2472 and the wild-type strain Rt24.2. (DOCX 25 kb)
Autoaggregation of the R. leguminosarum bv. trifolii wild-type Rt24.2 strain, the rosR mutant Rt2472, and the Rt2472 strain after complementation by rosR. (A) The autoaggregation ability of the analyzed strains determined after 24 h; (B) Cell morphology of the analyzed strains tagged with a gusA reporter gene (pHC60 plasmid)  in fluorescence microscopy (magnification, 200× and 630×). (PNG 547 kb)
Motility assay of R. leguminosarum bv. trifolii wild-type Rt24.2, rosR mutant Rt2472, and Rt2472(pRC24) cells on 0.3 % and 0.7 % 79CA agar plates in 3-day and 3-week experiments. (DOCX 13 kb)
The cells of R. leguminosarum bv. trifolii wild-type Rt24.2, the rosR mutant Rt2472, and Rt2472(pRC24) visualized in transmission electron microscopy. The strain Rt24.2 growing on 79CA (A) and TY medium (B); The strain Rt2472 growing on 79CA (C) and TY medium (D); The strain Rt2472(pRC24) growing on TY medium (E and F). Flagella are marked with red arrows and the extracellular layer surrounding bacterial cells is marked with a blue arrow. (PNG 1167 kb)
Sequence motifs identified in regulatory regions of the genes differentially expressed in the rosR mutant and the wild-type strain, which show a significant similarity to the RosR-box sequence. (XLSX 26 kb)
Bacterial strains and plasmids used in this study. (DOCX 16 kb)
Oligomeric primers used in this study. (DOCX 14 kb)
About this article
Cite this article
Rachwał, K., Matczyńska, E. & Janczarek, M. Transcriptome profiling of a Rhizobium leguminosarum bv. trifolii rosR mutant reveals the role of the transcriptional regulator RosR in motility, synthesis of cell-surface components, and other cellular processes. BMC Genomics 16, 1111 (2015). https://doi.org/10.1186/s12864-015-2332-4
- Transcriptional regulator RosR
- Rhizobium leguminosarum
- RNA sequencing
- Transcriptome analysis
- Cell-surface components
- Nitrogen-fixing symbiosis