A consolidated analysis of the physiologic and molecular responses induced under acid stress in the legume-symbiont model-soil bacterium Sinorhizobium meliloti.

Abiotic stresses in general and extracellular acidity in particular disturb and limit nitrogen-fixing symbioses between rhizobia and their host legumes. Except for valuable molecular-biological studies on different rhizobia, no consolidated models have been formulated to describe the central physiologic changes that occur in acid-stressed bacteria. We present here an integrated analysis entailing the main cultural, metabolic, and molecular responses of the model bacterium Sinorhizobium meliloti growing under controlled acid stress in a chemostat. A stepwise extracellular acidification of the culture medium had indicated that S. meliloti stopped growing at ca. pH 6.0-6.1. Under such stress the rhizobia increased the O2 consumption per cell by more than 5-fold. This phenotype, together with an increase in the transcripts for several membrane cytochromes, entails a higher aerobic-respiration rate in the acid-stressed rhizobia. Multivariate analysis of global metabolome data served to unequivocally correlate specific-metabolite profiles with the extracellular pH, showing that at low pH the pentose-phosphate pathway exhibited increases in several transcripts, enzymes, and metabolites. Further analyses should be focused on the time course of the observed changes, its associated intracellular signaling, and on the comparison with the changes that operate during the sub lethal acid-adaptive response (ATR) in rhizobia.

Rhizobia are Gram-negative α- and β-proteobacteria that have the ability to fix atmospheric nitrogen in symbiotic association with legumes12 as well as with nonleguminous plants in a few known examples3. That capacity gives these associations a central role in the N cycle, thus making rhizobia also a subject of interest in agricultural-production systems4567. The establishment of symbiosis in natural soils, however, is frequently limited by abiotic-stress conditions which affect that mutualistic interaction in different ways8. Thus, the acidity of soils limits nodulation and N2 fixation in many legume-rhizobia symbioses of agronomic interest9. Soil acidity is an extended edaphic condition in culturable lands over the entire globe10, and is known to be a major limiting circumstance for legume productivity. The increased concentration of hydrogen ions negatively affects at once the host plant1112, the rhizobia13, and the symbiotic relationship141516. Different stages in the interaction between rhizobia and the host roots have been reported to be affected by acidity, including the production of nodulation factors1417181920, the attachment of rhizobia to roots14, the number of root nodules15, and the nitrogenase activity18.

Sinorhizobium meliloti and Sinorhizobium medicae are rhizobia with the ability to establish symbiosis with legumes of the genera Medicago, Melilotus, and Trigonella. These rhizobia are particularly sensitive to low pH21, growing either very slowly or not at all at pHs between 6.0–5.5 depending on the medium and the cultivation condition122223. For this reason, acid tolerance has been considered a positive trait to perform well—in both survival and symbiosis—under acid conditions24. The screening for acid-tolerant alfalfa-nodulating isolates that could colonize and/or persist in acidic media and/or soils thus resulted in a selection of novel strains with enhanced saprophytic competence11 and/or improved symbiosis under moderately-acid conditions232526272829. On the basis of these and other previous results, the genetic analysis of acid tolerance in rhizobia in general, and in alfalfa symbionts in particular, emerged as a potential powerful approach for a rational improvement of the rhizobial performance at low pH. As a result of the screening of Tn5-mutant libraries, several genes for acid tolerance (namely act, among others) could be identified in S. medicae—including actA30 and actR/S3132, actP33, exoH34, and exoR9 along with some other genes that were induced at low pH, such as phrR35 and lpiA3637. The collected genetic evidence indicates that tolerance to acidity in these rhizobia is a multigenic phenotype. Such a view is consistent with the failure to convert acid-sensitive rhizobia into acid-tolerant variants simply by the transfer of single cosmid clones generated from a more tolerant genotype (results from our laboratory). A more extensive understanding at the molecular level of the metabolism and responses of rhizobia to low pH remain essential for making an effective use of genetic manipulations in order to enhance bacterial acid tolerance.

Aiming at exploring other alternatives for improving acid tolerance in rhizobia, we attempted to induce transient phenotypes of acid tolerance. As previously demonstrated in E. coli and other bacteria3839, an acid-tolerance response (ATR) could be induced in different rhizobia, including S. medicae40 and S. meliloti41 among others42. An increase in the tolerance to severe acid shocks was observed when rhizobia had been previously cultivated in batch under sublethal acidic conditions4142. Tiwari et al.31 had demonstrated that the genes actSR were necessary for the induction of an ATR in S. medicae. Later on, Draghi et al.41 demonstrated that the ATR in S. meliloti was coupled to an increased symbiotic competitiveness for nodule occupancy, opening the possibility of exploring the stabilization of the ATR in rhizobia to be inoculated into acid soils. Unfortunately, no data were available on the molecular changes associated with that acid adaptation.

The subsequent approaches to improve the basic understanding of the rhizobial responses to extracellular acidity were based on different experimental configurations, all directed at characterizing the genetic expression and protein profile under either acid growth or acid shock. Transcriptional fusions generated by Tn5 transposition demonstrated that acid-induced promoters were associated with cytochrome synthesis, potassium-ion cycling, and lipid biosynthesis and transport among other functions43. After performing a complementary study using proteomic tools, Reeve et al.43 suggested that the folding, proteolysis, and transport processes were key activities in S. medicae growing in acidity. While all these studies were performed on low-pH batch cultures, Hellweg et al.44 undertook an in-depth time-course analysis of the transcriptomic response in S. meliloti after shifting the bacterial cells from neutral to acid pH. The work revealed that a short-term exposure of the bacteria to low pH was sufficient to induce significant transcriptional changes in diverse rhizobial genes. All these studies provided useful—albeit fragmented—information on the phenotypes, and in certain examples, on the molecular components displayed by the rhizobia after challenge with a high concentration of extracellular hydrogen ions.

To date, most physiological and genetic studies on acid-stress in rhizobia have used batch cultures fraught with the consequent heterogeneity in cell physiology during cultivation, and also lacking an external control on the growth rate and other culture variables. The alternative use of continuous cultivation in a chemostat, already employed for the analysis of acid stress in rhizobia45, is an appropriate, and indeed preferable, choice for a more robust experimental design through the reduction of variables down to only the extracellular pH, thus leading to physiologically homogeneous bacterial populations and subsequent datasets46. Having used both classical and various omic tools, in this work we present a consolidated analysis of the responses of S. meliloti growing under a controlled acid stress. Changes in the transcriptome, proteome, and metabolome were integrated into a model to describe the cultural responses in the chemostat, together with the modifications in the enzymes and compounds of the central metabolic pathways.

Results

Cultivation parameters of S. meliloti 2011 growing in continuous culture in chemostat either at pH 7.0 or at nearly growth-limiting acidity

A continuous culture of S. meliloti 2011 was established at pH 7.0 in Evans minimal medium with glucose and ammonium as the respective carbon and nitrogen sources. The cells were cultured under the conditions indicated in Materials and Methods (i.e., D = 0.07.h−1, 28 °C, dissolved oxygen at saturation). Because of the C:N ratio in Evans defined medium the bacterial culture became N-limited at pH 7.0. Under steady-state conditions at pH 7.0 the culture reached a cell density of ca. 4.7 × 109 c.f.u.ml−1 (OD600nm = 4.6) at a biomass of 1.66 g.l−1 dry weight (Table 1). In order to determine the lowest extracellular pH at which S. meliloti was still able to grow in the chemostat, the pH of the extracellular medium was lowered stepwise at intervals of 0.2 pH units to achieve discrete steady states at each new pH (between pHs of 7.0 and 6.0). When the extracellular pH reached the value 6.0 (the acidic pHlimit), the culture washed out, indicating that the rhizobia either stopped growing or duplicated at a rate that was lower than the one imposed by the chemostat-dilution rate. The continuous cultivation provided an optimal design for determining the lowest extracellular pH that was still compatible with the duplication rate imposed by the experimental setting. The acidic pHlimit in the chemostat proved to be comparable to those previously obtained in batch cultivation for S. meliloti and S. medicae (i.e., between 5.8–5.6)2347. This information on the acidic pHlimit was used in a subsequent experiment to set up a continuous culture at pH 6.1 in order to have rhizobial cells growing just 0.1 pH unit above the condition where they had reached their limit of acid tolerance. Table 1 lists the cultivation parameters that characterized the growth of S. meliloti in the chemostat at pH 7.0 and at pH 6.1, including the data on the biomass production, number of culturable cells, and respiration indices. Under acidity a clear change in the limiting substrate occurred (an excess of near 7 mM ammonium was measured in the supernatant) accompanied by a 53% decrease in the biomass. The acid stress in the rhizobial cells at pH 6.1 was clearly reflected in that only 26% of the total bacteria counted in the Petroff-Hausser chamber were culturable (the number of living cells was estimated by plating). The more than 5-fold increase in the O2 consumption per living cell was observed in the culture at pH 6.1 compared to the same parameter in the culture at pH 7.0 indicated a higher aerobic respiration in the acid-stressed rhizobia. Finally, while cells produced significantly more EPS under acidity (YEPS/s; cf. Table 1), no major change in the amount of polyhydroxybutyrate per cell was detected. The connection between EPS production and acid tolerance in rhizobia has been suspected since long. In S. meliloti, however, there is no evidence that EPS could have a positive effect on the bacterial acid tolerance. It has been previously reported that an exoY S. meliloti mutant, which lacks EPS, did not display any significant difference in its death rate when exposed at low pH40.

