- Open Access
Evolutionary engineering of E. coli MG1655 for tolerance against isoprenol
Biotechnology for Biofuels volume 13, Article number: 183 (2020)
Isoprenol is the basis for industrial flavor and vitamin synthesis and also a promising biofuel. Biotechnological production of isoprenol with E. coli is currently limited by the high toxicity of the final product. Adaptive laboratory evolution (ALE) is a promising method to address complex biological problems such as toxicity.
Here we applied this method successfully to evolve E. coli towards higher tolerance against isoprenol, increasing growth at the half-maximal inhibitory concentration by 47%. Whole-genome re-sequencing of strains isolated from three replicate evolutions at seven time-points identified four major target genes for isoprenol tolerance: fabF, marC, yghB, and rob. We could show that knock-out of marC and expression of mutated Rob H(48) → frameshift increased tolerance against isoprenol and butanol. RNA-sequencing showed that the deletion identified upstream of yghB correlated with a strong overexpression of the gene. The knock-out of yghB demonstrated that it was essential for isoprenol tolerance. The mutated Rob protein and yghB deletion also lead to increased vanillin tolerance.
Through ALE, novel targets for strain optimization in isoprenol production and also the production of other fuels, such as butanol, could be obtained. Their effectiveness could be shown through re-engineering. This paves the way for further optimization of E. coli for biofuel production.
3-Methyl-3-buten-1-ol, or isoprenol, is a hemiterpene belonging to the class of naturally occurring terpenoid compounds . It is the basis for the chemical synthesis of flavor compounds, such as menthol, citral, vitamin A, E, and several carotenoids . It has also been discussed for several other applications such as biofuel , as an anti-knocking additive in gasoline , and as a lead nutraceutical for longevity .
The biotechnological production of terpenoid compounds in microorganisms relies on the natural precursor isopentenyl diphosphate (IPP) from which isoprenol can be obtained by simple dephosphorylation. So far, strain engineering has focused on increasing the intracellular concentration of IPP. In Escherichia coli this has been achieved by introducing an additional metabolic pathway that produces IPP, the DXP pathway, resulting in a product titer of 61 mg/L . However, the intermediate IPP has been identified as a major obstacle in those processes due to its high toxicity . Using a specially designed pathway that by-passes the intermediate IPP, a product titer of 3.7 g/L was achieved. Applying the resulting strain in a fed-batch two-phase system, the highest published titer of 10.8 g/L isoprenol could be obtained . Present research projects try to develop integrated processes where isopentenol (a mixture of prenol and isoprenol) is obtained from hydrolyzed polysaccharides originating from biomass .
Product toxicity towards microorganisms is a key issue that most economically viable bioprocesses have to face. This is also true for the biotechnological production of terpenoids and biofuels [10, 11]. Kang et al. found that indeed also isoprenol production was limited by product toxicity. A way to alleviate this toxicity is process engineering, for instance through the introduction of in situ product extraction into an organic phase. Isoprenol toxicity could be alleviated by using such a two-phase system , an approach that also works for other isoprenoids . But two-phase systems rely on the partitioning coefficient of the target molecule between the aqueous and the organic phase. This means that the hydrophobic parts of the cell (e.g., the membrane) will most likely still be saturated with the target product and potentially the extractant, affecting structure and function. Improving the inherent resistance of the microorganism towards the target product through metabolic engineering is hence complementary to introducing in situ extraction.
Improving tolerance through metabolic engineering is possible through overexpression of export pumps [13, 14], but the prediction of engineering targets not related to transport is difficult. An approach to overcome this problem is to use the power of evolution to come up with solutions to alleviate toxicity, this can best be done by using adaptive laboratory evolution (ALE) [10, 11, 15].
In this study, we used ALE to adapt E. coli MG 1655 against isoprenol. We isolated highly tolerant strains that showed up to a 47% increase in growth at 50 mM isoprenol (4.3 g/L) and could grow in the presence of 80 mM isoprenol (6.8 g/L), a concentration at which the parental strain ceased to grow. With a combination of DNA and RNA-sequencing, we identified 4 target mutations for high tolerance. To confirm the role of the mutations, we re-engineered those into the parental strain and quantified the effects on isoprenol tolerance.
Experimental evolution results in increased isoprenol tolerance
To design the adaptive evolution experiment with the appropriate isoprenol stress, we first quantified its toxic effect on E. coli K12 MG1655 growth. We experimentally determined the half-maximal inhibitory concentration IC50 to be 53 mM (Fig. 1a). The growth rate decreased with increasing isoprenol concentration and there was no growth at the highest tested concentration of 80 mM.
We evolved E. coli by serial passaging of three independent populations for approximately 200 generations on isoprenol-spiked M9 minimal medium supplemented with 5 g/L glucose. The initial isoprenol concentration was set to 60 mM (5.2 g/L) corresponding to approximately 70% growth reduction in all three replicates. This meant that cells could grow for 24 h without reaching the stationary phase at the chosen seed-densities. During the course of the evolution, samples were taken and conserved with 50% v/v glycerol at − 80 °C. The isoprenol concentration was gradually increased to approximating a constant growth rate and reaching 80 mM after approximately 100 generations (Fig. 1a "Isolate B 100 gen"). At this time-point we already obtained tolerant strains, able to grow at isoprenol concentrations were no growth in the parental strain was observed (Fig. 1a). However, an attempt to further raise the isoprenol concentration to 90 mM after 140 generations resulted in the extinction of the cultures. The evolution was restarted from the prior cryo-stocks and from then onwards 80 mM isoprenol was used. Finally, at seven time-points during the evolution strains were isolated by selecting the five largest colonies on isoprenol-spiked LB-agar plates. The isoprenol amounts depended on the respective concentration in the evolution experiments at the time the samples were taken and were 64, 72, and 80 mM, respectively (Additional file 1: Table S1). Of those five largest colonies the fastest growing strain was determined in a falcon tube experiment with isoprenol (60 mM) (Additional file 1: Table S2) and used for genome re-sequencing.
An exact determination of the fitness benefit was performed for three final isolates of the biological replicates (Fig. 1b). At the half-maximal toxic concentration of 50 mM isoprenol, we observed an increase of up to 47% in growth rate compared to the initial MG1655. The isolated strains did not differ in their growth behavior in medium without isoprenol (Fig. 1a). Contrary to the initial strain, the final isolates could grow on LB-agar plates supplemented with 80 mM isoprenol. In summary, we successfully obtained E. coli strains that showed an increased tolerance against isoprenol using ALE.
Identification of mutations by whole-genome re-sequencing
To identify the genetic basis for the increased tolerance of the adapted strain and to uncover the temporal occurrence of target mutations, genome re-sequencing was applied to the three parallel evolution experiments. There are two important points to be raised before looking at the actual results. On the one hand, the commercial isoprenol used for the evolution had 3% impurities. We found that formaldehyde was present as a major impurity, which most likely caused the mutations in the formaldehyde repressor frmR. We tested the final isolates against further purified isoprenol to ensure the observed phenotypes were also stable without formaldehyde. We could not observe differences in their fitness patterns. In addition, the re-sequencing revealed two strains from two parallel evolution cultures at a single time-point (culture A and C at T6) that carry the same set of four mutations. Around this time-point, all three clones isolated from the three replicates showed a great deal of exact mutational convergence, whereas most mutations previously seen in these replicates, besides yghB, were no longer detected. This comes after the evolution cultures had to be restarted from cryo-stocks due to extinction (see above). While this could be coincidental, it could also point to cross-contamination from another evolution. In this context, it is important to emphasize that we only sequenced a single clone per time-point and do not know the genetic heterogeneity of the cultures, but colony size on isoprenol-spiked agar plates varied (data not shown). We analyzed this computationally and were able to show that a simple model can produce our observed mutations without the assumption of cross-contamination. Adding the assumption of 5% cross-contamination modestly increases the log likelihood from − 6.1 in the model without contamination to − 2.8 in a model with contamination after T5 (which is the highest likelihood of different contamination scenarios). However, a likelihood ratio test reveals, that the differences are not statistically significant (Chi2-test with df = 10, p-value = 0.78). On the basis of the available data, we cannot conclude whether a cross-contamination occurred (see Additional file 1: Additional text and Tables S3–S5).
The parental MG1655, the final isolates, and isolates from the intermittent time-points of the evolution course were sequenced. All mutations identified in the re-sequencing of the isolated E. coli strains are listed below (Table 1) with their time-point of occurrence and their respective frequency in the whole sample of sequenced strains. In an analysis of the GO terms linked to the mutated genes, we found that target genes are associated with the membrane and fatty acid or phospholipid biosynthesis. But there were also genes acting as transcriptional regulators (Additional file 1: Table S6).
Occurrence and persistence of mutations in the evolution experiment
To visualize the evolution of the three cultures and identify persistent mutations, we plotted the mutations in the different cultures in their temporal order in Fig. 2. In the isolates at the first time-point approximately after 32 generations, we could not identify any mutations. In the second isolates, the first mutations occur, but none of those persist. In the isolates after 108 generations the FabF F(74) → C and MarC M(35) → stop mutation appear and in culture B persist until the end of the evolution experiment. At the next time-point T4, a mutation upstream of yghB appears that is also present in the final isolates. Finally, the first mutation in the rob gene appears at T5, in culture C. After T6 all isolated strains carry the fabF mutation, the marC mutation, and the yghB promoter mutation. Also, all strains isolated at T6 and T7 have one of two mutations of the rob gene. In two of three strains isolated at the final time-point, there are also two distinct mutations in the frmR gene. All other mutations that occurred in the intermediate isolates were not found in the final isolates, i.e., they did not persist. Interestingly in culture A, we isolated apparently the same genotype at T4 and T5 twice carrying mutations in the genes gltA, plsX, marC, and rraA. Of those target genes, marC is also found in the final genotype, also distinct plsX and rraA mutation can be identified in other isolates at other time-points making them additional targets for the tolerance phenotype.
Transcriptional response of adapted strains
To understand how the genetic adaptations influence the cell's phenotype, we performed an RNA-Seq experiment with the parental strain and strains isolated at T7. All strains were put under isoprenol stress (50 mM) and subsequently, mRNA levels were quantified. We were primarily interested in the transcriptome changes that are consistent among all three mutant strains since the final isolates share three exact mutations and one target gene.
The top ten up- and down-regulated transcripts compared to the parental strain are listed below (Table 2). Similar to the target genes of the mutations, a gene ontology analysis revealed that highly differentially regulated transcripts were associated with the cell membrane and integral membrane components (Additional file 1: Table S7). Mutations in three genes were observed, that could directly affect the regulation of gene expression: FrmR is a transcriptional repressor, the mutation at the upstream of yghB could influence its expression and finally, Rob is a global transcriptional regulator. For the yghB gene, we found in all samples a strong up-regulation in the isolated strains compared to the parental strain. For both strains with a mutation in the frmR gene, we observed an up-regulation of the frmRAB operon compared to the parental strain (data not shown). Most of the strongly differentially regulated transcripts were not directly linked to genetic mutations found in the isolates. For instance, a very strong down-regulation was also observed for hslU (Table 2). HslU is part of the protease system responsible for the degradation of the Arc repressor, usually induced by heat-shock.
Characterization of mutations
Finally, we wanted to analyze the contribution of each highly persistent target mutation to the tolerance phenotype, except for the highly prevalent fabF mutation. This has already been discovered and analyzed previously in a mutagenesis study for butanol tolerance . Because of the similarity of isoprenol and butanol in terms of structure and polarity, we did not test the fabF mutation under our conditions.
The mutation found upstream of the yghB gene deleted a 15-bp portion that overlaps with the − 35 region of the promoter. Interestingly there appears to be a 9-bp DNA motif that lies in the deletion region and is repeated inverted after the − 35 region (see Fig. 3a). Those repeated regions might be the binding site of an unknown transcriptional regulator or might form a secondary structure. In the case of an activator, deletion would reduce yghB expression in case of a repressor the opposite would be the case. In fact, the RNA-Seq data showed that yghB was 14-fold upregulated (Additional file 1: Figure S1) compared to the parental strain suggesting the deletion of a repressor binding site.
Knock-out of yghB (strain ΔyghB) clearly increased the sensitivity against isoprenol (Fig. 3b), while introducing a yghB carrying plasmid into the knock-out strain (strain ΔyghB yghB) complemented the effect. Due to the leakiness of the promotor on the high-copy number plasmid, induction was not even necessary.
The Rob protein is a transcriptional activator that regulates together with SoxRS and MarA a stress response regulon . In the evolution experiment, we identified three variant mutations G(273) → stop, Y(103) → stop and H(48) → frameshift which appeared in two final strains and in total with the highest frequency. The fact that with increasing frequency of the rob mutations less of the protein should remain functional led us to hypothesize that the fitness benefit originated from a loss of function. Consequently, we tested the fitness of a rob knock-out (strain Δrob) with isoprenol stress, but found no increase in fitness (Fig. 4 light gray). The mutated Rob H(48) → fs protein might retain some of its N-terminal DNA-binding domain and thereby still have a function as a transcription factor with a different or no sensory response since its C-terminal domain cannot have its original function due to the truncation. Plasmid-based introduction of the mutated RobH(48) → fs protein (strain Δrob robH) increased the fitness approx. 14%. At a low induction level of 10 µM isopropyl-β-d-thiogalactopyranoside (IPTG) there was still a positive fitness effect; the fitness of the strain became negative at high induction levels of 100 µM IPTG. Since the knock-out of rob did not have a fitness benefit and expression of the mutant protein did, it appeared that the mutated version had an altered function in the cell that benefitted isoprenol tolerance. The RNA-sequencing analysis showed that indeed only 3 out of 21 genes of the rob-regulon (acrZ, yhbW, and inaA) were significantly down-regulated in all mutant strains (Additional file 1: Figure S2), while other genes of the regulon remained unchanged or even slightly upregulated. AcrZ is a membrane protein associated with an efflux pump . Rob-controlled expression of inaA was observed during dipyridyl stress  and yhbW has been predicted to be a monooxygenase  and was shown to be strongly upregulated during ALE experiments on vanillin but without activity against vanillin .
MarC is a transmembrane protein  and a positive effect of gene deletion has been found in isobutanol tolerance . Minty et al. found a disruption of the marC locus by a transposon insertion and hypothesized that this insertion would have a positive effect. Indeed a knock-out of marC (strain ΔmarC) resulted in increased isobutanol tolerance. In our experimental setup knock-out of marC also had a positive effect on isoprenol tolerance (Fig. 4). Since the highest frequency mutation of marC introduced a stop codon after the 35th methionine residue, we wondered whether the protein fragment had an additional fitness benefit. Expression of the MarC M(35) → stop mutant in the marC knock-out background (strain ΔmarC marC35) did not significantly change the fitness compared to the knock-out (Fig. 4) and appeared to only lead to a loss of MarC function.
The extent of the tolerance mechanism to additional compounds
In our analysis, we found that some of the target genes were identified in evolution experiments against C-4 alcohols. Therefore, we tested the novel mutations for their fitness benefit against butanol (Fig. 5). Knock-out of yghB resulted in a slight fitness decrease. Analogously to isobutanol tolerance a marC knock-out also increased tolerance against butanol. In contrast to isoprenol tolerance, knocking out the regulator rob already had a strong fitness benefit. Introducing the plasmid expressing the rob H(48) mutant only raised fitness slightly. For butanol tolerance, the mutations found in the isoprenol evolution behaved mostly similar, except for rob where a knock-out did not improve tolerance against isoprenol.
Finally, we wanted to test a compound with different physicochemical properties and evaluate if still the same mutations benefit tolerance. In earlier experiments, we observed no improved tolerance (data not shown) of our evolved strains towards geraniol (logP 2.5, solubility in water 0.68 g/L), another terpenoid of commercial relevance. We chose vanillin (logP 1.37, solubility in water 10 g/L) as a compound with a logP value between geraniol and isoprenol to test our target mutations (Fig. 5). We found that the yghB gene still played a significant role for tolerance, however, its role was the opposite compared to isoprenol and butanol tolerance since the knock-out had a positive fitness effect. Knock-outs of marC and rob did not increase fitness strongly. But a strong fitness benefit could be obtained by expression of the mutant rob H(48) protein, a similar effect as observed on isoprenol tolerance.
Using an ALE approach, we were able to identify strains that are highly tolerant against isoprenol compared to the parental E. coli MG1655 strain. At the IC50 the initial MG1655 strain was 50% inhibited, while in the isolated strains exhibited growth at 75% of the unchallenged maximum. At the same time, these strains were capable of growing at 80 mM isoprenol a concentration completely inhibiting the parental MG1655. By re-sequencing multiple strains isolated throughout the evolution, we could identify approximate time-points of the appearance of mutation and calculate their persistence. Performing the same analysis on three biological replicates further increased our ability to judge the importance of certain mutations. The combination of genome re-sequencing with RNA-sequencing made it possible to determine the regulatory effects of mutations that are in the intergenic region.
The temporal order of those key mutations could hint at strong epistatic effects, i.e., that fitness benefits of a late mutation such as PyghB are much stronger in a marC background. In addition, the temporal order could be the consequence of increasing isoprenol concentration and thereby increased stress during evolution. If this would be the case, late mutations would exhibit a much larger fitness increase under higher isoprenol concentrations. However, it was beyond the scope of this study to test epistatic and stressor concentration-dependent fitness effects systematically. We find that at the half-maximal isoprenol concentration of 50 mM, the mutations tested in this study (PyghB, ΔmarC, and Rob H(48) → fs) would have an additive relative fitness that overlaps with the fitness benefit of the final isolates that have the additional fabF mutation. This is an indication that at this isoprenol concentration epistatic interactions must be negative, otherwise, at least one mutation would not have persisted in the evolution [23, 24].
We observed no mutations at the first time-point of isolation (32 generations) although our growth data indicated already an improvement of growth on isoprenol (Additional file 1: Table S1). This means that in addition to genetic adaptation, a phenotypic or possibly epigenetic adaptation occurred.
Contribution of single mutations to tolerance
Two of the target genes in the set of highly persistent mutations have been identified in previous studies. The MarC membrane protein had initially been assumed to play a role in multiple antibiotic resistance . In an adaptive evolution against isobutanol, Minty et al. discovered an insertion element in the marC gene and confirmed that a knock-out of marC provided increased tolerance against isobutanol. Similarly, Atsumi et al. found the deletion of the complete mar operon to provide isobutanol tolerance . We showed in our study that knock-out of marC also increases tolerance against isoprenol and expression of the mutated MarC M(35) → stop does not have a strong additional fitness effect.
A mutation in the fabF gene resulting in an amino acid exchange F(74) → C occurred with a high frequency in our re-sequencing dataset. The identical mutation has previously been identified in a mutagenesis study of E. coli for butanol stress . In a subsequent study, the investigators showed that the amino acid change led to a change in fatty acid composition increasing the major unsaturated fatty acid, cis-vaccenic acid .
Following the fabF and the marC mutation, the deletion upstream of yghB occurred with the second-highest frequency. YghB belongs to the DedA protein family of membrane proteins and is necessary for temperature tolerance . It was speculated that yghB and its family members are proton-dependent transporters . The analysis of the deleted sequence and its context revealed that upstream of yghB an unknown 9 bp sequence motif is present that is repeated as an inverted repeat. Interestingly, part of the motif overlaps with the -35 region possibly obstructing binding of the sigma-factor; the motif could be recognized by a repressor protein or inhibit transcription by cruciform extrusion . Knock-out of yghB decreases isoprenol tolerance showing its relevance for the tolerance mechanism. The knock-out could be complemented with yghB on an inducible plasmid. The yghB promoter mutation also showed significant effects on downstream genes. In fact, the RNA-seq data showed that especially yqhC and yqhD were also significantly upregulated (Additional file 1: Figure S2).
The third most abundant gene target rob belongs to the AraC/XylS protein family of transcription factors . It has been implicated in organic solvent tolerance [31, 32] and together with MarA and SoxS it regulates a large regulon . Interestingly, complementary tolerance engineering and functional genomic approaches yielded gene targets such as soxS and marA [33, 34] that share a regulon with rob via the marbox DNA motif . Rob consists of two major domains, an N-terminal DNA-binding domain comprised two HTH motifs and a C-terminal domain that binds to its signal molecules and is able to form protein complexes. The G(273) → stop mutation shortened the C-terminal domain, Y(103) → stop truncated it completely. These mutations might lead to a signal independent activation of the Rob-regulon . The most prevalent mutation, H(48) → fs, however, resulted in an altered protein after the first of the two HTH motifs of the N-terminal DNA-binding domain . Structural data shows that both HTH motifs of Rob bind to DNA , the Rob H(48) → fs might therefore not be functional. However, we could show that a rob deletion mutant did not have a fitness benefit, instead a plasmid-based expression of the mutated version of rob lead to a significant fitness increase. rob has been shown to trigger the expression of the AcrAB multidrug efflux pump and confer tolerance to multiple substances . However, our RNA-sequencing analysis showed that acrAB expression was not significantly changed in the evolved strains and that the multidrug efflux pump accessory protein gene acrZ was even down-regulated (Additional file 1: Figure S2). The marR regulator is also part of the rob operon (Additional file 1: Figure S2) and was slightly downregulated (FC around − 0.8). This did, however, not lead to a significant change in the expression of marABC (data not shown). How exactly this mutant Rob protein binds to DNA and how the altered structure results in an altered function remains unclear and needs to be addressed in future studies.
Both rraA and plsX are targets of two distinct mutations and occur in conjunction in an intermediate genotype of culture A. rraA is an inhibitor of RNase E and globally influences RNA stability . The more prevalent mutation V(96) → E introduced an amino acid change in the highly conserved V(96) residue  and possibly had a deleterious effect. Such a mutation would consequently globally increase mRNA levels. Possibly isoprenol also targets the intracellular translation process which could be compensated for by increasing mRNA levels. Interestingly, rraA was the transcript that was downregulated the strongest (up to 2–14-fold) in the isolated final mutants compared to the initial strain.
The fatty acid/phospholipid synthesis protein PlsX has also been identified as a target in an evolution study against isobutanol . In our study, two mutations that locate to the c-terminal α-helices were identified . The C-terminal parts of plsX are rather conserved and important for dimerization , amino acid changes might be deleterious or change the substrate specificity thereby changing the lipid composition in the cell membrane.
Tolerance against additional compounds
The mutations investigated in this study could also be used to engineer tolerance against butanol stress. Especially, yghB overexpression and expression of mutant rob lead to significant fitness increase. The fragrance compound vanillin appeared to have different toxic effects on MG1655, here surprisingly deletion of yghB increased the growth rate by approximately 20%. A knock-out of marC did not increase fitness, however, the mutant rob had a strong fitness effect. We found that similar gene targets play a role in tolerance against structurally different compounds such as vanillin although in a qualitatively different manner.
ALE is a powerful tool to be used alongside strain development, especially for the development of chemical resistance. It has been shown previously that the toxicity of isoprenol is the current limitation for high-yield production . With our experimental evolution approach, we could identify four target mutations that could be implemented in isoprenol producing E. coli strains in the future to alleviate the toxic effects of the product. The same mutations also showed effects in E. coli during challenges with butanol or vanillin, albeit to a different extent.
Materials and methods
Chemicals and reagents
All chemicals used in this study were at least of analytical grade and purchased from Sigma-Aldrich. Aqueous stock solutions were made using ultra-pure water (resistance > 18 MΩ). Isoprenol was also purchased from Sigma with 97% purity. BASF identified formaldehyde as an impurity. To ensure strains were not selected on formaldehyde resistance only, the final isolates were also challenged with isoprenol depleted of formaldehyde (custom chemical provided by BASF SE). No effect on resistance could be observed.
The E. coli MG1655 (CGSC6300) strain  used in the evolution (GenBank: GCA_000005845.2) exhibited a reconstituted gatC gene, a functional glrR glycerol 3-phosphate repressor  and variation in the repeat REP321j. The following strains were used and constructed in this work (Tables 3 and 4):
For plate cultivation, E. coli MG1655 was streaked out on M9 medium supplemented with 5 g/L Glucose and grown at 37 °C. When colony formation was observed, 10 mL M9 medium was inoculated with a colony and incubated in a 100 mL baffled flask in a shaking incubator (Multitron or Ecotron, both 25 mm shaking throw at 200 rpm; Infors HT, Bottmingen, Switzerland,). From this culture, an over-night culture of 25 mL medium in a 250 mL baffled flask was inoculated and incubated for 16 h. Inoculation volume was calculated such that the finished culture could be harvested the next morning in the mid-exponential phase. This was used to inoculate 25 mL of medium in 250 mL baffled flasks with Teflon® liner screw cabs to an optical density (OD) of 0.2 and incubated in a shaking incubator as described above. The toxic compound was added to the specified concentration (for isoprenol see below). To characterize the relative fitness of reconstructed strains, growth experiments with butanol (5 g/L) and vanillin (1.5 g/L) were performed separately. OD was measured every hour at 600 nm against water and the growth rate was calculated by a linear fit to log-transformed OD-values. Each experiment was performed in triplicates.
Adaptive evolution was carried out in biological triplicates. 25 mL of M9 medium supplemented with 5 g/l glucose were inoculated with E. coli MG1655. Before the cell-culture reached the stationary phase, part of the culture was transferred to fresh medium in a fresh flask with a toxic compound. The mean growth rate of the culture was determined by comparing the initial OD and the culture OD before each passaging. Upon passaging the cell culture, 600 µL was withdrawn and mixed with 600 µL 50% v/v glycerol solution. The samples were stored at – 80 °C.
The initial isoprenol concentration was 60 mM (selected after an initial toxicity screen) and with increasing mean growth rates the isoprenol concentration was increased stepwise to 80 mM isoprenol after approx. 80 generations where it was kept until the end of the experiment.
For easier comparison, the growth rate of evolved or reconstructed strains was divided by the growth rate of the initial MG1655 under the tested conditions and termed ‘relative fitness’.
Selected strains were grown overnight in 5 mL lysogeny broth (LB) medium supplemented with 60 mM isoprenol (for mutant strains). Genomic DNA was isolated by LGC Genomics GmbH (Berlin, Germany) using the DNeasy® UltraClean® Microbial Kit (Qiagen, Netherlands). Genomic DNA was fragmented using Covaris (300 bp) (Covaris, MA, USA) and samples were subsequently purified using MinElute columns (Qiagen, Netherlands). Libraries were prepared using Ovation Rapid DR Multiplex System 1–96 (Tecan, Switzerland) and amplified for 13 cycles using MyTaq (Bioline, UK) and standard Illumina primers. Size selection was done on the Pippin Prep system (Sage Science, MA, USA) selecting a range between 300 and 500 bp. Final library purification and quality control of DNA libraries were carried out on BioAnalyzer (Agilent, CA, USA) and Qubit (Thermo Fisher Scientific, MA, USA). Sequencing was done on an Illumina NextSeq 500/550 (Illumina, CA, USA) with 2 × 150 read length following manufacturer’s instructions, achieving about 100× coverage.
Libraries were demultiplexed using Illumina bcl2fastq 184.108.40.206 software  followed by clipping of adapters and trimming to achieve a minimum average Phred quality score of 20 over a window of ten bases. Alignment against E. coli MG1655 genome was carried out using BWA-MEM version 0.7.12  and variant discovery and genotyping were performed using Freebayes v1.0.2-16 .
The parental MG1655 strain and the three final mutant strains were grown in biological triplicates with 50 mM isoprenol until an OD of 1.0 as described above (25 mL sealed flasks). Then 10 mL of cell-culture was vacuum filtered using a Supor® 800 Grid filter (Pall, NY, USA) with 0.8 µM pore size. The filter containing the cells was put in a 15 mL falcon tube containing 700 µL PGTX solution and immediately frozen in liquid nitrogen. The samples were stored at − 80 °C until further processing.
RNA was extracted using standard methods  and RNA quality was determined using Bioanalyzer (Agilent, CA, USA). Ribo-Zero rRNA Removal Kit for Bacteria (Illumina, CA, USA) was used to deplete rRNA. First-strand cDNA synthesis and second-strand synthesis was carried out using NEBNext RNA First Stand and Second Strand Synthesis Module (New England Biolabs, MA, USA). cDNA was then purified and concentrated using MiniElute Columns (Qiagen, Netherlands) and the Encore Rapid DR Multiplex System (Tecan, Switzerland) was used for library preparation. Libraries were amplified for 12 cycles using MyTaq (Bioline, UK) and standard Illumina primers. Size selection was done using a preparative Agarose Gel selecting fragments between 300 and 500 bp and quality control of libraries was performed with Bioanalyzer and Qubit. Finally, sequencing was carried out on Illumina NextSeq500/550 (Illumina, CA, USA) with 1 × 75 bp read length following the manufacturer’s instructions.
The sequencing data were demultiplexed using Illumina bcl2fastq 220.127.116.11 , adapters were clipped from the reads, and rRNA was filtered using RiboPicker 0.4.3 . The sequences were aligned using STAR 2.4 and rRNA or tRNA reads were filtered . TopHat-aligned reads were counted using htseq-count . Differential expression was determined using edgeR 3.2.3 , DESeq 1.12.0  and cuffdiff 2.1.1 . Raw p-values from the statistical test were adjusted for multiple testing by the Benjamini–Hochberg false discovery rate method .
Expression plasmids and knock-out strains
All plasmids and primer sequences used in this study are provided in the supplementary information (Additional file 1: Tables S8 and S9). Knock-out strains were constructed by amplification of resistance cassette with 25-bp overlap from corresponding Keio strains  (Primers 3 + 4, 5 + 6, and 7 + 8). The PCR products carrying a homologous 25-bp sequence and a kanamycin resistance were used to transform E. coli MG1655 using standard procedures . For over-expression plasmids, target genes were amplified with a 25-bp homology to the pAH030 overexpression plasmid. The plasmid was linearized using the SpeI restriction site and the PCR-product containing the gene of interest was inserted using Gibson assembly . Recombination was carried out using a standard RED/ET kit .
Availability of data and materials
The evolved strains can be obtained for non-commercial research purposes upon the signature of a materials transfer agreement. Genome re-sequencing data (PRJNA665813) and RNA-sequencing data (GSE158959) has been deposited to the Sequence Read Archive (SRA).
Adaptive laboratory evolution
Octanol–water partitioning coefficient
Mineral minimal medium
Polymerase chain reaction
Withers ST, Keasling JD. Biosynthesis and engineering of isoprenoid small molecules. Appl Microbiol Biotechnol. 2007;73:980–90.
Parker GL, Smith LK, Baxendale IR. Development of the industrial synthesis of vitamin A. Tetrahedron. 2016;72:1645–52.
Beller HR, Lee TS, Katz L. Natural products as biofuels and bio-based chemicals: fatty acids and isoprenoids. Nat Prod Rep. 2015;32:1508–26.
Mack JH, Rapp VH, Broeckelmann M, Lee TS, Dibble RW. Investigation of biofuels from microorganism metabolism for use as anti-knock additives. Fuel. 2014;117:939–43.
Pandey S, Phulara SC, Jha A, Chauhan PS, Gupta P, Shukla V. 3-Methyl-3-buten-1-ol (isoprenol) confers longevity and stress tolerance in Caenorhabditis elegans. Int J Food Sci Nutr. 2019;70:1–8.
Liu H, Wang Y, Tang Q, Kong W, Chung W-J, Lu T. MEP pathway-mediated isopentenol production in metabolically engineered Escherichia coli. Microb Cell Fact. 2014;13:135.
George KW, Thompson M, Kim J, Baidoo EEK, Wang G, Benites VT, et al. Integrated analysis of isopentenyl pyrophosphate (IPP) toxicity in isoprenoid-producing Escherichia coli. Metab Eng. 2018;47:60–72.
Kang A, Mendez-Perez D, Goh E-B, Baidoo EEK, Benites VT, Beller HR, et al. Optimization of the IPP-bypass mevalonate pathway and fed-batch fermentation for the production of isoprenol in Escherichia coli. Metab Eng. 2019;56:85–96.
Wang S, Cheng G, Dong J, Tian T, Lee TS, Mukhopadhyay A, et al. NaCl enhances Escherichia coli growth and isoprenol production in the presence of imidazolium-based ionic liquids. Bioresour Technol Rep. 2019;6:1–5.
Brennan TCR, Williams TC, Schulz BL, Palfreyman RW, Krömer JO, Nielsen LK. Evolutionary engineering improves tolerance for replacement jet fuels in Saccharomyces cerevisiae. Appl Environ Microbiol. 2015;81:3316–25.
Minty JJ, Lesnefsky AA, Lin F, Chen Y, Zaroff TA, Veloso AB, et al. Evolution combined with genomic study elucidates genetic bases of isobutanol tolerance in Escherichia coli. Microb Cell Fact. 2011;10:18.
Brennan TCR, Turner CD, Krömer JO, Nielsen LK. Alleviating monoterpene toxicity using a two-phase extractive fermentation for the bioproduction of jet fuel mixtures in Saccharomyces cerevisiae. Biotechnol Bioeng. 2012;109:2513–22.
Dunlop MJ, Dossani ZY, Szmidt HL, Chu HC, Lee TS, Keasling JD, et al. Engineering microbial biofuel tolerance and export using efflux pumps. Mol Syst Biol. 2011;7:487.
Basler G, Thompson M, Tullman-Ercek D, Keasling J. A Pseudomonas putida efflux pump acts on short-chain alcohols. Biotechnol Biofuels. 2018;11:136.
Dragosits M, Mattanovich D. Adaptive laboratory evolution—principles and applications for biotechnology. Microb Cell Fact. 2013;12:64.
Haeyoung J, Jihee H. Enhancing 1-butanol tolerance in Escherichia coli through repetitive proton beam irradiation. J Korean Phys Soc. 2010;56:2041.
Griffith KL, Fitzpatrick MM, Keen EF, Wolf RE. Two functions of the C-terminal domain of Escherichia coli Rob: mediating “sequestration-dispersal” as a novel off-on switch for regulating Rob’s activity as a transcription activator and preventing degradation of Rob by Lon protease. J Mol Biol. 2009;388:415–30.
Hobbs EC, Yin X, Paul BJ, Astarita JL, Storz G. Conserved small protein associates with the multidrug efflux pump AcrB and differentially affects antibiotic resistance. Proc Natl Acad Sci USA. 2012;109:16696–701.
Rosner JL, Dangi B, Gronenborn AM, Martin RG. Posttranscriptional activation of the transcriptional activator Rob by dipyridyl in Escherichia coli. J Bacteriol. 2002;184:1407–16.
Hayden JD, Ades SE. The extracytoplasmic stress factor, σE, is required to maintain cell envelope integrity in Escherichia coli. PLoS ONE. 2008;3:e1573.
Pattrick CA, Webb JP, Green J, Chaudhuri RR, Collins MO, Kelly DJ. Proteomic profiling, transcription factor modeling, and genomics of evolved tolerant strains elucidate mechanisms of vanillin toxicity in Escherichia coli. mSystems. 2019;4(4):e00163-19.
McDermott PF, McMurry LM, Podglajen I, Dzink-Fox JL, Schneiders T, Draper MP, et al. The marC gene of Escherichia coli is not involved in multiple antibiotic resistance. Antimicrob Agents Chemother. 2008;52:382–3.
Khan AI, Dinh DM, Schneider D, Lenski RE, Cooper TF. Negative epistasis between beneficial mutations in an evolving bacterial population. Science. 2011;332:1193–6.
Chou HH, Chiu HC, Delaney NF, Segrè D, Marx CJ. Diminishing returns epistasis among beneficial mutations decelerates adaptation. Science. 2011;332:1190–2.
Atsumi S, Wu TY, MacHado IMP, Huang WC, Chen PY, Pellegrini M, et al. Evolution, genomic analysis, and reconstruction of isobutanol tolerance in Escherichia coli. Mol Syst Biol. 2010;6:1–11.
Jeong H, Kim SH, Han SS, Kim MH, Lee KC. Changes in membrane fatty acid composition through proton-induced fabF mutation enhancing 1-butanol tolerance in E. coli. J Korean Phys Soc. 2012;61:227–33.
Thompkins K, Chattopadhyay B, Xiao Y, Henk MC, Doerrler WT. Temperature sensitivity and cell division defects in an Escherichia coli strain with mutations in yghB and yqjA, encoding related and conserved inner membrane proteins. J Bacteriol. 2008;190:4489–500.
Kumar S, Doerrler WT. Members of the conserved DedA family are likely membrane transporters and are required for drug resistance in Escherichia coli. Antimicrob Agents Chemother. 2014;58:923–30.
Bikard D, Loot C, Baharoglu Z, Mazel D. Folded DNA in action: hairpin formation and biological functions in prokaryotes. Microbiol Mol Biol Rev. 2010;74:570–88.
Kwon HJ, Bennik MHJ, Demple B, Ellenberger T. Crystal structure of the Escherichia coli Rob transcription factor in complex with DNA. Nat Struct Biol. 2000;7:424–30.
Nakajima H, Kobayashi K, Kobayashi M, Asako H, Aono R. Overexpression of the robA gene increases organic solvent tolerance and multiple antibiotic and heavy metal ion resistance in Escherichia coli. Appl Environ Microbiol. 1995;61:2302–7.
Aono R. Improvement of organic solvent tolerance level of Escherichia coli by overexpression of stress-responsive genes. Extremophiles. 1998;2:239–48.
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:e01932.
Rutherford BJ, Dahl RH, Price RE, Szmidt HL, Benke PI, Mukhopadhyay A, et al. Functional genomic study of exogenous n-butanol stress in Escherichia coli. Appl Environ Microbiol. 2010;76:1935–45.
Martin RG, Bartlett ES, Rosner JL, Wall ME. Activation of the Escherichia coli marA/soxS/rob regulon in response to transcriptional activator concentration. J Mol Biol. 2008;380:278–84.
Rosenberg EY, Bertenthal D, Nilles ML, Bertrand KP, Nikaido H. Bile salts and fatty acids induce the expression of Escherichia coli AcrAB multidrug efflux pump through their interaction with Rob regulatory protein. Mol Microbiol. 2003;48:1609–19.
Lee K, Zhan X, Gao J, Qiu J, Feng Y, Meganathan R, et al. RraA: a protein inhibitor of RNase E activity that globally modulates RNA abundance in E. coli. Cell. 2003;114:623–34.
Monzingo AF, Gao J, Qiu J, Georgiou G, Robertus JD. The X-ray structure of Escherichia coli RraA (MenG), a protein inhibitor of RNA processing. J Mol Biol. 2003;332:1015–24.
Kim Y, Li H, Binkowski TA, Holzle D, Joachimiak A. Crystal structure of fatty acid/phospholipid synthesis protein PlsX from Enterococcus faecalis. J Struct Funct Genomics. 2009;10:157–63.
Veit A, Polen T, Wendisch VF. Global gene expression analysis of glucose overflow metabolism in Escherichia coli and reduction of aerobic acetate formation. Appl Microbiol Biotechnol. 2007;74:406–21.
Freddolino PL, Amini S, Tavazoie S. Newly identified genetic variations in common Escherichia coli MG1655 stock cultures. J Bacteriol. 2012;194:303–6.
Baba T, Ara T, Hasegawa M, Takai Y, Okumura Y, Baba M, et al. Construction of Escherichia coli K-12 in-frame, single-gene knockout mutants: the Keio collection. Mol Syst Biol. 2006;2006(2):0008.
bcl2fastq and bcl2fastq2 Conversion Software. Illumina. 2019. https://emea.support.illumina.com/sequencing/sequencing_software/bcl2fastq-conversion-software.html. Accessed 19 Nov 2019.
Li H, Durbin R. Fast and accurate short read alignment with Burrows–Wheeler transform. Bioinformatics. 2009;25:1754–60.
Garrison E, Marth G. Haplotype-based variant detection from short-read sequencing. arXiv preprint, arXiv:1008.0576. 2012.
Pinto F, Thapper A, Sontheim W, Lindblad P. Analysis of current and alternative phenol based RNA extraction methodologies for cyanobacteria. BMC Mol Biol. 2009;10:79.
Schmieder R, Lim YW, Edwards R. Identification and removal of ribosomal RNA sequences from metatranscriptomes. Bioinformatics. 2012;28:433–5.
Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, et al. STAR: Ultrafast universal RNA-seq aligner. Bioinformatics. 2013;29:15–21.
Anders S, Pyl PT, Huber W. HTSeq-A Python framework to work with high-throughput sequencing data. Bioinformatics. 2015;31:166–9.
Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26:139–40.
Anders S, Huber W. Differential expression analysis for sequence count data. Genome Biol. 2010;11:R106.
Trapnell C, Hendrickson DG, Sauvageau M, Goff L, Rinn JL, Pachter L. Differential analysis of gene regulation at transcript resolution with RNA-seq. Nat Biotechnol. 2013;31:46–53.
Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B. 1995;57:289–300.
Daniel G, Gibson DG, Young L, Chuang R-Y, Venter JC, Hutchison CA, et al. Enzymatic assembly of DNA molecules up to several hundred kilobases. Nat Methods. 2009;6:343–5.
Genebridges Red/ET Kit. 2019. https://www.genebridges.com/products/redet-kits. Accessed 9 May 2019.
The authors also would like to thank BASF partners Julia Ena Hartig, Michael Breuer, Heike Brüser, Roland Minges, Robert Thummer, Michael Braun, and Oliver Oswald for fruitful discussions. The authors also acknowledge UFZ technical support from Kristin Lindstaedt, Ron Stauder, Martin Wegner, Danielle Heiland and thank Anna Hoschek for kindly providing plasmid pAH030.
Open Access funding enabled and organized by Projekt DEAL. This work was financed by BASF SE Ludwigshafen, Germany.
Ethics approval and consent to participate
Consent for publication
The authors have transferred intellectual property rights to the commercial partner BASF SE, who has filed results for a patent. The authors have no further financial interest in the results.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Mean Growth Rates, Isoprenol Concentration. Table S2. Relative Fitness of fast-growing plate colonies. Figure S1. Gene expression changes in response to PyghB mutation. Figure S2. Gene expression changes in the rob regulon. Figure S3. Clustering of observed genotypes based on Euclidian distance and farthestneighbor (complete linkage clustering).(A) Labeling according to time-point and replicate (B) Labeling with 9 selected mutations. Table S3. Genotype labels. Figure S4. Proposed model for the evolution of cultures. Figure S5. Schematic presentation of implementation of ALE Model. Figure S6. Simulation of evolution three adaptive evolution experiments from ALE model. Figure S7. Sensitivity analysis of the adapted fitness parameters for the ALE model. Figure S8. Sensitivity analysis of the adapted initial proportion parameter Pi for the ALE model. Table S4. Likelihood ratio of maximal likelihood estimates of models with contamination compared to reference model without contamination. Table S5. Likelihood ratio of maximal likelihood estimates of models with contamination compared to reference model without contamination (different contamination proportions). Table S6. GO IDs (by protein) of target mutations. Table S7. GO IDs (by protein) of top 20 highly differentially expressed genes in adaptedstrains. Table S8. Plasmids. Table S9. Primers.
About this article
Cite this article
Babel, H., Krömer, J.O. Evolutionary engineering of E. coli MG1655 for tolerance against isoprenol. Biotechnol Biofuels 13, 183 (2020). https://doi.org/10.1186/s13068-020-01825-6
- Adaptive laboratory evolution
- E. coli