Skip to main content

Rapid acid treatment of Escherichia coli: transcriptomic response and recovery



Many E. coli genes show pH-dependent expression during logarithmic growth in acid (pH 5–6) or in base (pH 8–9). The effect of rapid pH change, however, has rarely been tested. Rapid acid treatment could distinguish between genes responding to external pH, and genes responding to cytoplasmic acidification, which occurs transiently following rapid external acidification. It could reveal previously unknown acid-stress genes whose effects are transient, as well as show which acid-stress genes have a delayed response.


Microarray hybridization was employed to observe the global gene expression of E. coli K-12 W3110 following rapid acidification of the external medium, from pH 7.6 to pH 5.5. Fluorimetric observation of pH-dependent tetR-YFP showed that rapid external acidification led to a half-unit drop in cytoplasmic pH (from pH 7.6 to pH 6.4) which began to recover within 20 s. Following acid treatment, 630 genes were up-regulated and 586 genes were down-regulated. Up-regulated genes included amino-acid decarboxylases (cadA, adiY, gadA), succinate dehydrogenase (sdhABCD), biofilm-associated genes (bdm, gatAB, and ymgABC), and the Gad, Fur and Rcs regulons. Genes with response patterns consistent with cytoplasmic acid stress were revealed by addition of benzoate, a membrane-permeant acid that permanently depresses cytoplasmic pH without affecting external pH. Several genes (yagU, ygiN, yjeI, and yneI) were up-regulated specifically by external acidification, while other genes (fimB, ygaC, yhcN, yhjX, ymgABC, yodA) presented a benzoate response consistent with cytoplasmic pH stress. Other genes (the nuo operon for NADH dehydrogenase I, and the HslUV protease) showed delayed up-regulation by acid, with expression rising by 10 min following the acid shift.


Transcriptomic profiling of E. coli K-12 distinguished three different classes of change in gene expression following rapid acid treatment: up-regulation with or without recovery, and delayed response to acid. For eight genes showing acid response and recovery (fimB, ygaC, yhcN, yhjX, ymgABC, yodA), responses to the permeant acid benzoate revealed expression patterns consistent with sensing of cytoplasmic pH. The delayed acid response of nuo genes shows that NADH dehydrogenase I is probably induced as a secondary result of acid-associated metabolism, not as a direct response to cytoplasmic acidification.


Extreme-acid survival is an important virulence factor for human pathogens such as E. coli strain O157:H7 [1]. In order to colonize the gastrointestinal tract, Escherichia coli and other enteric bacteria must be able to grow in environments at extreme pH such as the duodenum (pH 9–10) and the stomach (pH 2–4) [24]. In the gastrointestinal tract, enteric bacteria are subjected to acid stress from strong acid (HCl) as well as bacterial fermentation products such as acetic, propionic, and butyric acids, which are membrane-permeant weak acids [5]. E. coli preserves the integrity of proteins and nucleic acids present in the cytoplasm by maintaining cytoplasmic pH at approximately pH 7.6, over a wide range of external pH [6, 7]. Low pH up-regulates genes required for survival under more extreme acid conditions, such as the arginine- and glutamate-dependent acid resistance systems [810].

In E. coli K-12, many studies show that adaptation to low or high pH stress involves regulation of gene expression and protein synthesis, as well as post-translational and regulation of protein function [1114]. For example, acid induces the lysine decarboxylase operon cadBA under control of the signal-transducing regulator CadC [13, 14]. Lysine decarboxylase (CadA) removes CO2 from lysine, releasing cadaverine, a base that counteracts acidity. Other kinds of genes up-regulated by acid include periplasmic chaperones, inner-membrane and outer-membrane proteins, acetate-stress proteins, and systems utilizing the proton gradient, such as motility and chemotaxis [11]. Base-induced proteins include fermentation pathways generating acidic products, and metabolic complexes that import protons or minimize proton export, such as the F1Fo ATPase and cytochrome d oxidoreductase (cydAB) [11, 12].

While gene expression as a function of steady-state external pH has received extensive study, less is known about the dynamic response of E. coli to a sudden pH change. Rapid acidification (within 10 s) of the external medium causes the cytoplasmic pH to fall, followed by recovery within minutes to close to the original value [15, 16]. For the cadBA operon, expression is upregulated within 4 min of rapid acidification, and downregulated within 4 min after pH neutralization [17]. But other genes showing pH-dependent expression ratios under growth for several generations at different pH values may be dependent on secondary effects of growth with pH stress. The gene products required to recover from rapid change in external pH (with a transient failure of pH homeostasis) may differ substantially from those elevated or repressed during steady-state growth at low or high pH, where pH homeostasis is maintained. Furthermore, few studies distinguish the effects of cytoplasmic versus external pH on gene expression; for example, yfiD expression responds to membrane-permeant weak acids that depress cytoplasmic pH without affecting external pH [18].

In this study we used microarray hybridization to investigate the transcriptomic response to a sudden acid shift. This approach enabled us to distinguish between genes up-regulated immediately upon acid shift, versus those whose expression increases after some time in the presence of acid. Genes of interest were followed up with quantitative real-time PCR (qPCR) using a membrane-permeant weak acid to assess whether the effect of acid treatment could be associated with the decrease of cytoplasmic pH.

Results and Discussion

Cytoplasmic pH change

The effect of external acidification on the E. coli cytoplasmic pH under our experimental conditions was tested, in order to assess the state of the cytoplasm during our transcriptional profiling. The probe used was pH-dependent fluorescence of tetR-YFP [16].

The external pH was acidified by adding HCl (8.5 mM) to cell cultures of E. coli W3110 containing a tetR-YFP plasmid, suspended in media adjusted to pH 7.6. Upon HCl addition, the cytoplasmic pH fell within 5 s to approximately pH 6.4 (Fig. 1A). Recovery began approximately 4 s after the cytoplasmic pH reached the lowest point (~pH 6.4). After 10 min, the cytoplasmic pH ranged between pH 7.5 and pH 7.7. Each replicate culture exhibited a biphasic recovery, which is consistent with the recovery pattern previously described by ref [16]. As a control, KCl was added instead of HCl (8.5 mM KCl, pH 7.6) in order to detect changes in fluorescence signal that were independent of pH change. Upon KCl addition, there was no significant change in fluorescence signal.

