- Open Access
Identification of functional butanol-tolerant genes from Escherichia coli mutants derived from error-prone PCR-based whole-genome shuffling
- Xueting He†1, 2,
- Tingli Xue†1, 2,
- Yuanyuan Ma1, 3, 4, 5Email author,
- Junyan Zhang1,
- Zhiquan Wang1, 2,
- Jiefang Hong1,
- Lanfeng Hui6,
- Jianjun Qiao4, 5,
- Hao Song4, 5 and
- Minhua Zhang1, 3
© The Author(s) 2019
- Received: 22 August 2018
- Accepted: 12 March 2019
- Published: 1 April 2019
Butanol is an important biofuel and chemical. The development of butanol-tolerant strains and the identification of functional butanol-tolerant genes is essential for high-yield bio-butanol production due to the toxicity of butanol.
Escherichia coli BW25113 was subjected for the first time to error-prone PCR-based whole-genome shuffling. The resulting mutants BW1847 and BW1857 were found to tolerate 2% (v/v) butanol and short-chain alcohols, including ethanol, isobutanol, and 1-pentanol. The mutants exhibited good stability under butanol stress, indicating that they are potential host strains for the construction of butanol pathways. BW1847 had better butanol tolerance than BW1857 under 0–0.75% (v/v) butanol stress, but showed a lower tolerance than BW1857 under 1.25–2% (v/v) butanol stress. Genome resequencing and PCR confirmation revealed that BW1847 and BW1857 had nine and seven single nucleotide polymorphisms, respectively, and a common 14-kb deletion. Functional complementation experiments of the SNPs and deleted genes demonstrated that the mutations of acrB and rob gene and the deletion of TqsA increased the tolerance of the two mutants to butanol. Genome-wide site-specific mutated strains DT385 (acrB C1198T) and DT900 (rob AT686–7) also showed significant tolerance to butanol and had higher butanol efflux ability than the control, further demonstrating that their mutations yield an inactive protein that enhances butanol resistance characteristics.
Stable E. coli mutants with enhanced short alcohols and high concentrations of butanol tolerance were obtained through a rapid and effective method. The key genes of butanol tolerance in the two mutants were identified by comparative functional genomic analysis.
- Escherichia coli
- Butanol-tolerant mutants
- Comparative functional genome
- Transcription regulator factor
- Efflux pump
- Quorum-sensing signal transporter
Butanol is a potential superior alternative to ethanol as biofuel, with a higher energy density, lower hygroscopicity and volatility, and less corrosive. Butanol has also been used in a wide range of fields, including the food, plastics, and pharmaceutical industries . Engineering a butanol biosynthetic pathway into Escherichia coli (E. coli) could produce n-butanol at the level of grams per liter and compete with the traditional industrial producer, genus Clostridium [1–3]. However, the poor butanol tolerance of bacteria leads to the inhibition of cell growth and a decrease in butanol yield, which has been a bottleneck in biobutanol production. Thus, the economic production of butanol relies on the improvement in butanol tolerance of the bacterial producers [4, 5].
To improve the robustness of chassis strains, rational metabolic engineering strategies have been used to generate butanol or other fuel-tolerant mutants of E. coli [4, 6]. Genes in relation to membrane function and ion transport system have been reported to respond to butanol stress, the tolerance to which could be improved by the overexpression of genes encoding the efflux pump [7–9], ion transport proteins entC and feoA [10, 11], and membrane-targeted metallothionein . Heat shock protein (HSP) genes also respond to butanol stress; as such, their upregulation or overexpression could also lead to an increase in tolerance .
The mechanism of butanol tolerance, thus, involves multiple physiological processes mediated by a gene network. Butanol-tolerant E. coli strains are also obtained using random engineering strategies, such as global transcription machinery engineering (GTME) and evolution engineering, among others [14–17]. These mutants developed by random mutation and selection generally show a greater tolerance than those with functional changes in one gene, and could tolerate 1–2.0% (v/v) butanol [16, 18]. This may due to the fact that multiple mutated genes endow strains with a greater tolerance. When membrane-related functional genes are overexpressed in butanol-tolerant E. coli through evolution, the final engineered strain can grow under 2% butanol stress , indicating that a combinational strategy could effectively improve tolerance. If the key butanol-tolerant genes could be identified and clarified, these favorable mutations could be introduced to chassis strain, also referred to as inverse metabolic engineering . Nevertheless, the genes of many mutants have not yet been functionally identified, which limits our understanding of the butanol tolerance mechanism as well as the potential to obtain butanol-tolerant characteristics.
It is still challenging to develop high butanol-tolerant strains using metabolic engineering or combinatorial approach because of the current limited understanding of the tolerance mechanisms; it is first necessary to establish a simple and rapid evolution strategy for the acquisition of resistant mutants, to thereafter identify the functional stress-tolerant genes and elucidate the mechanism of butanol tolerance. In this study, error-prone PCR (epPCR)-based whole-genome shuffling was used to improve the butanol tolerance of E. coli BW25113. As a result, two mutants tolerating 2% (v/v) butanol were successfully obtained. The mutation sites of the two strains were then identified by whole-genome resequencing, and comparative functional genomic analysis was performed to identify the key butanol-tolerant genes.
Isolation and characterization of butanol-tolerant mutants by error-prone PCR-based whole-genome shuffling
Tolerance of BW1847 and BW1857 to various short-chain alcohols
Growth and stability under varying butanol stress
BW1857 showed slower growth and lower cell density than BW25113 and BW1847 both under 0.75% butanol stress conditions and without butanol, demonstrating that growth of BW1857 was inhibited (Fig. 3a, b). However, BW1857 had a higher cell density than BW1847 in the presence of 1–1.5% butanol, indicating that BW1857 exhibited a better tolerance in high butanol concentrations (Fig. 3c–e). The cell densities of BW1847 and BW1857 were 68% and 204% higher, respectively, than that of the control BW25113 at 1.5% butanol, demonstrating their tolerance to high concentrations of butanol (Fig. 3e). BW1857 showed an inhibited and enhanced growth under low and high concentrations of butanol, respectively, implying that it contains multiple gene mutations, some of which may lead to improved butanol tolerance, with others resulting in growth inhibition.
Growth under high concentrations of butanol was evaluated (Fig. 3f). BW1847 and BW1857 showed 1.5- and 2.7-fold, and 3.5- and 5.9-fold increased growth compared to the wild-type in the presence of 1.6 and 1.8% butanol, respectively. The cell growth of BW1847 and BW1857 was 1.9–3.7% and 16.5–17.3%, respectively, of the control strain under high butanol stress (1.9–2.0%). These data demonstrate that the butanol tolerance of the two mutants was significantly higher than that of BW25113, and they could tolerate 2% butanol. BW1857 showed a relatively higher tolerance than BW1847 in concentrations of butanol above 1.25% (Fig. 3c–f).
Genome-wide identification of mutations by comparative genomics
Mutations identified in the BW1847 and BW1857 strains
BW1847 and BW1857
G (251) C
V (149) G
Formate hydrogenlyase subunit
A (300) E
Anaerobic sulfatase maturase AslB
V (408) E
Multidrug efflux RND transporter permease subunit AcrB
L (400) F
Bifunctional (p)ppGpp synthetase II/guanosine-3′,5′-bis pyrophosphate 3′-pyrophosphohydrolase
nusG transcription termination/antitermination
G (146) C
Right oriC-binding transcriptional activator, AraC family
50S ribosomal protein L2
I (164) N
Translation initiation factor IF-2
225 bp deletionc
Functional identification of mutated genes in BW1847 and BW1857
Strains used in functional complementary study
E. coli Top10F′
Host for plasmid construction
E. coli BW25113
Host for plasmid expression, F− Δ(araD-araB)567 ΔlacZ4787(::rrnB-3)λ−rph-1 Δ(rhaD-rhaB)568 hsdR514
pBAD30 is a arabinose-inducible expression vector, with ori 15A replicon, ApR
BW25113 (pBAD30 carries wild/mutated pgsA)
BW25113 (pBAD30 carries wild/mutated yehQ)
BW25113 (pBAD30 carries wild/mutated hycD)
BW25113 (pBAD30 carries wild/mutated aslB)
BW25113 (pBAD30 carries wild/mutated cdsA)
BW25113 (pBAD30 carries wild/mutated acrB)
BW25113 (pBAD30 carries wild/mutated nusG)
BW25113 (pBAD30 carries wild/mutated rob)
BW25113 (pBAD30 carries wild/mutated spoT)
BW25113 (pBAD30 carries wild/mutated infB)
BW25113 (pBAD30 carries wild/mutated rplB)
BW25113 ΔpgsA::cmr (pBAD30 carries wild/mutated pgsA)
BW25113 ΔyehQ::cmr (pBAD30 carries wild/mutated yehQ)
BW25113 ΔhycD::cmr (pBAD30 carries wild/mutated hycD)
BW25113 ΔaslB::cmr (pBAD30 carries wild/mutated aslB)
BW25113, ΔcdsA::cmr (pBAD30 carries wild/mutated cdsA)
BW25113, ΔacrB::cmr (pBAD30 carries wild/mutated acrB)
BW25113, ΔnusG::cmr (pBAD30 carries wild/mutated nusG)
BW25113, Δrob::cmr (pBAD30 carries wild/mutated rob)
BW25113, ΔspoT::cmr (pBAD30 carries wild/mutated spoT)
BW25113, ΔinfB::cmr (pBAD30 carries wild/mutated infB)
The deletion and overexpression of the rob gene (RS22900) also had no influence on growth (Fig. 4c). The rob deletion mutant D900 (pBAD30) showed 1.6-fold higher maximum cell density than BW25113 (pBAD30). Nevertheless, the overexpression of the rob gene resulted in a 37.5–57.0% decrease of growth under 0.7% butanol stress, indicating that the deficiency of rob led to an enhanced butanol tolerance. The complementation strain D900 (pM900) showed an inhibited growth similar to that of D900 (pBAD30) under butanol stress (Fig. 4d), indicating that the overexpression of the truncated rob gene was unable to restore the function of rob. This also demonstrates that the deletion of the AT686–7 base in the rob gene in BW1847 causes an inactivation of Rob function, which improves the butanol tolerance of BW1847. Another 9 mutated genes, including yhcD (RS14165) and nusG (RS20660), were identified but had no distinct influence on growth or butanol tolerance.
Functional identification of deleted genes
The deleted 14-kb fragment contained 17 genes, in which the pntA (RS08390) and pntB (RS08385) belonged to one operon and encoded the alpha and beta subunits of pyridine nucleotide transhydrogenase. The deletion mutant of each gene and mutant D123 with simultaneous deletion of the 17 genes (Additional file 1: Table S1) were constructed to test the butanol-tolerant function of these genes (Additional file 1: Table S5, Additional data). The large fragment deletion caused a slight growth inhibition but resulted in a slight improvement of butanol tolerance under 1–1.5% (v/v) butanol stress (Additional file 1: Figure S2).
Growth improvement of the TqsA deletion mutant D380
Butanol % (V/V)
6.0 ± 0.6
3.3 ± 1.0
Multiple gene-deleted strains, D1 (with deletion of genes RS08320–RS08335), D2 (with deletion of genes RS08340–RS08375), D3 (with deletion of genes RS08380–RS08395), D12, D13, and D23, were constructed to investigate the effect of combined deletion on growth and tolerance. The growth and tolerance of the D2, D12, D23, and D13 strains were lower than that of the other strains, indicating that the deletion of RS08340–RS08375 leads to a decrease in growth (Additional file 1: Figure S5). The deleted fragment contained dgsA and mdtJ, further indicating that the simultaneous deletion of these two genes causes a slight inhibition of growth and butanol resistance. Therefore, the deletion of TqsA is a key factor for improving the butanol tolerance characteristics of D123, BW1847, and BW1857 grown under 1–1.5% butanol stress.
Butanol tolerance function of rob and acrB genes
The rob gene encodes OriC-binding transcriptional activator, and AT686–7 deletion of rob in genome did not affect growth (Fig. 5c) but caused a significant increase in growth under 0.75% (v/v) butanol stress, indicating that deficiency of AT686–7 in rob can improve butanol tolerance(Fig. 5d). The rob-deleted strain D900 also showed higher tolerance, demonstrating that the truncated rob leads to its inactivation, causing an enhanced butanol tolerance. Moreover, the butanol efflux ability of the two mutants was also significantly improved (Fig. 5d). The inactivation of Rob could change the gene expression profile and corresponding physiological/biochemical processes under butanol stress. Some of the response genes regulated by Rob may function in the efflux of potentially toxic molecules to the outside of the cells.
Bacterial genomes can be conveniently and rapidly modified using epPCR whole-genome shuffling with a highly efficient recombination system. The recombination efficiency of E. coli recombinase recA is very low . Therefore, BW25113 (pKD46) containing a high-efficient recombinase system was chosen as the host. This resulted in two mutants, BW1847 and BW1857, which tolerated 2% butanol stress. Even though the growth of BW25113 decreased by 50% under 0.8% butanol stress , the two mutants in the study had a fourfold higher maximum cell density than BW25113 under 0.75% butanol stress, indicating a great improvement in their butanol tolerance. Until now, most E. coli strains with single gene modifications could only tolerate 0.75–1.5% butanol stress [7, 13]. However, the two mutants obtained in this study showed better tolerance than most reported strains. A few mutants tolerating 2% butanol have been previously obtained via combined strategies or multi-gene mutations [16, 18], by which the function of multiple genes involved in various physiology events is changed to respond to stress. Even though multiple mutations can easily change physiological and biochemical functions, thereby yielding a better tolerance, they can lead to a decrease in tolerance function. Consequently, the identification of key resistant genes is necessary.
Many butanol-tolerant strains are also reported to show relevant tolerance to other short-chain alcohols, such as ethanol, isobutanol, 1-pentanol, 1,2,4-butanetriol, and n-propanol [25, 26]. These short-chain alcohols possess similar molecular structures and exhibit toxicity. The fatty acid composition, structure, and fluidity of the cell membrane are similarly changed in response to stress by these short-chain alcohols. Therefore, the tolerance mechanism towards these alcohols is similar [27, 28].
BW1847 and BW1857 have 5 common SNPs, which also should be the mutations from BW184. However, functional complementation experiments carried out on each of these five genes did not confirm their butanol tolerance. As such, the tolerance of BW184 to butanol stress may result from the functional accumulation of all five mutated genes. The maximum cell density of BW1857 was lower than that of BW25113 and BW1847, and BW1857and BW1847 contained two different mutations, infBdeletion331–555 and rplBT491A. Nevertheless, the functional complementation experiments of the two mutated genes did not elucidate growth inhibition. Moreover, the deletion mutant and the point mutant strain of rplB could not be successfully obtained due to the deletion being lethal, neither is it easily edited by Cas9-sgRNA. However, the deletion of a 14-kb fragment in both strains leads to slight growth inhibition, suggesting that the combined actions of the two SNPs and the deleted genes cause growth inhibition of BW1857.
Interestingly, a 14-kb deletion was observed in both mutants, but large DNA fragment deletion was not reported to be present in other mutants obtained by epPCR-based genome shuffling [21, 29], which is because only a few mutants, obtained using this technology, are analyzed by genome resequencing . Consequently, an attempt was made to demonstrate the origin of the 14-kb deletion by experimentation and inference. If the deletion could be obtained by the lambda-Red recombination system of epPCR-based shuffling, then the deletion could have resulted from crossover recombination of a certain random epPCR product with coincident left and right homologous arm sequences of the target deletion fragment. The lambda-Red recombination system (pKD46), used in the shuffling, was used to knockout the target region, and twenty-nine transformants with 14-kb deletions were obtained (Additional file 1: Additional data, Figure S9, Table S11), indicating the possibility of deriving targeted deletions using this technology. Besides, final mutation rates of the two mutants were 1000-fold higher than the spontaneous mutation rates (Additional file 1: Additional data). Therefore, the 14-kb deletion was more likely due to homologous recombination of random epPCR products via lambda-Red recombinase used in the epPCR shuffling technology, than from spontaneous mutations during the recovery and screening phases after electroporation.
Three important tolerance genes were identified in this study: mutant acrB, mutant rob, and the deletion of TqsA, all increased butanol tolerance. Mutated AcrBI466F is reported to pump butanol out of the cell, leading to an improvement in butanol tolerance . AcrBL400F in this study is also found to exhibit butanol-efflux ability, and the mutated I466F and L400F are both located on the transmembrane domain, and would change the spatial structure. AcrB forms a trimer composed of loose (L), tight (T), and open (O) conformational states. Toxic compounds were bound to the L monomer and then transported to the T monomer, before finally being excreted out of cells by the O monomer  (Additional file 1: Figure S6 A–C). The mutation in AcrB, thus, allowed for the binding and transportation of butanol, resulting in improved efflux function. Interestingly, acrB deficiency was reported to yield an improvement in butanol tolerance for the first time. This deletion resulted in an open funnel-shaped compound formed by AcrA and TolC in the periplasmic space (Additional file 1: Figure S6D), such that butanol was freely and directly transported out of cells by the uncompleted efflux pump system.
The transcription factor Rob responds to many stimuli, including antibiotics, organic solvents, and oxidative stressors, and its function in resistance to the stress is shown by overexpression of Rob through a plasmid [32–34]; while, the rob-deleted strains exhibit a decreased tolerance to the stresses [33, 34]. Nevertheless, the inactivation of Rob in this study enhanced the butanol tolerance (Figs. 4, 5). The Rob-inactive strain would have a different expression profile than the wild-type under butanol stress. The response genes, and thus the physiological roles regulated by Rob would be interesting and complex, and should be investigated to provide insights into the mechanism of butanol tolerance in future studies.
The deletion of TqsA, encoding an AI-2 transporter, was found to improve the butanol tolerance. Transmembrane protein TqsA transports quorum-sensing molecule AI-2 out of the cells, and its deletion increases the concentration of intracellular AI-2, which enhances the motility of cells, allowing cells to escape from toxic compounds [35, 36], and improves the biofilm thickness by up-regulating the expression of lipopolysaccharide (LPS)-synthesis-related genes [37, 38]. Therefore, the deletion of TqsA may lead to the improvement of intracellular AI-2, which subsequently increases cell motility and biofilm thickness, resulting in enhanced butanol tolerance.
Both the mutation and the deletion strain of acrB and rob improved growth by about 1.7- to 2.0-fold, compared to the control BW25113 under 0.75% butanol stress (Additional file 1: Table S10). The deletion of TqsA led to an improved tolerance (1.7-fold) in high concentration butanol media (1–1.25%). These three genes result in features for butanol tolerance, which could provide target genes for the engineering of a butanol-tolerant strain. The co-expression of these three genes is speculated to be able to result in an accumulated improvement of tolerance, which should be interesting and would be investigated further in future studies.
Among the 17 deleted and 11 mutated genes from the two mutants, only three genes were identified to exhibit significant butanol tolerance in this study. Correspondingly, among 16 mutated genes in mutant E72, only three mutations were beneficial . These results demonstrate that only a few of multiple mutations may be valuable for functional improvement. Some mutations may produce an unobvious objective function or cause decreased target function. Therefore, improved characteristics in a mutant with random mutations may result from the combined actions of several key mutated genes. Consequently, the identification of key functional genes is important for revealing a corresponding mechanism and providing genetic resources for the construction of optimal chassis strains using rational metabolic engineering strategies.
Two stable and high-concentration butanol-tolerant E. coli mutants were obtained by epPCR-based genome shuffling, a rapid and effective method. The two mutants were found to also tolerate short alcohols. Their mutation genes were revealed through genome resequencing. A comparative functional genomic study identified the key genes involved in butanol tolerance in both mutants, including rob, encoding transcript regulator factor Rob, acrB, encoding efflux pump transporter, and tqsA, encoding a quorum-sensing molecule transporter. The mutation or deletion of these three genes resulted in an improvement in butanol tolerance. Interestingly, the present data found for the first time that the deletion of acrB and rob yielded a significantly increased butanol-efflux ability, and the deletion of tqsA caused slight growth inhibition and high-concentration butanol-tolerance characteristics. These surprising and novel gene functions could provide a new point of view for the investigation of tolerance mechanism. This study, thus, provides valuable information for the construction of butanol-resistant strains and insights into the mechanism of butanol tolerance.
Strains and culture conditions
E. coli BW25113 (pKD46) was used as a starting strain. Strain BW1847 was deposited in the China Center for Type Culture Collection (CCTCC; Accession NO. M2015365). All E. coli strains were cultured in Luria–Bertani (LB) medium (5 g/L yeast extract, 10 g/L Casein tryptone, and 10 g/L NaCl) at 30 °C or 37 °C with 180–200 rpm shaking for growth.
Error-prone PCR-based genome shuffling of E. coli
Genomic DNA of E. coli BW25113 was extracted from 2 mL of overnight cell culture and then used for the epPCR templates. EpPCR was conducted according to the protocols developed by Ye et al.  with some modifications  (Fig. 1a). The amplified PCR products were then concentrated 5- to 10-fold by ethanol-precipitated for electro-transformation (Additional file 1: Additional data).
Escherichia coli BW25113 (pKD46) were cultured in LB media with 100 μg/mL ampicillin at 30 °C until optical density at 600 nm (OD600) reached 0.15–0.2, and then L-arabinose was added to the media at a final concentration of 20 mM to induce expression of lambda-Red recombinase. When the cultures reached an OD600 of 0.45–0.5, electro-competent cells were prepared. Competent cells (40 μL) were mixed with 1–5 μL of ethanol-precipitated PCR products, and the mixture was incubated in ice bath for 5 min and then transferred to a pre-chilled 1-mm-gap electroporation cuvette (Eppendorf, Germany). The cells were electroporated with a pulse of 1800 v at 5 ms using Electroporator 2510 (Eppendorf, Germany), and the suspensions were immediately mixed with 1 mL of LB broth and transferred to a microcentrifuge tube for culture at 37 °C with 100–110 rpm shaking for 1 h. About 200 μL of culture was then spread on LB agar plates containing an initial 1.2%, 1.4%, 1.6% and 1.8% (v/v) butanol, respectively. After 16–24 h of 37 °C incubation, the transformants were picked up from these plates and inoculated into 4-mL LB containing 0.7–0.9% butanol for growth evaluation. The best butanol-tolerant mutant was used for subsequent round of shuffling. The pKD46 plasmid was extracted and then transformed to the target strain to provide recombinase system for next round of shuffling (Fig. 1a).
Butanol tolerance assay
Determination of short-chain alcohol tolerance of BW1847 and BW1857
Seed cultures of BW1847, BW1857 and BW25113 were prepared as aforementioned and then inoculated into 4-mL LB with butanol (0–1.5%), isobutanol (0–1.5%), and ethanol (0–7%), 1-pentanol (0–0.5%), acetic acid (0–0.3%), respectively. Strains were cultured at 37 °C with 200-rpm shaking, and initial OD600 was controlled at 0.1. The OD600 values at 12 h were determined, and the relative OD600 ratios of transformants to control BW25113 cultured under same conditions were calculated to evaluate the butanol tolerance.
Growth and butanol concentration detection
All growth detection experiments in butanol stress were performed using screw-cap test tubes or conical flasks as aforementioned. Cell density was determined by OD600 with a Cary 50 spectrophotometer (Varian Inc.). Butanol concentration in media during culture was detected by high-performance liquid chromatography (HPLC) using Organic Acid Analysis Column (Aminex HPX-87H Ion Exclusion Column, 300 mm × 7.8 mm) at 45 °C and sulphuric acid (4 mM) as mobile phrase at 0.8 mL/min. The butanol concentration at 12 h was subtracted from the initial butanol concentration in the media, and then divided by the OD600 value of strain grown at 12 h to obtain the relative butanol content in one cell from each strain, which is used to measure the butanol-efflux ability of cells. All experiments for growth evaluation were performed in triplicates, and statistical differences were analyzed using ANOVA, Single star indicated p ≤ 0.05. The heights of the bars and the error bars in figures indicated the means ± standard deviations.
Evaluation of butanol-tolerance stability of BW1847 and BW1857
Strain BW25113, BW1847 and BW1857 were cultured in LB medium until OD600 value reached 4. These strains were then serially transferred ten times (about 50–60 generation) into fresh LB medium with an initial inoculation concentration of 1% (v/v). The serial cultures and the cultures stored in glycerol were all inoculated into 3-mL LB media, and pre-cultured for 12–15 h with 200-rpm shaking. All the pre-cultures were then transferred into 10-mL LB media for the second pre-culture. After 12 h of incubation, the second pre-cultures were inoculated into 250-mL screw-cap conical flask with 50-mL LB containing 0, 0.75, 1.0, 1.25 and 1.5% (v/v) butanol for growth evaluation, respectively. The initial OD600 was controlled at 0.2, and the OD600 value was determined every 4 h for comparison of growth curves.
Genome resequencing of mutants and identification of Single Nucleotide Polymorphism (SNP)
Genomic DNA of BW1847 and BW1857 was submitted to Meiji Company (Shanghai) and Novogene (Beijing) for genome resequencing, respectively. Genome sequence of BW25113 (NZ-CP009273.1) was used as a reference sequence for mutation analysis. The sequencing data were deposited in Genebank (SRA accession No. PRJNA397315 and No. PRJNA399799 for BW1847 and BW1857, respectively). The mutated genes were amplified by PCR with 2×Pfu Mastermix (Kangwei Biotech Co. Beijing, China) for sequencing confirmation (Additional file 1: Additional data).
To confirm a deletion of a large gene fragment in genome of BW1857, PCR identification was performed using the genomic DNA of BW1857 as template and primer sets QS-F and QS-R (Additional file 1: Table S1, Additional data), which was located in upstream and downstream of the deleted gene fragment, respectively. The PCR products were then sequencing verified (Sangong Biotech, Shanghai, China).
Construction of deletion and overexpression strains
Knockout of each mutated gene was performed according to Datsenko et al.’s reports  with a minor modification (Additional file 1: Additional data) for functional identification of the mutated genes in BW1847 and BW1857. Overexpression vectors were constructed by ligating the target wild-type gene and its mutated gene into multiple clone sites of pBAD30 vector (as described in Additional file 1: Additional data, Tables S1–S5, Figure S7), yielding the pBAD30-Wgene and the pBAD30-Mgene plasmid, respectively (Table 2). The plasmids were then transferred into strain BW25113 and a target-gene-deleted strain Dgene (Additional file 1: Table S5), and the serial strains were obtained (Table 2), including BW25113 (pBAD30-Wgene), BW25113 (pBAD30-Mgene), BW25113 (pBAD30), Dgene (pBAD30-Wgene), Dgene (pBAD30-Mgene) and Dgene (pBAD30).
Functional identification of deleted genes in BW1847 and BW1857
Knockout of each gene in the large deleted fragment of BW1857 (Table 1) was performed according to the above description (Additional file 1: Additional data). These mutants (Additional file 1: Table S5) and control BW25113 were inoculated into 3-mL LB liquid media for 12 h pre-culture. The subcultures were then transferred into 25 mL of tubes containing 4-ml LB liquid medium with 0, 0.5%, and 0.75% (v/v) butanol, respectively. The initial OD600 was controlled at 0.1, and the OD600 were measured at 0 and 12 h to preliminarily evaluate their growth difference. The strains that showed higher cell densities than BW25113 were selected to measure growth curves. These transformants and BW25113 were subcultured as described above, and then transferred into 250-mL conical flasks which contained 50-mL LB liquid medium with 0, 0.5% and 0.75% (v/v) butanol, respectively. The initial OD600 was controlled at 0.1, and the OD600 value was measured every 1.5 or 2 h to observe the growth difference between transformants and control.
Site-specific genomic mutation and deletion in E. coli
Strains with the large fragment deletion (D123) and partially deleted genes (D1, D2, D3, D12, D13 and D23) were constructed using CRISPR/Cas9  system according to the order of deleted genes on genome (Additional file 1: Tables S5, S7 and S8, Additional data) to detect the combined impact of several deleted genes on 14-kb deletion fragment on growth. Genome-wide site-specific mutant DT385 (mutation of C1198T in acrB gene) and strain DT900 (AT686–7 deletion in rob gene) were constructed using scarless gene modification strategy according to the previous reports [41, 42] with some modification (Additional file 1: Additional data, Table S9, Figure S8).
YM designed the experiments, analyzed the data and wrote the manuscript. XH, TX, JZ, JH carried out experiments and drafted the manuscript. XH, TX and JZ carried out all experiments. LH, JQ and HS designed the experiments and drafted the manuscript. YM, JQ, HS and MZ supervised the study and revised the manuscript. All authors read and approved the final manuscript.
The authors appreciate Prof. Sheng Yang for kindly providing CRISPR-Cas9 system, and thank Prof. Tao Chen and Prof. Zhiwen Wang for kindly providing site-specific mutagenesis tool. The author would like to thank Zheng Hu and Jingxiao Ji for data analysis of genome resequencing. The author would also like to thank the members of Novogene for their valuable informations on the sequence submission to NCBI.
The authors declare that they have no competing interests.
Availability of data and materials
Consent for publication
Authors agreed to publish this article.
Ethics approval and consent to participate
This research was supported by the National Natural Science Foundation of China (NSFC 30900033), and the Natural Science Foundation of Tianjin (18JCYBJC24200).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
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.
- Jiang Y, Liu J, Jiang W, Yang Y, Yang S. Current status and prospects of industrial bio-production of n-butanol in China. Biotechnol Adv. 2015;33(7):1493–501.View ArticleGoogle Scholar
- Dong H, Zhao C, Zhang T, Zhu H, Lin Z, Tao W, Zhang Y, Li Y. A systematically chromosomally engineered Escherichia coli efficiently produces butanol. Metab Eng. 2017;44:284–92.View ArticleGoogle Scholar
- Shen CR, Lan EI, Dekishima Y, Baez A, Cho KM, Liao JC. Driving forces enable high-titer anaerobic 1-butanol synthesis in Escherichia coli. Appl Environ Microbiol. 2011;77(9):2905–15.View ArticleGoogle Scholar
- Dong H, Zhao C, Zhang T, Lin Z, Li Y, Zhang Y. Engineering Escherichia coli cell factories for n-butanol production. Adv Biochem Eng Biotechnol. 2016;155:141–63.PubMedGoogle Scholar
- Nanda S, Golemi-Kotra D, McDermott JC, Dalai AK, Gökalp I, Kozinski JA. Fermentative production of butanol: perspectives on synthetic biology. Nat Biotechnol. 2017;37(Pt B):210–21.Google Scholar
- Siu Y, Fenno J, Lindle JM, Dunlop MJ. Design and selection of a synthetic feedback loop for optimizing biofuel tolerance. ACS Synth Biol. 2018;7(1):16–23.View ArticleGoogle Scholar
- Boyarskiy S, Davis López S, Kong N, Tullman-Ercek D. Transcriptional feedback regulation of efflux protein expression for increased tolerance to and production of n-butanol. Metab Eng. 2016;33:130–7.View ArticleGoogle Scholar
- Dunlop MJ, Dossani ZY, Szmidt HL, Chu HC, Lee TS, Keasling JD, Hadi MZ, Mukhopadhyay A. Engineering microbial biofuel tolerance and export using efflux pumps. Mol Syst Biol. 2011;7:487.View ArticleGoogle Scholar
- Fisher MA, Boyarskiy S, Yamada MR, Kong N, Bauer S, Tullman-Ercek D. Enhancing tolerance to short-chain alcohols by engineering the Escherichia coli AcrB efflux pump to secrete the non-native substrate n-butanol. ACS Synth Biol. 2014;3(1):30–40.View ArticleGoogle Scholar
- Reyes LH, Abdelaal AS, Kao KC. Genetic determinants for n-butanol tolerance in evolved Escherichia coli mutants: cross adaptation and antagonistic pleiotropy between n-butanol and other stressors. Appl Environ Microbiol. 2013;79(17):5313–20.View ArticleGoogle Scholar
- Reyes LH, Almario MP, Kao KC. Genomic library screens for genes involved in n-butanol tolerance in Escherichia coli. PLoS ONE. 2011;6(3):e17678.View ArticleGoogle Scholar
- Chin WC, Lin KH, Chang JJ, Huang CC. Improvement of n-butanol tolerance in Escherichia coli by membrane-targeted tilapia metallothionein. Biotechnol Biofuels. 2013;6(1):130.View ArticleGoogle Scholar
- Abdelaal AS, Ageez AM, Abd El-Hadi AE, Abdallah NA. Genetic improvement of n-butanol tolerance in Escherichia coli by heterologous overexpression of groESL operon from Clostridium acetobutylicum. 3 Biotech. 2015;5(4):401–10.View ArticleGoogle Scholar
- Lee JY, Yang KS, Jang SA, Sung BH, Kim SC. Engineering butanol-tolerance in Escherichia coli with artificial transcription factor libraries. Biotechnol Bioeng. 2011;108(4):742–9.View ArticleGoogle Scholar
- Reyes LH, Almario MP, Winkler J, Orozco MM, Kao KC. Visualizing evolution in real time to determine the molecular mechanisms of n-butanol tolerance in Escherichia coli. Metab Eng. 2012;14(5):579–90.View ArticleGoogle Scholar
- Si HM, Zhang F, Wu AN, Han RZ, Xu GC, Ni Y. DNA microarray of global transcription factor mutant reveals membrane-related proteins involved in n-butanol tolerance in Escherichia coli. Biotechnol Biofuels. 2016;9:114.View ArticleGoogle Scholar
- Zhu L, Cai Z, Zhang Y, Li Y. Engineering stress tolerance of Escherichia coli by stress-induced mutagenesis (SIM)-based adaptive evolution. Biotechnol J. 2014;9(1):120–7.View ArticleGoogle Scholar
- Bui LM, Lee JY, Geraldi A, Rahman A, Lee JH, Kim SC. Improved n-butanol tolerance in Escherichia coli by controlling membrane related functions. J Biotechnol. 2015;204:33–44.View ArticleGoogle Scholar
- Zhao XM, Chen T, Wang ZW. Metabolic engineering. Beijing: Higher Education Press; 2015.Google Scholar
- Huang S, Xue T, Wang Z, Ma Y, He X, Hong J, Zou S, Song H, Zhang M. Furfural-tolerant Zymomonas mobilis derived from error-prone PCR-based whole genome shuffling and their tolerant mechanism. Appl Microbiol Biotechnol. 2018;102(7):3337–47.View ArticleGoogle Scholar
- Ye L, Zhao H, Li Z, Wu JC. Improved acid tolerance of Lactobacillus pentosus by error-prone whole genome amplification. Bioresour Technol. 2013;135:459–63.View ArticleGoogle Scholar
- Seeger MA, Schiefner A, Eicher T, Verrey F, Diederichs K, Pos KM. Structural asymmetry of AcrB trimer suggests a peristaltic pump mechanism. Science. 2006;313(5791):1295–8.View ArticleGoogle Scholar
- Datsenko KA, Wanner BL. One-step inactivation of chromosomal genes in Escherichia coli K-12 using PCR products. Proc Natl Acad Sci. 2000;97(12):6640–5.View ArticleGoogle Scholar
- Rutherford RJ, Dahl RH, Price RE, Szmidt HL, Benke PI, Mukhopadhyay A, Keasling JD. Functional genomic study of exogenous n-butanol stress in Escherichia coli. Appl Environ Microbiol. 2010;76(6):1935–45.View ArticleGoogle Scholar
- Royce LA, Yoon JM, Chen Y, Rickenbach E, Shanks JV, Jarboe LR. Evolution for exogenous octanoic acid tolerance improves carboxylic acid production and membrane integrity. Metab Eng. 2015;29:180–8.View ArticleGoogle Scholar
- Zingaro KA, Papoutsakis ET. GroESL overexpression imparts Escherichia coli tolerance to i-,n-, and 2-butanol, 1,2,4-butanetriol and ethanol with complex and unpredictable patterns. Metab Eng. 2013;15:196–205.View ArticleGoogle Scholar
- Li J, Zhao JB, Zhao M, Yang YL, Jiang WH, Yang S. Screening and characterization of butanol-tolerant microorganisms. Lett Appl Microbiol. 2010;50(4):373–9.View ArticleGoogle Scholar
- Liu S, Qureshi N. How microbes tolerate ethanol and butanol. Nat Biotechnol. 2009;26(3–4):117–21.Google Scholar
- Luhe AL, Tan L, Wu J, Zhao H. Increase of ethanol tolerance of Saccharomyces cerevisiae by error-prone whole genome amplification. Biotechnol Lett. 2011;33(5):1007–11.View ArticleGoogle Scholar
- Duval V, Lister IM. MarA, SoxS and Rob of Escherichia coli—global regulators of multidrug resistance, virulence and stress response. Int J Biotechnol Wellness Ind. 2013;2(3):101–24.PubMedPubMed CentralGoogle Scholar
- 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(6):2302–7.PubMedPubMed CentralGoogle Scholar
- Sharma P, Haycocks JRJ, Middlemiss AD, Kettles RA, Sellars LE, Ricci V, Piddock LJV, Grainger DC. The multiple antibiotic resistance operon of enteric bacteria controls DNA repair and outer membrane integrity. Nat Commun. 2017;8(1):1444.View ArticleGoogle Scholar
- Bennik MH, Pomposiello PJ, Thorne DF, Demple B. Defining a rob regulon in Escherichia coli by using transposon mutagenesis. J Bacteriol. 2000;182(13):3794–801.View ArticleGoogle Scholar
- White DG, Goldman JD, Demple B, Levy SB. Role of the acrAB locus in organic solvent tolerance mediated by expression of marA, soxS, or robA in Escherichia coli. J Bacteriol. 1997;179(19):6122–6.View ArticleGoogle Scholar
- Herzberg M, Kaye IK, Peti W, Wood TK. YdgG (TqsA) controls biofilm formation in Escherichia coli K-12 through autoinducer 2 transport. J Bacteriol. 2006;188(2):587–98.View ArticleGoogle Scholar
- Kim YJ, Im SY, Lee JO, Kim OB. Potential swimming motility variation by AcrR in Escherichia coli. J Microbiol Biotechnol. 2016;26(10):1824–8.View ArticleGoogle Scholar
- De Araujo C, Balestrino D, Roth L, Charbonnel N, Forestier C. Quorum sensing affects biofilm formation through lipopolysaccharide synthesis in Klebsiella pneumoniae. Res Microniol. 2010;161(7):595–603.View ArticleGoogle Scholar
- Shao C, Shang W, Yang Z, Sun Z, Li Y, Guo J, Wang X, Zou D, Wang S, Lei H, et al. LuxS-dependent AI-2 regulates versatile functions in Enterococcus faecalis V583. J Proteome Res. 2012;11(9):4465–75.View ArticleGoogle Scholar
- Zhang B, Li N, Wang ZW, Tang YJ, Chen T, Zhao XM. Inverse metabolic engineering of Bacillus subtilis for xylose utilization based on adaptive evolution and whole-genome sequencing. Appl Microbiol Biotechnol. 2015;99:885–96.View ArticleGoogle Scholar
- Jiang Y, Chen B, Duan C, Sun B, Yang J, Yang S. Multigene editing in the Escherichia coli genome via the CRISPR-Cas9 system. Appl Environ Microbiol. 2015;81(7):2506–14.View ArticleGoogle Scholar
- Lin Z, Xu Z, Li Y, Wang Z, Chen T, Zhao X. Metabolic engineering of Escherichia coli for the production of riboflavin. Microb Cell Fact. 2014;13:104.PubMedPubMed CentralGoogle Scholar
- Zhang Y, Lin Z, Liu Q, Li Y, Wang Z, Ma H, Chen T, Zhao X. Engineering of serine-deamination pathway, Entner-Doudoroff pathway and pyruvate dehydrogenase complex to improve poly (3-hydroxybutyrate) production in Escherichia coli. Microb Cell Fact. 2014;13:172.View ArticleGoogle Scholar