Skip to main content
  • Research article
  • Open access
  • Published:

The time course of the transcriptomic response of Sinorhizobium meliloti1021 following a shift to acidic pH

Abstract

Background

The symbiotic soil bacterium Sinorhizobium meliloti often has to face low pH in its natural habitats. To identify genes responding to pH stress a global transcriptional analysis of S. meliloti strain 1021 following a pH shift from pH 7.0 to pH 5.75 was carried out. In detail, oligo-based whole genome microarrays were used in a time course experiment. The monitoring period covered a time span of about one hour after the pH shift. The obtained microarray data was filtered and grouped by K-means clustering in order to obtain groups of genes behaving similarly concerning their expression levels throughout the time course.

Results

The results display a versatile response of S. meliloti 1021 represented by distinct expression profiles of subsets of genes with functional relation. The eight generated clusters could be subdivided into a group of four clusters containing genes that were up-regulated and another group of four clusters containing genes that were down-regulated in response to the acidic pH shift. The respective mean expression progression of the four up-regulated clusters could be described as (i) permanently and strong, (ii) permanently and intermediate, (iii) permanently and progressive, and (iv) transiently up-regulated. The expression profile of the four down-regulated clusters could be characterized as (i) permanently, (ii) permanently and progressive, (iii) transiently, and (iv) ultra short down-regulated. Genes coding for proteins with functional relation were mostly cumulated in the same cluster, pointing to a characteristic expression profile for distinct cellular functions. Among the strongest up-regulated genes lpiA, degP1, cah, exoV and exoH were found. The most striking functional groups responding to the shift to acidic pH were genes of the exopolysaccharide I biosynthesis as well as flagellar and chemotaxis genes. While the genes of the exopolysaccharide I biosynthesis (exoY, exoQ, exoW, exoV, exoT, exoH, exoK exoL, exoO, exoN, exoP) were up-regulated, the expression level of the flagellar and chemotaxis genes (visR, motA, flgF, flgB, flgC, fliE, flgG, flgE, flgL, flbT, mcpU) simultaneously decreased in response to acidic pH. Other responding functional groups of genes mainly belonged to nitrogen uptake and metabolism (amtB, nrtB, nirB, nirD), methionine metabolism (metA, metF, metH, metK, bmt and ahcY) as well as ion transport systems (sitABCD, phoCD). It is noteworthy, that several genes coding for hypothetical proteins of unknown function could be identified as up-regulated in response to the pH shift.

Conclusion

It was shown that the short term response to acidic pH stress does not result in a simple induction or repression of genes, but in a sequence of responses varying in their intensity over time. Obviously, the response to acidic pH is not based on a few specific genes, but involves whole sets of genes associated with various cellular functions.

Background

The symbiotic interaction between rhizobia and leguminous plants plays an important role in global nitrogen fixation. During symbiosis rhizobia colonize the root nodules and induce nodule formation. Rhizobia in turn differentiate into bacteroids and live as endosymbionts inside plant cells. They fix atmospheric nitrogen and provide the fixed nitrogen to the host plant. The efficiency of this symbiosis is constrained by several factors relating to the soil and the rate of nodulation and nitrogen fixation is diminished. The most commonly observed factors are water deficiency, high temperature, high salt content and low pH (for review see [1]). At acidic pH conditions the bacterial partner is limited in survival and persistence and the nodulation efficiency is reduced [2–4]. Another situation where rhizobia are commonly facing a low pH environment is the rhizoplane of their leguminous host plants, where the pH is decreased by protons and organic acids excreted by the plants [5]. Once a symbiosis has been established the symbiosome has been postulated to form an acidic and lytic compartment [6]. Several research groups have been trying to identify pH tolerant strains [3, 7] and to reveal the genetic mechanisms enabling those strains to outperform other strains in low pH soils, however up until now the basis of the rhizobial pH tolerance remains unknown.

Since the genome of S. meliloti 1021 is well characterised [8–11]S. meliloti 1021 is considered to represent an ideal candidate to analyse its behaviour under environmental conditions. The response to several stresses has been studied for individual genes of S. meliloti on a proteomic as well as a transcriptomic scale [12–15]. But the cellular response of S. meliloti to acid stress has so far not been investigated on a genome-wide level. pH stress can affect cells in several ways and therefore different responses exist. Acid tolerance in general is a mechanism of the cell to face an unfavourable acidic condition, whereas an adaptive acid tolerance (ATR) is defined as increased tolerance against low pH after growing cells in moderately low pH media [16] (for review see [17]). For rhizobia most studies about genes involving the acid stress response have been conducted with S. medicae (formerly classified as S. meliloti WSM 419). By using a transposon mutagenesis system [18] a functionally diverse set of pH responsive and acid tolerance related genes could be identified [19]. Gene products required for acid tolerance in S. medicae are for example ActP, a CPx heavy metal transporting ATPase [20], and ActA, an apolipoprotein acyl transferase [21]. A gene coding for a regulatory protein known to be required for the acid tolerance in S. medicae is actR [22]. The encoded response regulator ActR is activated by its corresponding sensor histidine kinase ActS, whose loss also leads to sensitivity to low pH. The cbbS gene involved in CO2 fixation and the narB gene involved in nitrate assimilation as well as the nitrogen fixation regulator genes fixK and nifA could be identified as target genes for the regulator ActR [23]. Along with the genes required for low pH tolerance some further genes up-regulated by low pH were identified for S. medicae [19, 24]. Among these was lpiA, a gene found to be necessary for the adaptive acid tolerance (ATR). In Rhizobium tropici, the bacterial symbiont of Phaseolus vulgaris, this gene was also up-regulated by low pH and was found to be necessary for an increased nodulation competitiveness [25].

In this study the transcriptional response of S. meliloti strain 1021 following a pH shift from pH 7.0 to pH 5.75 was analysed on a genome wide level. Using whole-genome Sm6kOligo microarrays [15] the expression of S. meliloti genes responding to this environmental change was monitored over a period of one hour. The data obtained was filtered and clustered to obtain groups of genes with a similar behaviour.

Results and Discussion

Growth analysis of S. meliloti1021 cultures exposed to neutral and acidic pH

The aim of this study was to analyse the transcriptional response of S. meliloti 1021 following a shift from a neutral to an acidic pH. Since adaptation to new environmental conditions means passing through an evolving process of cellular responses until reaching a steady state balance, it was decided to monitor the transcriptional response over a certain period of time. One critical point concerns the correct choice of parameters for the pH shift. The pH stress should be applied to S. meliloti in a manner that it is not lethal for the cells, the central metabolism should also not be turned down on a broad scale, otherwise a global stress and not a pH-specific stress response is to be expected. In order to define appropriate experimental conditions for the pH shift, growth tests in Vincent minimal medium were carried out by varying the pH from 5.5 to 7.0 in 0.25 increments. It turned out that S. meliloti 1021 is not able to grow at pH 5.5 while above pH 6.0 only minor deviations from the growth curve at pH 7.0 occurred (data not shown). At pH 5.75 S. meliloti 1021 showed a reduced growth rate, but the cell titer counts documented that this pH was not yet lethal (data not shown).

The aim of this study was to identify genes of S. meliloti that directly respond to changes of the environmental pH, the transcriptional short term response within the first hour after a pH change was therefore the focus of our interest. In a time course experiment the global gene expression of S. meliloti cells exposed to a pH change from 7.0 to 5.75 was compared to the gene expression of untreated cells. To ensure identical conditions and treatment S. meliloti 1021 cells were grown in VMM at pH 7.0 until an o.D.580 of 0.8 was reached (Fig. 1), subsequently the culture was split in two and centrifuged. After centrifugation of the split cultures, the used growth medium was decanted and exchanged by fresh VMM adjusted to pH 5.75 (as testing condition) and to pH 7.0 (as reference), respectively. All manipulation steps were carried out very gently by using pre-warmed equipment and material to avoid any unwanted influences on the cells. The growth curves show the effect of the lowered pH on the growth of the S. meliloti 1021 culture (Fig. 1). The culture that was shifted to pH 5.75 grew slower than the pH 7.0 culture. For the duration of the time course experiment, the pH value of both cultures did not change. At later time points an alkalisation of the growth medium could be observed for the low pH culture (data not shown).