Figure 1
figure 1

Effects of external acid shift and benzoate addition on cytoplasmic pH. (A) Cultures were suspended at pH 7.6 in supplemented M63 minimal media (20 mM HOMOPIPES). At time zero, 8.5 mM HCl was added to shift the external pH from pH 7.6 to 5.5. For the control, KCl (8.5 mM at pH 7.6) was added instead of HCl (B) The permeant weak acid sodium benzoate (30 mM) was added at time zero to cultures suspended at pH 7.0 in supplemented M63 minimal meida (50 mM HOMOPIPES). As a control, HOMOPIPES (pH 50 mM at 7.0) was added instead of sodium benzoate. For both experimental conditions, fluorescence was converted to pH units using the standard curves as described in Materials and Methods. Under each condition, three independent cultures were tested.

For comparison, the effect of a permeant-acid stress was measured under our experimental conditions [16]. Addition of a permeant acid such as benzoate at high concentration depresses cytoplasmic pH with little or no recovery, without affecting external pH [15]. Benzoate (30 mM) was added to cell cultures of W3110 containing p-tetR-YFP, which were suspended in media pH adjusted to pH 7.0. After addition of 30 mM benzoate, the cytoplasmic pH fell within 8 s to that of the external medium (pH 7.0, Fig. 1B). Over the period of 5 min during which spectra were recorded, the cytoplasmic pH made a minimal recovery. (A period of only 5 min is shown because the YFP protein begins to destabilize in the presence of benzoate.) Buffer (50 mM HOMOPIPES, pH 7.0) was added as a control, and did not cause a significant change in fluorescence signal.

RNA hybridization following rapid acidification

For acid treatment and RNA hybridization, E. coli strain W3110 was cultured to early log phase (OD600 = 0.2) at pH 7.6 (+/- 0.2 pH units). pHext 7.6 is equivalent to the cytoplasmic pH, and thus cells suspended at pHext 7.6 exhibit minimal pH regulated gene expression relative to rapidly acidified cultures. At time zero (OD600 0.2), HCl was added rapidly, resulting in the acidification of the external medium to pH 5.5. At time points 0, 1, 5, and 10 min, 10 mL of culture were removed. The pH of the media remained at pH 5.5 (+/- 0.2 pH units) after the 10 min sample was taken.

At time points 0, 1, 5, and 10 min, RNA was purified and cDNA generated as described under Methods. The cDNA from five independent cultures of each time-point were hybridized to Affymetrix antisense E. coli arrays. Model-based expression indices were calculated from the probe-level data from Affymetrix DAT files using dChip software [19, 12]. Array data were deposited at the NCBI Gene Expression Omnibus (accession GSE4778).

Time-dependent variation

To determine the effect of time on differential gene expression, the log2 transforms of normalized model-based expression values from the microarray probe level data were analyzed during the 10-min interval. A mixed-effects linear model for longitudinal data was applied using the software SAS Mixed Procedure (Proc Mixed). For each gene, the mixed-effects model was applied at a significance level of 0.001 to determine whether or not the gene expression profile deviates significantly from a flat line over time, which indicates no time effect. The model takes into account time (fixed effect), the random deviation among biological replicates (random effects), and random error to determine significant time-dependence for gene variation [2022].

For array hybridization, our experimental design and analysis were consistent with the "consensus" recommendations of Allison [23] in that we included an ample number of biological replicates (five independent cultures for each time point), assuring a low false-positive rate (0.001, or one false-positive in 1,000 genes) as well as a high power of detection. The power of detection of changes in gene expression over time was estimated based on the observed biological variation among flask cultures (90% of the genes showed a standard deviation of less then 0.08) and the technical error term (90% of the genes showed a standard deviation of less than 0.28). The power of detection was estimated for two cases of expression change using a Monte Carlo simulation based on 5,000 trials. For a gene whose expression increases two-fold over 10 min, the power of detection was 96%. For a gene whose expression increases two-fold by 5 min, with partial recovery by 10 min, the power of detection was 97%.

Of the 4,377 genes on the array, 1,582 (36% of the E. coli genome) showed significant time-dependent variation in expression among the four time points. Their expression indices and log2 expression ratios relative to time zero are tabulated in Additional File 1. The percent recovery of expression level was calculated as follows from the values of the log2 expression indices [log2(E)] for each gene:

Percent recovery = | ( p e a k log 2 ( E ) t i m e 10 log 2 ( E ) ) ( p e a k log 2 ( E ) t i m e 0 log 2 ( E ) ) | × 100 MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGacaGaaiaabeqaaeqabiWaaaGcbaGaeeiuaaLaeeyzauMaeeOCaiNaee4yamMaeeyzauMaeeOBa4MaeeiDaqNaeeiiaaIaeeOCaiNaeeyzauMaee4yamMaee4Ba8MaeeODayNaeeyzauMaeeOCaiNaeeyEaKNaeyypa0ZaaqWaaKqbagaadaWcaaqaamaabmaabaGaemiCaaNaemyzauMaemyyaeMaem4AaSMagiiBaWMaei4Ba8Maei4zaCMaeGOmaiJaeiikaGIaeeyrauKaeiykaKIaeyOeI0IaemiDaqNaemyAaKMaemyBa0MaemyzauMaeeiiaaIaeGymaeJaeGimaaJagiiBaWMaei4Ba8Maei4zaCMaeGOmaiJaeiikaGIaeeyrauKaeiykaKcacaGLOaGaayzkaaaabaWaaeWaaeaacqWGWbaCcqWGLbqzcqWGHbqycqWGRbWAcyGGSbaBcqGGVbWBcqGGNbWzcqaIYaGmcqGGOaakcqqGfbqrcqGGPaqkcqGHsislcqWG0baDcqWGPbqAcqWGTbqBcqWGLbqzcqqGGaaicqaIWaamcyGGSbaBcqGGVbWBcqGGNbWzcqaIYaGmcqGGOaakcqqGfbqrcqGGPaqkaiaawIcacaGLPaaaaaaakiaawEa7caGLiWoacqGHxdaTcqaIXaqmcqaIWaamcqaIWaamaaa@8966@