Table 1

Cultivation parameters of S. meliloti 2011 grown in the chemostat at extracellular pH 7.0 and pH 6.1.

Cultivation parameters	Chemostat pH 7.0	Chemostat pH 6.1
Dilution rate (h−1)	0.073 ± 0.002	0.070 ± 0.003
Optical Density(600nm)	4.57 ± 0.02	1.82 ± 0.12
Bacterial cell number (cells.ml−1)a	4.70 × 109 ± 0.9 × 109	1.79 × 109 ± 1 × 109
Culturable cells/total bacteria (%) (c.f.u. × 100 total bacteria−1)b	100%	26%
Dry weight (mg.ml−1)	1.66 ± 0.11	0.78 ± 0.05 (0.20 g of these being able to grow)b
Cells/Dry weight (cells.mg−1)	2.8 × 109 ± 0.6 × 109	2.3 × 109 ± 1.3 × 109
Dry weight/cell (pg.cell−1)	0.35 ± 0.02	0.43 ± 0.02
EPS production (mg.mg dry weight−1)	0.04 ± 0.22	0.24 ± 0.06
PHB content (% w.w−1)	6.8 ± 0.12	5.3 ± 0.04
Glucose consumption (g.l−1)	3.04 ± 0.15	1.46 ± 0.08
Ammonia (limiting nutrient at pH 7.0) (mM)	ND	6.9
rCO2 (mmol.l−1.h−1)	3.58 ± 0.07	1.66 ± 0.01
(mmol.C-mol−1 h−1)	56 ± 2	55 ± 2
CO2 production/culturable cells (fmol.cfu−1.h−1)	0.76 ± 0.01	3.56 ± 0.06 [↑ ca. 4.7 times]
rO2 (mmol.l−1.h−1)	3.97 ± 0.1	2.19 ± 0.1
(mmol.C-mol−1 h−1)	62 ± 1	72 ± 1
O2 consumption/culturable cells (fmol.cfu−1.h−1)	0.8 ± 0.01	4.7 ± 0.2 [↑ ca. 5.8 times]
RQ (rCO2.rO2−1)	0.9 ± 0.04	0.76 ± 0.04
YX/S (C-mol.C-mol−1)	0.64 ± 0.04	0.62 ± 0.06
(mol.C-mol−1)	0.48 ± 0.01	0.49 ± 0.01
YEPS/S (YP/S estimation) (C-mol.C-mol−1)	0.02 ± 0.01	0.13 ± 0.03
C balance ( + YX/S + YEPS/S)	1.14 ± 0.01	1.24 ± 0.1

All values are expressed with their corresponding standard deviations.

rO2: rate of O2 consumption (mmol.l−1.h−1).

qO2/X: mmol O2 consumption/C-mol of biomass/h.

rCO2: rate of CO2 production (mmol.l−1.h−1).

qCO2/X: mmol of CO2 production/C-mol of biomass/h.

Yx/s: biomass-growth yield (C-mol of biomass produced/C-mol of consumed glucose).

YCO2/S: CO2 production (molCO2/C-mol of consumed glucose).

YEPS/S: exopolysaccharide production (C-mol EPS/C-mol of consumed glucose).

RQ: Respiratory Quotient.

ND: Non-detected.

aNumber of bacteria counted in Petroff-Hausser chamber.

bProportion of bacteria counted by plating with respect to those counted in the Petroff-Hausser chamber (%).

Identification of protein markers induced during the continuous cultivation of S. meliloti 2011 at the acid pH 6.1 compared to neutrality

The cells collected from the continuous cultures were used to investigate the presence of molecular markers specifically associated with growth at either an acid pH or neutrality. In contrast to the limitations imposed by studies on batch cultures, bacterial cultivation in the chemostat at a constant dilution rate (D) provided a suitable experimental protocol for keeping rhizobia growing at the same population-doubling time, at an extracellular pH of either 6.1 or 7.0. The rhizobia in the culture at pH 6.1 continued to grow at 0.1 pH unit above their acidic pHlimit (i.e., pH 6.0, cf. the previous section).

The cytoplasmic proteomes from cells harvested from each chemostat culture were isolated as described in Materials and Methods and analyzed by 2-D–gel electrophoresis. Figure 1A,B show representative gels corresponding to the rhizobia grown at pH 6.1 and at 7.0, respectively. The identification of polypeptides was achieved by UV-MALDI-TOF™ peptide-mass fingerprinting (cf. Materials and Methods). The gel position of most identified products corresponded to the expected molecular mass and isoelectric point of each protein species as inferred from the S. meliloti genomic sequence (Fig. S1). With the aid of the ImageMaster 2DTM software we automatically recognized 382 protein spots (the red-outlined spots in Figs S2 and S3), with those representing more than 6% of the total translation products predicted from the genomic sequence48. Such a proportion is consistent with the 13% previously reported by Djordjevic et al.49 when they analyzed the S. meliloti proteome under five different physiologic conditions. In some instances, different molecular forms of a same protein were found at slightly discrepant positions in the gel, an observation that was more frequent in the proteome from rhizobia grown at pH 6.1 (i.e. five molecular forms for Pnp, three for DegP1 and FusA, and two for GyrB and Tig; Fig. 1A) than in the rhizobia grown at pH 7.0 (i.e. four molecular forms for DppA2 and two for LivK; Fig. 1B).

Figure 1

2-D–gel–electrophoretic profile of S. meliloti 2011 cytosolic proteins from cells grown in the chemostat at pH 6.1 and at pH 7.0.

Isoelectric focussing of 100 μg of cytosolic proteins (first dimension, horizontal) was carried out on 24 cm IPGTM (GE) strips with a final 4 to 7 pH gradient. The gel shown is representative of 4 independent technical replicates. The labels in the figure are the names of the polypeptides identified by UV-MALDI-TOF peptide-mass fingerprinting (Materials and Methods). Panel A Markers overexpressed at pH 6.1 compared to their relative expression at pH 7.0 (gel in B) were detected with the aid of the ImageMaster 2DTM software (Fig. S2) and by visual inspection, and are all listed in Table 2. Panel B Markers overexpressed at pH 7.0 compared to their relative expression at pH 6.1 (gel in A) were also detected with the aid of the ImageMaster 2DTM software (Fig. S3) and by visual inspection, and are all listed in Table 3.

After a comparative analysis of the 2-D gels, both by visual inspection and by using the ImageMaster 2DTM software, we detected 43 and 28 protein markers specifically overexpressed under the conditions of acid-stress and neutrality, respectively (cf. the list of differentially expressed markers at pH 6.1 and pH 7.0 in Tables 2 and 3, respectively). Two types of such markers were defined: class-I, corresponding to polypeptides detected under only one of the two conditions analyzed (i.e., either in cells from the acidic or the neutral culture: that is, markers with an on-off expression pattern); and class-II, corresponding to polypeptides present under both conditions but expressed at different relative intensities. Out of the 71 differentially induced polypeptides, only 9 were class-I markers (on-off). Of these 9, 8 were found in the rhizobia grown at pH 7.0 (SMc00242, SMc01827, IlvE1, AglE, ExoN2, SMa1507, SMa2361, and SMb20605), and only one was found in the rhizobia grown at pH 6.1 (GyrB). Most of the 62 class-II markers corresponded to proteins overexpressed at pH 6.1, thus revealing a higher set of induced proteins under the acid-stress conditions.

Table 2

S. meliloti 2011 protein markers differentially overexpressed under continuous cultivation of the rhizobia at extracellular pH 6.1 compared with the expression of the same markers at pH 7.0.

