Characterization of acetic acid-detoxifying Escherichia coli evolved under phosphate starvation conditions
© Moreau and Loiseau. 2016
Received: 21 September 2015
Accepted: 7 February 2016
Published: 19 February 2016
During prolonged incubation of Escherichia coli K-12 in batch culture under aerobic, phosphate (Pi) starvation conditions, excess glucose is converted into acetic acid, which may trigger cell death. Following serial cultures, we isolated five evolved strains in two populations that survived prolonged incubation.
We sequenced the genomes of the ancestral and evolved strains, and determined the effects of the genetic changes, tested alone and in combination, on characteristic phenotypes in pure and in mixed cultures.
Evolved strains used two main strategies: (1) the constitutive expression of the Trk- and Kdp-dependent K+ transport systems, and (2) the inactivation of the ArcA global regulator. Both processes helped to maintain a residual activity of the tricarboxylic acid cycle, which decreased the production of acetic acid and eventually allowed its re-consumption. Evolved strains acquired a few additional genetic changes besides the trkH, kdpD and arcA mutations, which might increase the scavenging of organophosphates (phnE +, lapB, and rseP) and the resistance to oxidative (rsxC) and acetic acid stresses (e14−/icd +).
Evolved strains rapidly acquired mutations (phnE + lapB rpoS trkH and phnE + rseP kdpD) that were globally beneficial to growth on glucose and organophosphates, but detrimental to long-term viability. The spread of these mutant strains might give the ancestral strain time to accumulate up to five genetic changes (phnE + arcA rsxC crfC e14−/icd +), which allowed growth on glucose and organophosphates, and provided a long-term survival. The latter strain, which expressed several mechanisms of protection against endogenous and exogenous stresses, might provide a platform for producing toxic recombinant proteins and chemicals during prolonged incubation under aerobic, Pi starvation conditions.
In nature, bacteria are frequently starved for essential nutrients such as phosphate (HPO4 2−, Pi) [1, 2]. When the levels of Pi in the medium decrease below 4 μM, Escherichia coli induces the PhoBR regulon, which helps to scavenge low levels of Pi and secondary sources of Pi such as organophosphates and phosphonates . When the cells cannot find a sufficient source of Pi to maintain growth, they enter the so-called stationary phase . The E. coli K-12 strain MG1655 incubated in batch culture under aerobic, Pi-limiting conditions (0.1 mM K2HPO4 and 40 mM glucose) enters stationary phase after 10 h of incubation when Pi levels drop below 1 μM .
Whereas E. coli is normally doomed to die when glutamate is not added in Pi-limiting medium, evolved strains can sweep populations after generally three serial batch cultures. Variants that survive prolonged incubation exhibit different mutant phenotypes refereed to as RpoS− (colonies grown on LB medium produce low levels of glycogen and exhibit a reduced catalase activity) and Glg− (colonies produce reduced levels of glycogen) . From one evolved population, we isolated three strains exhibiting RpoS− (ENZ1901), Glg− (ENZ1902), and Glg+ (ENZ1903) phenotypes. In a parallel-evolved population, which contained a majority of Glg− mutants and eventually Glg+ cells, we isolated two Glg+ evolved strains (ENZ1904 and ENZ1905) .
Compared to the ancestral strain, the evolved strain ENZ1901 (RpoS−) harbors a single-nucleotide deletion (rpoS-ΔG214) that inactivates RpoS, and an 8-bp deletion in phnE that activates the PhnE permease . The lack of RpoS activity allows a substantial activity of the tricarboxylic acid cycle in Pi-starved cells, which decreases the production of acetic acid and allows its re-consumption (Ace+ phenotype) [1, 5, 11, 12]. The PhnE permease can scavenge organophosphates excreted into the medium . The combination of the phnE and rpoS evolved alleles triggers a novel growth-under-Pi-starvation (GPS) phenotype, which refers to the ability of a mutant strain—incubated for 1 day in Pi-limiting medium and added in minority into a 1-day-old culture of the ancestral strain in Pi-limiting medium—to resume growth between days 1 and 4 of incubation by using primarily organophosphates released into the mixed culture by the ancestral strain . The GPS phenotype may account at least in part for the spread of the evolved strain ENZ1901 in populations starved for Pi.
If the presence of an rpoS-null mutation may help to explain the behavior of the evolved strain ENZ1901 (RpoS−), the majority of the evolved strains did not exhibit an RpoS-negative phenotype . This prompted us to determine the complete genotypes, the metabolic patterns and the growth capacities in Pi-limiting medium of the evolved and reconstructed mutant strains. We show here that mutants could overtake populations and detoxify acetic acid by using two novel strategies: (1) K+ accumulation (triggered by trkH or kdpD mutations) combined with changes in outer membrane permeability (lapB or rseP mutations), and (2) inactivation of ArcA activity (arcA mutation) combined with an increased resistance to oxidative stress (rsxC mutation) and to acetic acid stress (crfC and e14 − / icd + mutations).
Genotype of the evolved strains
E. coli K-12 strains
F− λ− rfb-50 rph-1 ilvG phnE ylbE + icd::e14
As 535 but ΔlacIZ
As 1734 (535 ΔlacIZ) but poxB176::lacZ.Cm r
As 535 but ΔlacY::kan
As 535 but lacY::Tn10
ΔlacX74 (λ rpoHP3::lacZ)
ΔlacU169 (λ cydA::lacZ-km)
ΔlacU (λ PsoxS::lacZ)
As 535 but phnE + trkH80 rpoS214 lapB43
Evolved 535 RpoS− 
As 535 but phnE + kdpD460 rseP98
Evolved 535 Glg− 
As 535 but phnE + arcA172 rsxC525 crfC693 e14−/icd +
Evolved 535 Glg+ 
As 535 but phnE +
Evolved 535 Glg+ 
As 535 but phnE +
Evolved 535 Glg+ 
As 535 but ΔrpoS::kan
As 1734 (535 ΔlacIZ) but ΔrpoS::kan
As 1944 (535 ΔrpoS) but Kans
As 1945 (535 ΔlacIZ ΔrpoS) but Kans
As 1903 but ΔlacY::kan
As 1905 but ΔlacY::kan
As 1901 but ΔlacY::kan
As 1902 but ΔlacY::kan
As 1904 but ΔlacY::kan
As 1944 (535 ΔrpoS::kan) but cysC95::Tn10
As 1946 (535 ΔrpoS) but ΔlacY::kan
As 1901 but ΔlacIZ
As 1902 but ΔlacIZ
As 1903 but ΔlacIZ
As 535 but rpoS214
As 1901 but rpoS + (phnE + trkH80 lapB43)
As 1902 (phnE + kdpD460 rseP98)
As 1734 (535 ΔlacIZ) but poxB176::lacZ.Cm r
As 1947 (535 ΔlacIZ ΔrpoS) but poxB176::lacZ.Cm r
As 2041 (535 rpoS214) but ΔlacY::kan
As 2043 (phnE + trkH80 lapB43) but ΔlacY::kan
As 2044 (phnE + kdpD460 rseP98) but ΔlacY::kan
As 2044 but ΔybfH::kan kdpD + (phnE + rseP98)
As 2043 (1901 rpoS +) but ΔyihL::kan (phnE + lapB43 trkH80)
As 2043 (1901 rpoS +) but ΔyihL::kan trkH + (phnE + lapB43)
As 1902 but lacY::Tn10
As 2044 (1902) but ΔybfH::kan
As 535 but ΔyihL::kan trkH80
As 535 but ΔyihL::kan trkH +
As 535 but ΔybfH::kan kdpD460
As 535 but ΔybfH::kan kdpD +
As 535 but ΔproP::kan phnE + rpoS214
As 1901 but lacY::Tn10
As 535 but ΔtrkH::kan
As 535 but phnE +
As 2309 (535 ΔtrkH) but phnE +
As 1902 but ΔyaeH::kan
As 1734 (535 ΔlacIZ) but phnE +
As 1832 (rpoHP3::lacZ) but phnE +
As 1842 (Pkdp::lacZ) but phnE +
As 2045 (poxB176::lacZ.Cm r ) but phnE +
As 2329 (535 ΔlacIZ phnE +) but poxB176::lacZ.Cm r
As 2032 (1901 ΔlacIZ) but poxB176::lacZ.Cm r
As 2033 (1902 ΔlacIZ) but poxB176::lacZ.Cm r
As 2034 (1903 ΔlacIZ) but cydA::lacZ-kan
As 2315 (535 phnE +) but ΔlacY::kan
As 2315 (535 phnE +) but lacY::Tn10
As 2337 (Pkdp::lacZ) but ΔybfH::kan kdpD460
As 2337 (Pkdp::lacZ) but ΔybfH::kan kdpD +
As 2163 (535 ΔyihL::kan trkH80) but phnE +
As 2315 (535 phnE +) but ΔybfH::kan kdpD +
As 2315 (535 phnE +) but ΔybfH::kan kdpD460
As 2315 (535 phnE +) but ΔyihL::kan trkH80
As 535 but ΔarcA::kan
As 535 but ΔrsxC::kan
As 1903 but ΔarcA::kan
As 1903 but ΔrsxC::kan
As 2034 (1903 ΔlacIZ) but cydA::lacZ-kan
As 1843 (PsoxS::lacZ) but ΔarcA::kan
As 1843 (PsoxS::lacZ) but ΔrsxC::kan
As 2413 (535 ΔarcA) but Kans
As 2420 (PsoxS::lacZ ΔarcA) but Kans
As 2424 (535 ΔarcA) but ΔrsxC::kan
As 2430 (PsoxS::lacZ ΔarcA) but ΔrsxC::kan
As 2434 (535 ΔarcA ΔrsxC::kan) but phnE +
As 2413 (535 ΔarcA::kan) but phnE +
As 2413 (535 ΔrsxC::kan) but phnE +
As 2333 (rpoHP3::lacZ phnE +) but ΔrseA::kan
As 2333 (rpoHP3::lacZ phnE +) but ΔyaeH::kan rseP +
As 2333 (rpoHP3::lacZ phnE +) but ΔyaeH::kan rseP98
As 2341 (poxB176::lacZ.Cm r ) but ΔyihL::kan trkH +
As 2341 (poxB176::lacZ.Cm r ) but ΔyihL::kan trkH80
As 2341 (poxB176::lacZ.Cm r ) but ΔybfH::kan kdpD +
As 2341 (poxB176::lacZ.Cm r ) but ΔybfH::kan kdpD460
In the ancestral strain ENZ535, the sequence of the ylbE gene differed from the reference sequence by a single-nucleotide polymorphism and a single-nucleotide insertion that removed a stop codon. The so-called ylbE + allele in ENZ535 is also present in the MG1655 strain referred to as ATCC47076 . The function of YlbE is not known.
In the five evolved strains (PhnE-expressing), the sequence of the phnE gene differed from the reference sequence by a deletion of 8 bp, which removed a frameshift inactivating phnE .
The evolved strains ENZ1904 and ENZ1905 (Glg+) acquired no others mutations besides the 8-bp deletion in phnE (phnE + ).
- 4.The evolved strains ENZ1901 (RpoS−), ENZ1902 (Glg−), and ENZ1903 (Glg+) harbored unique genetic changes scattered around the chromosome (Fig. 2): rpoS-ΔG214, trkH-L80Q and lapB-V43G in ENZ1901; kdpD-D460V and rseP-A98V in ENZ1902; and arcA-F172Y, rsxC-I525F, crfC-E693Stop and e14− /icd + in ENZ1903. In the latter strain, excision of the lambdoid prophage e14 restored the icd + sequence that differed from the reference sequence, icdA, by 12 single-nucleotide polymorphisms [18, 19] (Fig. 3).
Considering the phenotypic and genetic diversity that occurred in evolved populations, a number of questions arise regarding (1) the minimum number of mutations required by a variant to sweep a population [20–23], (2) the types of mutations (loss- or gain-of-function; global regulatory or specific changes) [20–27], and (3) the diversity or convergence of the metabolic changes occurring in evolved strains [28, 29].
In ENZ1901 (phnE + rpoS214 trkH80 lapB43), the trkH80 and rpoS214 mutations similarly affect the metabolic flux
To analyze the effects of the trkH80 and lapB43 mutations present in the evolved strain ENZ1901 (phnE + rpoS214 trkH80 lapB43), we exchanged the rpoS214 allele for the rpoS + allele from the ancestral strain ENZ535, which gave rise to the strain ENZ2043 (phnE + trkH80 lapB43). This strain exhibited a Glg − phenotype, whereas strains that harbored a single phnE + allele (e.g. ENZ1904) exhibited a Glg + phenotype. Thus, besides the rpoS214 mutation, which decreased the cellular levels of glycogen and of catalase, the trkH80 and/or lapB43 mutation(s) somehow decreased the cellular levels of glycogen.
LapB is an inner-membrane protein that plays a key role in the biogenesis of the glycolipid lipopolysaccharide leaflet in the outer membrane [30, 31]; at first glance, it seemed unlikely that a lapB mutation might affect the glycogen content, a storage compound. In contrast, changes in the activity of TrkH, the primary potassium (K+) transporter, might possibly change the levels of glycogen because the deregulation of K+ homeostasis, which normally occurs in response to an osmotic stress, affects many metabolic processes .
To determine whether mutations in trkH could change the Glg phenotype, we transduced the ΔtrkH::kan deletion and the trkH80 evolved mutation (co-transduced with ΔyihL::kan) into the ancestral strain ENZ535. Compared to the parental strain, the ΔtrkH (ENZ2317) and trkH80 (ENZ2163) mutant strains produced more and less glycogen, respectively. Considering that these opposite effects on glycogen contents might reflect changes in K+ transport, and that the rate of transport of K+ is decreased in ΔtrkH mutants , we concluded that the rate of transport of K+ could be increased in trkH80 mutants. The TrkH-L80Q mutant protein might therefore mimic the form of TrkH (or of the TrkAEGH complex) that is normally present following an osmotic stress, when the rate of transport of K+ through TrkH transiently increases .
RpoS activity measured with the levels of expression of the poxB176::lacZ fusion in strains starved for N or P
24 h in N-limiting medium
24 h in P-limiting medium
2045: rpoS +
104 ± 1 (n = 3)
439 ± 25 (n = 3)
1 ± 0.2 (n = 2)
4 ± 1 (n = 2)
2343: rpoS214 trkH80 lapB43
1 ± 0.2 (n = 2)
3 ± 0.4 (n = 2)
2493: yihL::kan trkH +
108 ± 5 (n = 3)
437 ± 2 (n = 3)
2494: yihL::kan trkH80
104 ± 4 (n = 2)
404 ± 64 (n = 2)
2345: kdpD460 rseP98
109 ± 10 (n = 2)
435 ± 5 (n = 2)
2495: ybfH::kan kdpD +
112 ± 2 (n = 3)
430 ± 63 (n = 3)
2496: ybfH::kan kdpD460
113 ± 6 (n = 3)
469 ± 88 (n = 3)
Effects of the trkH80, lapB43 and rpoS214 mutations in mixed cultures
The efficiency of the GPS phenotype (determined on day 5 of incubation) decreased in the order: ≈5 × 107 CFU/ml for phnE + trkH80 lapB43 rpoS214 mutants (ENZ2000) (Fig. 5a) ; ≈2 × 107 CFU/ml for phnE + trkH80 lapB43 mutants (ENZ2067 and ENZ2136) (Fig. 5c, d); ≈2 × 106 CFU/ml for phnE + trkH80 (ENZ2395) (Fig. 5e); ≈7 × 104 CFU/ml for phnE + lapB43 mutants (ENZ2137) (Fig. 5g); ≈ 3 × 104 CFU/ml for single trkH80 mutants (ENZ2163) (Fig. 5f); and ≈104 CFU/ml for single phnE + (ENZ1984 and ENZ2003) (Fig. 5h) and parental strains (ENZ1791 and ENZ2166) (Fig. 5b) .
When the co-cultures were re-diluted into fresh medium, Kanr strains were generally overtook by the ancestral strain (Tcr) and both strains lost viability owing to the high levels of acetic acid excreted (up to 30 mM acetic acid at pH 4.8, on day 34 of incubation). This notably occurred in the case of the evolved strain ENZ2000 (phnE + trkH80 lapB43 rpoS214) (Fig. 5a). In sharp contrast, triple mutants phnE + trkH80 lapB43 (ENZ2067 and ENZ2136) overtook the ancestral strain, survived prolonged incubation, and afforded a robust cross-protection to the ancestral strain as a result of the detoxification of acetic acid in the incubation medium (≤0.02 mM glucose and acetate at pH 6.8, on day 35 of incubation) (Fig. 5c, d). The phnE + trkH80 double mutants (ENZ2395) sometimes failed to overtake the ancestral strain, which then triggered the death of the whole population during prolonged incubation (Fig. 5e). The phnE + lapB43 strain (ENZ2137) was overtook by the ancestral strain and died, whereas it rapidly evolved RpoS− variants (Fig. 5g).
Therefore, the lapB43 mutation, which had no effect on the metabolic pattern in monoculture (Fig. 4) and only a weak effect on the GPS phenotype, might improve the GPS phenotype and the long-term viability of phnE + trkH80 lapB43 triple mutants in mixed culture. In contrast, the addition of the rpoS214 mutation into the phnE + trkH80 lapB43 strain slightly improved the GPS phenotype but dramatically decreased the long-term viability, which indicates that the presence of the rpoS214 mutation in the evolved strain ENZ1901 (phnE + trkH80 lapB43 rpoS214) was eventually deleterious.
The kdpD-D460V mutation in ENZ1902 triggers a constitutive expression of the kdpFABC operon
The rseP-A98V mutation triggers a GPS phenotype
Compared to the evolved strain ENZ1902 (phnE + kdpD460 rseP98), the reconstructed strain ENZ2402 (phnE + kdpD460), which lacked the rseP98 mutation, generally re-consumed acetic acid more efficiently in monoculture (Fig. 6g, h), and in mixed culture (Fig. 7a, b). These results suggest that in the evolved strain ENZ1902, the kdpD460 mutation played a beneficial role, whereas the rseP98 mutation could play a detrimental role in metabolism.
To further characterize the role of the rseP98 mutation in metabolism, we exchanged the kdpD460 mutant allele for the kdpD + allele in ENZ1902, and tested the new strain, ENZ2112 (phnE + rseP98), in monoculture and in mixed culture. In monoculture, the strain ENZ2112 excreted somewhat less acetic acid and consumed less glucose than a phnE + strain (Fig. 6f, h), which confirmed that the rseP98 mutation might cause a slight defect in metabolism during prolonged incubation. Surprisingly, the strain ENZ2112 (phnE + rseP98) exhibited a strong GPS phenotype but did not survive prolonged incubation in mixed culture (Fig. 7c).
The RpoE activity measured with the levels of expression of the rpoHP3::lacZ fusion is not affected by the rseP98 mutation
Exponential growth phase
(7 h in P-limiting medium)
(24 h in P-limiting medium)
37 ± 3 (n = 2)
277 ± 21 (n = 9)
2480: 2333 ΔrseA::kan
142 ± 1 (n = 2)
1043 ± 6 (n = 2)
2489: 2333 yaeH::kan rseP +
39 ± 1 (n = 2)
282 ± 19 (n = 4)
2490: 2333 yaeH::kan rseP98
37 ± 1 (n = 2)
292 ± 3 (n = 4)
Synergy of arcA and rsxC mutations on acetic acid detoxification in ENZ1903
The dual regulator ArcA primarily inhibits the expression of genes implicated in the activity of the tricarboxylic acid cycle (e.g. acnB and fumA) and of the aerobic respiratory chain when the concentration of dissolved oxygen decreases. ArcA also increases the synthesis of CydAB, an enzyme of the aerobic respiratory chain that scavenges low levels of oxygen, which increases the energetic metabolism and extents the protection of cytoplasmic proteins (e.g. FNR) against oxygen at low pO2 [42, 43]. Inactivation of arcA increases the metabolic flux into the tricarboxylic acid cycle in exponentially growing cells .
RsxC belongs to the Rsx NADPH-dependent reducing system that inactivates SoxR and thus shunts off the SoxRS response to oxidative stress . The SoxRS response comprises acnA and fumC, which encode proteins of the tricarboxylic acid cycle that are resistant to oxidative stress . In exponentially growing cells, a ΔrsxC mutation increases the basal level of expression of the SoxRS regulon primarily because the [2Fe-2S] clusters in SoxR that are spontaneously oxidized cannot be re-reduced . In addition, in arcA mutants, the oxygen-sensitive transcription factor FNR should be inactive, which might increase the cellular levels of SoxR and SoxS [42, 47].
During growth on acetate, the activity of Icd (isocitrate dehydrogenase) is decreased as a result of its phosphorylation, which helps to redirect part of the metabolic flux from the tricarboxylic acid cycle towards the glyoxylate shunt . In E. coli K-12, the phage e14 is normally integrated in icd (Fig. 3). Although the sequences of icd in lysogenic (e14+; icdA) and in non-lysogenic strains (e14−; icd +) differ by 12 single-nucleotide polymorphisms, 10 synonymous and two non-synonymous, both forms of the Icd protein exhibit the same activity in vitro. Recent data suggest that excision of the prophage e14 and/or restoration of the icd + gene could somehow decrease or increase the resistance of E. coli K-12 to different stresses [18, 19].
We found evidence that the arcA-F172Y allele encoded an inactive protein. First, we transduced a cydA::lacZ fusion, which is induced by ArcA , into the evolved strain ENZ1903 and into the ancestral strain ENZ535. After 24 h of incubation under aerobic, Pi starvation conditions, the levels of β-galactosidase were 3.4-fold lower in ENZ1903 (1787 ± 303 units, n = 2) than in ENZ535 derivatives (6062 ± 401 units, n = 2), which suggests that the evolved strain ENZ1903 might exhibit little, if any, ArcA activity . Second, we exchanged the arcA172 allele in ENZ1903 for a ΔarcA::kan deletion (ENZ2415), which did not change significantly the metabolic patterns of the strains in monocultures (Fig. 11a, c, e, g). Similarly, exchange of the rsxC525 allele in ENZ1903 for a ΔrsxC::kan deletion (ENZ2416) had no significant effect on the metabolic patterns, which suggests that the rsxC-I525F mutation was null or neutral (Fig. 11a, c, e, g).
We transduced the ΔarcA::kan mutation (the kanamycin resistance marker could be eventually excised) and/or the ΔrsxC::kan mutation into the ancestral strain ENZ535, and determined the metabolic patterns of the different strains during incubation in monocultures (Fig. 11). Whereas the ΔrsxC mutation alone (ENZ2414) had no effect on the cellular viability and on the kinetics of production of acetic acid, the ΔarcA mutation alone (ENZ2413) had moderate effects: values were in-between the values obtained for the ancestral strain and the evolved strain ENZ1903 (Fig. 11b, d, f, h). Surprisingly, transduction of the ΔrsxC::kan mutation into the ΔarcA mutant strain (ENZ2434: ENZ535 ΔarcA ΔrsxC) dramatically improved the viability of the strain and its ability to consume the low levels of acetic acid that were previously excreted (Fig. 11b, d, f, h).
A ΔrsxC mutation increases the SoxRS oxidative stress response in ΔarcA mutants
Could the synergy of the ΔarcA and ΔrsxC mutations be explained by an effect of the rsxC mutation on the SoxRS oxidative stress response? To answer this question, we transduced the ΔarcA and/or ΔrsxC mutations into the strain ENZ1843, which carries a soxS::lacZ fusion that is induced by SoxR in its active (oxidized) form , and we measured the β-galactosidase activity in exponentially growing cells (6 h of incubation; OD600 = 0.1) and in non-growing cells (24 h of incubation in P-limiting medium).
The arcA and rsxC mutations are not sufficient for a long-term viability in mixed culture
The constructed strain ENZ2439 (ENZ535 phnE + ΔarcA ΔrsxC::kan), added in minority in a mixed culture, grew vigorously but died more rapidly than the evolved strain ENZ1982 (ENZ1903 ΔlacY::kan) between days 6 and 9 of incubation (Fig. 12a, b). Moreover, following the dilution of the mixed cultures into fresh medium on day 9 of incubation, the reconstructed strain ENZ2439 grew slowly and died off (Fig. 12b), whereas the evolved strain ENZ1982 overtook the population and survived prolonged incubation (Fig. 12a). A simple interpretation of these results is that the icd + and/or crfC693 alleles, present in the evolved strain (phnE + arcA172 rsxC525 e14− /icd + crfC693), but not in the reconstructed strain ENZ2439 (phnE + ΔarcA ΔrsxC icd::e14), might help cells to tolerate a stress that built up in mixed culture.
It is generally thought that cells that evolve under nutrient limitation conditions primarily accumulate mutations that decrease the activity of the RpoS (σs) factor [18, 20–22, 25, 26, 28, 29, 49]. We show here that in strains that evolved in serial batch cultures under aerobic, Pi starvation conditions, rpoS-null mutations, which were first beneficial, were eventually detrimental. In contrast, we identified novel combinations of mutations that were beneficial to the long-term viability of the cells.
K+ homeostasis and metabolic flux
The trkH80 and kdpD460 mutations played a key role in the survival of the evolved strains ENZ1901 (phnE + rpoS-ΔG214 trkH-L80Q lapB-V43G) and ENZ1902 (phnE + kdpD-D460V rseP-A98V), respectively. The trkH80 and kdpD460 mutations triggered a constitutive expression of the low-affinity high-rate Trk system and of the high-affinity low-rate Kdp system, respectively.
The trkH-L80Q mutation was previously identified in strains that evolved in Pi-limited chemostat cultures and in serial batch cultures containing the antibiotic gentamicin sulfate [25, 50]. In the strain evolved in Pi-limited chemostat, the trkH-L80Q mutation was not further characterized because of the presence of numerous mutations . In strains resistant to gentamicin, it was suggested that excess K+ could prevent the transport of the antibiotic inside the cells [50, 51]. The kdpD-D460V constitutive mutation was not previously described but the KdpD-D460V evolved protein is reminiscent of the KdpD-D474A engineered protein, which cannot sense K+ in the periplasm and thus constitutively activates the KdpE regulator of the kdpFABC operon .
How excess K+ could change the metabolism of glucose and of acetic acid in evolved strains? Potassium glutamate constitutes the most abundant ions inside the cell; K+ transport is linked to Pi transport and to internal pH (K+ transport inside the cell triggers a net H+ excretion) [32, 33, 36, 38, 52]. The activity of the K+ transporters can rapidly increase when the osmotic pressure of the medium increases and when the concentration of K+ in the medium decreases [32, 36], but excess K+ may be eventually toxic . Interestingly, IlvBN, the first enzyme in the pathway of biosynthesis of branched-chain amino acids is especially sensitive to allosteric inhibition by K+ . Because ilvB evolved through duplication/deletion events from poxB , it is tempting to speculate that both IlvB and PoxB might be inhibited by excess K+, which would direct the metabolic flux from pyruvate towards the tricarboxylic acid cycle rather than towards the synthesis of branched-chain amino acids and especially acetic acid in Pi-starved cells (Fig. 1).
Changes in outer membrane protein activities
The rseP98 mutation, present in ENZ1902, triggered growth on glucose and organophosphates but hampered growth on acetic acid and Pi, which threatened the long-term viability of the cells in mixed culture. Although the primary role of RseP is to cleave the anti-sigma E factor, RseA , the rseP98 mutation had no effect on the level of expression of the RpoE (σE) regulon, which was normally induced in Pi-starved cells [5, 7]. Because inner membrane metalloproteins such as FtsH and RseP may exhibit dual functions in membrane protein biogenesis, as protease and as chaperone [56–58], it is tempting to speculate that the RseP-A98V protein might affect the folding of outer membrane proteins, which might increase the PhnE-dependent scavenging of organophosphates and decrease the PhoE-dependent scavenging of Pi (Fig. 14). Such trade-off would be reminiscent of the apparent changes in activities of the porins OmpC and PhoE, which occur in chemostat cultures limited in glucose .
Resistance to oxidative and acetic acid stresses
The arcA172 mutation might play a key role in the metabolic changes that occurred in the evolved strain ENZ1903. The global regulator ArcA progressively inhibits the transcription of genes of the tricarboxylic acid cycle when the concentration of dissolved oxygen decreases, which redirects the metabolic flux towards the production of weak acids such as acetic acid : the inactivation of ArcA might thus maintain the activity of the tricarboxylic acid cycle in Pi-starved cells, which might decrease the production of acetic acid and allow its re-consumption (Fig. 1). Although different arcA mutations occurred in strains that evolved under nutrient stress conditions, the arcA172 mutation, which affects the DNA-binding domain of ArcA, was not previously described [21, 22, 26].
A single ΔarcA mutation moderately decreased the production of acetic acid under aerobic, Pi starvation conditions, which suggests that the tricarboxylic acid cycle was poorly expressed. It is generally thought that the limited activity of the tricarboxylic acid cycle in arcA mutants might result from an insufficient activity of the aerobic respiratory chain, which would limit the re-oxidation of NADH; excess NADH may inhibit the activity of the citrate synthase, GltA, which controls the entry into the tricarboxylic acid cycle [11, 59, 60] (Fig. 1). As our results revealed novel synergistic effects between the ΔarcA and ΔrsxC mutations on the metabolism of acetate and on the expression of the SoxRS oxidative stress response, we would like to suggest that metabolic defects in ΔarcA mutants would result not from the inhibition of enzyme activities by excess NADH but rather from the degradation of enzymes of the tricarboxylic acid cycle (e.g. AcnB and FumA) by reactive oxygen species adventitiously produced by NADH dehydrogenases  (Fig. 1).
In mixed culture, the evolved strain ENZ1903 (phnE + arcA172 rsxC525 e14−/icd + crfC693) survived better than the ENZ535 phnE +ΔarcA ΔrsxC reconstructed strain, which suggests that the e14−/icd + and/or crfC693 genetic changes might protect the strain ENZ1903. Interestingly, strains evolving at pH 4.8 accumulate mutations in icd, including an icd-D398E mutation . We realized that this mutation was identical to a single-nucleotide polymorphism induced by the excision of the prophage e14 from icd: the icdA (icd::e14) and icd + (e14−) alleles differed by 12 single-nucleotide polymorphisms, 10 synonymous and 2 non-synonymous including the D398E change (Fig. 3). Although the enzymes encoded by the icdA and icd + alleles exhibit the same activity in vitro, it is likely that active enzymes were produced at different levels in vivo as a result of changes in codon usage, which may affect the rate of synthesis and the folding of the proteins [62, 63]. We propose that the isocitrate dehydrogenase activity might increase as a result of the e14− /icd + genetic changes, which might increase the levels of α-ketoglutarate, of glutamate (synthesized from α-ketoglutarate), and of NADH (produced by the NAD-dependent α-ketoglutarate dehydrogenase): glutamate may be used by the RpoS-dependent GadB decarboxylase , which may consume protons inside the cell, and the re-oxidation of NADH by the aerobic respiratory chain may trigger the excretion of protons (Fig. 1). Both processes might help the evolved strain ENZ1903 to resist internal acidification by acetic acid (pH 4.8) that was excreted into the medium by the ancestral strain present in majority in mixed culture between days 1 and 9 of incubation. In this light, it is tempting to speculate about how the truncation of CrfC (crfC-E693Stop), which normally tethers clamp-replication forks complexes , could favor the e14−/icd + genetic changes: the CrfC693 mutant protein might destabilize the replication forks , generate excess single-stranded DNA, activate the RecA protein, and trigger the excision of the lambdoid prophage e14 (e14−) [65, 66].
Parallel evolutions that change acetate producers in acetate consumers
Collectively, our results revealed that to survive during prolonged incubation under aerobic, Pi starvation conditions, evolved strains acquired several types of mutations that helped them to grow on glucose and organophosphates (GPS phenotype), to produce low levels of acetic acid, and to combat endogenous oxidative stress and exogenous acid stress. Then, evolved strains that overtook the population could grow on acetic acid and Pi previously excreted (Ace+ phenotype), thereby detoxifying the medium .
Based on the dynamics of evolving populations (generally, RpoS− strains spread first, followed by Glg− and eventually Glg+ strains) , we would like to suggest a possible scenario for the time course of spreading of evolving strains. It is likely that the phnE + change was primarily selected because the 8-bp deletion in phnE, which results from a slippage between direct repeats during replication, occurs with a high frequency . Then, acquisition of the lapB43 and especially rpoS214 mutations, which triggered growth and thus DNA replication under Pi starvation conditions, might increase the probability of acquiring the trkH80 constitutive mutation, thereby generating the strain ENZ1901 that exhibited a high metabolic rate . A similar evolutionary process might occur with the acquisition of the phnE + and rseP98 mutations, which triggered a GPS phenotype, and ultimately of the kdpD460 constitutive mutation, thereby generating the strain ENZ1902. Although the evolved strains ENZ1901 (RpoS−) and ENZ1902 (Glg−) were eventually doomed to die, these strains could provide a transient protection to the population by consuming glucose while producing low levels of acetic acid, which could provide sufficient time for the progressive accumulation of the arcA172 and rsxC525 mutations, and then the crfC93 and e14− /icd + genetic changes, thereby generating the strain ENZ1903. Eventually, during prolonged incubation in detoxified spent medium at pH 7, phnE + single mutant strains (ENZ1904 and ENZ1905) could sweep the population, thereby restoring the ancestral metabolism.
Considering the characteristics of the evolved (ENZ1903) and reconstructed strains (trkH80 lapB43 and kdpD460), we suggest that these strains could be useful in human health and in industrial production of bacterial products.
As previously suggested, E. coli strains that evolved under aerobic, Pi starvation conditions might help to colonize the human intestinal tract, in which both Pi and glutamate are reabsorbed . In this light, it would be interesting to determine whether evolved strains could exhibit probiotic effects, alone or in consortium .
The consumption of energy and amino acids for growth and the generation of metabolic stresses (e.g. envelope stress) can limit the production and excretion of proteins and chemicals. To resolve these problems, it has been suggested to stop growth but not protein synthesis, and to use engineered strains (e.g. pta ackA poxB and arcA iclR mutants) that produce less acetic acid and are more tolerant to stress [68–71]. We would like to suggest that the evolved strain ENZ1903 (phnE + arcA172 rsxC525 crfC693 e14− /icd +) incubated in batch culture under aerobic, Pi starvation conditions could afford a platform for producing toxic proteins and chemicals because of two major characteristics: (1) the sustained activity of the tricarboxylic acid cycle in stationary phase, which helps to produce energy and building blocks rather than acetic acid, and (2) the spontaneous induction of several stress responses, i.e. the SoxRS oxidative stress response and a novel e14− /icd + acid resistance system, which are specifically induced in ENZ1903, in addition to the Cpx and RpoE envelope stress responses and the RpoS general stress response, which are normally induced in Pi-starved cells .
The E. coli K-12 laboratory strain MG1655 is very close to commensal strains . Among several MG1655 isolates , we chose the strain CF1648, therein called ENZ535 (Table 1), because it exhibited a high growth rate and a high survival rate during prolonged incubation in rich medium, in glycerol-minimal medium under aerobic conditions, and in glucose-minimal medium under aerobic and anaerobic conditions [1, 10, 73]. Strains of the Keio collection  were used as donors to transduce single-gene in-frame deletions replaced with a kanamycin-resistance cassette (Kanr). Transductions were performed with phage P1vir and cells grown in LB5 medium containing 5 g/l NaCl [73, 74]. To construct Kans derivatives of deletion mutants, the kan cassette flanked by the FLP-recombination-target sites was removed by introducing the FLP recombinase-expressing plasmid pCP20 at 32 °C and purifying clones at 42 °C . We co-transduced evolved mutations with Kanr markers from the Keio collection: trkH with ΔyihL::kan, kdpD with ΔybfH::kan, and rseP with ΔyaeH::kan. The ΔlacIZ mutants were constructed as previously described . Spontaneous PhnE+ strains were selected on MOPS minimal medium plates containing 0.05 mM methyl phosphonate as the sole source of P .
Before we sequenced the genomes of the evolved strains, we exchanged rpoS alleles for rpoS + in ENZ1901 and ENZ1902. Practically, we transduced the cysC::Tn10 tetracycline-resistance marker (Tcr; 7 kb apart from rpoS, frequency of co-transduction of 35 %) from ENZ2005 (ENZ535 cysC::Tn10 ΔrpoS::kan) into evolved strains, selected cysC::Tn10 ΔrpoS::kan (Tcr Kanr) strains, and transduced the cysC + rpoS + region from the ancestral strain ENZ535 (Cys+ Kans). The evolved strains ENZ1901 (RpoS−) gave rise to ENZ2043 (Glg−), and the evolved strains ENZ1902 (Glg−) to ENZ2044 (Glg−). In the latter case, the strains ENZ1902 and ENZ2044 were indeed identical (rpoS +) as revealed by sequencing of the genomes.
Media and culture conditions
The MOPS minimal medium used for the pre-cultures contained notably 40 mM MOPS (pH 7.4), 86 mM NaCl, 9.5 mM NH4Cl, 20 mM glucose, and 5 mM K2HPO4. The P-limiting MOPS minimal medium contained 40 mM MOPS (pH 7.4), 86 mM NaCl, 9.8 mM KCl, 9.5 mM NH4Cl, 40 mM glucose, and 0.1 mM K2HPO4. The N-limiting MOPS minimal medium contained 40 mM MOPS (pH 7.4), 86 mM NaCl, 2 mM NH4Cl, 40 mM glucose, and 5 mM K2HPO4 .
To perform monocultures, strains grown for 24 h in MOPS medium were inoculated (0.1 ml) into 50 ml of P- or N-limiting medium in 500-ml Erlenmeyer flasks (time zero) and incubated with agitation at 150 rpm in a covered water bath rotary shaker at 37 °C. To perform mixed cultures, the strains tested as a minority were grown as monocultures in P-limiting medium for 1 day, diluted 103-fold in MOPS0 buffer (MOPS medium without P, N and C source), added (0.5 ml) into 50 ml of 1-day-old cultures of the ancestral strain in P-limiting medium, and incubated for 8 days. On day 9 of incubation, 1 ml of mixed culture was added into 50 ml of fresh P-limiting medium and incubated further for up to 30 days. In order to distinguish the different strains in mixed cultures, we transduced the lacY::Tn10 (Tcr) and ΔlacY::kan (Kanr) mutations .
Determination of RpoS, Glg and PhnE phenotypes
The RpoS− (Glg− KatE−) and Glg− (Glg− KatE+) phenotypes were determined on isolated colonies grown on LB10 medium plates containing 10 g/l NaCl as previously described [1, 74]. The PhnE+ phenotype was determined on isolated colonies grown on glucose MOPS minimal medium plates containing 0.05 mM methyl phosphonate as the sole source of P .
Measurement of cell viability
To assess cell viability, serial dilutions were prepared in M9 buffer and aliquots (10 μl) were spotted in triplicate onto LB10 medium plates, which were spread with 2000 units catalase and might contain 30 μg/ml kanamycin or 12 μg/ml tetracycline [1, 74]. In the figures, the values of 10 CFU/ml in parentheses indicate that no CFU were detected when 5 × 20-μl portions of the cultures were directly plated.
Levels of glucose and of acetic acid
The pH of the culture supernatants were determined at 25 °C, adjusted to pH 7 and the concentrations of glucose and acetate were determined by enzymatic tests (R-Biopharm) .
Measurement of β-galactosidase activity
Whole-genome sequencing was performed with genomic DNA isolated with the “Wizard Genomic DNA Purification Kit” (Promega). Genomic DNA library preparation (Insert size: ~330 bp; Tiles: 120), sequencing (“Kit v5”, Illumina) with the “Genome Analyzer GA-IIx Illumina” (74 cycles), and data analysis (CASAVA-1.8.2) were performed by Imagif (CNRS). Tablet (18.104.22.168) was used for the visualization of sequence assemblies .
Sequencing of the genes of interest in the evolved and constructed strains was performed from PCR fragments generated with the “PCR Master Mix” (Promega) from colonies. PCR products were purified using a “QIAquick” PCR purification kit (Qiagen) and sequencing was performed by Beckman Coulter Genomics.
LL carried out the molecular genetic studies. PLM conceived of the study, performed the experiments and wrote the manuscript. Both authors read and approved the final manuscript.
We thank A. Battesti for comments on the manuscript. The Centre National de la Recherche Scientifique and Aix-Marseille Université provided funding.
Both authors declare that they have no competing interests.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Guillemet ML, Moreau PL. Activation of the cryptic PhnE permease promotes rapid adaptive evolution in a population of Escherichia coli K-12 starved for phosphate. J Bacteriol. 2012;194:253–60.View ArticleGoogle Scholar
- Moreau PL. Protective role of the RpoE (σE) and Cpx envelope stress responses against gentamicin killing on nongrowing Escherichia coli incubated under aerobic, phosphate starvation conditions. FEMS Microbiol Lett. 2014;357:151–6.Google Scholar
- Hsieh Y-J, Wanner BL. Global regulation by the seven-component Pi signaling system. Curr Opin Microbiol. 2013;13:198–203.View ArticleGoogle Scholar
- Moreau PL. Diversion of the metabolic flux from pyruvate dehydrogenase to pyruvate oxidase decreases oxidative stress during glucose metabolism in nongrowing Escherichia coli cells incubated under aerobic, phosphate starvation conditions. J Bacteriol. 2004;186:7364–8.View ArticleGoogle Scholar
- Guillemet ML, Moreau PL. Fur-dependent detoxification of organic acids by rpoS mutants during prolonged incubation under aerobic, phosphate starvation conditions. J Bacteriol. 2008;190:5567–75.View ArticleGoogle Scholar
- Battesti A, Majdalani N, Gottesman S. The RpoS-mediated general stress response in Escherichia coli. Ann Rev Microbiol. 2011;65:189–213.View ArticleGoogle Scholar
- Gopalkrishnan S, Nicoloff H, Ades SE. Co-ordinated regulation of the extracytoplasmic stress factor, sigma E, with other Escherichia coli sigma factors by (p)ppGpp and DksA may be achieved by specific regulation of individual holoenzymes. Mol Microbiol. 2014;93:479–93.View ArticleGoogle Scholar
- Weber H, Polen T, Heuveling J, Wendisch VF, Hengge R. Genome-wide analysis of the general stress response network in Escherichia coli: σS-dependent genes, promoters, and sigma factor selectivity. J Bacteriol. 2005;187:1591–603.View ArticleGoogle Scholar
- Korshunov S, Imlay JA. Detection and quantification of superoxide formed within the periplasm of Escherichia coli. J Bacteriol. 2006;188:6326–34.View ArticleGoogle Scholar
- Moreau PL. The lysine decarboxylase CadA protects Escherichia coli starved of phosphate against fermentation acids. J Bacteriol. 2007;189:2249–61.View ArticleGoogle Scholar
- Cronan JE, Laporte D. Tricarboxylic acid cycle and glyoxylate bypass. EcoSal Plus. 2013. doi:10.1128/ecosalplus.3.5.2.Google Scholar
- Wolfe AJ. The acetate switch. Microb Mol Biol Rev. 2005;69:12–50.View ArticleGoogle Scholar
- Bekker M, Alexeeva S, Laan W, Sawers G, Teixeira de Mattos J, Hellingwerf K. The ArcBA two-component system of Escherichia coli is regulated by the redox state of both the ubiquinone and the menaquinone pool. J Bacteriol. 2010;192:746–54.View ArticleGoogle Scholar
- Fried L, Lassak J, Jung K. A comprehensive toolbox for the rapid construction of lacZ fusion reporters. J Microbiol Methods. 2012;91:537–43.View ArticleGoogle Scholar
- Nunoshiba T, Hidalgo E, Li Z, Demple B. Negative autoregulation by the Escherichia coli SoxS protein: a dampening mechanism for the soxRS redox stress response. J Bacteriol. 1993;175:7492–4.Google Scholar
- Baba T, Ara T, Hasegawa M, Takai Y, Okumura Y, Baba M, Datsenko KA, et al. Construction of Escherichia coli K-12 in-frame, single-gene knockout mutants: the Keio collection. Mol Syst Biol. 2006;2006(2):0008.Google Scholar
- Freddolino PL, Amini S, Tavazoie S. Newly identified genetic variations in common Escherichia coli MG1655 stock cultures. J Bacteriol. 2012;194:303–6.View ArticleGoogle Scholar
- Charusanti P, Conrad TM, Knight EM, Venkataraman K, Fong NL, Xie B, Gao Y, Palsson BØ. Genetic basis of growth adaptation of Escherichia coli after deletion of pgi, a major metabolic gene. PLoS Genet. 2010;6(11):e1001186.View ArticleGoogle Scholar
- Wang X, Kim Y, Ma Q, Hong SH, Pokusaeva K, Sturino JM, Wood TK. Cryptic prophages help bacteria cope with adverse environments. Nat Commun. 2010;1:147. doi:10.1038/ncomms1146.View ArticleGoogle Scholar
- Zinser ER, Schneider D, Blot M, Kolter R. Bacterial evolution through the selective loss of beneficial genes: trade-offs in expression involving two loci. Genetics. 2003;164:1271–7.Google Scholar
- Plucain J, Hindré T, Le Gac M, Tenaillon O, Cruveiller S, Médigue C, et al. Epistasis and allele specificity in the emergence of a stable polymorphism in Escherichia coli. Science. 2014;343:1366–9.View ArticleGoogle Scholar
- Saxer G, Krepps MD, Merkley ED, Ansong C, Deatherage Kaiser BL, Valovska M-T, et al. Mutations in global regulators lead to metabolic selection during adaptation to complex environments. PLoS Genet. 2014;10(12):e1004872.View ArticleGoogle Scholar
- Lang GI, Rice DP, Hickman MJ, Sodergren E, Weinstock GM, Botstein D, Desai MM. Pervasive genetic hitchhiking and clonal interference in forty evolving yeast populations. Nature. 2013;500:571–4.View ArticleGoogle Scholar
- Johnson MD, Bell J, Clarke K, Chandler R, Pathak P, Xia Y, et al. Characterization of mutations in the PAS domain of the EvgS sensor kinase selected by laboratory evolution for acid resistance in Escherichia coli. Mol Microbiol. 2014;93:911–27.View ArticleGoogle Scholar
- Wang L, Spira B, Zhou Z, Feng L, Maharjan RP, Li X, et al. Divergence involving global regulatory gene mutations in an Escherichia coli population evolving under phosphate limitation. Genome Biol Evol. 2010;2:478–87.View ArticleGoogle Scholar
- Herron MD, Doebeli M. Parallel evolutionary dynamics of adaptive diversification in Escherichia coli. PLoS Biol. 2013;11(2):e1001490.View ArticleGoogle Scholar
- Hottes AK, Freddolino PL, Khare A, Donnell ZN, Liu JC, Tavazoie S. Bacterial adaptation through loss of function. PLoS Genet. 2013;9(7):e1003617.View ArticleGoogle Scholar
- Maharjan R, Seeto S, Notley-McRobb L, Ferenci T. Clonal adaptive radiation in a constant environment. Science. 2006;313:514–7.View ArticleGoogle Scholar
- Tenaillon O, Rodriguez-Verdugo A, Gaut RL, McDonald P, Bennett AF, Long AD, Gaut BS. The molecular diversity of adaptive convergence. Science. 2012;335:457–61.View ArticleGoogle Scholar
- Klein G, Kobylak N, Lindner B, Stupak A, Raina S. Assembly of lipopolysaccharide in Escherichia coli requires the essential LapB heat shock protein. J Biol Chem. 2014;289:14829–53.View ArticleGoogle Scholar
- Mahalakshmi S, Sunayana MR, SaiSree L, Reddy M. yciM is an essential gene required for regulation of lipopolysaccharide synthesis in Escherichia coli. Mol Microbiol. 2014;91:145–57.View ArticleGoogle Scholar
- Altendorf K, Booth IR, Gralla J, Greie J-C, Rosenthal AZ, Wood JM. Osmotic stress. EcoSal Plus. 2013. doi:10.1128/ecosalplus.5.4.5.Google Scholar
- Lee C-R, Cho S-H, Kim H-J, Kim M, Peterkofsky A, Seok Y-J. Potassium mediates Escherichia coli enzyme IIANtr-dependent regulation of sigma factor selectivity. Mol Microbiol. 2010;78:1468–83.View ArticleGoogle Scholar
- Zafar MA, Carabetta VJ, Mandel MJ, Silhavy TJ. Transcriptional occlusion caused by overlapping promoters. Proc Natl Acad Sci USA. 2014;111:1557–61.View ArticleGoogle Scholar
- Heermann R, Weber A, Mayer B, Ott M, Hauser E, Gabriel G, et al. The universal stress protein UspC scaffolds the KdpD/KdpE signaling cascade of Escherichia coli under salt stress. J Mol Biol. 2008;386:134–48.View ArticleGoogle Scholar
- Pflüger-Grau K, Görke B. Regulatory roles of the bacterial nitrogen-related phosphotransferase system. Trends Microbiol. 2010;18:205–14.View ArticleGoogle Scholar
- Reaves ML, Rabinowitz JD. Characteristic phenotypes associated with ptsN-null mutants in Escherichia coli K-12 are absent in strains with functional ilvG. J Bacteriol. 2011;193:4576–81.View ArticleGoogle Scholar
- Laermann V, Cudic E, Kipschull K, Zimmann P, Altendorf K. The sensor kinase KdpD of Escherichia coli senses external K+. Mol Microbiol. 2013;88:1194–204.View ArticleGoogle Scholar
- Ades SE. Regulation by destruction: design of the σE envelope stress response. Curr Opin Microbiol. 2008;11:535–40.View ArticleGoogle Scholar
- Ades SE, Grigorova IL, Gross CA. Regulation of the alternative sigma factor σE during initiation, adaptation, and shutoff of the extracytoplasmic heat shock response in Escherichia coli. J Bacteriol. 2003;185:2512–9.View ArticleGoogle Scholar
- Costanzo A, Nicoloff H, Barchinger SE, Banta AB, Gourse RL, Ades SE. ppGpp and DksA likely regulate the activity of the extracytoplasmic stress factor σE in Escherichia coli by both direct and indirect mechanisms. Mol Microbiol. 2008;67:619–32.View ArticleGoogle Scholar
- Rolfe MD, Ocone A, Stapleton MR, Hall S, Trotter EW, Poole RK, et al. Systems analysis of transcription factor activities in environments with stable and dynamic oxygen concentrations. Open Biol. 2012;2:120091.View ArticleGoogle Scholar
- Park DM, Akhtar MS, Ansari AZ, Landick R, Kiley PJ. The bacterial response regulator ArcA uses a diverse binding site architecture to regulate carbon oxidation globally. PLoS Genet. 2013;9(10):e1003839.View ArticleGoogle Scholar
- Perrenoud A, Sauer U. Impact of global transcriptional regulation by ArcA, ArcB, Cra, Crp, Cya, Fnr, and Mlc on glucose catabolism in Escherichia coli. J Bacteriol. 2005;187:3171–9.View ArticleGoogle Scholar
- Koo M-S, Lee J-H, Rah S-Y, W-s Yeo, Lee J-W, Lee K-L, et al. A reducing system of the superoxide sensor SoxR in Escherichia coli. EMBO J. 2003;22:2614–22.View ArticleGoogle Scholar
- Imlay JA. Cellular defenses against superoxide and hydrogen peroxide. Ann Rev Biochem. 2008;77:755–76.View ArticleGoogle Scholar
- Kumar R, Shimizu K. Transcriptional regulation of main metabolic pathways of cyoA, cydB, fnr, and fur gene knockout Escherichia coli in C-limited and N-limited aerobic continuous cultures. Microb Cell Fact. 2011;10:3.View ArticleGoogle Scholar
- Lee PE, Demple B, Barton JK. DNA-mediated redox signaling for transcriptional activation of SoxR. Proc Natl Acad Sci USA. 2009;106:13164–8.View ArticleGoogle Scholar
- Maharjan R, Zhou Z, Ren Y, Li Y, Gaffé J, Schneider D, et al. Genomic identification of a novel mutation in hfq that provides multiple benefits in evolving glucose-limited populations of Escherichia coli. J Bacteriol. 2010;192:4517–21.View ArticleGoogle Scholar
- Lázár V, Singh GP, Spohn R, Nagy I, Horvath B, Hrtyan M, et al. Bacterial evolution of antibiotic hypersensitivity. Mol Syst Biol. 2013;9:700.View ArticleGoogle Scholar
- Taber HW, Mueller JP, Miller PF, Arrow AS. Bacterial uptake of aminoglycoside antibiotics. Microbiol Rev. 1987;51:439–57.Google Scholar
- Lüttmann D, Göpel Y, Görke B. The phosphotransferase protein EIIANtr modulates the phosphate starvation response through interaction with histidine kinase PhoR in Escherichia coli. Mol Microbiol. 2012;86:96–110.View ArticleGoogle Scholar
- Jensen RA. Evolution of metabolic pathways in enteric bacteria. In: Neidhardt FC, Curtiss III R, Ingraham JL, Lin ECC, Low KB, Magasanick B, Resnikoff WS, Riley M, Schaechter M, Umbarger HE, editors. Escherichia coli and Salmonella: cellular and molecular biology. Washington, DC: ASM Press; 1996. p. 2649–62.Google Scholar
- Kato A, Chen HD, Latifi T, Groisman EA. Reciprocal control between a bacterium’s regulatory system and the modification status of its lipopolysaccharide. Mol Cell. 2012;47:897–908.View ArticleGoogle Scholar
- Price SB, Wright JC, DeGraves FJ, Castanie-Cornet MP, Foster JW. Acid resistance systems required for survival of Escherichia coli O157:H7 in the bovine gastrointestinal tract and in apple cider are different. Appl Environ Microbiol. 2004;70:4792–9.View ArticleGoogle Scholar
- Douchin V, Bohn C, Bouloc P. Down-regulation of porins by a small RNA bypasses the essentiality of the regulated intramembrane proteolysis protease RseP in Escherichia coli. J Biol Chem. 2006;281:12253–9.View ArticleGoogle Scholar
- Hizukuri Y, Akiyama Y. PDZ domains of RseP are not essential for sequential cleavage of RseA or stress-induced σE activation in vivo. Mol Microbiol. 2012;86:1232–45.View ArticleGoogle Scholar
- Li W, Rao DK, Kaur P. Dual role of the metalloprotease FtsH in biogenesis of the DrrAB drug transporter. J Biol Chem. 2013;288:11854–64.View ArticleGoogle Scholar
- Prüss BM, Nelms JM, Park C, Wolfe AJ. Mutations in NADH: ubiquinone oxidoreductase of Escherichia coli affect growth on mixed amino acids. J Bacteriol. 1994;176:2143–50.Google Scholar
- Vemuri GN, Altman E, Sangurdekar DP, Khodursky AB, Eiteman MA. Overflow metabolism in Escherichia coli during steady-state growth: transcriptional regulation and effect of the redox ratio. Appl Environ Microbiol. 2006;72:3653–61.View ArticleGoogle Scholar
- Harden MM, He A, Creamer K, Clark MW, Hamdallah I, Martinez KA, et al. Acid-adapted strains of Escherichia coli K-12 obtained by experimental evolution. Appl Environ Microbiol. 2015;81:1932–41.View ArticleGoogle Scholar
- Goodman DB, Church GM, Kosuri S. Causes and effects of N-terminal codon bias in bacterial genes. Science. 2013;342:475–9.View ArticleGoogle Scholar
- Kim SJ, Yoon JS, Shishido H, Yang Z, Rooney LA, Barral JM, Skach WR. Translational tuning optimizes nascent protein folding in cells. Science. 2015;348:444–8.View ArticleGoogle Scholar
- Ozaki S, Matsuda Y, Keyamura K, Kawakami H, Noguchi Y, Kasho K, et al. A replicase clamp-binding dynamin-like protein promotes colocalization of nascent DNA strands and equipartitioning of chromosomes in Escherichia coli. Cell Rep. 2013;4:985–95.View ArticleGoogle Scholar
- Moreau PL. Effects of overproduction of single-stranded DNA-binding protein on RecA-dependent processes in Escherichia coli. J Mol Biol. 1987;194:621–34.View ArticleGoogle Scholar
- Moreau PL, Fanica M, Devoret R. Cleavage of λ repressor and induction of RecA protein synthesis elicited by aflatoxin B1 metabolites in Escherichia coli. Carcinogenesis. 1980;1:837–48.View ArticleGoogle Scholar
- Bron PA, van Baarlen P, Kleerebezem M. Emerging molecular insights into the interaction between probiotics and the host intestinal mucosa. Nat Rev Microbiol. 2012;10:66–78.Google Scholar
- Causey TB, Shanmugam KT, Yomano LP, Ingram LO. Engineering Escherichia coli for efficient conversion of glucose to pyruvate. Proc Natl Acad Sci USA. 2004;101:2235–40.View ArticleGoogle Scholar
- Waegeman H, Beauprez J, Moens H, Maertens J, DeMey M, Foulquié-Moreno M, et al. Effect of iclR and arcA knockouts on biomass formation and metabolic fluxes in Escherichia coli K12 and its implications on understanding the metabolism of Escherichia coli BL21 (DE3). BMC Microbiol. 2011;11:70.View ArticleGoogle Scholar
- Foo JL, Jensen HM, Dahl RH, George K, Keasling JD, Lee TS, et al. Improving microbial biogasoline production in Escherichia coli using tolerance engineering. mbio. 2014;5(6):e01932-14.View ArticleGoogle Scholar
- Mahalik S, Sharma AK, Mukherjee KJ. Genome engineering for improved recombinant protein expression in Escherichia coli. Microb Cell Fact. 2014;13:177.View ArticleGoogle Scholar
- Vital M, Chai B, Ostman B, Cole J, Konstantinidis KT, Tiedje JM. Gene expression analysis of Escherichia coli strains provides insights into the role of gene regulation in diversification. ISME J. 2015;9:1130–40.View ArticleGoogle Scholar
- Gérard F, Dri A-M, Moreau PL. Role of Escherichia coli RpoS, LexA and H-NS global regulators in metabolism and survival under aerobic, phosphate-starvation conditions. Microbiology. 1999;145:1547–62.View ArticleGoogle Scholar
- Miller JH. Experiments in molecular genetics. New York: Cold Spring Harbor Lab; 1972.Google Scholar
- Milne I, Bayer M, Cardle L, Shaw P, Stephen G, Wright F, Marshall D. Tablet-next generation sequence assembly visualization. Bioinformatics. 2010;26:401–2.View ArticleGoogle Scholar