The time course of the transcriptomic response of Sinorhizobium meliloti1021 following a shift to acidic pH
© Hellweg et al; licensee BioMed Central Ltd. 2009
Received: 02 August 2008
Accepted: 15 February 2009
Published: 15 February 2009
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.
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.
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.
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 ). 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 . Once a symbiosis has been established the symbiosome has been postulated to form an acidic and lytic compartment . 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  (for review see ). 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  a functionally diverse set of pH responsive and acid tolerance related genes could be identified . Gene products required for acid tolerance in S. medicae are for example ActP, a CPx heavy metal transporting ATPase , and ActA, an apolipoprotein acyl transferase . A gene coding for a regulatory protein known to be required for the acid tolerance in S. medicae is actR . 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 . 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 .
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  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).
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  (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  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 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 . 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 . 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 . 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  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 . 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 .
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 . Altogether 15 out of 41 described genes being rpoE2 dependent regulated under heat stress  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 . 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 . For Brucella abortus it was demonstrated that an interruption in the orthologue of the qxtAB operon, named cydAB, caused high acid sensitivity . 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  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  is well characterised and a broad substrate range has been identified . 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 . 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 . 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 . A phoB dependent control of phoCDET could be observed in S. meliloti  and in E. coli it could be shown that phoB is involved in the acid shock response .
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 .
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 .
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) , 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 , 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 . 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.
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 , 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 ) 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.
Media and growth conditions
S. meliloti strain 1021 was cultivated in Erlenmeyer flasks at 30°C in Vincent minimal medium (VMM)  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  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 was isolated according to the protocol published by Rüberg et al. . 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  was employed for transcriptional profiling. For each preparation of Cy3 and Cy5 labelled cDNAs 10 μg of total RNA were used . 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 . 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.
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.
- Zahran HH: Rhizobium-legume symbiosis and nitrogen fixation under severe conditions and in an arid climate. Microbiol Mol Biol Rev. 1999, 63: 968-89.PubMed CentralPubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- Marschner H: Mineral nutrition of higher plants. 2006, Academic Press, LondonGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Foster JW, Hall HK: Adaptive acidification tolerance response of Salmonella typhimurium. Journal of Bacteriology. 1990, 172: 771-778.PubMed CentralPubMedGoogle Scholar
- Bearson S, Bearson B, Foster JW: Acid stress responses in enterobacteria. FEMS Microbiol Lett. 1997, 147: 173-180.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Merlin C, Masters M, McAteer S, Coulson A: Why is carbonic anhydrase essential to Escherichia coli ?. Journal of Bacteriology. 2003, 185: 6415-6424.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- Oke V, Long SR: Bacterial genes induced within the nodule during the Rhizobium-legume symbiosis. Mol Microbiol. 1999, 32: 837-849.PubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedGoogle Scholar
- 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.PubMed CentralPubMedGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- Dubin DT, Rosenthal SM: The acetylation of polyamines in Escherichia coli. J Biol Chem. 1960, 235: 776-782.PubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedGoogle Scholar
- 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.PubMed CentralPubMedGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- Glenn AR, Reeve WG, Tiwari RP, Dilworth MJ: Acid tolerance in root nodule bacteria. Novartis Found Symp. 1999, 221: 112-26.PubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMed CentralPubMedGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- Bowra BJ, Dilworth MJ: Motility and chemotaxis towards sugars in Rhizobium leguminosarum. Journal of General Microbiology. 1981, 126: 231-235.Google Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedGoogle Scholar
- Khan S, Macnab RM: Proton chemical potential, proton electrical potential and bacterial motility. J Mol Biol. 1980, 138: 599-614.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMed CentralPubMedView ArticleGoogle Scholar
- Foster JW: Escherichia coli acid resistance: tales of an amateur acidophile. Nat Rev Microbiol. 2004, 2: 898-907.PubMedView ArticleGoogle Scholar
- 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-PubMed CentralPubMedView ArticleGoogle Scholar
- Vincent JM: A manual for the practical study of root nodule bacteria. Oxford: IBP Handbook No 15. 1970Google Scholar
- Beringer JE: R factor transfer in Rhizobium leguminosarum. J Gen Microbiol. 1974, 84: 188-98.PubMedGoogle Scholar
- Derisi JL, Iyer VR, Brown PO: Exploring the metabolic and genetic control of gene expression on a genomic scale. Science. 1997, 278: 680-686.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.