- Open Access
DNA microarray of global transcription factor mutant reveals membrane-related proteins involved in n-butanol tolerance in Escherichia coli
Biotechnology for Biofuelsvolume 9, Article number: 114 (2016)
Escherichia coli has been explored as a platform host strain for biofuels production such as butanol. However, the severe toxicity of butanol is considered to be one major limitation for butanol production from E. coli. The goal of this study is therefore to construct butanol-tolerant E. coli strains and clarify the tolerance mechanisms.
A recombinant E. coli strain harboring σ70 mutation capable of tolerating 2 % (v/v) butanol was isolated by the global transcription machinery engineering (gTME) approach. DNA microarrays were employed to assess the transcriptome profile of butanol-tolerant strain B8. Compared with the wild-type strain, 329 differentially expressed genes (197 up-regulated and 132 down-regulated) (p < 0.05; FC ≥ 2) were identified. These genes are involved in carbohydrate metabolism, energy metabolism, two-component signal transduction system, oxidative stress response, lipid and cell envelope biogenesis and efflux pump.
Several membrane-related proteins were proved to be involved in butanol tolerance of E. coli. Two down-regulated genes, yibT and yghW, were identified to be capable of affecting butanol tolerance by regulating membrane fatty acid composition. Another down-regulated gene ybjC encodes a predicted inner membrane protein. In addition, a number of up-regulated genes, such as gcl and glcF, contribute to supplement metabolic intermediates for glyoxylate and TCA cycles to enhance energy supply. Our results could serve as a practical strategy for the construction of platform E. coli strains as biofuel producer.
Concerns on global energy crisis and environmental problems have prompted the development of renewable biofuels as potential alternatives for replacing traditional fossil fuels. Among biofuels, butanol has attracted much attention due to its higher energy density, miscibility with gasoline and lower corrosivity . Escherichia coli, as an important platform microorganism, has been widely engineered as an alternative host for the production of various biofuels due to its advantages of fast growth and easy genetic manipulation [2, 3]. However, most biofuels are toxic to E. coli, which barely tolerate organic solvents with LogP values lower than 3.4–3.8 . For example, cell growth of E. coli is completely inhibited in the presence of 1 % (v/v) n-butanol . The poor butanol tolerance of E. coli has been a major limitation in the development of butanol-producing strains. Therefore, it is urgently necessary to improve the butanol tolerance of E. coli.
Most industrial biofuel-producing strains were obtained through solvent stress adaptation, genetic and metabolic engineering and traditional mutagenesis [6–9]. However, long-term adaptive evolution and traditional mutagenesis are often time-consuming. Significantly improved microbial tolerance requires a complex and multigenic phenotype. With recent development of direct mapping between the transcriptome and phenotype, strain improvement efforts have been focused on the manipulation of transcription factors. Global transcription machinery engineering (gTME) emerged as a promising strategy and has been widely used to evolve the desired phenotypes in recent years [10, 11]. Several transcription factors, including sigma factor, CRP, Spt15, H-NS and Hha, have been successfully engineered to improve organic solvent tolerance of various microbial strains [12–15]. σ70 is a subunit of RNA polymerase encoded by rpoD, which regulates over 1000 genes in E. coli. Studies showed that σ70 mutations could alter the promoter preferences of RNA polymerase and therefore affect transcriptome at a global level [16–18]. Our group previously isolated an E. coli harboring σ70 mutant C9, which could grow in the presence of 69 % (v/v) of cyclohexane .
Many solvent tolerance-related genes and their mechanisms have been reported in E. coli strains. MarA (encoded by marA) was confirmed to be a transcription factor capable of inducing the expression of mar-sox regulon genes [20, 21]. Disruption of proV and marR genes could increase the n-hexane tolerance of E. coli cells, possibly due to their function in regulation of osmotic pressure . Rutherford and coworkers reported that n-butanol stress response genes are also involved in other stress responses, such as oxidative stress (sodA, sodC and yqhD), heat shock and cell envelope stress (rpoE, clpB, htpG, cpxR and cpxP) . In our previous study, E. coli strain overexpressing mmsB (encoding 3-hydroxyisobutyrate dehydrogenase) exhibited high solvent tolerance by generating more energy to pump out intracellular organic solvent . These studies demonstrate the complexity and diversity of solvent tolerance mechanisms in E. coli.
It has been reported that E. coli strain capable of growing in the presence of 1.2 % (v/v) butanol was obtained by engineering the global transcription factor cyclic AMP receptor protein . Recently, n-butanol tolerance of E. coli was expanded to 2 % (v/v) by engineering an artificial transcription factor combining with controlling membrane-related functions . In this work, we aimed to improve the butanol tolerance of E. coli by a feasible and efficient approach gTME and understand the mechanism between rpoD mutagenesis and the evolved phenotype. An E. coli strain carrying σ70 mutant capable of tolerating as high as 2 % (v/v) butanol was isolated, which is close to the highest butanol tolerance level of E. coli reported so far. DNA microarrays were employed to identify critical genes related to the n-butanol tolerance based on the transcriptome profile of mutant B8 compared with its wild type (WT). Several membrane-related genes (such as yibT, yghW and ybjC) were recognized to affect membrane fatty acid composition or function as a membrane protein. Other genes (such as gcl and glcF) involved in enhancing metabolic intermediates levels for glyoxylate and TCA cycles were also identified.
Results and discussion
Isolation of σ70 mutants
To obtain σ70 mutants with high n-butanol tolerance, random mutagenesis was performed to construct a rpoD mutant library of around 106. In the first round of screening, 483 mutants were selected from the mutagenesis library on agar plate containing 0.5 % (v/v) butanol. They were further inoculated into 24-well plate under butanol pressure. Among them, ten mutants with OD660 of 1.3–2.0 were selected due to their advantageous growth under 0.5 % (v/v) butanol than other mutants (OD660 < 1) (Fig. 1a). Then these mutants were cultured in medium containing 0.5–1.5 % (v/v) n-butanol (0.1 % gradient increasing). One mutant exhibited the highest tolerance (0.76 OD660) at 1.2 % (v/v) n-butanol and was designated as D3, which was confirmed with two amino acid substitutions (I41L, P97Q).
To further improve the butanol tolerance, the second round of random mutagenesis was performed to construct a variant library using D3 as a template. In the second round of screening, eight mutants exhibited higher butanol tolerance (OD660 > 2.0) than others (Fig. 1b). In further screening under 0.5–1.5 % (v/v) butanol (0.1 % gradient increasing), one mutant B8 with three substitutions (I41L, E57D, P97Q) exhibited excellent growth in the presence of 1.2 % (v/v) butanol, reaching 1.432 OD660 compared with 0.778 of D3 and 0.247 of WT (Additional file 1: Figure S1). Then, mutants B8, D3 and WT were further cultured under higher n-butanol concentration (1.2–2.2 %) to assess their maximum butanol tolerance (Fig. 2a). The cell growth of WT was completely inhibited when butanol reached over 1.2 % (v/v), and mutant D3 could hardly grow in the presence of 1.6 % (v/v) butanol. Nevertheless, for mutant B8, OD660 reached 0.372 after 8 h incubation with 2.0 % (v/v) butanol, suggesting a much higher butanol tolerance than D3 and WT, while no cell growth was observed for B8 at over 2.0 % (v/v) n-butanol. Finally, B8 was selected as the best candidate for further study.
Characterization of butanol-tolerant mutants
Solvent shock experiments
To further characterize the tolerance of B8, solvent shock experiments were performed in medium containing higher concentration of n-butanol. As shown in Fig. 2b, the growth of B8 and WT has no significant difference in the absence of butanol. After treatment with 3 % (v/v) butanol, WT showed growth only at tenfold dilution, while B8 could grow well at 103-fold dilution on LB agar plate, indicating B8 has significantly higher butanol tolerance than WT.
The cell morphology of microorganisms can adapt to the environmental changes . Both WT and B8 grew in the absence or presence of butanol. In transmission electron microscope observation (Fig. 3), mutant B8 cells were significantly longer than WT cells after butanol treatment. In the absence of n-butanol, the cell size was (1.72 ± 0.17) × 0.73 μm for WT (Fig. 3a) and (1.89 ± 0.22) × 0.87 μm for B8 (Fig. 3c), respectively. After 0.8 % (v/v) n-butanol treatment, the cell size of WT and B8 shifted to (2.18 ± 0.25) × 0.53 μm (Fig. 3b) and (2.94 ± 0.39) × 0.52 μm, respectively (Fig. 3d), based on the measurement of around 100 cells. The changes in cell size were confirmed to be statistically significant (p < 0.05). Although the cell morphology was observed based on 100 cells per strain, the cell sizes could also be affected by reasons unrelated to the tolerance effects. Notably, B8 cells were longer than WT when subjected to n-butanol stress. A similar observation on the cell size of ethanol-adapted Saccharomyces cerevisiae has also been reported . Neumann and coworkers reported that the microorganism cells of larger size could be more advantageous over the smaller ones under stress conditions, because the ratio of surface area to volume (S:V = 4πr2:4/3πr3) of larger cells is relatively lower than that of the smaller ones . Similarly, B8 cells with smaller S:V value exhibited greater butanol tolerance than WT. Meanwhile, it was observed that the cytoplasm of B8 shrank in the presence of butanol. Occasionally, invaginated bodies (Fig. 3d) appeared in B8 and the inner membrane of B8 cells was not broken or leaky under butanol treatment. The huddling cytoplasm in the inner membrane could be a self-protection mechanism which protects the cells from damage due to the solvent . Aono and coworkers also reported a similar phenomenon in E. coli cell membrane or cytoplasm in the presence of n-hexane or cyclohexane . Overall, butanol-tolerant mutant B8 demonstrated a stronger stress response capability. Our results suggest that microbial cells could adapt to solvent stress via morphology change.
DNA microarrays and data analysis
Organic solvent tolerance-related genes are usually correlated in a precise regulatory network. The global gene expression profile of mutant B8 and WT was analyzed by DNA microarray. Data were analyzed by Gene Spring Software (Santa Clara, CA, USA) to identify differentially expressed genes. For correlation of gene expression difference under solvent challenge, an over twofold change in gene expression was required with a p value <0.05 . The results show that 329 genes (including 197 up-regulated and 132 down-regulated) exhibited differential expression (p < 0.05; FC ≥ 2) between B8 and WT under 0.8 % (v/v) butanol stress (Additional file 2: Figure S2). Supporting information on the detailed description of differentially expressed genes is provided in Additional file 3: Table S3. Furthermore, genes showing significant difference in transcription level were selected for clustering analysis, which helps to understand the relationships and discrepancy of samples more comprehensively and intuitively (Additional file 4: Figure S3). The same types of genes were gathered in a cluster with similar biological functions.
To identify the function of differentially expressed genes, their biological pathways corresponding to their functional class were summarized by KEGG (Table 1). These genes mainly involve in cell envelope biogenesis, carbohydrate metabolism, energy metabolism, nucleotide metabolism and two-component signal transduction system. As shown in Table 1, most genes (such as ybbQ, glxR, hyi, gcl, glcB, glcD, glcF) involved in glyoxylate and dicarboxylate metabolism pathway were up-regulated for over fivefold, and other genes were also enhanced by more than twofold. The glyoxylate and dicarboxylate metabolism pathway could generate ATP, and its intermediate metabolite NADH could assist in regulating proton gradient and membrane potential. This might be a key factor for improving butanol tolerance, because energy supply is essential for overcoming butanol pressure . In addition, genes related to the ABC transporter pathway were also significantly up-regulated by 2.18- to 5.63-fold. ABC transporter systems consist of different transmembrane protein components and share a common ATP-binding site. ABC transporters play roles in bacterial virulence, cell growth and development, and survival under various environments [33, 34]. Based on our microarray data, ABC transporters annotated as molybdate/arginine/ferric/phosphate type transporters were identified (Table 1), which are responsible for transporting molybdate/arginine/ferric/phosphate substrates across biological membranes. Ferric is important for the regulation of intracellular redox respiratory system, while molybdate/arginine/phosphate substrates are essential for the cellular anabolism and pH homeostasis [35–38]. Therefore, these ABC transporters could affect the adaptation of E. coli strains under severe environments, such as solvent challenge.
Identification of genes associated with butanol tolerance
Analysis of membrane-related down-regulated genes
Many genes have been confirmed to be related to the organic solvent tolerance of the E. coli strain. For example, overexpression of marA could enhance the function of AcrAB-TolC efflux pump, so that the toxic substances could be extruded more efficiently [21, 39]. Escherichia coli mutant △lon (cell envelop related gene) showed enhanced solvent tolerance level . An E. coli mutant with gene disruptions in both proV and marR showed increased solvent tolerance due to their functions in regulation of osmotic pressure . In this study, seven down-regulated genes (yibT, yghW, ymgI, yhcN, yrbL, ECs4086, ybjC) were selected for further study. They were rarely investigated previously and exhibited over sixfold changes in microarray analysis. All seven genes were annotated as predicted protein or hypothetical protein. The down-regulation folds of these genes are 46.72, 13.00, 16.02, 14.74, 14.07, 12.91 and 8.73, respectively. The knockout strains were constructed for further investigation. Among them, △yibT, △yghW and △ybjC exhibited higher n-butanol tolerance, while other knockouts showed similar growth to the control (Fig. 4). So far, there has been no report on the solvent tolerance-related functions of yghW and yibT. It was, however, noticed that these two predicted proteins are closely related to membrane proteins in stitch networks (http://www.stitch.embl.de/). Stitch is a database which contains interaction information for connected proteins and chemicals. It allows querying by genes name and metabolic pathways . In stitch networks, gnsA and gnsB are predicted regulators of phosphatidylethanolamine synthesis and unsaturated fatty acids regulatory proteins, respectively, which are both linked to yibT. Gene yghW is linked to the inner membrane protein encoding genes ybjO and ybjM, as well as genes encoding the predicted lipoprotein or conserved protein in stitch networks. It is therefore presumed that modulation of membrane fatty acid compositions is one possible defense mechanism of yibT and yghW. It has been reported that the properties of membrane fatty acids, such as fatty acid chain length, branching pattern and unsaturation degree of fatty acids, could change when exposed to environment stress [42–45]. In this study, phospholipids of △yghW and △yibT were extracted and analyzed. Our results show that the main components of membrane fatty acids are C16:1, C16:0, C18:0 and C18:1, which account for over 70 % of the total fatty acids in all strains (Table 2) and are responsible for the integrity and fluidity of the membrane . It was noticed that the proportion of unsaturated fatty acid (UFA) in total fatty acids was increased in both knockouts, especially palmitoleic acid (C16:1) and oleic acid (C18:1). Oleic acid (C18:1) has been regarded as the most important UFA in counteracting the toxic effects of the solvent by modulating plasma membrane fluidity . Additionally, palmitoleic acid (C16:1) could influence the rigidity and integrity of membrane lipid bilayer [46, 48]. It is speculated that a higher proportion of oleic acid and palmitoleic acid in △yghW and △yibT contributes to a lower membrane fluidity, higher rigidity and membrane integrity, which might be a compensatory advantage of the membrane challenged by butanol. As a result of changes in membrane fatty acid composition, other physicochemical properties of the membrane, such as proton permeability and lipid–protein interactions could also be affected [49, 50]. Our findings suggest that the butanol tolerance mechanism of △yghW and △yibT are related to membrane fatty acid composition.
Moreover, the hydrophobicity, acidity and alkalinity of cell surface are important properties related to the solvent tolerance [51–53]. Using the MATS method, the adsorptions of three strains (△yghW, △yibT and JM109) to organic solvents with different hydrophobicities such as chloroform, hexadecane, ethyl acetate and decane were determined (Additional file 5: Figure S4). Our results show that the adsorptions of △yghW, △yibT and JM109 to chloroform were 1.89-, 1.86- and 1.61-fold of those to hexadecane, indicating that the alkaline strength of the knockout strains was higher than that of WT. The alkaline strength of the cell surface has been reported to be proportional to the adsorption ratio of chloroform adsorption to hexadecane . Similarly, the adsorptions of △yghW, △yibT and JM109 to ethyl acetate were 1.81-, 1.79- and 1.65-fold of those to decane. It has also been reported that the acidity strength of cell surface was proportional to the adsorption ratio of ethyl acetate to decane . It is supposed that the increase of surface acidity and alkalinity is due to the presence of proteins and charged chemical groups on the cell surface, such as PO4 2− and COO− , which may assist strains to counteract extracellular toxic substances. Additionally, the adhesion to hexadecane and decane reflects the hydrophobicity of the cell surface. In Additional file 5: Figure S4, the adhesion of △yghW and △yibT cells to hexadecane and decane was weaker compared with that of the control, suggesting that the surface hydrophobicity of the knockouts was higher than that of the control. The contact angle measurement (CAM) was also conducted to determine the surface hydrophobicity. Our results show that the contact angle of the control, △yghW and △yibT were 21.17 ± 1.78°, 36.25 ± 2.13° and 34.65 ± 2.04°, respectively (Additional file 6: Figure S5). Consistent with the results of MATS, CAM assay represents a higher surface hydrophobicity of △yghW and △yibT. Other studies suggest that the expansion of surface hydrophobic region could promote the interactions between the phospholipids and embedded proteins and bonds between cation and electronegative phospholipids in the membrane [55, 56]. It is therefore rational to conjecture that the cell surface of △yghW and △yibT is less permeable, which could help to prevent the intruding of toxicity compounds.
The gene (ybjC) encoding a predicted inner membrane protein was also investigated. To verify the location of YbjC and its tolerance-related mechanism, YbjC–GFP (green fluorescent protein) fusion protein was constructed. As shown in Fig. 5a, a clear fluorescence signal surrounding the membrane region was detected in E. coli cells expressing fusion protein YbjC–GFP. Differently, the entire cells were filled with green fluorescence when GFP was expressed alone (Fig. 5b). Here, a known membrane protein YidC fused with GFP was expressed as the positive control (Fig. 5c), and the microscopy pattern of YbjC–GFP fusion was similar to YidC–GFP. This further suggests that YbjC is a membrane protein. Unfortunately, the function of YidC related to the butanol tolerance mechanism is not clear yet.
Analysis of carboxylic acid metabolism-related up-regulated genes
Gene cluster glc, encoding enzymes involved in glyoxylate and dicarboxylate metabolism, was up-regulated significantly in this study. For example, genes glcF, glcA, glcD, glcG, glcB and gcl were up-regulated by 22.60-, 18.17-, 12.78-, 12.75- 7.86- and 10.60-fold, respectively. The above genes and the glc gene cluster were selected for overexpression to investigate their functional relevance to butanol tolerance. As shown in Additional file 7: Figure S6, glcA overexpression strain could barely grow even without butanol. For glcB, glcD and glcG overexpression strains, similar growth with control strain (JM109 carrying empty plasmid pQE80L) was observed in the presence and absence of butanol. Strains overexpressing glcF and gcl showed improved cell growth compared with the control. Furthermore, glcDEFGBA cluster co-expression strains also exhibited better growth than the control with or without butanol. The possible explanation is that glcA encodes glycolate transporter, and its overexpression could result in metabolism imbalance of carbon source and energy in cells . glcB encodes malate synthase G which is not essential for growth. glcD and glcG encode different subunits of glycolate oxidase, which have no significant effect on the oxidase activity . Gene glcF encodes iron–sulfur subunit of glycolate oxidase that catalyzes the conversion of glycolate to glyoxylate, and gcl encodes glyoxylate carboligase that transforms glyoxylate into phosphoglycerate [59, 60]. One possible explanation is that these two enzymes, glycolate oxidase and glyoxylate carboligase, directly improve the TCA cycle via modulating glyoxylate and pyruvate, which replenish metabolic intermediates for the TCA cycle (Additional file 8: Figure S7). In this study, rpoD mutagenesis approach conduced to enhanced OST phenotype by regulating the expression level of hundreds of genes. Our results also suggest that the co-expression of glc gene cluster is more favorable for cell growth than individual gene overexpression, which might have minor or negative effect on butanol tolerance. To elucidate the possible mechanism, the concentrations of carboxylic acids in TCA cycle were monitored in glcF and gcl overexpression strains. For both strains, pyruvate concentrations were increased in the presence of butanol, whereas higher pyruvate levels were also observed without butanol (Additional file 9: Table S4). Moreover, the fumarate concentration was also increased in the overexpression strains without butanol, but declined slightly when 0.8 % (v/v) butanol was added. This phenomenon may be attributed to the inhibition of certain TCA cycle enzymes (such as succinate dehydrogenase) that are sensitive to the oxidative stress induced by butanol. A similar phenomenon was reported by Fu and coworkers . Although the elevated pyruvate level in gcl and glcF overexpression strains does not seem to be directly related to butanol stress, it is assumed that the up-regulation of these two genes is favorable for the cell growth by replenishing TCA metabolic intermediates.
In summary, we successfully isolated an E. coli strain harboring rpoD mutant B8 with 2 % (v/v) butanol tolerance using global transcriptional machinery engineering approach. Based on DNA microarrays results, 329 genes (including 197 up-regulated and 132 down-regulated) showed over twofold difference in the expression level compared with WT after butanol treatment. These genes are mainly involved in the metabolic pathways including ABC pump, ascorbate and aldarate metabolism, energy metabolism, two-component signal transduction system and amino acid metabolism. Tolerance mechanisms of several critical genes have been elucidated. Among them, down-regulated genes yghW and yibT were proved to improve n-butanol tolerance due to their regulatory roles in membrane fatty acids composition. YbjC was confirmed to be a membrane protein, while up-regulated genes gcl and glcF could replenish TCA cycle metabolic intermediates to improve cell growth and metabolism. These results could provide the potential approach for the construction of E. coli strain as a bio-butanol producer.
Strains, plasmids and culture conditions
Escherichia coli JM109 was used as the host strain. Gene deletion strains were generated by Red-mediated recombination approach and overexpression strains were generated using pQE80L as an expression vector. Strains and plasmids are listed in Additional file 10: Table S1. Primers used in this study are listed in Additional file 11: Table S2. Plasmid pQE80L was purchased from Qiagen GmbH (Hilden, Germany). pHACM-rpoD WT was presented as a kind gift of Dr. Huimin Yu (Tsinghua University, China). Restriction enzymes and PrimeSTAR®HS DNA Polymerase were purchased from Takara (Tokyo, Japan). The dam-methylated DNA-specific restriction enzyme DpnI was purchased from New England Biolabs (Ipswich, MA, USA). GeneMorph II Random Mutagenesis Kit was obtained from Stratagene (La Jolla, CA, USA).
All strains were grown in Luria–Bertani (LB) medium (tryptone 10 g/L, yeast extract 5 g/L, NaCl 10 g/L) at 37 °C, 120 rpm. Butanol was added as specified in each experiment. When necessary, antibiotics chloramphenicol (34 μg/mL), ampicillin (100 μg/mL) and kanamycin (50 μg/mL) were added to the media. For gene overexpression, 0.2 mM IPTG was added to the medium at around 0.3 OD660.
Construction of random mutagenesis library
Using plasmid pHACM-rpoD WT as the template, random mutagenesis was performed by the GenemorphII Random Mutagenesis kit (Stratagene) with a mutation rate of approximately 4.5–9 mutations/kb. The error-prone PCR program was set as follows: 5 min at 95 °C, 30 cycles of 95 °C for 30 s, 57 °C for 1 min, followed by 72 °C for 2 min, and 10 min at 72 °C. After running the whole plasmid PCR, the PCR product mixture was digested with DpnI and then transformed into E. coli JM109. Escherichia coli transformants were plated on LB agar plates containing 34 μg/mL chloramphenicol and incubated at 37 °C. Then the colonies were scraped off to create a σ70 mutant library for further butanol-tolerant phenotype selection. The total mutant library size was approximately 106.
Phenotype selection of n-butanol-tolerant mutants
The first round of screening
First, 483 colonies were selected from the σ70 mutant library on agar plate containing 0.5 % butanol (v/v). Then, these colonies were inoculated into 24-well plates with LB/Cm liquid medium. Butanol (0.5 %, v/v) was added into the culture when OD660 reached 0.2, and the cells were further incubated for 8 h at 120 rpm and 37 °C. Cell density was measured, and ten mutants with OD660 above 1.0 were selected. Then, the selected mutants were cultured in a medium containing higher concentrations (0.5–1.5 %, v/v) of n-butanol by 0.1 % (v/v) gradient, until a strain exhibiting the highest tolerance was obtained. The plasmid harboring the rpoD mutant was then sequenced and designated as D3.
The second round of screening
To further improve the n-butanol tolerance, σ70 mutant D3 was used as a template to construct a random mutagenesis library for a second round of screening. The screening method was the same as that in the first round. The best mutant B8 was selected and evaluated under higher n-butanol concentrations (0.1 % gradient increasing from 1.2 to 2.2 %) to determine its maximum tolerant level.
Solvent shock treatment
The overnight cell culture was inoculated (1 %) into a fresh LB medium. When OD660 reached 0.8, 3 % (v/v) n-butanol was added into the culture. After incubation at 37 °C for 1.5 h, the cultures were serially diluted for 105, 104, 103, 102 and 10-fold with aseptic water. Then, 10 μL of the each diluted culture was spotted onto LB/agar plates and further incubated at 37 °C overnight.
The σ70 mutant B8 and WT strains were cultured overnight and inoculated (1 %) into fresh LB/Cm liquid medium for incubation at 37 °C and 120 rpm for 8 h with or without 0.8 % (v/v) n-butanol. The cells were diluted and spread on LB/Cm agar plates. Single colonies were treated as described in the literature and observed using transmission electron microscope . Briefly, single colonies were picked and fixed by immersion in 2.5 % glutaraldehyde at 4 °C for 3 h. The cell suspension was mixed once every half-hour. Then, cells were washed four times with 0.1 M phosphate buffer (pH 7.2), and the samples were diluted for cell morphology observation under Hitachi-H7650 transmission electron microscopy (Japan). The average cellular size of WT and B8 was counted on the electron microscope based on around 100 cells.
Escherichia coli strains harboring σ70 mutant B8 and WT were cultured overnight and inoculated (1 %) into fresh medium. n-Butanol (0.8 %, v/v) was added at 0.8 OD660 and further incubated for 1.5 h. Cells were harvested by centrifugation (8800 g, 4 °C). Total RNA was extracted using Qiagen RNeasy kit (Hilden, Germany) following the manufacturer’s instructions. Qualified total RNA was further purified by Qiagen RNeasy mini kit and RNase-Free DNase Set. Three biological replicates of RNA samples were stored in dry ice and subjected to further DNA microarray analysis. The microarray service was provided by Shanghai Biotechnology Co., Ltd. (Shanghai, China) using Agilent SurePrint E. coli 8 × 15 K slides, and the quality and integrity of RNA were examined before analysis.
Slides were scanned by Agilent Microarray Scanner (Santa Clara, CA, USA), and data were extracted with Agilent Feature Extraction software 10.7. Raw data were normalized by Quantile algorithm, Gene Spring Software 11.0 (Santa Clara, CA, USA). Differentially expressed genes were identified using the rank product method . The MeV (TM4) software was used for clustering and other expression profile analysis . The related metabolic pathway of differentially expressed genes was analyzed using Kyoto Encyclopedia of Gene and Genomics (KEGG) database . The microarray data have been deposited at the gene expression omnibus (GEO) under the accession number GSE79305.
Assay of n-butanol sensitivity of knockout and overexpression strains
Gene knockout (for down-regulated genes) and overexpression (for up-regulated genes) strains were constructed based on the microarrays results. Spot assay as a tolerance confirmation method was conducted for knockout strains. n-Butanol (0.8 %, v/v) was added when the cell density reached 0.8 OD660, then the cells were further cultured for 1.5 h. Cell culture was serially diluted at a tenfold gradient, and 10 μL of the diluted culture was spotted onto LB/agar plates for incubation at 37 °C overnight. For overexpression strains, the overnight culture was inoculated (1 %) into fresh medium and cultured at 37 °C and 120 rpm, and 0.2 mM IPTG was added when OD660 reached 0.3. Then, 0.8 % (v/v) butanol was added at 0.8 OD660 for further incubation at 30 °C for 8 h. The growth of the strains was monitored by measuring cell density at OD660.
Lipid extraction was performed as described by Bligh . Briefly, the stationary phase cells of E. coli JM109, △yibT and △yghW were collected by centrifugation and washed with 10 mM, pH 7.4 sodium phosphate buffer. The cell pellets (0.3 g wet weight) were suspended in a mixture consisting of buffer, methanol and chloroform in the proportion of 0.8:2:1 (v/v) and incubated at room temperature for 2 h, with brief shaking every half-hour. Then, the extract was centrifuged at 4 °C after diluting with methanol and chloroform (1:1, v/v). The lower chloroform phase (containing lipids) was added into the methanol solution containing 5 % (v/v) H2SO4. After 2 h of methylation at 70 °C, the mixture was cooled to room temperature and extracted three times with pentane. Then, the samples were air dried to remove pentane. Finally, the samples were reconstituted with n-hexane for GC–MS analysis.
Determination of cell surface hydrophobicity
MATS (microbial adhesion to solvents) analysis was performed following the MATH method developed by Bellon-Fontaine with slight modification [53, 65]. Briefly, the strains were cultured in LB medium, and butanol (0.8 %, v/v) was added when OD660 reached 0.8. Cells were further incubated for 12 h at 37 °C and were collected by centrifugation (8800 g and 10 min at 4 °C), then washed with 100 mM, pH 6.0 potassium phosphate buffer and centrifuged again. The cell concentration was then adjusted to around 1.0 OD400 (A0) using the same potassium phosphate buffer. The cell suspension (4.8 mL) was mixed with 0.8 mL of chloroform, hexadecane, ethyl acetate and decane, respectively. The two-phase mixture was mixed by vortexing for 90 s and then incubated statically at room temperature for 15 min. The aqueous phase was removed and OD400 was measured (A) to calculate the adhesion ratio [Adhesion % = (1 − A/A0) × 100 %]. The above experiment was repeated for three times.
Contact angle measurements (CAM) was performed to measure the surface hydrophobicity . Briefly, strains were cultured in the LB medium, and butanol (0.8 %, v/v) was added when OD660 reached 0.8. Cells were further incubated for 12 h at 37 °C and were collected by centrifugation at 8800 g for 10 min at 4 °C and washed twice by 0.85 % saline. Then, the concentration of cell suspension was adjusted to 50 mg wet cells/mL, and 10 mL cells suspension was filtered through polyvinylidene difluoride membrane (0.22 μm, 50 mm in diameter) under vacuum. For each strain, three biological replicates were performed and measured independently. The contact angle was measured in three phases: the bacterial lawn, n-tetradecane and a droplet of distilled water using a contact angle meter (Dataphysics, Germany).
Fatty acid component analysis
The fatty acid component analysis was performed as previously described [66, 67]. The fatty acid profile was assessed with Trace1310 GC equipped with TSQ8000 MSD (Thermo Fisher Scientific, Massachusetts, USA) and an HP-5MS methylpolysiloxane phase column (30 m × 0.25 mm × 0.25 μm). The GC conditions are as follows: the initial temperature of 60 °C for 2 min, 8 °C/min to 150 °C, then to 250 °C at a rate of 3 °C/min, finally at 10 °C/min to 280 °C and held for 5 min. Fatty acids and other volatile compounds were identified by mass spectral library search.
Quantification of pyruvate and fumarate
Overexpression strains and control (JM109 carrying empty plasmid pQE80L) were cultured overnight and inoculated into fresh LBG medium (1 %, v/v). IPTG (0.2 mM) was added at 0.3 OD660, cells were further cultured for 8 h with or without 0.8 % n-butanol. The cell culture was centrifuged and the supernatant was subjected to HPLC analysis using Agilent 1260 (Palo Alto, USA) equipped with an Agilent Hi-Plex H PL1170-6830 column (7.7 × 300 mm, 8 µm) at 55 °C. 5 mM H2SO4 was used as an eluent at a flow rate of 0.6 mL/min.
global transcription machinery engineering
Kyoto Encyclopedia of Gene and Genomics
unsaturated fatty acid
saturated fatty acid
microbial adhesion to solvents
contact angle measurement
green fluorescent protein
Bruno TJ, Wolk A, Naydich A. Composition-explicit distillation curves for mixtures of gasoline with four-carbon alcohols (butanols). Energy Fuels. 2009;23:2295–306.
Atsumi S, Cann AF, Connor MR, Shen CR, Smith KM, Brynildsen MP, Chou KJ, Hanai T, Liao JC. Metabolic engineering of Escherichia coli for 1-butanol production. Metab Eng. 2008;10(6):305–11.
Baez A, Cho KM, Liao JC. High-flux isobutanol production using engineered Escherichia coli: a bioreactor study with in situ product removal. Appl Microbiol Biotechnol. 2011;90(5):1681–90.
Inoue A, Horikoshi K. Estimation of solvent-tolerance of Bacteria by the solvent parameter Log-P. J Ferment Bioeng. 1991;71(3):194–6.
Atsumi S, Hanai T, Liao JC. Non-fermentative pathways for synthesis of branched-chain higher alcohols as biofuels. Nature. 2008;451(7174):86–9.
Bajwa PK, Shireen T, D’Aoust F, Pinel D, Martin VJ, Trevors JT, Lee H. Mutants of the pentose-fermenting yeast Pichia stipitis with improved tolerance to inhibitors in hardwood spent sulfite liquor. Biotechnol Bioeng. 2009;104(5):892–900.
Geddes CC, Mullinnix MT, Nieves IU, Peterson JJ, Hoffman RW, York SW, Yomano LP, Miller EN, Shanmugam KT, Ingram LO. Simplified process for ethanol production from sugarcane bagasse using hydrolysate-resistant Escherichia coli strain MM160. Bioresour Technol. 2011;102(3):2702–11.
Shi DJ, Wang CL, Wang KM. Genome shuffling to improve thermotolerance, ethanol tolerance and ethanol productivity of Saccharomyces cerevisiae. J Ind Microbiol Biotechnol. 2009;36(1):139–47.
Tomas CA, Beamish J, Papoutsakis ET. Transcriptional analysis of butanol stress and tolerance in Clostridium acetobutylicum. J Bacteriol. 2004;186(7):2006–18.
Alper H, Moxley J, Nevoigt E, Fink GR, Stephanopoulos G. Engineering yeast transcription machinery for improved ethanol tolerance and production. Science. 2006;314(5805):1565–8.
Alper H, Stephanopoulos G. Global transcription machinery engineering: a new approach for improving cellular phenotype. Metab Eng. 2007;9(3):258–67.
Klein-Marcuschamer D, Stephanopoulos G. Assessing the potential of mutational strategies to elicit new phenotypes in industrial strains. Proc Natl Acad Sci USA. 2008;105(7):2319–24.
Zhang H, Chong H, Ching CB, Song H, Jiang R. Engineering global transcription factor cyclic AMP receptor protein of Escherichia coli for improved 1-butanol tolerance. Appl Microbiol Biotechnol. 2012;94(4):1107–17.
Liu H, Yan M, Lai C, Xu L, Ouyang P. gTME for improved xylose fermentation of Saccharomyces cerevisiae. Appl Biochem Biotechnol. 2010;160(2):574–82.
Hong SH, Lee J, Wood TK. Engineering global regulator Hha of Escherichia coli to control biofilm dispersal. Microb Biotechnol. 2010;3(6):717–28.
Malhotra A, Severinova E, Darst SA. Crystal structure of a sigma 70 subunit fragment from E. coli RNA polymerase. Cell. 1996;87(1):127–36.
Owens JT, Miyake R, Murakami K, Chmura AJ, Fujita N, Ishihama A, Meares CF. Mapping the sigma70 subunit contact sites on Escherichia coli RNA polymerase with a sigma70-conjugated chemical protease. Proc Natl Acad Sci USA. 1998;95(11):6021–6.
Yang X, Ma C, Lewis P. A vector system that allows simple generation of mutant Escherichia coli RNA polymerase. Plasmid. 2014;75:37–41.
Zhang F, Qian X, Si H, Xu G, Han R, Ni Y. Significantly improved solvent tolerance of Escherichia coli by global transcription machinery engineering. Microb Cell Fact. 2015;14(1):175.
Hayashi S, Aono R, Hanai T, Mori H, Kobayashi T, Honda H. Analysis of organic solvent tolerance in Escherichia coli using gene expression profiles from DNA microarrays. J Biosci Bioeng. 2003;95(4):379–83.
Keeney D, Ruzin A, McAleese F, Murphy E, Bradford PA. MarA-mediated overexpression of the AcrAB efflux pump results in decreased susceptibility to tigecycline in Escherichia coli. J Antimicrob Chemother. 2008;61(1):46–53.
Doukyu N, Ishikawa K, Watanabe R, Ogino H. Improvement in organic solvent tolerance by double disruptions of proV and marR genes in Escherichia coli. J Appl Microbiol. 2012;112(3):464–74.
Rutherford BJ, 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.
Ni Y, Song L, Qian X, Sun Z. Proteomic analysis of Pseudomonas putida reveals an organic solvent tolerance-related gene mmsB. PLoS ONE. 2013;8(2):e55858.
le Bui M, Lee JY, Geraldi A, Rahman Z, Lee JH, Kim SC. Improved n-butanol tolerance in Escherichia coli by controlling membrane related functions. J Biotechnol. 2015;204:33–44.
Dinh TN, Nagahisa K, Hirasawa T, Furusawa C, Shimizu H. Adaptation of Saccharomyces cerevisiae cells to high ethanol concentration and changes in fatty acid composition of membrane and cell size. PLoS ONE. 2008;3(7):e2623.
Kubota S, Takeo I, Kume K, Kanai M, Shitamukai A, Mizunuma M, Miyakawa T, Shimoi H, Iefuji H, Hirata D. Effect of ethanol on cell growth of budding yeast: genes that are important for cell growth in the presence of ethanol. Biosci Biotechnol Biochem. 2004;68(4):968–72.
Neumann G, Veeranagouda Y, Karegoudar TB, Sahin O, Mausezahl I, Kabelitz N, Kappelmeyer U, Heipieper HJ. Cells of Pseudomonas putida and Enterobacter sp. adapt to toxic organic compounds by increasing their size. Extremophiles. 2005;9(2):163–8.
Aono R. Improvement of organic solvent tolerance level of Escherichia coli by overexpression of stress-responsive genes. Extremophiles. 1998;2(3):239–48.
Aono R, Kobayashi H, Joblin KN, Horikoshi K. Effects of organic solvents on growth of Escherichia coli K-12. Biosci Biotech Bioch. 1994;58(11):2009–14.
Yang J, Hong Y, Yuan C, Fu Z, Shi Y, Zhang M, Shen L, Han Y, Zhu C, Li H, Lu K, Liu J, Feng X, Lin J. Microarray analysis of gene expression profiles of Schistosoma japonicum derived from less-susceptible host water buffalo and susceptible host goat. PLoS ONE. 2013;8(8):e70367.
Fu F, Cheng VW, Wu Y, Tang Y, Weiner JH, Li L. Comparative proteomic and metabolomic analysis of Staphylococcus warneri SG1 cultured in the presence and absence of butanol. J Proteome Res. 2013;12(10):4478–89.
Fath MJ, Kolter R. ABC transporters: bacterial exporters. Microbiol Rev. 1993;57(4):995–1017.
Linton KJ, Higgins CF. The Escherichia coli ATP-binding cassette (ABC) proteins. Mol Microbiol. 1998;28(1):5–13.
Self WT, Grunden AM, Hasona A, et al. Molybdate transport. Res Microbiol. 2001;152(3):311–21.
Wissenbach U, Six S, Bongaerts J, Ternes D, Steinwachs S, Unden G. A third periplasmic transport system for l-arginine in Escherichia coli: molecular characterization of the artPIQMJ genes, arginine binding and transport. Mol Microbiol. 1995;17(4):675–86.
Chenault SS, Earhart CF. Organization of genes encoding membrane proteins of the Escherichia coli ferrienterobactin permease. Mol Microbiol. 1991;5(6):1405–13.
Horinouchi T, Tamaoka K, Furusawa C, Ono N, Suzuki S, Hirasawa T, Yomo T, Shimizu H. Transcriptome analysis of parallel-evolved Escherichia coli strains under ethanol stress. BMC Genom. 2010;11:579.
Mingardon F, Clement C, Hirano K, Nhan M, Luning EG, Chanal A, Mukhopadhyay A. Improving olefin tolerance and production in E. coli using native and evolved AcrB. Biotechnol Bioeng. 2015;112(5):879–88.
Watanabe R, Doukyu N. Improvement of organic solvent tolerance by disruption of the lon gene in Escherichia coli. J Biosci Bioeng. 2014;118(2):139–44.
Kuhn M, von Mering C, Campillos M, Jensen LJ, Bork P. STITCH: interaction networks of chemicals and proteins. Nucleic Acids Res. 2008;36:D684–8.
Archana KM, Ravi R, Anu-Appaiah KA. Correlation between ethanol stress and cellular fatty acid composition of alcohol producing non-Saccharomyces in comparison with Saccharomyces cerevisiae by multivariate techniques. J Food Sci Technol. 2015;52(10):6770–6.
Fozo EM, Quivey RG. Shifts in the membrane fatty acid profile of Streptococcus mutans enhance survival in acidic environments. Appl Environ Microbiol. 2004;70(2):929–36.
Nurlinawati, Vanoirbeek K, Aertsen A, Michiels CW. Role of 1-acyl-sn-glycerol-3-phosphate acyltransferase in psychrotrophy and stress tolerance of Serratia plymuthica RVH1. Res Microbiol. 2015;166(1):28–37.
Ramos JL, Duque E, Gallegos MT, Godoy P, Ramos-Gonzalez MI, Rojas A, Teran W, Segura A. Mechanisms of solvent tolerance in gram-negative bacteria. Annu Rev Microbiol. 2002;56:743–68.
Ramos JL, Duque E, Rodriguez-Herva JJ, Godoy P, Haidour A, Reyes F, Fernandez-Barrero A. Mechanisms for solvent tolerance in bacteria. J Biol Chem. 1997;272(7):3887–90.
You KM, Rosenfield CL, Knipple DC. Ethanol tolerance in the yeast Saccharomyces cerevisiae is dependent on cellular oleic acid content. Appl Environ Microbiol. 2003;69(3):1499–503.
Mannazzu I, Angelozzi D, Belviso S, Budroni M, Farris GA, Goffrini P, Lodi T, Marzona M, Bardi L. Behaviour of Saccharomyces cerevisiae wine strains during adaptation to unfavourable conditions of fermentation on synthetic medium: cell lipid composition, membrane integrity, viability and fermentative activity. Int J Food Microbiol. 2008;121(1):84–91.
Quivey RG Jr, Faustoferri R, Monahan K, Marquis R. Shifts in membrane fatty acid profiles associated with acid adaptation of Streptococcus mutans. FEMS Microbiol Lett. 2000;189(1):89–92.
Sikkema J, De Bont JA, Poolman B. Mechanisms of membrane toxicity of hydrocarbons. Microbiol Rev. 1995;59(2):201–22.
Hamadi F, Latrache H. Comparison of contact angle measurement and microbial adhesion to solvents for assaying electron donor-electron acceptor (acid-base) properties of bacterial surface. Colloids Surf B Biointerfaces. 2008;65(1):134–9.
Aono R, Kobayashi H. Cell surface properties of organic solvent-tolerant mutants of Escherichia coli K-12. Appl Environ Microbiol. 1997;63(9):3637–42.
BellonFontaine MN, Rault J, VanOss CJ. Microbial adhesion to solvents: a novel method to determine the electron-donor/electron-acceptor or Lewis acid-base properties of microbial cells. Colloid Surface B. 1996;7(1–2):47–53.
Merritt K, An YH. Factors influencing bacterial adhesion. Handbook of bacterial adhesion: principles, methods, and applications. Totowa: Humana Press; 2000. p. 53–72.
Guan L, Smirnova IN, Verner G, Nagamoni S, Kaback HR. Manipulating phospholipids for crystallization of a membrane transport protein. Proc Natl Acad Sci USA. 2006;103(6):1723–6.
de Planque MR, Killian JA. Protein-lipid interactions studied with designed transmembrane peptides: role of hydrophobic matching and interfacial anchoring. Mol Membr Biol. 2003;20(4):271–84.
Núñez MF, Kwon O, Wilson TH, Aguilar J, Baldoma L, Lin EC. Transport of L-Lactate, D-Lactate, and glycolate by the LldP and GlcA membrane carriers of Escherichia coli. Biochem Biophys Res Commun. 2002;290(2):824–9.
Pellicer MT, Badia J, Aguilar J, Baldoma L. glc locus of Escherichia coli: characterization of genes encoding the subunits of glycolate oxidase and the glc regulator protein. J Bacteriol. 1996;178(7):2051–9.
Clark DP, Cronan JE. Two-carbon compounds and fatty acids as carbon sources. EcoSal Plus. 2005. doi:10.1128/ecosalplus.3.4.4.
Nemeria N, Binshtein E, Patel H, Balakrishnan A, Vered I, Shaanan B, Barak Z, Chipman D, Jordan F. Glyoxylate carboligase: a unique thiamin diphosphate-dependent enzyme that can cycle between the 4′-aminopyrimidinium and 1′,4′-iminopyrimidine tautomeric forms in the absence of the conserved glutamate. Biochemistry. 2012;51(40):7940–52.
Breitling R, Armengaud P, Amtmann A, Herzyk P. Rank products: a simple, yet powerful, new method to detect differentially regulated genes in replicated microarray experiments. FEBS Lett. 2004;573(1–3):83–92.
Saeed AI, Sharov V, White J, Li J, Liang W, Bhagabati N, Braisted J, Klapa M, Currier T, Thiagarajan M, Sturn A, Snuffin M, Rezantsev A, Popov D, Ryltsov A, Kostukovich E, Borisovsky I, Liu Z, Vinsavich A, Trush V, Quackenbush J. TM4: a free, open-source system for microarray data management and analysis. Biotechniques. 2003;34(2):374–8.
Kanehisa M, Goto S, Furumichi M, Tanabe M, Hirakawa M. KEGG for representation and analysis of molecular networks involving diseases and drugs. Nucleic Acids Res. 2010;38:D355–60.
Bligh EG, Dyer WJ. A rapid method of total lipid extraction and purification. Can J Biochem Physiol. 1959;37(8):911–7.
Rosenberg M. Microbial adhesion to hydrocarbons: twenty-five years of doing MATH. FEMS Microbiol Lett. 2006;262(2):129–34.
Dembitsky VM, Řezanková H, Řezanka T, Hanuš LO. Variability of the fatty acids of the marine green algae belonging to the genus Codium. Biochem Syst Ecol. 2003;31(10):1125–45.
Ni Y, Reye J, Chen RR. lpp deletion as a permeabilization method. Biotechnol Bioeng. 2007;97(6):1347–56.
HS carried out the DNA microarray data analysis, solvent tolerant, membrane lipid analysis, cell surface hydrophobicity and quantification of metabolites, and drafted the manuscript. FZ carried out screening of σ70 mutants, DNA microarray sample preparation, transmission electron microscopy and participated in drafting the manuscript. AW carried out gene knockout/overexpression experiments and analyzed the data. RH participated in the design of the study, analyzed the data and helped to revise the manuscript. GX helped to design the experiments and revise the manuscript. YN conceived the study, designed and coordinated the experiments and critically revised the manuscript. All authors read and approved the final manuscript.
We thank Mr. Kang Wu and Prof. Xiaodong Gao for help with the transmission electron micrographs and fluorescent micrograph analysis.
Availability of supporting data
The microarray data have been deposited at Gene Expression Omnibus (http://www.ncbi.nlm.nih.gov/geo/info/submission.html) under the accession number GSE79305.
The authors declare that they have no competing interests.
We are grateful to the National Natural Science Foundation of China (21276112, 31401634), Natural Science Foundation of Jiangsu Province (BK20150003, BK20140135), Fundamental Research Funds for the Central Universities (JUSRP51409B), the Program of Introducing Talents of Discipline to Universities (111-2-06) and a project funded by the Priority Academic Program Development of Jiangsu Higher Education Institutions for the financial support of this research.