Gel Code (a)	ORF	Protein name (a)	Accession number	Sequence coverage (%)	Score (MASCOT)	Ratio of intensities (pH 6.1/pH 7.0) (b)	M Value (h) (transcriptome)	Predicted/putative function	COG (i)
1	SMc02365	DegP1 (c)	Q52894	41	154	4.80	1.97	Probable serine protease do-like precursor (EC 3.4.21.-).	O
2	SMc02365	DegP1 (c)		34	148
3	SMc02365	DegP1 (c)		38	181
4	SMc02050	Tig (d)	Q92Q12	46	181	1.49	1.55	Trigger factor (TF).	O
5	SMc02050	Tig (d)		47	193
6	SMc02782	GyrB (e)	Q92TE4	24	136	ON	0.59	DNA gyrase subunit B (EC 5.99.1.3).	L
7	SMc02782	GyrB (e)		10	81
8	SMc03157	–	Q92LX5	39	155	2.78	1.23	Outer membrane lipoprotein 3 precursor (PLP3).	P
9	SMb20724	–	Q92TR9	30	84	2.00	0.65	Hypothetical protein.	S
10	SMc01312	FusA (f)	Q92QH2	25	116	1.43	0.61	Elongation factor G (EF-G).	J
11	SMc01312	FusA (f)		41	192
12	SMc01312	FusA (f)		17	109
13	SMc00324	Pnp (g)	Q92SW0	35	180	1.37	1.97	Polyribonucleotide nucleotidyltransferase (EC 2.7.7.8).	J
14	SMc00324	Pnp (g)		18	111
15	SMc00324	Pnp (g)		20	129
16	SMc00324	Pnp (g)		30	210	4.98	5.36
17	SMc00324	Pnp (g)		20	111	4.98	5.40
18	SMc03978	Tkt2	Q92M82	18	118	2.25	1.11	Transketolase (EC 2.2.1.1).	G
19	SMc00043	SodB	Q9XD74	48	82	1.60	0.28	Superoxide dismutase [Mn] (EC 1.15.1.1).	P
20	SMc02560	ChvI	P50350	32	68	1.86	1.23	Transcriptional regulatory protein ChvI.	T-K
21	SMc01126	Tme	O30808	27	132	2.28	−0.27	NADP-dependent malic enzyme (EC 1.1.1.40).	C
22	SMc00861	–	Q92KK7	14	62	2.00	0.40	Putative signal peptide protein.	R
23	SMc01326	TufB	Q925Y6	30	1,82	1.82	0.88	Elongation factor Tu (EF-Tu).	J
24	SMc02692	RplY	Q92N68	55	200	1.73	1.87	Probable 50S ribosomal protein L25.	J
25	SMc00565	RplI	Q92QZ9	29	68	1.74	NA	Probable 50S ribosomal protein L9.	J
26	SMc01333	PrfB	Q92QJ2	16	59	3.47	1.21	Peptide chain release factor 2 (RF-2).	J
27	SMc00357	Efp	Q92ST6	33	65	1.90	0.40	Elongation factor P (EF-P).	J
28	SMc01917	NuoE1	P56909	29	80	1.59	0.97	NADH-quinone oxidoreductase chain E 1 (EC 1.6.99.5).	C
29	SMc01845	–	Q92NK3	18	86	1.51	1.23	Membrane-bound lytic mureintransglycosylase B precursor (EC 3.2.1.-).	M
–	SMc01285	RpoA	Q925Z2	31	93	SVNA	0.55	DNA-directed RNA polymerase alpha chain (EC 2.7.7.6).	K
–	SMc01318	RplL	Q92QH8	66	88	SVNA	0.47	50S ribosomal protein L7/L12.	J
–	SMc00335	RpsA	P14129	43	233	SVNA	1.86	30S ribosomal protein S1.	J
–	SMc01934	ProS	Q92QN2	53	168	SVNA	0.77	Putativeprolyl-tRNAsynthetase (EC 6.1.1.15).	J
–	SMc00480	Icd	Q92PG6	37	148	SVNA	0.49	Isocitrate dehydrogenase [NADP] (EC 1.1.1.42).	C
–	SMc01270	AdhC1	Q92QD7	37	144	SVNA	0.48	Alcohol dehydrogenase class III (EC 1.1.1.1).	C
–	SMc03826	ArgG	Q92L73	46	189	1.24	1.46	Argininosuccinatesynthase (EC 6.3.4.5).	E
–	SMc00643	PurA	Q92MA5	40	175	SVNA	1	Adenylosuccinatesynthetase (EC 6.3.4.4).	F
–	SMc02064	SerS	Q92Q22	41	150	SVNA	0.48	Seryl-tRNAsynthetase (EC 6.1.1.11).	J
–	SMc00641	SerA	Q92MA3	25	86	SVNA	0.71	D-3-phosphoglycerate dehydrogenase (EC 1.1.1.95) (PGDH).	E
–	SMc01192	MetS	Q92PX0	28	109	SVNA	0.26	Methionyl-tRNAsynthetase (EC 6.1.1.10).	J
–	SMc04461	TolB	Q926C2	42	162	SVNA	0.75	TolB protein precursor.	U
–	SMc02686	PrsA	Q92N73	24	55	SVNA	NA	Probable ribose-phosphate pyrophosphokinase.	F
–	SMc00153	—	Q92PB9	56	95	SVNA	1.80	Probable ATP synthase delta chain (EC 3.6.3.14).	S
–	SMc00595	NdK	Q92QX9	25	82	SVNA	0.88	Nucleoside diphosphatekinase (EC 2.7.4.6).	F
–	SMc01288	Adk	Q93FE6	55	124	SVNA	NA	Probable adenylatekinase.	F
–	SMc03242	TypA	Q92LG8	27	140	SVNA	3.39	GTP-binding protein typA/BipA.	T
–	SMc03925	Pgm	Q92M12	31	81	SVNA	1.39	Phosphoglucomutase (EC 5.4.2.2).	G
–	SMc02100	Tsf	Q92Q54	62	174	SVNA	−0.19	Elongation factor Ts (EF-Ts).	J
–	SMc02163	Pgi	Q92SC4	28	101	SVNA	1.43	Glucose-6-phosphate isomerase (EC 5.3.1.9).	G
–	SMc02495	Tal	Q92LK3	38	80	SVNA	1.70	Probable transaldolase (EC 2.2.1.2).	G
–	SMc00040	—	Q92RC9	27	67	SVNA	1.04	Hypothetical protein.	O
–	SMc01233	Ssb	P56898	64	189	SVNA	NA	Probable single-strand binding protein.	L
–	SMc01345	AccC	Q92QK1	30	117	SVNA	0.25	Biotin carboxylase (EC 6.3.4.14).	I

a: The code numbers (first column) and protein names (third column) correspond to those as in Fig. S2 and Fig. 1A, respectively. The polypeptides from each spot were identified by means of UV-MALDI-TOF MS and the peptide-mass fingerprint as described in Materials and Methods. The bold data corresponds to polypeptides detected as overexpressed both by the ImageMaster 2D™ (GE) software (using 50% volume difference compared to the spot a pH 7.0, and p test lower than 0.05 in a t-test), and by (a priori) comparative visual inspection of the 2D gels (i.e. by the comparison of spots from gels shown in Fig. 1A,B).

b: Ratios of spot intensities (spot intensity at pH 6.1 vs. spot intensity of the homolog spot at pH 7.0) as estimated by the spot volume function of the ImageMaster 2D™ software.

c–g: The same letters correspond to different isoforms and/or/modifications or degradation products of a same protein detected at different positions in the 2D gel from Fig. 1A (or its replicates).

h: M values correspond to the log2 of the relative transcriptional activity of the indicated gene at pH 7.0 compared to the corresponding activity at pH 6.1, as listed in Supplementary Table S1.

i: COG classes as previously defined elsewhere103. Cf. consolidated COG analysis from the differential proteome and transcriptome in Fig. 2.

NA: Transcriptome data not available.

SVNA: “Spot volume not available” when using the ImageMaster sofware.

ON: Polypeptides that are present only in the cytoplasm of rhizobia grown at pH 6.1 (absent in the homologous sample from rhizobia grown at pH 7.0).

Table 3

S. meliloti 2011 protein markers differentially overexpressed under continuous cultivation of the rhizobia at extracellular pH 7.0 compared with the expression of the same markers at pH 6.1.