where "peak log2(E)" refers to the highest mean log2(E) of an up-regulated gene (at 1 min or at 5 min), or the lowest mean log2(E) of a down-regulated gene. For example, a gene whose expression increase peaks at 100-fold, then declines to 10-fold, would show 50% recovery with a remaining elevation of 10-fold; thus it would be "half recovered" in that the peak/recovered expression ratio equals the recovered/original expression ratio.

Genes that showed a significant increase or decrease in mean expression at 1 min and/or 5 min, and exhibited less than a 200% recovery (that is, over 100% past the baseline expression level), were defined as up-regulated or down-regulated, respectively (Additional File 1). Following the rapid acid shift, 630 genes (14% of the genome) were up-regulated and 586 genes (13%) were down-regulated. Their expression ratios relative to time zero, and their percent recovery, are tabulated. Of the genes defined as up-regulated or down-regulated, those whose expression recovered fully by 10 min (a percent recovery of between 100%–200%) were designated "transient." Thus defined, 52 genes (1.2% of the genome) showed transient up-regulation, and 45 genes (1%) showed transient down-regulation.

Genes that showed up-regulation or down-regulation by 5 min, and reached expression ratios of 2-fold or greater, are tabulated in Table 1 and Additional file 2, respectively. Many of these also show pH-dependent expression under steady-state conditions, in a study by Maurer et al in which the bacteria were cultured for several generations at pH 5.0, pH 7.0, and pH 8.7 [11]. The percentage of genes with elevated steady-state expression in acid culture [11] was 19%, for genes up-regulated at 1 min; 24%, at 5 min; and 34%, at 10 min (Table 1). For genes down-regulated after acid shift, the percentage down-regulated in steady-state acid culture [11] was 17% at 1 min, 17% at 5 min, and 27% at 10 min (Additional file 2). Thus, the longer a gene showed up-regulation (or down-regulation) in our study, the more likely it was to show elevated (or depressed) expression under steady-state growth at low pH.

Table 1 Genes upregulated at least 2-fold after acid shift.a

Genes whose up-regulation (or down-regulation) was transient would not be expected to show steady-state expression ratios. In fact, of the 52 genes showing transient up-regulation to levels 2-fold or greater, none are reported in Ref. [11]; and of the 45 genes showing transient down-regulation, only one is reported. Thus, our study revealed transiently acid-responsive genes that were not observed under steady-state growth conditions.

Up-regulated genes following an acid shift

yagU and other acid stress genes of unknown function

The gene encoding inner-membrane protein YagU [24] contributes to acid resistance [12] and is one of the genes most strongly up-regulated (16-fold) at low pH [11, 12]. Following acid shift, yagU was one of the top genes up-regulated at 1 min, with peak expression at 5 min and a 10% recovery at 10 min (Table 1).

The up-regulation of yagU was confirmed by real-time PCR (Fig. 2). The time course of yagU expression relative to time 0 was similar to that of the microarray data (Table 1). In order to determine whether up-regulation was related to external pH or cytoplasmic pH, a second time course was conducted using the permeant weak acid benzoate. Addition of benzoate to cultures suspended at pHext 7.0 acidifies the cytoplasmic pH while the external pH remains constant. The yagU gene showed no significant up-regulation following benzoate addition. Therefore, yagU response is more consistent with direct sensing of external pH, not cytoplasmic pH.

Figure 2
figure 2

Real-Time PCR expression ratios of selected genes. Bacteria were cultured and real-time PCR was conducted as described in Materials and Methods. For each experimental condition, the 1, 5, and 10 min expression values were normalized to time 0. For each gene, log2 (expression ratios) are shown for amplified message after the rapid addition of either HCl (15 mM, solid line) or sodium benzoate (30 mM, dotted line).

Three other genes observed by real-time PCR (ygiN, yjeI, and yneI) showed approximately 16-fold increase in expression after rapid acid shift, but no effect of benzoate addition (Fig. 2). Expression of yjeI recovered after acid shift, but the benzoate experiment ruled out a direct response to cytoplasmic pH. Eight other genes tested (fimB, ygaC, yhcN, yhjX, ymgA,ymgB,ymgC, and yodA) increased expression more than 4-fold with full or partial recovery, and were up-regulated after benzoate addition. This pattern is consistent with response to intracellular pH change. The gene fimB encodes the flagellar phase switching mediator FimB [25]; the others have unknown functions. None of these genes were previously reported to be up-regulated by acid.

Acid-inducible amino-acid decarboxylases

The most well-documented class of response to external acid is the degradative amino-acid decarboxylases, which convert specific amino acids to amines while exporting a proton [14, 26]. All of these were up-regulated within 5 min following the acid shift, with little or no recovery. These include the lysine decarboxylase and transporter cadBA [27] and the arginine decarboxylase adiY [28] (see Table 1 and Additional File 1) as well as the glutamate decarboxylase and transporter (gadA and gadBC) which are part of the gad acid response regulon [9, 29]. Other members of the gad regulon up-regulated after the acid shift included gadE, gadWX, and evgAS.

Metabolism and transport

A number of metabolic enzymes and transporters were up-regulated following acid shift. Several of these are associated with biofilm formation, such as isocitrate lyase (aceA) and malate synthase (aceB) [30]. The proteins of galactitol transport were upregulated strongly (gatA, B, Y, Z). The gat genes are involved in biofilm formation [30] and they show elevated expression in steady-state acid, which appears to favor substrates that produce fewer acidic fermentation products [11, 12]. Biofilm development genes ymgABC were also up-regulated after acid shift, as was the biofilm-dependent modulation protein bdm [31]. The expression levels of ymgABC began recovery by 10 min. These genes have not been reported previously to show elevated expression under steady-state acid conditions.

Other metabolic genes known for acid-enhanced expression [11, 12] were found to be up-regulated without recovery following acid shift. These include sdhABCD, encoding succinate dehydrogenase, [32]; the dppA, B, C, D, F dipeptide transport system [33]; and glycerol transport components glpD, F, Q, T, of which glpD encodes a multidrug-resistance component associated with persister cell formation [34].