Figure 1
figure 1

Growth of S. meliloti 1021 before and after a shift to low pH. An S. meliloti 1021 preculture has been grown in VMM buffered at pH 7.0 until it reached an o.D.580 of 0.8 (dotted line with triangles). Afterwards the pre-culture has been separated into even parts, centrifuged and re-suspended in VMM at pH 5.75 and VMM at pH 7.0, respectively. The growth of the pH 5.75 culture is given by lines with crosses and the growth of the pH 7.0 culture is given by lines with plus-symbols. The arrows in the diagram indicate the time points where cell culture probes were taken for transcriptional profiling. Remarks indicate the time in minutes passed after the splitting of the S. meliloti preculture.

Cluster analysis of expression profiles of S. melilotigenes following a shift to acidic pH

Cells were harvested from both cultures grown at pH 7.0 and pH 5.75 after 3, 8, 13, 18, 33 and 63 minutes (Fig. 1). Because both the sample (pH 5.75) and control (pH 7.0) were obtained from the same starting culture and were treated simultaneously, the differences in the transcriptomic profile should only be a result of the pH shift. For each time point the transcribed and labelled RNA of the pH 5.75 grown culture was hybridised together with the differently labelled RNA of the pH 7.0 reference culture to the Sm6kOligo microarray. The whole procedure was performed in three biological replicates to ensure the validity of the microarray data. The microarray images were analysed using the Imagene Software and EMMA [26] (For microarray data see: http://www.ebi.ac.uk/microarray-as/ae/).

As expected, the microarray analysis for the six successive time points revealed a high number of genes with different expression characteristics over the tested period. In order to identify genes that presumably play a significant role in the cellular response to acidic pH the following filtering criteria were applied. Only genes with a log2 fold difference in spot intensities on the microarray slides (M value) of ≥ 2 or ≤ -2 were considered. Because we were also interested in genes that were only transiently active, this limit of significance had to be achieved for at least one time point during the time series. In addition, it was of importance for clustering that each gene was represented with an evaluable expression value (R ≥ 1.5 for both channels) for at least 5 out of the 6 time points. 230 genes fulfilled these filtering criteria. To estimate the number of false positive genes after filtering the false discovery rate (FDR) control was applied for all expression data of these 230 genes. The FDR control revealed a proportion of less than 1% false positives. Additionally, the tendency of the microarray results was confirmed by qRT-PCR for two of the obtained genes (lpiA and phoC) (data not shown).

Since the S. meliloti genome is composed of three replicons with distinctive functional features [8] the distribution of the 230 genes fulfilling the filtering criteria was determined. The percentage of differentially expressed genes of the total number of genes was 3.95% for the chromosome, 2.48% for pSymA, and 4.20% for pSymB. Therefore, compared to the chromosome genes located on pSymB were slightly over represented whereas genes of pSymA were noticeably under represented in the time course experiment. A possible explanation is that pSymA carries mostly symbiosis related genes which are not responding, whereas pSymB and the chromosome contain housekeeping genes. The slight over representation of pSymB might be based on the up-regulation of exopolysaccharide biosynthesis genes (see below).

In the next step, a clustering of these genes was performed according to their expressional characteristics over time. By hierarchical clustering, a separation into eight different clusters was estimated. Since K-means clustering offers the possibility to group complex time course expression data into a pre selected number of distinct clusters according to their similarity, this method was applied. Eight K-means clusters (see additional files 1, 2, 3, 4, 5, 6, 7 and 8: Heat maps of the generated Clusters A to H; additional file 9: combined spread sheet of the clustered genes) was found to be the smallest number resulting in clusters with clearly distinguishable expression characteristics. These expression characteristics become apparent by calculating the arithmetic mean expression profile of all contained genes (Fig. 2).

Figure 2
figure 2

The eight clusters of the transcriptomic profiling of S. meliloti 1021 following a shift to acidic pH. The diagrams of clusters A to H show the mean M value (y-axis) obtained by the Sm6kOligo microarray analyses for each time point (x-axis) after pH shift. The standard deviations are represented by the vertical lines crossing each point of the graph. The dotted line divides the cellular response into two parts. The location of the dotted line was chosen according to the observation that most of the cellular response happened in the first 20 minutes.

The cluster analysis generated different groups with discrete expression profiles for up-regulated genes (cluster A to D in Fig. 2) and down-regulated genes (cluster E to H in Fig. 2), not only differing in the intensity of their expression level, but also in their time dependent expression behaviour. These time dependent behaviours can be roughly separated into clusters containing genes that were permanently differentially expressed and clusters containing genes which were only transiently differentially expressed (Fig. 3). An inspection of individual expression profiles (data not shown) indicated that for borderline cases the passage between clusters is fluent. The mean expression profiles of the clusters indicated that the main changes in response to low pH occurred approximately within the first 20 minutes (Fig. 2). After this period of time a constant differential expression level or a constantly changing differential expression level can be observed for most clusters. It is also noticeable that several clusters contain genes organised in operons and groups of genes belonging to related or similar cellular functions.

Figure 3
figure 3

Grouping of S. meliloti 1021 genes following a shift to acidic pH. The eight calculated gene clusters were characterised by their specific transcriptomic response. The figure shows the classification of the genes chosen for clustering by single attributes into the eight clusters calculated by K-means. The tables below give the gene names of the genes distributed to the corresponding clusters.

In cluster A, genes exhibiting a strong and permanent induction accumulated. Genes in this cluster remained up-regulated for the whole observation period. It therefore seems that these genes have a special impact for S. meliloti in facing low pH conditions.

Clusters B also contains genes that remained permanently up-regulated in response to the pH shift, but not as strong as those in cluster A.

Cluster C consists of genes that showed a progressive permanent induction in their mean expression behaviour. The expression level of these genes is continuously increasing for the duration of the experiment. Presumably the expression of these genes will reach a plateau phase at later time points.

The mean expression curve of cluster D shows that the genes therein were only transiently up-regulated during the first 10 to 30 minutes following the pH shift. This expressional characteristic suggests that the encoded functions of these genes were only needed for a short period of time.

Cluster E is composed of genes that showed a decrease in their expression up to 20 minutes after the pH shift and thereafter they remained permanently down-regulated.

In contrast, genes contained in cluster F had a progressive permanent repression for the whole duration of the experiment. Similar to cluster C a steady state can be expected for later points of time.

Cluster G consists of genes that were transiently down-regulated in their expression level with a minimum occurring within 20 minutes after the pH shift. In contrast to the genes of cluster D, the encoded functions of the genes in cluster G are likely to be temporarily not needed for the cell.

Cluster H represents the smallest cluster and consists of genes with an ultra short transient repression observed for the first time point of 3 minutes after pH shift.

Most of cluster A genes encode proteins carrying signal peptides

Cluster A contains 16 genes that showed an increasing induction in response to the pH shift up to 18 minutes and maintained constantly high expression values thereafter (Fig. 2A). It is of special interest that 9 of these 16 genes encode products carrying signal peptides for secretion. This is remarkable since this is by far the highest percentage of putatively secreted proteins for all clusters analysed. It is therefore possible that one of the most immediate and strongest responses of S. meliloti to face acidic pH is the secretion of proteins. Five of these genes code for hypothetical proteins (smb21025, smb21026, smb21440, smc01580 and smc01774). Two of the genes encode for putatively secreted lytic enzymes like a protease (degP1) and a putative lysozyme (smc01855). An orthologous protein of the degP1 gene product could also be identified in S. medicae following the growth for 5 days at pH 5.7 [27]. Another interesting gene pair in cluster A namely, smc02366 and smc02367 coding for a two component system, was found to be located downstream of degP1. Whether this gene pair is involved in the transcriptional regulation of degP1 has to be investigated. Due to its location next to the highly expressed degP1 gene a polar effect influencing the smc02366-smc02367 expression can not be excluded.

Among cluster A smc00611 represents one of the highest up-regulated genes during the time course. An orthologue in S. medicae termed lpiA, is already well known as a pH responsive gene and it could be shown that it is necessary for the adaptation to low pH [19, 28]. For R. tropici it was demonstrated that lpiA is pH responsive and symbiotically relevant [25]. Recently it was shown that lpiA is necessary for the lipid lysyl-phosphatidylglycerol formation in R. tropici in low pH minimal media and confers an increased resistance to the cationic peptide polymyxin B [29]. This points to a modification of the exterior cell wall by a change of the lipid-structure. In addition smc0612 located downstream of lpiA was also found to be highly expressed, but since its expression level was slightly lower it was included in cluster B. The two open reading frames smc00612 and smc00613 are obviously products of a frameshift mutation of the orthologous gene acvB [28] and are therefore probably not functional. It could be shown that a complementation of S. meliloti 1021 with the lpiA and acvB genes of A. tumefaciens resulted in an enhanced tolerance to acidic pH (C. Sohlenkamp, personal communication).

It has been proposed that a modulated or enhanced lipid biosynthesis, as indicated by the high induction of lpiA, can increase the biosynthetic need for bicarbonate [30]. A raised demand for bicarbonate can be associated with the strongly up-regulated expression of cah, also found in cluster A. The gene cah is coding for a carbonic anhydrase that catalyses the fixation of bicarbonate. Since this gene was also highly up-regulated in response to phosphate starvation of S. meliloti it seems not to be specific for low pH stress [15].

Another early induction was observed for exoV and exoH coding for proteins of the exopolysaccharide I biosynthesis (EPS I). The discussion of this and further genes involved in EPS I biosynthesis will be addressed in a later section.

The large cluster B contains some exo genes responsible for the biosynthesis of succinoglycan and several rpoE2dependently regulated genes

The expression level of the genes comprising cluster B increased to a medium level in the first 10–20 minutes after the pH shift, and remained at this level until the end of the time course experiment (Fig. 2B).

Cluster B represents the biggest cluster and includes 74 genes. This cluster mainly consists of genes coding for hypothetical or conserved hypothetical proteins (41 genes) predominantly located on pSymA or on the chromosome. For these genes no further functional prediction can be given. Besides these genes, eight exo genes were found whose products together with the three exo genes grouped in cluster A and C are involved in the synthesis of exopolysaccharide I, also termed succinoglycan. In addition, the gene chvI coding for a regulator is part of this cluster. The genes of the EPS I biosynthesis are discussed in more detail in a following section.

The gene katC present in cluster B was annotated as a catalase. The induction of a catalase in response to low pH seems reasonable to decompose hydrogen peroxide, since a lowered pH favours the generation of radicals by the Fenton reaction. The gene katC is known to be regulated in a heat dependent mechanism by rpoE2 in S. meliloti 1021 [31]. Altogether 15 out of 41 described genes being rpoE2 dependent regulated under heat stress [31] were found exclusively in cluster B. This is not only indicating a possible role of RpoE2 in the pH stress response but also a specific expression profile of the target genes. Besides katC, ndiA, glgA2 and glgX2 the remaining 11 genes are coding for hypothetical proteins. The rpoE2 gene itself was filtered for clustering with maximum log2 fold expression values of 1.36 and 1.07 at time points 18 minutes and 33 minutes, respectively.

Cluster C contains among others genes coding for a chaperone and a component of a low O2affinity oxidase

Cluster C contains 31 genes whose expression continuously increased during the time course experiment (Fig. 2C). With over 50% (16 of 31 genes) this cluster resembles cluster B composed of a large amount of genes coding for hypothetical proteins.

In this cluster groEL5 could be found, which was the only differentially expressed gene coding for a chaperone. This gene has recently been shown to be specialised for the S. meliloti stress response [32]. Besides the DegP1 protease encoding gene, this is the only quality control system found to be up-regulated after the pH shift. In contrast to degP1 the groEL5 gene was not immediately up-regulated after the pH shift, but slowly increasing in its expression level during the time course. With nex18 a gene with unknown function could be detected, which was already shown to be higher expressed during symbiosis and in response to nutrient deprivation stress [33, 34].

The gene cyoB of the cyoABC operon was also included in cluster C. The operon codes for a cytochrome o ubiquinol oxidase, a low O2 affinity oxidase with a high proton pumping activity. It is noteworthy that qxtA, a gene coding for part of the subunit of a high O2 affinity oxidase displayed an expression profile similar to genes of cluster C, but was filtered out for clustering analysis due to missing values for three time points. It is known that an increased ΔpH affects the expression of genes of the oxidative phosphorylation. In S. medicae the transcriptional induction of fixN, a symbiosis related high O2-affinity oxidase with a low proton pumping activity was observed after overnight growth at low pH [19]. For Brucella abortus it was demonstrated that an interruption in the orthologue of the qxtAB operon, named cydAB, caused high acid sensitivity [35]. In E. coli the gene expression of the orthologues of the low O2 affinity oxidase encoded by cyoABC and the qxtAB encoded high O2 affinity oxidase was dependent of the pH [36] with a preferred expression of the high O2 affinity oxidase at low pH. Since both, the cyoABC and the qxtAB systems of S. meliloti have so far not been further investigated, their specific role in the pH response cannot be defined.

Cluster D comprises carbon uptake and fatty acid degradation genes

The mean expression of the genes grouped into cluster D shows a transient induction during the time course (Fig. 2D). In cluster D the dctA gene coding for the DctA dicarboxylate import system was found. The DctA dicarboxylate import system [37] is well characterised and a broad substrate range has been identified [38]. This dicarboxylate import system is known to be essential for symbiosis since it is supposed to provide the cells in the bacteroid state with tricarbonic acid (TCA) cycle intermediates from the host plant, e.g. succinate, malate, and fumarate.

A group of genes in this cluster points to an induced fatty acid degradation. The gene smc00976 is coding for a putative enoyl CoA hydratase and smc00977 and smc02229 are coding for putative acyl CoA dehydrogenase proteins. With glpD, a gene coding for a glycerol-3-phosphate dehydrogenase involved in the glycerol degradation could also be found in cluster D. The transient induction of genes involved in fatty acid degradation might be related to a lack of energy or the modification of the membrane lipid composition.

Cluster E contains genes involved in nitrogen metabolism, ion transport and methionine metabolism

Cluster E consists of 22 genes whose expression was lowered in response to the pH shift. The expression was lowered up to 10 minutes after pH shift and then stayed constant until the end of the time course experiment (Fig. 2E).

Cluster E contains genes involved in nitrogen metabolism. The gene glnK codes for a PII nitrogen regulatory protein activated under nitrogen limiting conditions and forms together with amtB, which encodes a high affinity ammonium transport system, an operon. The GlnK protein could also be identified as lower expressed after a short exposure of S. medicae cells to low pH [27]. It was argued by Reeve et al. that this observation might be related to some crosstalk between nitrogen and pH sensing systems during the early pH adaptation [27]. With metF, metK, bmt, and ahcY four genes involved in the methionine metabolism were also grouped in this cluster, while two other met genes were grouped into cluster F (metA) and cluster G (metH), respectively. The distribution of these genes to two other clusters of down-regulated genes might be due to the fact the met genes are not organised in an operon, but dispersed over the chromosome. S-adenosylmethionine is formed from methionine by MetK and is the major methylation compound of the cell that is needed e.g. for polyamine- or phosphatidylcholine biosynthesis. The connection between the down-regulation of the methionine metabolism and the pH response is not clear. It was shown that various abiotic stresses result in a rapid change of cellular polyamine levels [39–41].

Several genes belonging to ion uptake systems were located in cluster E, like the complete sitABCD operon and phoC and phoD of the phoCDET operon. The sitABCD operon codes for a manganese/iron transport system [42, 43]. Since an acidic pH causes a higher solubility of ions, the repression of metal ion import systems is likely to be a response to their higher availability. The phoCDET operon, codes for a high affinity phosphate transport system [44]. A phoB dependent control of phoCDET could be observed in S. meliloti [15] and in E. coli it could be shown that phoB is involved in the acid shock response [45].

Cluster F is almost exclusively composed of genes playing a role in chemotaxis and motility

Cluster F consists of genes whose expression was continuously lowered during the time course experiment (Fig. 2F). It mainly comprises genes (10 of 22) belonging to chemotaxis and flagellar biosynthesis (flgB, flgG, flgL, flgF, flgC, flgE, fliE, flbT, motA, mcpU). This phenomenon will be discussed in more detail later.

Another gene in cluster F was lppB coding for a lipoprotein, which is a major outer membrane component was grouped in cluster F. A similar expression profile as those of the flagellar biosynthesis and chemotaxis genes confirms a possible co-regulation as was observed in Salmonella enterica [46].

Cluster G consists of several genes involved in nitrogen uptake and utilization

Cluster G consists of genes whose expression was transiently lowered between 8 and 18 minutes following the pH shift and afterwards returned nearly to the ground state (Fig. 2G).

The genes nirB, nirD and narB were distributed to this cluster and are coding for nitrite and nitrate reductases forming ammonia from nitrate. A homologue of narB was found to be regulated by the low pH and microaerobiosis regulator ActR in S. medicae. A gene coding for an element of a nitrate import system (nrtB) could also be found in this cluster, while the remaining two elements of this system encoded by nrtA and nrtC were not included in the analysis because their expression values were below the threshold for filtering.

Additionally, genes coding for an ABC transport system (smb21707, smb20602, smb20603, smb20604 and smb20605) sharing homologies with amino acid and urea/short chain amide transport systems are present in cluster G. In addition to this transport system genes (smb20141, smb20142) of an ABC transport system homologous to the Dpp system from E. coli were also grouped in this cluster. This system is known for the import of dipeptides to provide the cell with essential amino acids, nitrogen and energy [47].

Cluster H is formed by genes with distinct biological functions and a high variation in their expression levels

Cluster H consists of 13 genes that were transiently lower expressed on the very beginning of the time course (Fig. 2H). The proposed encoded functions of the genes in this cluster were found to be very diverse. A secreted peroxidase gene (sma1944) [48], a flagellar biosynthesis gene (fliP), a chemotaxis sensory gene (mcpW), a nodulation gene (nodP1) and several hypothetical protein encoding genes were identified to be in cluster H. Interestingly, nodP was described as up-regulated at acidic pH in Rhizobium tropici [49], but in the study of Moron et al. cells were grown for at least one day in low pH media.

The time resolved expression profile of the S. meliloti1021 exo genes and flagellar genes following a shift to acidic pH

Overall the number of differentially expressed genes belonging to the group of EPS I biosynthesis genes and to the group of genes involved in flagellar biosynthesis and motility is striking. Most exo genes were joined together in cluster B whereas most flagellar genes were grouped together in cluster F. Furthermore, it is noticeable that the expression of the two groups of genes displayed oppositional characteristics. The EPS I biosynthesis genes responded with a fast then constant induction for the duration of the time course, whereas the flagellar genes were increasingly down-regulated. For A. tumefaciens a similar response in succession to pH stress could be identified [50]. In case of A. tumefaciens the transcriptome profiling was performed after 7 hours of growth in low pH. Also in our experiment the expressional characteristics of the exo and flagellar genes indicated that their response to acidic pH conditions lasts longer than the monitored period of one hour.

The regulator coding gene chvI was with most of the exo genes distributed to cluster B. Like in A. tumefaciens the gene chvI was up-regulated together with several genes responsible for the succinoglycan biosynthesis [50], although it is believed that chvI is a negative regulator of the exo genes [51]. A closer view on the individual expression levels of the genes of the EPS I biosynthesis gene cluster on pSymB during the time course (Fig. 4) reveals the high induction levels for the majority of the exo genes. The maximum induction in the observation period was always reached at 63 minutes after pH shift. Besides the eight exo genes found in cluster B, three exo genes grouped in cluster A and C. The exo genes in cluster A (exoV and exoH) were among the strongest up-regulated genes in this experiment. The products of these genes are responsible for the final steps of the EPS I biosynthesis. They are involved in the succinylation and pyruvilation of EPS I. It could already be shown for S. meliloti that a mutant strain of exoH is sensitive to low pH [52], indicating a particular impact of exoH on the pH tolerance and of the EPS I biosynthesis genes on the pH tolerance in general. The higher expression value of exoH compared to other exo genes might also be caused by its position as the first gene in a large operon (exoHKLAMONP) [53]. The central genes of this operon (exoA and exoM) did not show a significant change in their expression level during the time course in contrast to the bordering genes. This might be caused by mRNA instability and degradation effects. Three predicted open reading frames of unknown function within the exo gene region (smb21673, smb20952 and smb20953) did not show a differential expression during the time course indicating that they are not protein coding or that they are not connected to the EPS I biosynthesis. Overall the observed induction of exo genes is in agreement with the mucoid phenotype observed for S. meliloti after growing on low pH plates (data not shown). In low pH soils this response could be a strategy of the cell to establish a more favourable microenvironment by secreting succinoglycan. It was shown that an EPS I overproduction results in a reduced nodulation efficiency [54], therefore the induction of EPS I biosynthesis genes could also be one of the reasons for the observed limited nodulation efficiency of rhizobia in low pH soils [2].

Figure 4
figure 4

Map of genes in the EPS I biosynthesis region on pSymB and their expression in response to acidic pH. The EPS I biosynthesis gene region on pSymB is schematically displayed with its genes given by open arrows coloured according to the K-means cluster distribution. Gene names are given below. Black arrows indicate known operon structures in this region. The graph above shows on the Y-axis the time after pH-shift and on the Z-axis for each time point the expression of the corresponding genes by the M value.

Whereas the exo gene expression was increased, several genes of chemotaxis and flagellar biosynthesis (flgB, flgG, flgL, flgF, flgC, flgE, fliE, flbT, motA, mcpU) were decreased in their expression levels. After 63 minutes of low pH treatment the genes have reached the highest level of repression. VisR is the main activator of the flagellar genes and forms together with VisN the top layer of a hierarchy of three expression classes. Since the visN gene expression was decreased early in the time course experiment (therefore visN was grouped into cluster E) the other flagellar genes follow the repression of their activator [55]. The gene coding for the subordinated regulator Rem [56] was also decreasingly expressed with time, but did not reach the threshold for clustering. A detailed consideration of the expression levels of the flagellar biosynthesis genes on the chromosome (Fig. 5) reveals a repression of the complete region, with some parts responding stronger than others. The decreased expression level of motA, flgF and flgE is likely to be a result of their first position in an operon. It is noticeable that among the 10 down-regulated and strongly responding flagellar genes in cluster F five are coding for parts of the rod (flgF, flgB, flgC, fliE and flgG) and two for parts of the hook (flgE and flgL) of the flagellum. The genes motA, fliM, fliN and fliG are proposed to form an operon [55]. While the expression of motA, which is coding for a transmembrane proton channel protein, was decreased in the time course experiment, the other three genes which encode flagellar switch proteins did not respond to the shift to acidic pH. If this behaviour is caused by a specific regulation or is due to mRNA degradation processes cannot be answered. It is known that the assembly of the flagellar apparatus follows a strict sequence of steps beginning with the MS- and C- rings followed by the export apparatus, rod, hook and ending with the flagellar filament and energizing complexes [57]. The results presented here indicate that the disassembly is also performed in a defined order. The loss of flagellar motility at low pH could already be shown for the closely related Rhizobium leguminosarum bv.viciae and A. tumefaciens [50, 58], whereas the more distantly related enterobacteria E. coli and Salmonella enterica serovar Thyphimurium showed an opposite response [59–61]. For cases of induced motility it was argued that at low pH the large ΔpH drives flagellar rotation [62]. Since there are also reports of E coli where it could be demonstrated that motility is lost at low pH [63] the picture is ambiguous. A turndown of the flagellar motility genes of S. meliloti was also observed for other stresses like osmotic stress [14, 64], heat shock and nutrient starvation [31]. It is therefore apparent that this response is a general stress response of S. meliloti 1021 and not an answer specific for pH stress. Since cell motility is very energy consumptive, the repression of the motility genes is likely to save energy which is needed to face the low pH e.g. by enhancing the EPS I biosynthesis.

Figure 5
figure 5

Map of genes of the flagellar biosynthesis region on the chromosome of S. meliloti 1021 and their expression in response to acidic pH. A part of the flagellar gene region is schematically displayed with its genes given by open arrows coloured according to the K-means cluster distribution. Gene names are given below. Black arrows indicate known operon structures. The graph above shows on the Y-axis the time after pH-shift and on the Z-axis for each time point the expression of the corresponding genes by the M-value. For clarity a region of 13 consecutive genes of the flagellar operon (flgA – fliK) has been omitted. The location of the omitted region is indicated by the orthogonal lines. The ending of a flagellar operon within the omitted region is depicted by a dotted black arrow.

Conclusion

This study demonstrates the complexity of the cellular response of S. meliloti to adapt to a new environmental conditions. The mechanism of the cell to face the low pH is a mixture of several distinct reactions which follow a particular order in time. By applying K-means clustering analysis the diversity of different responses of individual genes was reduced to 8 main expression profiles. By this method a reasonable distinction between differently behaving up-regulated and down-regulated genes could be performed. Furthermore, within the obtained clusters, groups of genes with functional relationship were often joined together. Additionally, this analysis revealed that within the first 20 minutes after the shift to acidic pH the cell appears to perform the main changes necessary to adapt to the new environmental circumstances on the transcriptional level.

The immediate response of S. meliloti to low pH in the main appears to be the modification of the cell surface, since several among the highest up-regulated genes (cluster A) encode for probably secreted proteins. The differential expression of some genes was obviously only of temporary need for the cell until about 20 minutes after pH shift (as indicated by clusters D and G). Possibly an increasing demand for energy causes the activation of the dicarboxylate transport system gene dctA and of several genes of the fatty acid degradation (cluster D) while at the same time genes for nitrogen uptake and utilization (cluster G) and amino acid biosynthesis were lower expressed. The latter was clearly indicated by the lowered expression of several methionine metabolism genes.

Several genes contributing to the EPS I biosynthesis were up-regulated in response to the acidic pH shift. The secretion of EPS I might be an attempt of the cell to ameliorate the environment. In parallel a decreasing expression of motility genes can be regarded as an attempt of the cell to save energy.

The transcriptional response of S. meliloti 1021 towards low pH showed several parallels to the response in A. tumefaciens [50], with the induction of the exo genes and the repression of motility genes. Mechanisms to actively compete against a lowered pH like e.g. in E. coli by decarboxylation of amino acids (for review see [65])[66] could not be identified. Possibly in oligotrophic soils a metabolisation of amino acids is inappropriate.

Overall this work showed that the short term response to acidic pH stress does not result in a simple induction or repression of genes, but in a sequence of responses varying in their intensity over time. This indicates that a comprehensive analysis of the transcriptional response of a cell confronted with a new environmental situation requires a monitoring over a longer period of time and not only the analysis of a snap shot. Obviously, the response to acidic pH is not based on a few specific genes, but involves several genes associated with various cellular functions. On the other hand, a considerable part of the responding genes belongs to the group of hypothetical genes. These genes represent promising objectives for future investigations.

Methods

Media and growth conditions

S. meliloti strain 1021 was cultivated in Erlenmeyer flasks at 30°C in Vincent minimal medium (VMM) [67] and shaken at 140 rpm. With exception of 37 μM iron(III) choride no additional metals have been added to the VMM. The pH of the VMM was adjusted by using either HCl or NaOH. Precultures were grown in tryptone yeast complex medium [68] with appropriate antibiotics (600 μg/ml streptomycin).

For pH shift experiments cells of three independent cultures were grown in 100 ml buffered VMM (20 mM BisTris) to an o.D.580 of 0.8. All of the following steps were carried out under gentle conditions using pre-warmed equipment. Cell cultures of each flask were halved in two even parts of 50 ml, centrifuged (10000 × g, 2 min, 30°C) and the supernatant was discarded. The cell pellets were resuspended in 50 ml VMM with pH 5.75 and 50 ml VMM with pH 7.0, respectively, and incubated at 30°C. At six time points cell suspension probes of 5 ml were harvested from each flask. Immediately centrifuged (10000 × g, 1 min, 4°C) the resulting pellets were instantly frozen by liquid nitrogen for later RNA preparation. Cell suspension probes were harvested at 3, 8, 13, 18, 33, and 63 minutes following the pH shift.

RNA isolation

RNA was isolated according to the protocol published by Rüberg et al. [14]. Total RNA was prepared using the RNeasy mini kit (QIAGEN, Hildesheim, Germany). By ribolysation (30 s; speed, 6.5; Hybaid, Heidelberg, Germany) cells were disrupted in the RLT buffer provided with the kit in Fast Protein Tubes (Qbiogene, Carlsbad, CA).

Transcriptional profiling using the SM6kOligo whole genome microarray

The well established Sm6kOligo microarray described by Krol and Becker [15] was employed for transcriptional profiling. For each preparation of Cy3 and Cy5 labelled cDNAs 10 μg of total RNA were used [69]. To each microarray the cDNA of the pH 7.0 and pH 5.75 grown cultures were mixed and hybridised. The microarray experiments were performed in three biological replicates.

The acquisition of the microarray images was performed as described previously [14, 15]. By using the ImaGene 5.0 software (Biodiscovery Inc., Los Angeles, CA, USA) the mean signal and mean local background intensities for each spot were identified and calculated. If R was ≤ 1.5 in both channels, spots were flagged as "empty", the remaining spots were used for further analysis. The log2 value of the intensity ratios (Mi) was calculated for each spot with Mi = log2(Ri/Gi). Ri = Ich1(i)-Bgch1(i) and Gi = Ich2(i)-Bgch2(i) with Ich1i and Ich2i being the intensities of a spot in channel l or channel 2 and Bgch1(i) and Bgch2(i) being the background intensity of a spot in channel 1 or channel 2, respectively. The mean intensity was calculated for each spot with Ai = log2(RiGi)0.5. Normalization and t-statistics were carried out using the Emma 1.1 microarray data analysis software [26]. It should be mentioned that in this work genes with a positive M value are addressed as "up-regulated" and genes with a negative M value are addressed as "down-regulated", although a positive value will also be calculated if a gene is less strong down-regulated under pH 5.75 than under pH 7.0 and vice versa.

The microarray results were verified for specific genes (lpiA and phoC) by quantitative reverse transcription-PCR using a QuantiTect SYBR Green reverse transcription-PCR kit (QIAGEN, Hildesheim, Germany) according to the manufacturer's instructions.

Filtering and clustering analysis of the microarray data

For clustering purposes only those genes were taken into account which had an evaluable expression value for at least 5 of 6 time points and for which at least one time point had an M value of ≥ 2 or ≤ 2. Expression values are evaluable if the value for R is ≥ 1.5 for both channels. 230 genes fulfilled these criteria. For these 230 genes 444 time points showed an M value of ≥ 2 or ≤ -2. In testing these time points for an FDR (False Discovery Rate) corrected P value of ≥ 0.05, only 4 results (≈ 0.9%) were above this value. These were: t3 smc01523 P = 0.07, t33 smc04173 P = 0.09, t63 smb21026 P = 0.06, and t63 sma1736 P = 0.22. For K means clustering analysis of the microarray experiment data the Genesis software was used (Sturn, 2001; http://genome.tugraz.at/genesisclient/genesisclient_description.shtml). The K means clustering was carried out in 8 groups.

References

  1. Zahran HH: Rhizobium-legume symbiosis and nitrogen fixation under severe conditions and in an arid climate. Microbiol Mol Biol Rev. 1999, 63: 968-89.

    CAS  PubMed Central  PubMed  Google Scholar 

  2. Ibekwe AM, Angle JS, Chaney RL, vanBerkum P: Enumeration and N2 fixation potential of Rhizobium leguminosarum biovar trifolii grown in soil with varying pH values and heavy metal concentrations. Agriculture Ecosystems & Environment. 1997, 61: 103-111.

    Article  CAS  Google Scholar 

  3. Graham PH, Viteri SE, Mackie F, Vargas AT, Palacios A: Variation in acid soil tolerance among strains of Rhizobium phaseoli. Field Crops Research. 1982, 5: 121-128.

    Article  Google Scholar 

  4. Brockwell J, Pilka A, Holliday RA: Soil-pH is a major determinant of the numbers of naturally-occurring Rhizobium meliloti in noncultivated soils in central New South Wales. Australian Journal of Experimental Agriculture. 1991, 31: 211-219.

    Article  Google Scholar 

  5. Marschner H: Mineral nutrition of higher plants. 2006, Academic Press, London

    Google Scholar 

  6. Mellor RB: Bacteroids in the Rhizobium-legume symbiosis inhabit a plant internal lytic compartment – implications for other microbial endosymbioses. Journal of Experimental Botany. 1989, 40: 831-839.

    Article  CAS  Google Scholar 

  7. Priefer UB, Aurag J, Boesten B, Bouhmouch I, Defez R, Filali-Maltouf A, et al: Characterisation of Phaseolus symbionts isolated from Mediterranean soils and analysis of genetic factors related to pH tolerance. Journal of Biotechnology. 2001, 91: 223-236.

    Article  CAS  PubMed  Google Scholar 

  8. Galibert F, Finan TM, Long SR, Pühler A, Abola P, Ampe F, et al: The composite genome of the legume symbiont Sinorhizobium meliloti. Science. 2001, 293: 668-72.

    Article  CAS  PubMed  Google Scholar 

  9. Capela D, Barloy-Hubler F, Gouzy J, Bothe G, Ampe F, Batut J, et al: Analysis of the chromosome sequence of the legume symbiont Sinorhizobium meliloti strain 1021. Proc Natl Acad Sci USA. 2001, 98: 9877-82.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  10. Barnett MJ, Fisher RF, Jones T, Komp C, Abola AP, Barloy-Hubler F, et al: Nucleotide sequence and predicted functions of the entire Sinorhizobium meliloti pSymA megaplasmid. Proc Natl Acad Sci USA. 2001, 98: 9883-9888.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  11. Finan TM, Weidner S, Wong K, Buhrmester J, Chain P, Vorhölter FJ, et al: The complete sequence of the 1,683-kb pSymB megaplasmid from the N2-fixing endosymbiont Sinorhizobium meliloti. Proc Natl Acad Sci USA. 2001, 98: 9889-9894.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  12. Becker A, Berges H, Krol E, Bruand C, Rüberg S, Capela D, et al: Global changes in gene expression in Sinorhizobium meliloti 1021 under microoxic and symbiotic conditions. Mol Plant Microbe Interact. 2004, 17: 292-303.

    Article  CAS  PubMed  Google Scholar 

  13. Djordjevic MA, Chen HC, Natera S, Van Noorden G, Menzel C, Taylor S, et al: A global analysis of protein expression profiles in Sinorhizobium meliloti: discovery of new genes for nodule occupancy and stress adaptation. Mol Plant Microbe Interact. 2003, 16: 508-24.

    Article  CAS  PubMed  Google Scholar 

  14. Rüberg S, Tian ZX, Krol E, Linke B, Meyer F, Wang Y, et al: Construction and validation of a Sinorhizobium meliloti whole genome DNA microarray: genome-wide profiling of osmoadaptive gene expression. J Biotechnol. 2003, 106: 255-68.

    Article  PubMed  Google Scholar 

  15. Krol E, Becker A: Global transcriptional analysis of the phosphate starvation response in Sinorhizobium meliloti strains 1021 and 2011. Mol Genet Genomics. 2004, 272: 1-17.

    Article  CAS  PubMed  Google Scholar 

  16. Foster JW, Hall HK: Adaptive acidification tolerance response of Salmonella typhimurium. Journal of Bacteriology. 1990, 172: 771-778.

    CAS  PubMed Central  PubMed  Google Scholar 

  17. Bearson S, Bearson B, Foster JW: Acid stress responses in enterobacteria. FEMS Microbiol Lett. 1997, 147: 173-180.

    Article  CAS  PubMed  Google Scholar 

  18. Reeve WG, Tiwari RP, Worsley PS, Dilworth MJ, Glenn AR, Howieson JG: Constructs for insertional mutagenesis, transcriptional signal localization and gene regulation studies in root nodule and other bacteria. Microbiology. 1999, 145 (Pt 6): 1307-16.

    Article  CAS  PubMed  Google Scholar 

  19. Tiwari RP, Reeve WG, Fenner BJ, Dilworth MJ, Glenn AR, Howieson JG: Probing for pH-regulated genes in Sinorhizobium medicae using transcriptional analysis. J Mol Microbiol Biotechnol. 2004, 7: 133-9.

    Article  CAS  PubMed  Google Scholar 

  20. Reeve WG, Tiwari RP, Kale NB, Dilworth MJ, Glenn AR: ActP controls copper homeostasis in Rhizobium leguminosarum bv. viciae and Sinorhizobium meliloti preventing low pH-induced copper toxicity. Mol Microbiol. 2002, 43: 981-91.

    Article  CAS  PubMed  Google Scholar 

  21. Tiwari RP, Reeve WG, Dilworth MJ, Glenn AR: An essential role for actA in acid tolerance of Rhizobium meliloti. Microbiology. 1996, 142 (Pt 3): 601-10.

    Article  CAS  PubMed  Google Scholar 

  22. Tiwari RP, Reeve WG, Dilworth MJ, Glenn AR: Acid tolerance in Rhizobium meliloti strain WSM419 involves a two-component sensor-regulator system. Microbiology. 1996, 142 (Pt 7): 1693-704.

    Article  CAS  PubMed  Google Scholar 

  23. Fenner BJ, Tiwari RP, Reeve WG, Dilworth MJ, Glenn AR: Sinorhizobium medicae genes whose regulation involves the ActS and/or ActR signal transduction proteins. FEMS Microbiol Lett. 2004, 236: 21-31.

    Article  CAS  PubMed  Google Scholar 

  24. Dilworth MJ, Howieson JG, Reeve WG, Tiwari RP, Glenn AR: Acid tolerance in legume root nodule bacteria and selecting for it. Australian Journal of Experimental Agriculture. 2001, 41: 435-446.

    Article  CAS  Google Scholar 

  25. Vinuesa P, Neumann-Silkow F, Pacios-Bras C, Spaink HP, Martinez-Romero E, Werner D: Genetic analysis of a pH-regulated operon from Rhizobium tropici CIAT899 involved in acid tolerance and nodulation competitiveness. Mol Plant Microbe Interact. 2003, 16: 159-168.

    Article  CAS  PubMed  Google Scholar 

  26. Dondrup M, Goesmann A, Bartels D, Kalinowski J, Krause L, Linke B, et al: EMMA: a platform for consistent storage and efficient analysis of microarray data. J Biotechnol. 2003, 106: 135-46.

    Article  CAS  PubMed  Google Scholar 

  27. Reeve WG, Tiwari RP, Guerreiro N, Stubbs J, Dilworth MJ, Glenn AR, et al: Probing for pH-regulated proteins in Sinorhizobium medicae using proteomic analysis. J Mol Microbiol Biotechnol. 2004, 7: 140-7.

    Article  CAS  PubMed  Google Scholar 

  28. Reeve WG, Bräu L, Castelli J, Garau G, Sohlenkamp C, Geiger O, et al: The Sinorhizobium medicae WSM419 lpiA gene is transcriptionally activated by FsrR and required to enhance survival in lethal acid conditions. Microbiology. 2006, 152: 3049-3059.

    Article  CAS  PubMed  Google Scholar 

  29. Sohlenkamp C, Galindo-Lagunas KA, Guan ZQ, Vinuesa P, Robinson S, Thomas-Oates J, et al: The lipid lysyl-phosphatidylglycerol is present in membranes of Rhizobium tropici CIAT899 and confers increased resistance to polymyxin B under acidic growth conditions. Molecular Plant-Microbe Interactions. 2007, 20: 1421-1430.

    Article  CAS  PubMed  Google Scholar 

  30. Merlin C, Masters M, McAteer S, Coulson A: Why is carbonic anhydrase essential to Escherichia coli ?. Journal of Bacteriology. 2003, 185: 6415-6424.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  31. Sauviac L, Philippe H, Phok K, Bruand C: An extracytoplasmic function sigma factor acts as a general stress response regulator in Sinorhizobium meliloti. Journal of Bacteriology. 2007, 189: 4204-4216.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  32. Bittner AN, Foltz A, Oke V: Only one of five groEL genes is required for viability and successful symbiosis in Sinorhizobium meliloti. J Bacteriol. 2007, 189: 1884-1889.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  33. Oke V, Long SR: Bacterial genes induced within the nodule during the Rhizobium-legume symbiosis. Mol Microbiol. 1999, 32: 837-849.

    Article  CAS  PubMed  Google Scholar 

  34. Davey ME, de Bruijn FJ: A homologue of the tryptophan-rich sensory protein TspO and FixL regulate a novel nutrient deprivation-induced Sinorhizobium meliloti locus. Appl Environ Microbiol. 2000, 66: 5353-5359.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  35. Endley S, McMurray D, Ficht TA: Interruption of the cydB locus in Brucella abortus attenuates intracellular survival and virulence in the mouse model of infection. J Bacteriol. 2001, 183: 2454-62.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  36. Cotter PA, Chepuri V, Gennis RB, Gunsalus RP: Cytochrome-O (cyoABCDE) and D (cydAB) oxidase gene-expression in Escherichia coli is regulated by oxygen, pH, and the fnr gene-product. Journal of Bacteriology. 1990, 172: 6333-6338.

    CAS  PubMed Central  PubMed  Google Scholar 

  37. Watson RJ, Chan YK, Wheatcroft R, Yang AF, Han SH: Rhizobium meliloti genes required for C-4-dicarboxylate transport and symbiotic nitrogen-fixation are located on a megaplasmid. Journal of Bacteriology. 1988, 170: 927-934.

    CAS  PubMed Central  PubMed  Google Scholar 

  38. Yurgel S, Mortimer MW, Rogers KN, Kahn ML: New substrates for the dicarboxylate transport system of Sinorhizobium meliloti. Journal of Bacteriology. 2000, 182: 4216-4221.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  39. Dubin DT, Rosenthal SM: The acetylation of polyamines in Escherichia coli. J Biol Chem. 1960, 235: 776-782.

    CAS  PubMed  Google Scholar 

  40. Munro GF, Hercules K, Morgan J, Sauerbier W: Dependence of the putrescine content of Escherichia coli on the osmotic strength of the medium. J Biol Chem. 1972, 247: 1272-1280.

    CAS  PubMed  Google Scholar 

  41. Yamamoto S, Yamasaki K, Takashina K, Katsu T, Shinoda S: Characterization of putrescine production in nongrowing Vibrio parahaemolyticus cells in response to external osmolality. Microbiol Immunol. 1989, 33: 11-21.

    Article  CAS  PubMed  Google Scholar 

  42. Chao TC, Becker A, Buhrmester J, Pühler A, Weidner S: The Sinorhizobium meliloti fur gene regulates, with dependence on Mn(II), transcription of the sitABCD operon, encoding a metal-type transporter. Journal of Bacteriology. 2004, 186: 3609-3620.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  43. Platero RA, Jaureguy M, Battistoni FJ, Fabiano ER: Mutations in sitB and sitD genes affect manganese-growth requirements in Sinorhizobium meliloti. Fems Microbiology Letters. 2003, 218: 65-70.

    Article  CAS  PubMed  Google Scholar 

  44. Bardin S, Dan S, Osteras M, Finan TM: A phosphate transport system is required for symbiotic nitrogen fixation by Rhizobium meliloti. Journal of Bacteriology. 1996, 178: 4540-4547.

    CAS  PubMed Central  PubMed  Google Scholar 

  45. Suziedeliene E, Suziedelis K, Garbenciute V, Normark S: The acid-inducible asr gene in Escherichia coli: Transcriptional control by the phoBR operon. Journal of Bacteriology. 1999, 181: 2084-2093.

    CAS  PubMed Central  PubMed  Google Scholar 

  46. Iyoda S, Kamidoi T, Hirose K, Kutsukake K, Watanabe H: A flagellar gene fliZ regulates the expression of invasion genes and virulence phenotype in Salmonella enterica serovar Typhimurium. Microbial Pathogenesis. 2001, 30: 81-90.

    Article  CAS  PubMed  Google Scholar 

  47. Olson ER, Dunyak DS, Jurss LM, Poorman RA: Identification and characterization of dppa, an Escherichia coli gene encoding a periplasmic dipeptide transport protein. Journal of Bacteriology. 1991, 173: 234-244.

    CAS  PubMed Central  PubMed  Google Scholar 

  48. Barloy-Hubler F, Cheron A, Hellegouarch A, Galibert F: Smc01944, a secreted peroxidase induced by oxidative stresses in Sinorhizobium meliloti 1021. Microbiology. 2004, 150: 657-64.

    Article  CAS  PubMed  Google Scholar 

  49. Moron B, Soria-Diaz ME, Ault J, Verroios G, Noreen S, Rodriguez-Navarro DN, et al: Low pH changes the profile of nodulation factors produced by Rhizobium tropici CIAT899. Chem Biol. 2005, 12: 1029-1040.

    Article  CAS  PubMed  Google Scholar 

  50. Yuan ZC, Liu P, Saenkham P, Kerr K, Nester EW: Transcriptome profiling and functional analysis of Agrobacterium tumefaciens reveals a general conserved response to acidic conditions (pH 5.5) and a complex acid-mediated signaling involved in Agrobacterium-plant interactions. J Bacteriol. 2008, 190: 494-507.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  51. Yao SY, Luo L, Har KJ, Becker A, Rüberg S, Yu GQ, et al: Sinorhizobium meliloti ExoR and ExoS proteins regulate both succinoglycan and flagellum production. Journal of Bacteriology. 2004, 186: 6042-6049.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  52. Glenn AR, Reeve WG, Tiwari RP, Dilworth MJ: Acid tolerance in root nodule bacteria. Novartis Found Symp. 1999, 221: 112-26.

    CAS  PubMed  Google Scholar 

  53. Becker A, Rüberg S, Baumgarth B, Bertram-Drogatz PA, Quester I, Pühler A: Regulation of succinoglycan and galactoglucan biosynthesis in Sinorhizobium meliloti. J Mol Microbiol Biotechnol. 2002, 4: 187-190.

    CAS  PubMed  Google Scholar 

  54. Cheng HP, Walker GC: Succinoglycan is required for initiation and elongation of infection threads during nodulation of alfalfa by Rhizobium meliloti. J Bacteriol. 1998, 180: 5183-5191.

    CAS  PubMed Central  PubMed  Google Scholar 

  55. Sourjik V, Muschler P, Scharf B, Schmitt R: VisN and VisR are global regulators of chemotaxis, flagellar, and motility genes in Sinorhizobium (Rhizobium) meliloti. Journal of Bacteriology. 2000, 182: 782-788.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  56. Rotter C, Mühlbacher S, Salamon D, Schmitt R, Scharf B: Rem, a new transcriptional activator of motility and chemotaxis in Sinorhizobium meliloti. J Bacteriol. 2006, 188: 6932-6942.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  57. Eggenhofer E, Rachel R, Haslbeck M, Scharf B: MotD of Sinorhizobium meliloti and related alpha-proteobacteria is the flagellar-hook-length regulator and therefore reassigned as FliK. J Bacteriol. 2006, 188: 2144-2153.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  58. Bowra BJ, Dilworth MJ: Motility and chemotaxis towards sugars in Rhizobium leguminosarum. Journal of General Microbiology. 1981, 126: 231-235.

    CAS  Google Scholar 

  59. Maurer LM, Yohannes E, Bondurant SS, Radmacher M, Slonczewski JL: pH regulates genes for flagellar motility, catabolism, and oxidative stress in Escherichia coli K-12. Journal of Bacteriology. 2005, 187: 304-319.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  60. de Jonge R, Ritmeester WS, van Leusden FM: Adaptive responses of Salmonella enterica serovar Typhimurium DT104 and other S. Typhimurium strains and Escherichia coli O157 to low pH environments. J Appl Microbiol. 2003, 94: 625-632.

    Article  CAS  PubMed  Google Scholar 

  61. Hickey EW, Hirshfield IN: Low-pH-induced effects on patterns of protein synthesis and on internal pH in Escherichia coli and Salmonella typhimurium. Appl Environ Microbiol. 1990, 56: 1038-1045.

    CAS  PubMed Central  PubMed  Google Scholar 

  62. Khan S, Macnab RM: Proton chemical potential, proton electrical potential and bacterial motility. J Mol Biol. 1980, 138: 599-614.

    Article  CAS  PubMed  Google Scholar 

  63. Soutourina OA, Krin E, Laurent-Winter C, Hommais F, Danchin A, Bertin PN: Regulation of bacterial motility in response to low pH in Escherichia coli: the role of H-NS protein. Microbiology. 2002, 148: 1543-1551.

    Article  CAS  PubMed  Google Scholar 

  64. Dominguez-Ferreras A, Perez-Arnedo R, Becker A, Olivares J, Soto MJ, Sanjuan J: Transcriptome profiling reveals the importance of plasmid pSymB for osmoadaptation of Sinorhizobium meliloti. J Bacteriol. 2006, 188: 7617-7625.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  65. Foster JW: Escherichia coli acid resistance: tales of an amateur acidophile. Nat Rev Microbiol. 2004, 2: 898-907.

    Article  CAS  PubMed  Google Scholar 

  66. Kannan G, Wilks JC, Fitzgerald DM, Jones BD, Bondurant SS, Slonczewski JL: Rapid acid treatment of Escherichia coli: transcriptomic response and recovery. BMC Microbiol. 2008, 8: 37-

    Article  PubMed Central  PubMed  Google Scholar 

  67. Vincent JM: A manual for the practical study of root nodule bacteria. Oxford: IBP Handbook No 15. 1970

    Google Scholar 

  68. Beringer JE: R factor transfer in Rhizobium leguminosarum. J Gen Microbiol. 1974, 84: 188-98.

    CAS  PubMed  Google Scholar 

  69. Derisi JL, Iyer VR, Brown PO: Exploring the metabolic and genetic control of gene expression on a genomic scale. Science. 1997, 278: 680-686.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

This work was performed in the framework of project QLK3-CT-2002-02097 funded by the commission of the European Communities. We thank Anke Becker for the possibility to use the Sm6kOligo microarrays and the analysis environment as well as Victoria Gödde and Manuela Meyer for the excellent technical support.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Stefan Weidner.

Additional information

Authors' contributions

CH and SW designed the study, CH performed all works. SW and AP provided critical expertise for the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

12866_2008_708_MOESM1_ESM.jpeg

Additional file 1: Heat map of cluster A. By K-means the transcriptional data obtained by microarray analysis of the S. meliloti 1021 pH shock time course experiment were grouped into eight clusters. In cluster A, genes exhibiting a strong and permanent induction were accumulated. Genes in this cluster remained up-regulated for the whole observation period. Presumably, these genes have a special impact for S. meliloti in facing low pH conditions. Each column of the heat map represents one time point after shift from pH 7.0 to pH 5.75 in the following order: 3, 8, 13, 18, 33, and 63 minutes. The values in the boxes are the M-values of a specific gene represented in a row. The background colour visualises the strength of the induction/lower expression (red/green) by the colour intensity. (JPEG 109 KB)

12866_2008_708_MOESM2_ESM.jpeg

Additional file 2: Heat map of cluster B of the eight clusters calculated by K-means clustering of the transcriptional data obtained by microarray analysis of the S. meliloti 1021 pH shock time course experiment. Cluster B is the largest cluster. The genes in this cluster are permanently up-regulated in response to the pH shift. It contains exo genes responsible for the biosynthesis of succinoglycan and several genes which are rpoE2 dependently regulated. Among the genes in cluster B several encode for hypothetical proteins. Each column of the heat map represents one time point after shift from pH 7.0 to pH 5.75 in the following order: 3, 8, 13, 18, 33, and 63 minutes. The values in the boxes are the M-values of a specific gene represented in a row. The background colour visualises the strength of the induction/lower expression (red/green) by the colour intensity. (JPEG 574 KB)

12866_2008_708_MOESM3_ESM.jpeg

Additional file 3: Heat map of cluster C of the eight clusters calculated by K-means clustering of the transcriptional data obtained by microarray analysis of the S. meliloti 1021 pH shock time course experiment. Cluster C contains over 50% hypothetical genes and a gene coding for a chaperone as well as for a component of a low O2 affinity oxidase. The genes in cluster C showed a progressive permanent induction in their mean expression behaviour. Each column of the heat map represents one time point after shift from pH 7.0 to pH 5.75 in the following order: 3, 8, 13, 18, 33, and 63 minutes. The values in the boxes are the M-values of a specific gene represented in a row. The background colour visualises the strength of the induction/lower expression (red/green) by the colour intensity. (JPEG 275 KB)

12866_2008_708_MOESM4_ESM.jpeg

Additional file 4: Heat map of cluster D of the eight clusters calculated by K-means clustering of the transcriptional data obtained by microarray analysis of the S. meliloti1021 pH shock time course experiment. Cluster D comprises carbon uptake and fatty acid degradation genes. The containing genes were transiently up-regulated during the first 10 to 30 minutes following the pH shift. Each column of the heat map represents one time point after shift from pH 7.0 to pH 5.75 in the following order: 3, 8, 13, 18, 33, and 63 minutes. The values in the boxes are the M-values of a specific gene represented in a row. The background colour visualises the strength of the induction/lower expression (red/green) by the colour intensity. (JPEG 210 KB)

12866_2008_708_MOESM5_ESM.jpeg

Additional file 5: Heat map of cluster E of the eight clusters calculated by K-means clustering of the transcriptional data obtained by microarray analysis of the S. meliloti1021 pH shock time course experiment. Cluster E contains genes involved in nitrogen metabolism, ion transport and amino acid biosynthesis. These genes were decreased in their expression value up to 20 minutes after pH shift and then stayed permanently down-regulated. Each column of the heat map represents one time point after shift from pH 7.0 to pH 5.75 in the following order: 3, 8, 13, 18, 33, and 63 minutes. The values in the boxes are the M-values of a specific gene represented in a row. The background colour visualises the strength of the induction/lower expression (red/green) by the colour intensity. (JPEG 236 KB)

12866_2008_708_MOESM6_ESM.jpeg

Additional file 6: Heat map of cluster F of the eight clusters calculated by K-means clustering of the transcriptional data obtained by microarray analysis of the S. meliloti1021 pH shock time course experiment. Cluster F is almost exclusively composed of genes playing a role in chemotaxis and motility. Genes in this cluster showed a progressive permanent repression for the duration of the time course. Each column of the heat map represents one time point after shift from pH 7.0 to pH 5.75 in the following order: 3, 8, 13, 18, 33, and 63 minutes. The values in the boxes are the M-values of a specific gene represented in a row. The background colour visualises the strength of the induction/lower expression (red/green) by the colour intensity. (JPEG 223 KB)

12866_2008_708_MOESM7_ESM.jpeg

Additional file 7: Heat map of cluster G of the eight clusters calculated by K-means clustering of the transcriptional data obtained by microarray analysis of the S. meliloti1021 pH shock time course experiment. Cluster G consists of several genes involved in nitrogen uptake and utilization. Genes in this cluster were transiently down-regulated with a minimum before 20 minutes after pH shift. Each column of the heat map represents one time point after shift from pH 7.0 to pH 5.75 in the following order: 3, 8, 13, 18, 33, and 63 minutes. The values in the boxes are the M-values of a specific gene represented in a row. The background colour visualises the strength of the induction/lower expression (red/green) by the colour intensity. (JPEG 292 KB)

12866_2008_708_MOESM8_ESM.jpeg

Additional file 8: Heat map of cluster H of the eight clusters calculated by K-means clustering of the transcriptional data obtained by microarray analysis of the S. meliloti1021 pH shock time course experiment. The small cluster H is formed by genes with distinct biological functions and a high variation in their expression levels. Genes in this cluster showed an ultra short transient repression for the first time point 3 minutes after pH shift. Each column of the heat map represents one time point after shift from pH 7.0 to pH 5.75 in the following order: 3, 8, 13, 18, 33, and 63 minutes. The values in the boxes are the M-values of a specific gene represented in a row. The background colour visualises the strength of the induction/lower expression (red/green) by the colour intensity. (JPEG 129 KB)

12866_2008_708_MOESM9_ESM.xls

Additional file 9: Spreadsheet of the 230 genes used for clustering analysis. Given is the name of each gene and its corresponding annotation, as well as the M-values calculated for the time course experiment. The last column indicates the cluster, in which the gene was distributed by K-means clustering. (XLS 62 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Hellweg, C., Pühler, A. & Weidner, S. The time course of the transcriptomic response of Sinorhizobium meliloti1021 following a shift to acidic pH. BMC Microbiol 9, 37 (2009). https://0-doi-org.brum.beds.ac.uk/10.1186/1471-2180-9-37

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/1471-2180-9-37

Keywords