Gel Code (a)	ORF	Protein name (a)	Accession number	Squence Coverage (%)	Score (MASCOT)	Ratio of intensities (pH 7.0/pH 6.1) (b)	M Value (h) (transcriptome)	Predicted/putative function	COG (i)
1	SMa2361	—	Q92XH5	43	162	ON	−3.05	Hypothetical protein.	R
2	SMc00786	DppA1	Q92RW0	18	87	2.21	−1.37	Periplasmic dipeptide transport protein precursor.	E
3	SMc01525	DppA2 (c)	Q92N26	41	201	5.46	−0.77	Periplasmic oligopeptide-binding protein precursor.	E
4		DppA2 (c)		46	156			Periplasmic oligopeptide-binding protein precursor.	E
5		DppA2 (c)		42	220			Periplasmic oligopeptide-binding protein precursor.	E
6		DppA2 (c)		41	168			Periplasmic oligopeptide-binding protein precursor.	E
7	SMb21196	OppA	Q9AKR0	25	72	SVNA	NA	ABC transporter of tetrapeptides and some tripeptides, periplasmic solute-binding protein.	E
8	SMc01946	LivK (d)	Q926C5	47	144	30.25	−2.58	Leucine- isoleucine- valine- threonine- and alanine-binding protein precursor.	E
9		LivK (d)		31	132			Leucine- isoleucine- valine- threonine- and alanine-binding protein precursor.	E
10	SMc00242	—	Q92PR4	27	87	ON	NA	Hypothetical signal peptide protein.	R
11	SMc01827	—	Q92K04	18	57	ON	−2.58	Hypothetical protein.	P
12	SMb20605	—	Q926E2	26	93	ON	−4.80	Aliphatic amidase expression-regulating protein.	E
13	SMb20745	GlnII	Q92TR3	20	58	5.00	−2.71	Glutamine synthetase II (EC 6.3.1.2).	E
14	SMb20964	—	Q9R9R2	38	101	1.44	NA	Peroxidase	P
15	SMc04023	ExoN2	Q92M48	25	110	ON	−0.40	UTP-glucose-1-phosphate uridylyltransferase (EC 2.7.7.9)	M
16	SMa1353	—	Q92YX3	27	70	4.55	−0.58	D-tagatose 3-epimerase (EC 5.3.1.-).	G
17	SMc02896	IlvE1	Q92SX8	23	59	ON	1.12	Probable branched-chain amino acid aminotransferase (EC 2.6.1.42).	E-H
18	SMc03061	AglE	Q9Z3R5	16	88	ON	0.53	Alpha-glucosides-binding periplasmic protein AglE precursor.	G
19	SMb20292	—	Q92WP7	51	143	1.51	−1.42	Hypothetical immunogenic protein.	R
20	SMa1507	—	Q92YP3	33	80	ON	−1.75	Transcriptional activator protein precursor.	R
–	SMc02356	LivA	Q92MNO	25	111	SVNA	−1.57	Putative branched chain amino acid binding periplasmic ABC transporter.	E-T
–	SMc02873	—	Q92SZ7	34	102	SVNA	−0.41	Multiple sugar-binding protein precursor.	G
–	SMc02118	AapJ	Q92Q71	48	187	SVNA	−2.82	General L-amino acid-binding periplasmic protein.	E-T
–	SMc03786	Bfr	Q92LA7	55	110	SVNA	−1.51	Bacterioferritin (BFR).	P
–	SMc02259	—	Q92S63	46	146	SVNA	−1.14	Lysine-arginine-ornithine-binding periplasmic protein precursor.	E-T
–	SMb21144	—	Q926G9	33	233	SVNA	−0.56	Putative choline uptake ABC transporter periplasmic solute-binding protein precursor.	M
–	SMc01605	—	Q92NI5	37	161	SVNA	−1.17	Putative periplasmic binding ABC transporter protein.	P
–	SMb20915	AslA1	Q92UC0	34	161	SVNA	−2.46	Putative arylsulfatase.	P
–	SMc00140	—	Q92PA9	55	125	SVNA	−1.30	Arginine-binding periplasmic protein 2 precursor.	E-T
–	SMc01642	PrbA	Q92NF1	27	97	SVNA	−1.38	Periplasmic dipeptide transportprotein precursor.	E
–	SMc00777	ThrC1	Q92RF5	38	116	SVNA	−0.45	Probable threonine synthase.	S

a: The code numbers (first column) and protein names (third column) correspond to those as in Fig. S3 and Fig. 1B, respectively. The polypeptides from each spot were identified by means of UV-MALDI-TOF MS and the peptide-mass fingerprint as described in Materials and Methods. Bold data corresponds to polypeptides detected as overexpressed both by the ImageMaster 2D™ (GE) software (using 50% volume difference compared to the spot a pH 6.1, and p test lower than 0.05 in a t-test), and by (a priori) comparative visual inspection of the 2D gels (i.e. by the comparison of spots from gels shown in Fig. 1A,B).

b: Ratios of spot intensities (spot intensity at pH 7.0 vs. spot intensity of the homolog spot at pH 6.1) as estimated by the spot volume function of the ImageMaster 2D™ software.

c–g: The same letters correspond to different isoforms and/or modifications or degradation products of a same protein detected at different positions in the 2D gel from Fig. 1B (or its replicates).

h: M values correspond to the log2 of the relative transcriptional activity of the indicated gene at pH 7.0 compared to the corresponding activity at pH 6.1, as listed in Supplementary Table S1.

i: COG classes as previously defined elsewhere103. Cf. consolidated COG analysis from the differential proteome and transcriptome in Fig. 2.

NA: Transcriptome data not available.

SVNA: “Spot volume not available” when using the ImageMaster sofware.

ON: Polypeptides that are present only in the cytoplasm of only rhizobia grown at pH 7.0 (absent in the homologous sample from rhizobia grown at pH 6.1).

Furthermore, clear differences were evident in the predicted functions associated with the induced markers under each pH condition (cf. Tables 2 and 3, last two columns). At neutral pH most overexpressed polypeptides (75% comprising class-I and -II markers) corresponded to genes belonging to clusters of orthologous groups (COGs) of transport and metabolism (Fig. 2, Panel A, blue bars). In the rhizobia growing at extracellular pH 6.1 most induced polypeptides corresponded to genes pertaining to COGs related to translation (28%), post-translational modifications (7%), and energy production and conversion (9%; Fig. 2, Panel A, red bars), thus indicating a change in the functional priorities under acid stress. While at pH 7.0 the transport and biosynthesis of cellular compounds emerge as quite active processes, under acid stress most of the overexpressed markers were associated with protein biosynthesis and energy metabolism. In the rhizobia growing at pH 6.1 the ribosomal polypeptides L25 (RplY) and L9 (RplI) were increased as well as the ATP-dependent chaperone Tig, a peptidylprolyl-cys-trans isomerase which interacts with short nascent polypeptide chains50. Likewise, related to translation at pH 6.1, a significant increase in the release factor PrfB and in the elongation factors EF-TU (TufB), EF-P (Efp), and EF-G (FusA) occurred; which proteins participate in the binding of aminoacyl-tRNAs to the ribosome and in the translocation of growing peptidyl-tRNAs51525354. These results, likely related to different needs for protein synthesis in the acid-stressed rhizobia, are consistent with recent publications where specific changes in ribosomal-protein composition were also observed in the responses of bacteria55 as well as in different eucaryotes5657 to stress. In our experiments, the observed changes appeared to be restricted to a set of translation-related proteins with no significant modifications in the cellular amounts of the 16S and 5S rRNAs (Fig. S4). All these results point to the higher demands on protein synthesis in the acid-stressed rhizobia that grow in the chemostat at the same population-doubling time as those at pH 7.0. The parallel increase in the protease DegP58 also suggests a higher rate of protein degradation and turnover under conditions of acidity. This protease had been previously reported to be necessary for the normal growth of S. meliloti under high-temperature stress59, and might likely be the consequence of an increased protein damage. Consistent with such a scenario, an overexpression of the Pnp ribonuclease60 was also observed at pH 6.1, likewise suggesting a parallel increase in RNA degradation.

Figure 2

Incidence of the main COGs among the differentially expressed markers at pH 6.1 and 7.0 as revealed by the proteomic (A) and transcriptomic (B) analyses. The abundance of the COG classes indicated on the abscissa is expressed as a percent (ordinate). pH 7.0: blue ; pH 6.1: red.

An analysis of Tables 2 and 3 also indicated that rhizobia in the acidic chemostat culture induced membrane-associated proteins such as the translation products of SMc03157 (a lipoprotein with homologs related to transport processes61), SMc01845 (a peptidoglycan-binding protein related to certain cell-wall–degradation enzymes6263), and the response regulator ChvI, which protein interacts with the membrane histidine-kinase receptor ExoS64. ChvG/ChvI, the Agrobacterium tumefaciens homolog of ExoS/ChvI, was found to be involved in the tolerance to acid stress65. In S. meliloti, this two-component system has been shown to be physically associated with the periplasmic ExoR as well as being involved in EPS biosynthesis and symbiosis66, with several molecular targets of this regulon having been recently identified67.

Transcriptomic analysis of S. meliloti growing in the chemostat at pH 6.1 compared to that of rhizobia growing at neutral pH