Osmotic shock response genes

Several genes up-regulated after acid shift are associated with osmotic shock, in particular the Rcs regulon [31, 35]. In the Rcs regulon, the transcriptional regulators RcsA and RcsB control expression of genes involved in colanic capsule biosynthesis. The Rcs regulon members up-regulated after acid shift included rcsABCD, bdm, and osmB [36]. Of these, only rcsA shows elevated expression in steady-state acid [11]. An osmotic shock gene outside the Rcs regulon, osmY [37], also was up-regulated after acid shift.

The Fur regulon

The ferric uptake regulator (Fur) governs the expression of various acid tolerance genes [38]. Of the genes that were strongly-up-regulated after rapid acidification, eight are genes previously shown to be controlled by the Fur regulon [39, 40, 32, 41, 42], Genes whose expression peaked at 5 min encode the cadmium induced metal binding protein (yodA) [39], ygaC [40], and the oxidative stress genes nrdH and nrdI [42]. The percent recovery of these genes ranged between 46 to 80% (Table 1). Fur also regulates the succinate dehydrogenase operon sdhABCD, which is up-regulated by acid under steady-state conditions [11] and was up-regulated after rapid acidification (Table 1).

Many genes of unknown or putative function exhibited expression profiles that paralleled the shift and recovery of cytoplasmic pH. Such genes include ygiN, whose gene product has recently been isolated and characterized as a quinol monooxygenase [43]. The gene ygiN exhibited peak up-regulation at 5 min and made a 63% recovery by 10 min (Table 1). Genes encoding the putative NAD+- dependent succinate semialdehyde dehydrogenase YneI [44] and the putative periplasmic protein YhcN, exhibited peak expression at 5 min and recovered by 57% and 23% respectively (Table 1).

Down-regulated genes following an acid shift

The most strongly down-regulated genes following acid shift are shown in Additional file 2. Many of these show steady-state elevated expression at high pH.

Genes encoding the high-affinity cytochrome oxidoreductase d (cyd operon) are favored under low-oxygen conditions as well as at high pH [11]. The genes cydA, B, C were down-regulated following acid shift. Similarly, the periplasmic extracytoplasmic response regulator CpxP [45, 46] is strongly up-regulated in base [11]. The cpxP gene was down-regulated immediately after acid shift, with little recovery (Additional file 2). The gene encoding the base-inducible putative membrane-bound redox modulator Alx [11, 47, 48] was also down-regulated, as were the genes encoding serine deaminase and transport, sdaABC (Additional Files 1 and 2). Amino-acid deaminases appear to be induced as a response to reverse pH increase [11, 18, 49].

Many genes involved in polyamine metabolism and transport show increased expression at high pH [11, 50]. Several of these were down-regulated following acid shift, including the spermidine metabolism genes speA, B, D, E [51] and polyamine transport genes potA, B, C, D [52].

Delayed response to acid shift

Some genes showed little or no response to acid shift until the 10-min time point (included in Additional File 1). Most striking was the delayed up-regulation of the nuo operon encoding NADH dehydrogenase I [53, 54]. Thirteen genes of the nuo operon were all up-regulated about 50% at 10 min, after little or no significant differences at 1 min or 5 min. The NADH dehydrogenase I is known to show elevated expression in steady-state acid [11] and under anoxic conditions; its expression declines with aeration [55]. Our new data showed that nuo expression does not increase immediately with acidity, but as a delayed consequence of the acid shift. The previously reported rise in NADH dehydrogenase could be a consequence of deenergization following acid stress. Another complex showing delayed up-regulation by acid was the HslVU protease [56]. The hsl genes are known to be up-regulated by acid under steady-state conditions [11]. Our new data is more consistent with a secondary response to acidification, such as a need to degrade acid-denatured proteins.


When cultures of E. coli K-12 W3110 were subjected to a rapid decrease in external pH, from pH 7.6 to pH 5.5, the cytoplasmic pH decreased but began to recover within less than a minute after the acid shift. Up-regulated genes included amino-acid decarboxylases (cadA, adiY, gadA), succinate dehydrogenase (sdhABCD), biofilm-associated genes (bdm, gatAB, and ymgABC), and components of the Gad, Fur and Rcs regulons. Many of the acid-up-regulated genes showed sustained elevation, such as cadA and sdhCD. On the other hand, the rapid acid exposure experiment revealed genes whose response to acid shift is transient, such as the biofilm development genes ymgABC and the flagellar phase variation modulator fimB. These genes were up-regulated after benzoate addition, a condition in which cytoplasmic pH is depressed permanently without change in external pH. Thus, these genes are candidates for direct response to cytoplasmic acid stress. The acid response of other genes, such as nuo and hsl, was shown to be delayed. Their delayed response suggests a secondary effect of growth in acid. Overall, our study of gene expression following rapid acid shift or rapid benzoate addition revealed several classes of acid stress response: genes showing sustained up-regulation (or down-regulation) following acid shift; genes showing transient up-regulation after acid shift, as well as up-regulation following permeant-acid treatment of the cytoplasm; and genes showing delayed up-regulation after acid shift.


Growth conditions

Rapid external acid stress

Escherichia coli K-12 strain W3110 [57] was cultured overnight in buffered LBK medium [11]. Overnight cultures were diluted 550-fold into 55 mLs of buffered LBK medium and aerated at 225 rpm in an orbital water bath. The medium was buffered with 20 mM homopiperazine-N, N'-bis-2-(ethane-sulfonic acid) (HOMOPIPES) (pKa 4.55), and was adjusted using KOH to pH 7.6. Cultures were grown at 37°C to an OD600 of 0.2 (approximately five generations). At OD600 = 0.2, the pH of the medium was lowered to pH 5.5 with the addition of 1 M HCl (840 μL). One 10-mL sample was taken from each of the five biological replicates immediately before the addition of HCl (0 min). After HCl addition, three 10-mL samples were taken from each of the five biological replicates, corresponding to 1, 5, and 10 min. The pH of the cultures was tested at the end of the acid shift to ensure they were maintained within 0.2 pH units of pH 5.5.

Permeant acid stress