With the aim at exploring the profiles of gene expression in acid-stressed and control rhizobia, we performed a global transcriptional analysis in microarrays. The relative gene expression (M value) could be estimated for 3,750 transcripts (p < 0.05), representing about 60% of the total number of genes. Of these loci, a higher proportion of those being induced occurred under acid stress. For example, we observed 393 transcripts induced at pH 6.1 (M ≥ 1) compared to 225 transcripts induced at pH 7.0 (see Table S1). The analysis of COGs in the induced transcripts was in striking parallel with the results obtained from the rhizobial proteomes (Fig. 2, Panel B compared to Fig. 2, Panel A), except for membrane associated markers. Several induced markers were associated with translation and amino-acid transport plus metabolism at pH 6.1 and 7.0, respectively. When the transcriptomic data became available for the differentially induced polypeptides listed in Tables 2 and 3 (columns 7 and 8), concordant results were observed with positive M values for the polypeptides overexpressed under acidity and with negative M values resulting for polypeptides overexpressed at neutral pH.

The transcriptomic analysis was particularly useful for investigating gene expression in envelope-associated markers that were excluded from the present cytosol-proteome analysis. Thus, genes of the motility apparatus and the chemotactic machine were observed to be repressed under the condition of acidity. The repressed genes included flaA, flaB, flaD (flagelins); flgB, flgF, flgG (flagellar basal-body rod proteins); flgL (flagellar hook-associated protein); fliE (flagellar-hook basal-body protein); fliM (flagellar motor-switch transmembrane protein); and cheE, motB, and motD (chemotaxis proteins). Conversely, at pH 6.1 we observed an induction of the genes murE (encoding an UDP-N-acetylmuramoylalanyl-D-glutamate-2,6-diaminopimelate ligase; M = 1.9) and murI (encoding a glutamate racemase; M = 2.1); which enzymes participate in peptidoglycan biosynthesis, an observation that suggests an intermembrane-space remodelling under acid stress. Likewise, the accompanying induction of genes such as accA (acetyl-CoA carboxylase carboxyltransferase subunit alpha; M = 2.3) and plsX (acyl-acyl carrier protein phosphate acyltransferase; M = 1.6)—those enzymes related to fatty-acid and phospholipid biosynthesis—points to modifications in the lipid composition of the membrane under acidity. Changes in membrane components related to energy metabolism were also evident from these transcriptome data, particularly with respect to the genes encoding electron transporters such as cycM (cytochrome c; SMc02897; M = 2.6), cycB2 (SMa1170; M = 2.0), azu1 (SMa1243; M = 2.3), fixN1 (SMa1220; M = 2.1), and fixG (SMa1211; M = 2.2), all of which loci increase their transcription at pH 6.1. The parallel increase in the transcriptional activity of genes encoding subunits of the membrane ATPase (atpA, the gene of the subunit alpha of the F0F1 ATP synthase; M = 1.6) and the NADH dehydrogenase (nuoA1, nuoB1, nuoC1, and nuoE1; genes encoding the subunits A, B, C, and E of the enzyme; M values from 2.5 to 1.4) strongly indicates the occurrence of a higher energy demand, electron-transport rate, and degree of oxygen consumption (Table 1) under acid stress.

Another key group of markers that frequently escape in-gel proteome analyses are those related to nucleic-acid metabolism and transcriptional regulation, owing to either their usually high isoelectric-point values—many are basic proteins—or their low cellular concentrations. In rhizobia growing at pH 6.1 we detected a higher expression of genes that encode transcription factors, and members of two-component systems such as: SMa0849 (syrM), SMc00929 (putative homolog of gcvA), SMa0760 (a gene encoding the antikinase FixT2), SMa1179 (nosR), SMc00458 (feuP, encoding a response regulator), SMa0815 (nifA), SMc02366 (ragA, encoding a response regulator), SMb20683 (encoding a transcriptional regulator of the LysR family), SMc00681 (lrp), SMc02560 (chvI, encoding a response regulator also detected in the proteome), SMc00109 (gene for a probable transcriptional regulator), SMa1705 (encoding a transcriptional regulator of the MucR family), and SMc01507 (encoding a sensor histidine kinase68). Of these genes, the S. meliloti feuP has been previously shown to be required for the expression of a wild-type hypoosmosis-tolerant phenotype through correct exportation of cyclic β-glucans69. These results therefore indicate that low pH induces central regulators which change the saccharidic composition of the bacterial envelope. Since in other bacteria gcvA —the SMc00929 homolog— and chvI have both been associated with responses to acid stress7071, these genes likely participate in an evolutionarily conserved response mechanism.

Finally, the transcriptome made evident as well the induction of several genes associated with the increase in general transcriptional and translational activities, supporting the results discussed in the previous sections. Accordingly, at low pH higher transcriptional activities of genes encoding ribosomal proteins, subunits of the RNA polymerase [SMc01317 (rpoB, gene for the DNA-directed–RNA-polymerase beta chain), SMc02408 (rnpO, gene for the DNA-directed–RNA-polymerase omega subunit)], and the Rho termination factor [SMc02796 (rho)] have been detected. An increase in RNA metabolism may likewise be inferred from the enhanced transcription of the ribonuclease genes SMc01365 (rnr, a probable exoribonuclease II) and SMc02652 (rnc, a probable ribonuclease III) and from the increase in different RNA helicases [e.g., SMc00522 (rhlE) and SMc03877 (mgpS)], which enzymes in prokaryotes have already been associated with housekeeping functions and with structural adaptations of RNAs and ribosomes to abiotic stress72.

Genome distribution of acid-induced markers

The proteome and transcriptome data were used to map, within the S. meliloti genome, the position of the markers with increased expression under the acidic condition (cf. Fig. S5). Consistent with the observation that many acid-induced markers are involved in central metabolic activities (Tables 1 and 2, and Supplementary Table S1), the S. meliloti chromosome proved to be the replicon with the highest density of acid-induced genes (on the average, ca. 8 genes with M ≥ 1 per 100 kb). With respect to the two symbiotic megaplasmids, pSymA encoded an average of 5 acid-induced genes per 100 kb, compared to the ca. 3 per 100 kb present in pSymB. Within this general context, the presence of a region of near 120 kb on the pSymA where we observed the striking value of ca. 21 acid-induced markers per 100 kb (indicated with asterisk on the pSymA in Fig. S5) is remarkable. The higher GC content within that region compared to that of the complete pSymA (61.36% vs. 60.4%) suggests that several of these genes might be part of a horizontally acquired region within the pSymA DNA. This region includes ORFs that codify for a transcriptional regulator (SMa1207), a cation transporter (SMa1087), a DegP-like protease (SMa1128), two Fix (SMa1216 and SMa1220, both encoding cytochrome-c subunits) and two non-Fix (SMa1170, SMa1243) electron transporters, proteins related to nitrous- and nitric-oxid metabolism (SMa1279, NorE, SMa1179, NosR), and many acid-induced hypothetical proteins among other markers.

Metabolome analysis of S. meliloti 2011 grown under either acidic or neutral extracellular conditions

In order to explore the metabolites that characterize the growth of S. meliloti under different pH conditions, we performed the metabolome analysis described in Materials and Methods. The results identified and quantified 66 different compounds (cf. the data in Supplementary Table S2), whose relative amounts per cell dry weight enabled the analysis as to whether or not rhizobia have specific metabolite-compositional patterns associated with their extracellular pH of growth. The MetaboAnalyst software revealed several compounds whose intracellular concentrations significantly changed when the extracellular pH was modified. Figure 3 shows a volcano plot where the blue dots correspond to those compounds with a fold change (FC) higher than 2 (increased at pH 6.1) or lower than 0.5 (increased at pH 7.0), and a t-test threshold with a p-value lower than 0.05 (see the numerical data in Supplementary Table S3). In rhizobia grown under the indicated acidic and neutral pH, 13 and 6 compounds, respectively, became increased. Of these metabolites, glucose-6-P, gluconate-1,5-lactone, fructose, galactose, and phosphoenolpyruvate are all members of the central carbon pathways in S. meliloti73.

Figure 3

Volcano plot showing metabolites with fold changes (FC) ≥ 2 or ≤0.5 (p ≤ 0.05), when the pH of the extracellular culture was modified from 7.0 to 6.1 (FC = amount at pH 6.1/amount at pH 7.0).

The plot shows the −log2 of the cytosolic amount of each metabolite at pH 6.1 with respect to the amount of the same species at pH 7.0. The signals for each metabolite were normalized to a constant cell dry weight and to the ribitol added as an internal standard (cf. Materials and Methods and Supplementary Table 2). The blue circles and names in blue indicate features (i.e., metabolites) with FC either ≥2 or ≤0.5 at a p ≤ 0.05. Calculations were performed by means of the MetaboAnalyst software. The data presented correspond to the statistical analysis of 4 technical replicates for each pH condition.

In order to have an unsupervised method directed at finding metabolome changes associated with pH variations, we performed a Principal-Component Analysis (PCA) using the compositional data obtained for the rhizobia grown under each pH condition. The results demonstrated that the data from rhizobia cultivated at pH 7.0 and 6.1 mapped at different regions in the two-dimensional space of components 1 and 2 (PC1 vs. PC2; cf. Fig. 4, Panel A and the vector-correlation plot in Supplementary Fig. S6). That the metabolome data from cells grown in an independent chemostat at pH 7.4 (which we used as a control sample of rhizobia grown under a slightly basic pH) could also be differentiated from the other two sets of data is indeed interesting, and shows that an extracellular pH difference as low as 0.4 units—and in the alkaline direction—is sufficient to make the metabolic profiles in S. meliloti distinctly different. The cluster analysis presented in Fig. 4, Panel B identified and quantified changes in which specific sets of metabolites (i.e., profiles, metabolic patterns) became positively or negatively correlated with particular values of the extracellular pH. Whereas the cytosolic concentrations of the metabolites from cluster I were positively correlated with the rhizobial growth in the acidic medium at pH 6.1 (with the several having FCs ≥ 2 also being shown in the volcano plot in Fig. 3), the levels of the metabolites from cluster II3 became predominantly increased only at pH 7.4 (i.e., O-acetyl-L-serine, phosphoenolpyruvate, D, L-cystathionine, glycerate, threalose, proline, alpha-ketoglutarate, and alpha-ketoaminobutyrate), thus providing a set of marker compounds for rhizobial growth in a more alkaline medium. Other metabolites, such as those from cluster II4a (i.e., N-acetylglutamate; glutamine; 4-aminobutyrate; ornithine, arginine, and citrulline; along with fructose, among others) increased at both pH 6.1 and pH 7.4 and thus constitute compounds marking culture-pH deviations from neutrality.

Figure 4

Multivariate chemometric and hierarchical-cluster analyses of the metabolome data from S. meliloti 2011 grown at pH 6.1, 7.0, and 7.4.

(A) PCA showing the metabolite distribution within the PC1 vs. PC2 space of variation. The PCA analysis was performed on the cytosolic concentration of each metabolite normalized to a constant cell dry weight and to the signal of a known amount of added ribitol as an internal standard (cf. Materials and Methods plus the data in Supplementary Table 2). The PC1 space (x-axis) served to clearly separate samples from the rhizobia grown under acidity (pH 6.1) from the corresponding ones from the rhizobia grown at pH 7.0 and 7.4, which samples could otherwise be discriminated through the PC2 space (y-axis). The colored areas surrounding the data points represent 95% confidence limits. (B) Cluster display of metabolome data from rhizobia grown at different pHs. Each metabolite is represented by a single row of colored boxes, while each rhizobial sample from cells grown at the indicated pH (with 4 technical replicates each) is represented by a single column. The dendrograms and color images were produced with the MetaboAnalyst software. The color scale indicates the fold change (FC) with respect to the average value over the three pH conditions for each specific metabolite (red: values higher than the average; blue: values lower than the average).

Discussion

Increased hydrogen-ion concentrations in the environment have been demonstrated to strongly affect viability and symbiosis in S. meliloti13, the related species S. medicae9, and several other rhizobia8. In view of these findings and on the worldwide distribution of acidity in agricultural soils, an understanding of the responses of rhizobia to low pH becomes a central issue in the attempt to improve symbiosis in moderately acidic soils. In the present work we investigated the physiologic and molecular responses to extracellular acidity of the model rhizobium S. meliloti. In order to characterize the bacterial changes induced by low pH, we set up S. meliloti steady-state cultures in a chemostat, an experimental system allowing a strict control of cultivation parameters. The use of a chemostat also enables, at the same time, the acquisition of reliable and homogeneous datasets, obviating the secondary growth effects intrinsically associated with the more heterogeneous batch cultures46. The results of the experiments summarized in Table 1 indicated a limiting growth in acidic medium close to pH 6.0–6.1 for S. meliloti 2011; a condition where we observed a significant proportion (ca. 74%) of bacteria that were unable to grow (results from the plating experiments), in agreement with the acidic stress imposed. Within such a context, the parallel increase in oxygen consumption by metabolizing bacterial cells is reflected in a higher energy demand for both repair and biosynthesis under conditions of acidity (i.e., for both maintenance and growth). The observed increase in several of the ribosome-associated proteins (e.g., RplI, RplY, RplL, RpsA, TufB, PrfB, and Efp; cf. Table 2) provided parallel evidence in support of a higher demand of protein synthesis during the growth at lower pH. Most notable is the difference between the two systems of cultivation in that, whereas in order to alleviate an acid stress, in batch cultures rhizobia decrease their duplication rate42; those same bacteria when confronted with an equally low pH in the chemostat are forced to maintain their replication (at a constant dilution rate) at the same rate as that used at pH 7.0.

As demonstrated here, the adaptive process included changes in transcripts associated with membrane components, such as the biosynthetic enzymes of phosphatidylglycerol derivatives (Smc00611, a homolog of lpiA from R. tropici36) and ornithine-containing lipids (Smc01116, namely olsA in S. meliloti74: note, too, that we also observed an increase in the cytosolic concentration of the ornithine substrate under acidity). Both of these lipids were relevant for acid tolerance in R. tropici through different proposed mechanisms (for a recent review see ref. 75). The correlation between the observations in R. tropici and S. meliloti suggests that membrane-lipid remodelling is a general phenomenon associated with acid adaptation in rhizobia, if not in other bacteria as well.

The remarkable increase in aerobic respiration (more than fivefold) under acid stress quantitatively reflects the magnitude of the impact of protons on the rhizobial energy demands. The omic analyses demonstrated an enhanced transcription of chromosomal genes encoding different cytochromes, part of the respiratory ATPase, and several subunits of the membrane-associated NADH dehydrogenase. A significant increase in the transcription of the pSymA-encoded cytochrome genes cycB2, fixN1, and fixG was likewise observed under acidity; in agreement with the previous evidence on the role of pSymA in stress responses76. The increased transcription of fixN1 in particular had also been reported in a strain of S. medicae incubated in batch at low pH77. The present results would indicate that the cellular role of certain of the fix genes exceeds their sole participation in nitrogen fixation during late symbiosis. The observed changes in several cytochrome-associated transcripts could be related to a reinforcement and/or switching-on of alternative destinations for the reducing power generated under acid stress. Other authors have suggested that alternative respiratory pathways may be associated with changes in the P/O ratios under specific circumstances, thus making the system more flexible to respond to environmental variations78.

In view of the known active pathways in S. meliloti73, we could infer that energy and reducing power for biosynthesis should derive from the Entner-Doudoroff pathway and from the pentose-phosphate pathway (PPP). According to our results, the higher expression under acidity (of the mRNAs and proteins) of glucose-6-phosphate dehydrogenase (zwf) together with the transketolase (Tkt2) and the transaldolase (tal) is strongly suggestive of an increased activity in both the oxidative and the nonoxidative branch of the PPP. The observed induction under acidity of phosphoglucose isomerase (Pgi), the first enzyme of the incomplete Embden-Meyerhof-Parnas pathway in S. meliloti, suggests an increased recycling of fructose-6-P synthesized in the nonoxidative branch of the PPP back into glucose-6-phosphate79. In agreement with this possibility, both the substrate fructose-6-P and the product glucose-6-P were found in higher concentrations in the acid-stressed rhizobia. Furthermore, the increased concentration of glucose-6-P and phosphoglucomutase are both consistent with the increased EPS biosynthesis also observed under acidity. The induction of the EPS-biosynthetic genes following a short exposure of rhizobia to low pH has already been reported44. In addition to these considerations and the complex hierarchical regulation of EPS biosynthesis in S. meliloti80, the increased production of EPS in the acidic continuous culture may also have been evoked by a metabolic overflow of carbon via the oxidative PPP, derived from higher demands for the NADPH used for biosynthesis and repair (cf. the schematic summary of these metabolic interactions in Fig. 5). The concept of metabolic overflow was originally introduced in microorganisms to describe the excretions of specific by-products as a consequence of their growth under glucose excess8182838485. The increased EPS production by rhizobia under acid stress resembles microorganism responses during metabolic overflow, where cells dispose of different compound(s) —those derived from central metabolic intermediates— in order to keep specific metabolic pathways active. Nonetheless, cellular needs that impel “classic” metabolic overflows are clearly different from those that operate under the acid stress studied here; where, in contrast to a repression in the respiratory chain, a strong increase in aerobic metabolism occurs. Unfortunately, we have no data on the exometabolome of the rhizobia in the chemostat in order to evaluate the presence of by-products that could be derived from the observed increase in the cytosolic amounts of three-carbon central metabolites (e.g., dihydroxyacetone-3-phosphate and pyruvate, Supplementary Table 3, Fig. 5).

Figure 5

Graphical representation of S. meliloti central carbon-metabolic pathways where significant changes in metabolites, transcripts, and/or specific enzymes were observed when rhizobia were grown under acid stress.