The permeant-acid time course was performed in order to generate a permanent intracellular acidification, thus revealing genes responsive to cytoplasmic pH stress. The growth conditions were the same as those used for HCl acidification, with the exception of the pH and buffering capacity of the medium. The medium was adjusted to pH 7.0 with KOH solution and was buffered with 50 mM HOMOPIPES (pKa 4.55). The external pH was set at 7.0 in order to ensure that the addition of weak acid led to a permanent decrease in cytoplasmic pH (from a cytoplasmic pH of approximately 7.6). Sodium benzoate (pH 7.2) was added once cultures reached an OD600 of 0.2, for a final concentration of 30 mM. As described above, samples were taken at time 0, 1, 5, and 10 min. The pH of the media remained at pH 7.0 ± 0.2. RNA was extracted to be analyzed by real-time PCR.

Cytoplasmic pH Measurement

Cytoplasmic pH was measured based on YFP fluorescence as reported [16]. Measurements were made on strain JLS0617, which consists of strain W3110 transformed with pSL38-YFP containing tetR-YFP [58]. Strain JLS0617 was cultured overnight in buffered M63 salts medium [20 mM homopiperazine-N,N-bis-2-(ethanesulfonic acid) (HOMOPIPES), pH 7.6] containing 1.5% casein hydrolysate, 0.8% glycerol, and 100 μg/ml ampicillin. The overnight culture was diluted 25-fold in the same medium and cultured to late log phase (OD600 = 0.8 to 0.9) at 37°C with aeration. The cultures were suspended to an OD600 of 0.4 in M63 medium containing 0.2% casein hydrolysate, 0.8% glycerol, and 20 mM HOMOPIPES (pH 7.5). The suspension media used for the permeant acid experiment contained 50 mM HOMOPIPES, in order to ensure that the external pH was maintained at pH 7.0.

Excitation spectra were recorded using a Fluoromax-3 spectrofluorimeter (Horiba Jobin Yvon). A cell suspension (3 ml) was placed into a Starna Spectrosil quartz cuvette with a path length of 10 mm. The temperature of the chamber was adjusted to 30°C, and aeration was provided by stirring. Spectra were recorded at 30°C because the YFP probe is temperature sensitive and exhibits the highest fluorescence signal at this temperature. YFP excitation was measured from 490 to 520 nm (slit width, 2 nm), using an emission wavelength of 550 nm (slit width, 20 nm). Continuous excitation spectra were obtained every 4 s for 1 min before aliquot addition (time zero) and then for 10 min after addition. In order to control for effects on fluorescence unrelated to pH change, KCl (8.5 mM, pH 7.6) was added to cell suspensions at time 0. For each condition, spectra were recorded for three biological replicates.

Data were fitted to a standard curve correlating cytoplasmic pH with fluorescence intensity. A linear interpolation was conducted between pH 7.6 and pH 5.5 signal intensities. The fluorescence corresponding to pH 5.5 was determined by adding an additional 20 mM benzoate to each test sample at the end of the time course (thus collapsing the cytoplasmic pH to pH 5.5.)

The pH curve was modified for experiments in which a permeant acid was added to cultures in medium at external pH 7.0. In this case, linear interpolation was conducted between pH 7.6 and pH 7.0. The fluorescence corresponding to pH 7.0 was determined by adding an additional 20 mM benzoate (50 mM benzoate was the final concentration) to each test sample at the end of the time course.

RNA isolation

Samples were transferred from growth culture to vials containing 1 mL ice-cold 10% phenol-ethanol stop solution [59] in order to stabilize the bacterial RNA. RNA was isolated as described previously [11, 12] using the RNeasy Kit with on-column DNA digestion (Qiagen). The quality of the isolated RNA was examined using the Agilent Bioanalyzer 2100.

cDNA preparation and array hybridization

Standard methods were used for cDNA synthesis, fragmentation, and end-terminus biotin labeling [11, 12]. Labeled cDNA samples were hybridized to Affymetrix GeneChip E. coli Antisense Genome Arrays. Hybridized arrays were stained with streptavidin-phycoerythrin using the Affymetrix Fluidic Station 450. After staining, arrays were scanned with a GC3000 scanner.

Analysis of Gene Expression

Log2 transformed model-based expression values were calculated from the probe level data using D-Chip software [19]. In order to determine what genes show significant time dependent variation in expression following the acid shift over 1, 5, and 10 min, a mixed effects linear model for longitudinal data was performed using Statistical Analysis System (SAS) [2022]. A p-value < 0.001 was considered a significant time dependence, which corresponds to a false positive rate of 1 in 1000 genes.

Log2 ratios for expression indices at times 1, 5, and 10 minutes post HCl addition were determined for each gene by subtracting the time 0 log2(E) from the time 1, 5, and 10 minute log2(E) respectively.

Real-Time Polymerase Chain Reaction (PCR)

Expression of mRNA for individual genes was quantified by real-time PCR using an ABI Prism 7500 DNA analyzer (Applied Biosystems) as performed previously [12]. Primer Express Software v2.0 (Applied Biosystems) was used for primer design. The primers chosen had minimal GC content and amplified 50–70 bp segments of the target genes (Table 2). The SYBR Green PCR One-Step RT-PCR protocol (Applied Biosystems) was used, in which cDNA reverse transcription and PCR amplification occur in the same well. Nucleic acid concentrations were: 0.1 nM forward primer, 0.1 nM reverse primer, and 50 ng target RNA. PCR cycling conditions were: reverse transcription at 48°C for 30 min, 95°C for 10 min, 40 cycles of denaturation at 92°C for 15s, and extension at 60°C for 1 min. For detection of primer dimerization or other artifacts of amplification, a dissociation curve was run immediately after completion of the real-time PCR. Individual gene expression profiles were normalized based on measurement of the original RNA sample amplified. All expression levels are presented relative to the expression at time 0.