Metabolites in red correspond to those with an increased concentration at pH 6.1 relative to neutrality according to the analysis shown in the volcano plot from Fig. 3 and/or in the cluster from Fig. 4, Panel B. The transcripts and enzymes in red color correspond to induced markers from the differential transcriptome and proteome shown in Supplementary Table 1 and Table 1, respectively. PPP: pentose-phosphate pathway, EDP: Entner-Doudoroff pathway, EMP: Embden-Meyerhof-Parnas pathway, TCA: tricarboxylic-acid cycle, arcA (Smc03091): arginine deaminase, zwf (Smc03070): glucose-6-phosphate dehydrogenase, tal (Smc02495): transaldolase, Tkt2 (Smc03978): transketolase, DHAP: dihydroxyacetone phosphate, PHB: polyhydroxybutyrate, EPS: exopolysaccharide, Pgi: phosphoglucose isomerase, Pgm: phosphoglucomutase.

Acid- and neutrality-grown rhizobia, though clearly recognizable through the differential-expression patterns observed in specific RNAs and proteins, could also be unequivocally distinguished through their cytosolic metabolomes. The possibility of associating metabolomes with extracellular pHs does not appear to be restricted to a comparison of widely different conditions (i.e., neutrality versus a growth-limiting acidity); and especially not so, since we could also clearly differentiate the metabolomes from cells grown at pHs 7.0 and 7.4 (Fig. 4). Thus, we observed a quite strong dependence of the metabolome of the rhizobia on (even limited) changes in the extracellular pH (<0.5 units).

The results presented in this report have demonstrated that S. meliloti has a complex process of adaptation to acidity in agreement with previous findings in support of the multigenic character of the response. The acidic condition of growth modified transcription (|M| ≥ 1, p ≤ 0.05) in nearly 10% of the genome (Supplementary Table 1). We have demonstrated here that the genetic information associated with the growth of S. meliloti under acidity is spread throughout the genome. The chromosomal location of most of the acid-regulated genes indicates that adaptation to acidity requires the readjustment of diverse central functions (e.g., the expression of membrane proteins related to respiration and lipid metabolism and changes in the central carbon metabolism). Besides such general mechanisms of adaptation—those likely existing since the presymbiotic state within the rhizobial evolution—an intriguing expression of specific fix and other pSym-encoded markers was also observed.

The induction of pSym markers by acid may be related to recent results revealing the presence of acidic environments along the infection pathway. Novel data obtained through the use of fluorescent probes demonstrated the existence of an acidic environment at the “bright spot” during the initiation of root-hair infections in M. sativa86 as well as later on during symbiosis in the peribacteroid space within zone III of N2-fixing nodules87. That some acid-responsive pSym markers detected in the acidic culture could also be induced in planta as a primary response against the reported extracellular low pH is thus conceivable. Recent transcriptome data from laser-dissected nodules showed that several pSymA genes induced within infected nodules were also overexpressed under acidic cultivation, thus confirming the versatile role of these pSym replicon88.

All this evidence taken together has indicated that the growth of S. meliloti under stressing acidic conditions in the chemostat resulted in clear biochemical modifications; whose main features are summarized in the scheme depicted in Fig. 5. The differential proteome and transcriptome revealed changes in proteins and genes associated with the bacterial surface (i.e., in lipid metabolism involving cell envelope), alterations in the central carbon metabolism and in respiration, and modifications in the amount of EPS secreted, among others. A new specific metabolome profile could be established, with changes in many compounds that correlated with the observed shifts in the proteome and transcriptome (i.e., intermediates of the PPP pathway). The quantitative rhizobial metabolome emerged as a substantially sensitive indicator of (even subtle) changes in the extracellular pH.

The question thus remains as to what could now constitute promising avenues in order to gain an understanding and improvement of the rhizobial behavior and symbiosis under moderate acidity. Unfortunately, we do not understand the relevance to acid tolerance of most of the differentially expressed markers identified in this present investigation. Thus, the phenotypic impact of specific mutants should now be investigated. Further analyses should also be focussed on the identification of the initial acid-induced molecular events that trigger the downstream biochemical changes described here. Additional approaches to improving the response of S. meliloti to low pH could also be addressed by the identification of the genes and processes that are part of the ATR, through which the bacteria transiently acquire both a higher tolerance to acidity, and an improved symbiotic competitiveness41. We need to learn, for example, how the molecular changes that take place during the ATR compare with the responses observed here. The connection between the transient tolerance to acidity and the concomitant improvement in symbiosis make the process of acid adaptation a suitable target for upcoming research towards an eventual use of the adapted phenotype in preconditioning rhizobia of inoculants.

Materials and Methods

Bacterial strains and growth conditions

Sinorhizobium meliloti Rm 2011, whose fully genomic sequence is available, was used in this work (GenBank accession numbers: NC_020528.1, NC_020527.1, and NC_020560.1). For starter batch cultures and nutrient-limited continuous cultures in chemostat Evans defined minimal medium was used89, containing glucose (10 g.l−1) and ammonium chloride (0.7 g.l−1) as the respective carbon and nitrogen sources. The use of glucose, which is metabolized by S. meliloti, facilitated the calculation of basic cultural parameters in the chemostat such as the yields YX/S, YCO2/S, and YEPS/S presented in Table 1.

Continuous-culture setup

Continuous cultures of S. meliloti were performed at 28 °C in a 2-L Bioflo IIe (New Brunswick Scientific Co., Edison, NJ, USA) reactor with a working volume of 1.5 l. The growth rate (i.e., dilution rate) was adjusted at 0.07 ± 0.01 h−1. The cultures were flushed with filtered sterile air (20 l.h−1) at flow rates measured with a bubblemeter. The dissolved-oxygen concentration was continuously measured with an Ingold (Wilmington, MA, USA) polarographic probe. Oxygen and CO2 concentrations were determined by a paramagnetic oxygen analyzer (Servomex 1100A; Norwood, MA, USA) and an infrared CO2 analyzer (Horiba PIR 2000; Japan), respectively. The rates of oxygen consumption and CO2 production were calculated by a mass-balance method according to Cooney90. The pH was automatically maintained at either 7.0 ± 0.05, 6.1 ± 0.05, or 7.4 ± 0.05 through the addition of 1 N NaOH. Glucose concentrations in the culture media and supernatants were determined with a glucose-oxidase enzymatic kit (Wiener Lab, Argentina) and ammonium concentrations in the media and supernatants were measured by the indophenol method91. The biomass dry weight and rates of oxygen consumption were calculated as previously described92. Bacterial cells were counted in a Petroff-Hausser chamber and the culturable cells determined by plating culture samples in agarized Evans medium. Cultures were considered to be under steady-state conditions when the biomass concentration and specific rate of oxygen consumption varied by less than 10% during three retention times.

Exopolysaccharide (EPS) and polyhydroxybutyrate determination

To isolate the EPS, the supernatant from continuous cultures was precipitated with three volumes of cold ethanol and incubated at −20 °C overnight. The mixture was centrifuged at 10,000 × g for 15 min. The pellet was dissolved and dialyzed against distilled water, lyophilized, and stored at −80 °C until use. The amount of EPS in these samples was determined by the Anthrone Method as previously described93.

The polyhydroxybutyrate content in the cell pellets was determined as chrotonic acid in H2SO4, as described by Law et al.94.

2-D Gel electrophoresis

One liter of culture was harvested and centrifuged at 10,000 × g for 20 min at 4 °C. The pellets were washed twice with cold phosphate-buffered saline; centrifuged as above; and the final pellets frozen immediately in liquid nitrogen, lyophilized, and stored at −80 °C until use.

To obtain the protein extracts, 100 mg of dry cells were suspended in 5 ml of 10 mM Tris-HCl pH 7.6 with Complete Protease Inhibitor Cocktail™ tablets (Roche) and sonicated for 1 min in 5–10 sec cycles in an ice-water bath at 50% power with a Sonifier 150™ (Brandson Ultrasonic Corporation; Connecticut, USA). The samples were then centrifuged at 10,000 × g, for 20 minutes at 4 °C and the supernatants centrifuged at 100,000 × g for 1 h at 4 °C to eliminate cell debris. The supernatants were incubated with 0.5 μg.ml Benzonase Nuclease™ (Sigma Aldrich; St. Louis MO, USA) for 30 min at 37 °C and the proteins finally precipitated with four volumes of cold acetone followed by an overnight incubation at −20 °C. After a 15-min centrifugation at 10,000 × g the pellets were washed twice with cold acetone. The air-dried pellets were finally suspended in 500 μl of rehydration buffer (RB) containing 8 M urea, 2% (w/v) 3-([3-chloramidopropyl]dimethyammonio)-1-propane-sulfonate (CHAPS), and 0.01% (w/v) bromophenol-blue. The total-protein contents were determined by the Bradford method (BioRad kit; Hercules, CA, USA).

Isoelectrofocussing and two-dimensional sodium-dodecyl-sulfate–polyacrylamide-gel electrophoresis (2-D-SDS-PAGE) of protein samples

Five hundred μg of total proteins dissolved in 450 μl of RB were mixed with 5 μl of IPG™ buffer for a strip-immobilized pH gradient from pH 4–7 (GE; Pittsburg, USA), and after the addition of 5 μl of 1 M dithiothreitol the sample was incubated on a shaker at 25 °C for 15 min. Then 450 μl were loaded onto a 24 cm, pH 4–7 Immobiline DryStrip™ (GE) and the proteins electrofocussed on the strips in an IPGphor™ unit (GE) for 1 h at 0 volts, 12 h at 30 volts, 2 h at 60 volts, 1 h at 500 volts, 1 h at 1,000 volts; <15 h at 8,000 volts. Isoelectric focusing was terminated after 75,000 volts h. After the electrofocussing, the strips were equilibrated in 5 ml of a solution containing 6 M urea, 50 mM Tris (pH 8.8), 30% (v/v) glycerol, 2% (w/v) SDS, and 2% (w/v) dithiothreitol on a tilt table for 15 min (the disulfide-reducing step). The solution was discarded and 5 ml of a second solution added for 15 min containing 6 M urea, 30% (v/v) glycerol, 2% (w/v) SDS, 2.5% (w/v) iodoacetamide, and 0.01% (w/v) bromophenol blue (SH alkylation step). The SDS-PAGE was performed on an EttanDalt™ (GE) electrophoresis unit. The strips were placed in a 1.5-mm thick 12.5% (w/v) polyacrylamide gel sealed with 0.1% (w/v) agarose in SDS-electrophoresis buffer (25 mM Tris, 192 mM glycine, 1% [w/v] SDS) containing 0.01% (w/v) bromophenol blue. The electrophoresis was run for 30 min at 4 watts per gel followed by a further run at 20 watts per gel until the bromophenol blue band reached the bottom of the gel. The gels were stained overnight with a solution of Coomassie Brilliant Blue (CBB: 5% [v/v] methanol, 42.5% [v/v] ethanol, 10% [v/v] acetic acid, 0.2% [w/v] CBB G250, 0.05% [w/v] CBB R250) and destained in the solution 5% (v/v) methanol, 42.5% (v/v) ethanol, 10% (v/v) acetic acid for 1 h then further in 7% (v/v) acetic acid until the background was sufficiently reduced. The gels were finally documented on an Imagescanner (GE).

Image Analysis

The gels were scanned through the use of a transmission Image Scanner (GE). Image analysis was performed with Image Master 2DTM Platinum (version 6.0) (GE) for at least three gels for each pH condition. Spot intensities were quantified according to the volume percent; and specific spots classified as up or down regulated whenever the average volume for a given spot at one pH condition was at least 50% different from the corresponding volume of the homologous spot under another pH condition, and presented simultaneously a p-value lower than 0.05 in a t-test.

In-gel tryptic digest of proteins and mass spectrometry

Protein spots were excised from the 2-D gels and placed into microtiter-plate wells that had been previously washed twice with a trifluoracetic acid: acetonitrile: water solution (0.1:60:40 [v/v]). The tryptic digestions were performed as described on the Keck home page (http://info.med.yale.edu/wmkeck/prochem/geldig3.htm). Samples containing digested proteins were mixed 1:1 with a solution of water:acetonitrile:trifluoracetic acid (67:33:0.1 [v/v]) saturated with α-cyano-4-hydroxycinnamic acid. The mass spectra were obtained on an Ultraflex MALDI-TOF-MS™ instrument (Brucker, Bremen). The identification of proteins from each peptide-mass fingerprint was carried out by means of the Mascot software (Matrix Sciences, London, UK) with the aid of a S. meliloti database. The following search parameters were set in the Mascot software: enzyme, trypsin; missed cleavages, 1; peptide tolerance, 150 ppm; MH+; and monoisotopic.

Transcriptomic analysis in microarrays in search of S. meliloti genes differentially expressed under acid and/or neutral extracellular conditions

For the microarray slides, RNA was isolated from the bacterial cells grown in the continuous culture stabilized at pH 7.0 or pH 6.1. The rhizobial cells were disrupted mechanically as described previously95. Cy3- and Cy5-labelled cDNAs were prepared by the method of De Risi et al. from 10 μg of total RNA96. Three slide hybridizations were performed with the labelled cDNA synthesized from both of the RNA preparations (i.e., pH 7.0, pH 6.1, 3 replicates each, technical replicates). Hybridization, image acquisition, and data analysis were carried out as previously described9597. The mean-signal and mean-local-background intensities were determined for each spot on the microarray images by means of the ImaGene 5.5 software for spot detection, image segmentation, and signal quantification (Biodiscovery Inc., Los Angeles, USA). M and A values for individual spots in the microarrays were calculated as previously described95. Here a normalization method was used based on a local regression that accounts for the intensity and spatial dependence on the dye biases98. Normalization and t-statistics were computed by means of the EMMA 1.1 microarray data-analysis software developed at the Bioinformatics Resource Facility, Center for Biotechnology, Bielefeld University (http://www.genetik.uni-bielefeld.de/EMMA/; ref. 99. Genes were considered differentially expressed if obtained for at least five of the nine replicate spots and if the confidence indicator (p) was ≤0.05 and the log2 of the expression ratio (M) ≥ 1 or ≤−1 (i.e., if at least a twofold difference was obtained between the two experimental conditions). The microarray results were verified for a specific acid-induced transcript (degP1 - SMc02365, MpH 6.1/pH 7.0 = 1.98; Supplementary Fig. 6) by real-time quantitative reverse-transcription PCR using a KAPA SYBR FAST One-Step qRT-PCR kit (Kapa Biosystems Inc. Wilmington, MA, USA) according to the manufacturer’s instructions. The arrays data have been deposited in the NCBI Gene Expression Omnibus database100 and are accessible through GEO Series accession number GSE74449.

Metabolite extraction and derivatization for the metabolome experiments

The extraction and sample derivatization were performed as described previously101. One ml of 80% (v/v) aqueous methanol containing 10 μM ribitol (internal standard) was added to 10–30 mg dry weight of S. meliloti cells in 1.5-ml screw-cap tubes with 0.5 g acid-washed glass beads (Sigma-Aldrich). Immediately after adding the methanol, the cells were disrupted three times in a FastPrep™ instrument (Qbiogene, Heidelberg, Germany) at 6.5 m/s for 45 s. For improved extraction, the mixtures were further incubated at 70 °C for 15 min at 1,400 rpm in a Thermomixer™ (Eppendorf, Hamburg, Germany). After centrifugation at 18,500 × g for 20 min at room temperature the clear supernatant was evaporated to dryness in a nitrogen stream. Methoximation of carbonyl moieties with 50 μl of a 20 mg.ml−1 solution of methoxylamine hydrochloride in pyridine was run at 37 °C for 90 min with constant stirring. Dissociable protons were protected with trimethylsilyl groups through reaction with 50 μl N-methyl-N-[trimethylsilyl]trifluoroacetamide at 37 °C for 30 min.

Gas chromatography–mass-spectrometry (GC-MS) analysis was performed with a TraceGC™ gas chromatograph and a PolarisQ™ ion-trap mass spectrometer equipped with an AS2000™ auto sampler (Thermo Finnigan, Dreieich, Germany). Sample volumes of 1 μl were injected at a 250 °C injector temperature without a split. The gas chromatograph contained a 30-m X 0.25-mm Equity-5™ column having a 0.25-μm 5% diphenyl- 95% dimethylsiloxane coating (Supelco, Bellfonte, Calif.). The interface temperature was adjusted to 250 °C and the ion source set at 200 °C. The helium carrier gas was set at a constant flow of 1 ml.min−1. After a 3 min constant heating at 80 °C the oven temperature was raised in steps of 3 °C.min−1 to 300 °C. In order to equilibrate the system for the next injection, the temperature was set at 80 °C for 5 min. Mass spectra were recorded at 2 scans.s−1 with a scanning range of 50–550 m/z. Evaluation of the chromatograms was performed with the Xcalibur 2.0 software (Thermo Finnigan). Metabolites were identified by comparison with the NIST 98 (NIST, Gaithersburg, MD) database and with purified standards. Automatic peak quantification of selected metabolites was implemented into the processing setup of the Xcalibur software. Virtual fragmentations were performed with the Mass Frontier 2.0 software (Thermo). The data were log-normalized and uni- and multivariate tests applied through the Metaboanalyst server102.

Additional Information

How to cite this article: Draghi, W. O. et al. A consolidated analysis of the physiologic and molecular responses induced under acid stress in the legume-symbiont model-soil bacterium Sinorhizobium meliloti. Sci. Rep. 6, 29278; doi: 10.1038/srep29278 (2016).