Table 2 Primers for real-time PCR


  1. Hayashi T, Makino K, Ohnishi M, Kurokawa K, Ishii K, Yokoyama K, Han CG, Ohtsubo E, Nakayama K, Murata T, Tanaka M, Tobe T, Iida T, Takami H, Honda T, Sasakawa C, Ogasawara N, Yasunaga T, Kuhara S, Shiba T, Hattori M, Shinagawa H: Complete genome sequence of enterohemorrhagic Escherichia coli O157:H7 and genomic comparison with a laboratory strain K-12. DNA Res. 2001, 8: 11-22. 10.1093/dnares/8.1.11.

    Article  CAS  PubMed  Google Scholar 

  2. de Jonge R, Takumi K, Ritmeester WS, Leusden FM: The adaptive response of Escherichia coli O157 in an environment with changing pH. J Appl Microbiol. 2003, 94: 555-560. 10.1046/j.1365-2672.2003.01865.x.

    Article  CAS  PubMed  Google Scholar 

  3. Russell TL, Berardi RR, Barnett JL, Dermentzoglou LC, Jarvenpaa KM, Schmaltz SP, Dressman JB: Upper gastrointestinal pH in seventy-nine healthy, elderly, North American men and women. Pharm Res. 1993, 10: 187-196. 10.1023/A:1018970323716.

    Article  CAS  PubMed  Google Scholar 

  4. Dressman JB, Berardi RR, Dermentzoglou LC, Russel TL, Schmaltz SP, Barnett JL, Jarvenpaa KM: Upper gastrointestinal (GI) pH in young, healthy men and women. Pharm Res. 1990, 7: 756-761. 10.1023/A:1015827908309.

    Article  CAS  PubMed  Google Scholar 

  5. Lin J, Smith MP, Chapin KC, Baik HS, Bennett GN, Foster JW: Mechanisms of acid resistance in enterohemorrhagic Escherichia coli . Appl Environ Microbiol. 1996, 62: 3094-3100.

    PubMed Central  CAS  PubMed  Google Scholar 

  6. USA, Slonczewski JL, Rosen BP, Alger JR, Macnab RM: pH homeostasis in Escherichia coli: measurement by 31P nuclear magnetic resonance of methylphosphonate and phosphate. Proc Natl Acad Sci. 1981, USA , 78: 6271-6275. 10.1073/pnas.78.10.6271.

    Google Scholar 

  7. Zilberstein D, Agmon V, Shuldiner S, Padan E: Escherichia coli intracellular pH, membrane potential, and cell growth. J Bacteriol. 1984, 158: 246-252.

    PubMed Central  CAS  PubMed  Google Scholar 

  8. Castanie-Cornet, Marie-Pierre, Penfound TA, Smith D, Elliott JF, Foster JW: Control of acid resistance in Escherichia coli. J Bacteriol. 1999, 181: 3525-3535.

    PubMed Central  CAS  PubMed  Google Scholar 

  9. Ma Z, Masuda N, Foster JW: Characterization of EvgAS-YdeO-GadE branched regulatory circuit governing glutamate-dependent acid resistance in Escherichia coli. J Bacteriol. 2004, 186 (21): 7378-7389. 10.1128/JB.186.21.7378-7389.2004.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  10. Masuda N, Church G: Regulatory network of acid resistance in Escherichia coli. Mol Microbiol. 2003, 48: 699-712. 10.1046/j.1365-2958.2003.03477.x.

    Article  CAS  PubMed  Google Scholar 

  11. 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. J Bacteriol. 2005, 187: 304-319. 10.1128/JB.187.1.304-319.2005.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  12. Hayes ET, Wilks JC, Sanfilippo P, Yohannes E, Tate DP, Jones BD, Radmacher MD, BonDurant SS, Slonczewski JL: Oxygen limitation modulates pH regulation of catabolism and hydrogenases, multidrug transports, and envelope composition in Escherichia coli K-12. BMC Microbiol. 2006, 6: 89-106. 10.1186/1471-2180-6-89.

    Article  PubMed Central  PubMed  Google Scholar 

  13. Auger EA, Redding KE, Plumb T, Childs LC, Meng SY, Bennett GN: Construction of lac fusions to the inducible arginine and lysine decarboxylase genes of Escherichia coli K-12. Mol Microbiol. 1989, 3: 609-620. 10.1111/j.1365-2958.1989.tb00208.x.

    Article  CAS  PubMed  Google Scholar 

  14. Neely MN, Dell CL, Olson ER: Roles of LysP and CadC in mediating the lysine requirement for acid induction of the Escherichia coli cad operon. J Bacteriol. 1994, 176: 3278-3285.

    PubMed Central  CAS  PubMed  Google Scholar 

  15. Slonczewski JL, Macnab RM, Alger JR, Castle AM: Effects of pH and repellent tactic stimuli on protein methylation levels in Escherichia coli. J Bacteriol. 1982, 384-399.

    Google Scholar 

  16. Wilks JC, Slonczewski JL: pH of the cytoplasm and periplasm of Escherichia coli: rapid measurement by green fluorescent protein fluorimetry. J Bacteriol. 2007, 189: 5601-5607. 10.1128/JB.00615-07.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  17. Neely MN, Olson ER: Kinetics of expression of the Escherichia coli cad operon as a function of pH and lysine. J Bacteriol. 1996, 178: 5522-5528.

    PubMed Central  CAS  PubMed  Google Scholar 

  18. Blankenhorn DJ, Phillips J, Slonczewski JL: Acid- and base-induced proteins during aerobic and anaerobic growth of Escherichia coli revealed by two-dimensional gel electrophoresis. J Bacteriol. 1999, 181: 2209-2216.

    PubMed Central  CAS  PubMed  Google Scholar 

  19. Li C, Wong WH: DNA-Chip Analyzer (dChip). In The analysis of gene expression data: methods and software G Parmigiani, E S Garrett, R Irizarry, and S L Zeger Springer-Verlag. 2003

    Google Scholar 

  20. Laird NM, Ware JH: Random-effects models for longitudinal data. Biometrics. 1982, 38: 963-974. 10.2307/2529876.

    Article  CAS  PubMed  Google Scholar 

  21. Ware JH: Linear models for the analysis of longitudinal studies. Am Stat. 1985, 39: 95-101. 10.2307/2682803.

    Google Scholar 

  22. Story JD, Xiao W, Leek JT, Tompkins RG, Davis RW: Significance analysis of time course microarray experiments. PNAS. 2005, 102: 12837-12842. 10.1073/pnas.0504609102.

    Article  Google Scholar 

  23. Allison DB, Cui X, Page GP, Sabripour M: Microarray data analysis: from disarray to consolidation and consensus. Nat Rev Genet. 2006, 7: 55-65. 10.1038/nrg1749.

    Article  CAS  PubMed  Google Scholar 

  24. Rapp M, Drew D, Daley DO, Nilsson J, Carvalho T, Melen K, Gier JWD, Heijne GV: Experimentally based topology models of E. coli inner membrane proteins. Protein Sci. 2004, 13: 937-945. 10.1110/ps.03553804.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  25. Dove SL, Smith SG, Dorman CJ: Control of Escherichia coli type 1 fimbrial gene expression in stationary phase: a negative role for RpoS. Mol Gen Genet. 1997, 254: 13-20. 10.1007/s004380050385.

    Article  CAS  PubMed  Google Scholar 

  26. Slonczewski JL, Foster JW: pH-regulated genes and survival at extreme pH. Escherichia coli and Salmonella: cellular and molecular biology. Edited by: Neidhardt, F C, III, Ingraham JL, Lin ECC, Low KB, Magasanik B, Reznikoff WS, Riley M, Schaechter M, Umbarger HE. 1996, Washington, DC , ASM Press, 1: 1539-1549. 2

    Google Scholar 

  27. Meng S, Bennett GN: Nucleotide sequence of the Escherichia coli cad operon: a system for neutralization of low extracellular pH. J Bacteriol. 1992, 174: 2659-2669.

    PubMed Central  CAS  PubMed  Google Scholar 

  28. Stim-Herndon KP, Flores TM, Bennett GM: Molecular characterization of adiY, a regulatory gene which affects expression of the biodegradative acid-induced arginine decarboxylase gene (adiA) of Escherichia coli. Microbiol. 1996, 142 (5): 1311-1320.

    Article  CAS  Google Scholar 

  29. Ma Z, Richard H, Tucker DL, Conway T, Foster JW: Collaborative regulation of Escherichia coli glutamate-dependent acid resistance by two AraC-like regulators, GadX and GadW (YhiW). J Bacteriol. 2002, 184: 7001-7012. 10.1128/JB.184.24.7001-7012.2002.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  30. Domka J, Lee J, Bansal T, Wood TK: Temporal gene-expression in Escherichia coli K-12 biofilms. Environ Microbiol. 2007, 9: 332-346. 10.1111/j.1462-2920.2006.01143.x.

    Article  CAS  PubMed  Google Scholar 

  31. Francez-Charlot A, Castanie-Cormet MP, Gutierrez C, Cam K: Osmotic regulatin of the Escherichia coli bdm (biofilm-dependent modulation) gene by the RcsCDB His-Asp phosphorelay. J Bacteriol. 2005, 187: 3873-3877. 10.1128/JB.187.11.3873-3877.2005.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  32. Nakamura K, Yamaki M, Sarada M, Nakayama S, Vibat CR, Gennis RB, Nakayashiki T, Inokuchi H, Kojima S, Kita K: Two hydrophobic subunits are essential for the heme b ligation and functional assembly of complex II (succinate-ubiquinone oxidoreductase) from Escherichia coli. J Biol Chem. 1996, 271: 521-527. 10.1074/jbc.271.26.15336.

    Article  CAS  PubMed  Google Scholar 

  33. Letoffe S, Delepelaire P, Wandersman C: The housekeeping dipeptide permease is the Escherichia coli heme transporter and functions with two optional peptide binding proteins. Proc Natl Acad Sci USA. 2006, 103: 12891-12896. 10.1073/pnas.0605440103.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  34. Spoering AL, Vulic M, Lewis K: GlpD and PlsB participate in persister cell formation in Escherichia coli . J Bacteriol. 2006, 188: 5136-5144. 10.1128/JB.00369-06.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  35. Hagiwara D, Sugiura M, Oshima T, Mori H, Aiba H, Yamashino T, Mizuno T: Genome-wide analysis revealing a signaling network of the RcsC-YojN-RcsB phosphorelay system in Escherichia coli. J Bacteriol. 2003, 185: 5735-5746. 10.1128/JB.185.19.5735-5746.2003.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  36. Nagahama H, Oshima T, Mori H, Matsumoto K, Hara H: Hyperexpression of the osmB gene in the acidic phospholipid-deficient Escherichia coli mutant. J Gen Appl Microbiol. 2007, 53: 143-151. 10.2323/jgam.53.143.

    Article  CAS  PubMed  Google Scholar 

  37. Weber A, Kogl SA, Jung K: Time-dependent proteome alterations under osmotic stress during aerobic and anaerobic growth in Escherichia coli. J Bacteriol. 2006, 188: 7165-7175. 10.1128/JB.00508-06.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  38. Bearson S, Bearson B, Foster JW: Acid stress response in enterobacteria. FEMS Microbiol Lett. 1997, 147: 173-180. 10.1111/j.1574-6968.1997.tb10238.x.

    Article  CAS  PubMed  Google Scholar 

  39. Puskarova A, Ferianc P, Kormanec J, Homerova D, Farewell A, Nystrom T: Regulation of yodA encoding a novel cadmium-induced protein in Escherichia coli. Microbiol. 2002, 148 (12): 3801-3811.

    Article  CAS  Google Scholar 

  40. Vassinova N, Kozyrev D: A method of direct cloning of Fur-regulated genes: identification of seven new Fur-regulated loci in Escherichia coli. Microbiol. 2000, 146 (12): 3171-3182.

    Article  CAS  Google Scholar 

  41. Monje-Casas F, Jurado J, Prieto-Alamo MJ, Holmgren A, Puevo C: Expression analysis of the nrdHIEF operon from Escherichia coli. Conditions that trigger the transcript level in vivo. J Biol Chem. 2001, 276: 18031-18037. 10.1074/jbc.M011728200.

    Article  CAS  PubMed  Google Scholar 

  42. Jordan A, Aslund F, Pontis E, Reichard P, Holmgren A: Characterization of NrdH. A glutaredoxin-like protein with a thioredoxin-like activity profile. J Biol Chem. 1997, 272: 18044-18050. 10.1074/jbc.272.29.18044.

    Article  CAS  PubMed  Google Scholar 

  43. Adams MA, Jia Z: Structural and biochemical evidence for an enzymatic quinone redox cycle in Escherichia coli. J Biol Chem. 2005, 280: 8358-8363. 10.1074/jbc.M412637200.

    Article  CAS  PubMed  Google Scholar 

  44. Chen L, Vitkup D: Predicting genes for orphan metabolic activities using phylogenetic profiles. Genome Biol. 2006, 7: R17-10.1186/gb-2006-7-2-r17.

    Article  PubMed Central  PubMed  Google Scholar 

  45. Danese PN, Silhavy TJ: CpxP, a stress-combative member of the Cpx regulon. J Bacteriol. 1998, 180: 831-839.

    PubMed Central  CAS  PubMed  Google Scholar 

  46. Raivio TL, Popkin DL, Silhavy TJ: The Cpx envelope stress response is controlled by amplification and feedback inhibition. J Bacteriol. 1999, 181: 5263-5272.

    PubMed Central  CAS  PubMed  Google Scholar 

  47. Bingham RJ, Hall KS, Slonczewski JL: Alkaline induction of a novel gene locus, alx, in Escherichia coli. J Bacteriol. 1990, 172: 2184-2186.

    PubMed Central  CAS  PubMed  Google Scholar 

  48. Stancik LM, Stancik DM, Schimidt B, Barnhart DM, Yoncheva YN, Slonczewski JL: pH-dependent expression of periplasmic proteins and amino acid catabolism in Escherichia coli. J Bacteriol. 2002, 184: 4246-4258. 10.1128/JB.184.15.4246-4258.2002.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  49. Yohannes E, Barnhart DM, Slonczewski JL: pH-dependent catabolic protein expression during anaerobic growth of Escherichia coli K-12. J Bacteriol. 2004, 186: 192-199. 10.1128/JB.186.1.192-199.2004.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  50. Yohannes E, Thurber AE, Wilks JC, Tate DP, Slonczewski JL: Polyamine stress at high pH in Escherichia coli K-12. BMC Microbiology. 2005, 5: 59-10.1186/1471-2180-5-59.

    Article  PubMed Central  PubMed  Google Scholar 

  51. Tabor CW, Tabor H: The speEspeD Operon of Escherichia coli: formation and processing of a proenzyme form of s-adenosylmenthionine decarboxylase. J Biol Chem. 1987, 262: 16037-16040.

    CAS  PubMed  Google Scholar 

  52. Igarashi K, Kashiwagi K: Polyamine transport in bacteria and yeast. Biochem J. 1999, 344 (3): 633-642. 10.1042/0264-6021:3440633.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  53. Leif H, Sled VD, Ohnishi T, Weiss H, Friedrich T: Isolation and characterization of the proton-translocating NADH: ubiquinone oxidoreductase from Escherichia coli. Eur J Biochem. 1995, 230: 538-548. 10.1111/j.1432-1033.1995.tb20594.x.

    Article  CAS  PubMed  Google Scholar 

  54. Friedrich T: The NADH:ubiquinone oxidoreductase (complex 1) from Escherichia coli. Biochim Biophys Acta. 1998, 1364 (2): 134-146. 10.1016/S0005-2728(98)00024-3.

    Article  CAS  PubMed  Google Scholar 

  55. Partridge JD, Scott C, Tang Y, Poole RK, Green J: Escherichia coli transcriptome dynamics during the transition from anaerobic to aerobic conditions. J Biol Chem. 2006, 281: 27806-27815. 10.1074/jbc.M603450200.

    Article  CAS  PubMed  Google Scholar 

  56. Rohrwild M, Coux O, Huang HC, Moerschell RP, Yoo SJ, Chung CH, Golderg AL: HslV-HslU: A novel ATP-dependent protease complex in Escherichia coli related to the eukaryotic proteasome. Proc Natl Acad Sci USA. 1996, 93: 5808-5813. 10.1073/pnas.93.12.5808.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  57. Smith MW, Neidhardt FC: 2-Oxoacid dehydrogenase complexes of Escherichia coli: cellular amounts and patterns of synthesis. J Bacteriol. 1983, 156: 81-88.

    PubMed Central  CAS  PubMed  Google Scholar 

  58. Lee PS, Grossman AD: The chromosome partitioning proteins Soj (ParA) and Spo0J (ParB) contribute to accurate chromosome partitioning, separation of replicated sister origins, and regulation of replication initiation in Bacillus subtilis. Mol Microbiol. 2006, 60 (4): 853-869. 10.1111/j.1365-2958.2006.05140.x.

    Article  CAS  PubMed  Google Scholar 

  59. Bhagwat AA: Regulation of the glutamate-dependent acid-resistance system of diarrheagenic Escherichia coli strains . FEMS Microbiol Lett. 2003, 227: 39-45. 10.1016/S0378-1097(03)00646-3.

    Article  CAS  PubMed  Google Scholar 

Download references


This work was supported by grant MCB-0644167 from the National Science Foundation. Fluorimetry and microarray analysis included support from an HHMI 2004 Undergraduate Science Education Award to Kenyon College.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Joan L Slonczewski.

Additional information

Authors' contributions

GK and JCW designed and conducted the experiments and authored the first draft manuscript. DMF contributed real-time PCR assays. BDJ devised and conducted the statistical analysis. SSB directed the array hybridization. JLS conceived the study, coordinated the project overall, and revised the final manuscript. All authors have read and approved the final manuscript.

Electronic supplementary material


Additional file 1: Expression indices, expression ratios and percent recovery of genes that show time-dependent expression in response to a rapid change in external pH. (XLS 394 KB)

Additional file 2: Table showing Genes downregulated at least 2-fold after acid shift. (DOC 214 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Rights and permissions

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

Reprints and permissions

About this article

Cite this article

Kannan, G., Wilks, J.C., Fitzgerald, D.M. et al. Rapid acid treatment of Escherichia coli: transcriptomic response and recovery. BMC Microbiol 8, 37 (2008).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: