- Research article
- Open Access
Improved n-butanol production via co-expression of membrane-targeted tilapia metallothionein and the clostridial metabolic pathway in Escherichia coli
BMC Biotechnology volume 17, Article number: 36 (2017)
N-Butanol has favorable characteristics for use as either an alternative fuel or platform chemical. Bio-based n-butanol production using microbes is an emerging technology that requires further development. Although bio-industrial microbes such as Escherichia coli have been engineered to produce n-butanol, reactive oxygen species (ROS)-mediated toxicity may limit productivity. Previously, we show that outer-membrane-targeted tilapia metallothionein (OmpC-TMT) is more effective as an ROS scavenger than human and mouse metallothioneins to reduce oxidative stress in the host cell.
The host strain (BUT1-DE) containing the clostridial n-butanol pathway displayed a decreased growth rate and limited n-butanol productivity, likely due to ROS accumulation. The clostridial n-butanol pathway was co-engineered with inducible OmpC-TMT in E. coli (BUT3-DE) for simultaneous ROS removal, and its effect on n-butanol productivity was examined. The ROS scavenging ability of cells overexpressing OmpC-TMT was examined and showed an approximately twofold increase in capacity. The modified strain improved n-butanol productivity to 320 mg/L, whereas the control strain produced only 95.1 mg/L. Transcriptomic analysis revealed three major KEGG pathways that were significantly differentially expressed in the BUT3-DE strain compared with their expression in the BUT1-DE strain, including genes involved in oxidative phosphorylation, fructose and mannose metabolism and glycolysis/gluconeogenesis.
These results indicate that OmpC-TMT can increase n-butanol production by scavenging ROS. The transcriptomic analysis suggested that n-butanol causes quinone malfunction, resulting in oxidative-phosphorylation-related nuo operon downregulation, which would diminish the ability to convert NADH to NAD+ and generate proton motive force. However, fructose and mannose metabolism-related genes (fucA, srlE and srlA) were upregulated, and glycolysis/gluconeogenesis-related genes (pfkB, pgm) were downregulated, which further assisted in regulating NADH/NAD+ redox and preventing additional ATP depletion. These results indicated that more NADH and ATP were required in the n-butanol synthetic pathway. Our study demonstrates a potential approach to increase the robustness of microorganisms and the production of toxic chemicals through the ability to reduce oxidative stress.
n-Butanol has many advantages over ethanol, including a higher energy density due to two extra carbons, and can be used in gasoline engines without modification. n-Butanol is less hygroscopic and volatile than ethanol and has been recently regarded as a more viable transportation biofuel than ethanol [1–3]. n-Butanol is also a versatile platform chemical that can be produced from a variety of biomass sources. In industry, many important chemicals derived from n-butanol are used extensively as solvents or as intermediates in the production of acrylates, n-butyl acetate, esters, and glycol ethers .
Several industrial microbes, including Escherichia coli, have been engineered to produce n-butanol. In the initial stage of engineering E. coli for n-butanol production, the whole n-butanol pathway from Clostridium was transferred to E. coli, including the thil, hbd, crt, bcd, etfA, etfB, and adhe (or adhel) genes, catalyzing the six-step conversion of two molecules of acetyl-CoA into one molecule of n-butanol. However, when the clostridial n-butanol pathway was first transferred to E. coli using plasmids in 2007, the engineered strain produced less than 1 g/L n-butanol (vs. a clostridial n-butanol titer of 10–20 g/L) [5, 6]. The results indicated that engineering efficient n-butanol-producing E. coli is not as simple as expressing several clostridial n-butanol pathway genes. Over the following years, the best heterologous n-butanol-producing strains derived from E. coli were able to produce 14–15 g/L n-butanol [7, 8], providing industrial advantages compared to clostridial strains . However, the toxicity of n-butanol to both natural producers and engineered hosts increases with accumulation [10, 11]. This toxicity complicates the economically efficient production of large titers of n-butanol; therefore, high-titer n-butanol (30 g/L) production currently relies on in situ product removal by continuous gas stripping .
Furthermore, although E. coli can convert sugars into n-butanol at relatively high levels, the cells cannot tolerate 2% (v/v) n-butanol  and produce n-butanol at insufficient levels. Considering the relationship between n-butanol tolerance and n-butanol production in Clostridial strains [12, 13], the toxicity of n-butanol to E. coli can be considered a production bottleneck. In fact, product toxicity is a widespread problem in the production of biorenewables. Therefore, the development of a stress-tolerant host strain for the bio-production of titer-dependent toxic chemicals is important.
To understand the effect of n-butanol toxicity on the host, cell-wide studies have been conducted to obtain a global view of the n-butanol stress-response at the transcript, protein, and metabolite levels. Transcriptomic analysis in Clostridium acetobutylicum indicated that the primary response involved the accumulation of transcripts encoding chaperones, proteases, and other heat shock-related proteins .
In E. coli, several transcriptional analyses have been performed to investigate the stress caused by alcohols, including ethanol, n-butanol, and isobutanol [15–17]. Additionally, observations from fluorescent dye-staining indicated a large increase in reactive oxygen species (ROS) during n-butanol stress . Increased oxidative stress is a common cellular response to extracellular xenobiotics, which may mediate macromolecular damage. These free radicals can directly attack the membrane via lipid peroxidation or cause DNA mutations, protein misfolding and fragmentation, and apoptosis [18–20].
To decrease ROS-induced oxidative damage, microorganisms synthesize many antioxidant enzymes, including catalases, superoxide dismutases and glutathione peroxidase [21, 22]. Recently, metallothioneins (MTs), beneficial antioxidant enzymes that occur widely in mammals, plants and fungi, have been identified . MTs are heat-stable, low-molecular-weight and cysteine-rich intracellular proteins [23–25]. MTs also act as a defense system against oxidative stress through their ROS scavenging ability . Furthermore, purified tilapia MT (TMT) has been shown to outperform glutathione (GSH) in scavenging both 2-diphenyl-1-picrylhydrazyl (DPPH●) and 2,2-azinobis (3-ethylbenzothiazoline- 6-sulfonic acid) diammonium salt (ABTS●+) . To reduce oxidative stress in the host cell, we previously engineered MTs into E. coli hosts for both cytosolic and outer-membrane-targeted (osmoregulatory membrane protein, OmpC, fused) expression. The ability of these engineered E. coli to scavenge intracellular or extracellular ROS was examined, and OmpC-fused TMT performed the best, growing in medium containing 1.5% n-butanol. n-Butanol tolerance was increased through the scavenging of intracellular and extracellular free radicals, and the fusion protein still contributed to osmotic tolerance via either glycine betaine or glucose uptake .
Intracellular ROS increase levels in E. coli after exposure to n-butanol . We therefore considered whether introducing OmpC-fused TMT into butanologenic E. coli would improve n-butanol tolerance and productivity. In the present study, we engineered the C. acetobutylicum n-butanol synthetic pathway into E. coli, along with OmpC-TMT co-expression to increase n-butanol tolerance. A synthetic biological technical platform (Ordered Gene Assembly in Bacillus subtilis – OGAB method) was employed for the genetic manipulation of the n-butanol production and tolerance gene clusters . Our results demonstrate that co-expressing OmpC-fused TMT in a butanologenic E. coli strain enhances n-butanol production.
Results and discussion
Construction of butanologenic E. coli strains
To establish butanologenic E. coli strains, we transferred only an essential set of genes (thil, crt, bcd, etfAB, hbd, adhe) for n-butanol production [5, 6]. In addition, our previous results indicated that the expression of the ompC-tmt gene enhances the tolerance for n-butanol through osmotic tolerance and extracellular radical scavenging capacity . These genes were cloned and constructed in the plasmid pGETS118 using the OGAB method . The plasmids pBUT, containing Pr-thil-crt-bcd-etfAB-hbd-adhe, pBUT-rO-tMT, containing Pr-thil-crt-bcd-etfAB-hbd-adhe-rbs-ompC-tmt, and pBUT-T7-rO-tMT, with Pr-thil-crt-bcd-etfAB-hbd-adhe-T7 promoter-rbs-ompC-tmt, were constructed. The constructed plasmids pBUT, pBUT-rO-tMT and pBUT-T7-rO-tMT (Fig. 1) were used to transform E. coli JM109 and E. coli JM109 (DE3) to obtain the butanologenic E. coli strains BUT1, BUT2, BUT3, BUT1-DE, BUT2-DE and BUT3-DE. Finally, pGETS118 was used to transform E. coli JM109 and E. coli JM109 (DE3) to obtain the control strains PGETS118 and PGETS118-DE (Table 1). Cell growth and n-butanol production were then investigated.
Growth profile and n-butanol production
Recently, many groups have reported successful n-butanol production in E. coli [5, 6, 8]; a similar approach was previously applied to n-butanol production in E. coli carrying the CoA-dependent synthetic pathway, yielding final n-butanol titers of 552 mg/L to 1.2 g/L [5, 6]. In addition, Shen et al.  further improved the strain by using NADH- and acetyl-CoA-driving forces; the best strain in their study showed high-titer production of n-butanol through in situ product removal. The production of n-butanol has been confirmed to affect host growth.
Compared with the control strains PGETS118 and PGETS118-DE, both recombinant strains BUT1 and BUT1-DE displayed increased n-butanol titer (both reached approximately 73 mg/L) after 20 h. Furthermore, 17-21% reductions in the growth of both strains in PYG medium were observed (Figs. 2 and 3). In addition, after 60 h of incubation, maximum n-butanol concentrations of 94.6 mg/L and 95.1 mg/L were recorded in BUT1 and BUT1-DE cultures, respectively. This result shows that cells at an OD600 value of 1 can convert 1 g of glucose to 1.75 mg/L of n-butanol in PYG medium, similar to the rate reported by Inui et al.  (1.5 mg/L n-butanol/1 g glucose). In the case of PGETS118 and PGETS118-DE, which were monitored as controls, no n-butanol was detected.
Improved tolerance and n-butanol production
A clostridial strain overexpressing GSH was reported to exhibit improved tolerance and increased n-butanol productivity . Boyarskiy et al.  reported increased n-butanol tolerance and production (approximately 35% increases) through transcriptional feedback regulation of efflux protein expression.
In the present study, we attempted to improve the n-butanol tolerance of butanologenic E. coli to further enhance n-butanol production by expressing OmpC-TMT. The concentrations of n-butanol found to be inhibitory to E. coli in n-butanol production strains (Fig. 2) were consistent with the trends noted in the literature for alcohol toxicity in E. coli . Hence, we examined the normal expression of OmpC-TMT (BUT2, BUT3) in butanologenic E. coli strains, and the results showed that cell growth was still affected by n-butanol. This result is similar to that obtained for the BUT1 strain (Fig. 2a).
However, the BUT2 and BUT3 strains showed similar n-butanol production (118 mg/L and 121 mg/L) after 20 h of incubation (Fig. 3a). After incubation for 60 h, the butanologenic E. coli strains BUT2 and BUT3 showed 1.6-fold and 2.7-fold higher n-butanol production, respectively, than BUT1 (94.6 mg/L). Notably, BUT3-DE achieved a 1.7-fold increase in n-butanol production compared to that achieved by BUT2-DE after 20 h of incubation (Fig. 3b). Furthermore, cell growth in the OmpC-TMT-overexpressing butanologenic E. coli (BUT3-DE) was nearly the same as in the control strain (PGETS118-DE) after 60 h of incubation (Fig. 2b). More importantly, n-butanol production increased nearly 3.4-fold compared with that of BUT1-DE (320 mg/L) and showed a significantly improved conversion rate from glucose to n-butanol (4.97 mg/L n-butanol/1 g glucose) (Fig. 3b). These results suggested that OmpC-TMT overexpression rescued the growth of the butanologenic E. coli strains, improved host robustness and increased n-butanol productivity. Our results are consistent with previous reports [28, 30, 31].
In addition, many genes were related to organic solvent tolerance in the E. coli strain, for example, mutant cyclic AMP receptor protein (CRP) , overexpression of groESL [33–35], mutant Δlon (cell envelope-related gene)  and control membrane-related functions (overexpression of ATF, fabD, feoA and srpABC) . Moreover, Rutherford et al.  reported that n-butanol stress-response genes are also involved in many stress responses, such as oxidative stress (sodA, sodC and yqhD), heat shock and cell envelope stress (rpoE, clpB, htpG, cpxR and cpxP). Recently, the tolerance mechanisms of several critical genes have been elucidated. Among them, the downregulated genes yghW and yibT were shown to improve n-butanol tolerance due to their regulatory roles in membrane fatty acid composition, while the upregulated genes gcl and glcF improved cell growth and metabolism by replenishing TCA cycle metabolic intermediates . However, these results reflect the addition of n-butanol to the medium rather than its production by E. coli (Table 2).
Aside from n-butanol, carboxylic acids are attractive biorenewable chemicals, such as fatty acid and aromatic carboxylic acids [39–41]. However, product toxicity is frequently encountered in metabolic engineering, and the toxicity of these carboxylic acids to the microbial biocatalyst appears to limit biocatalyst performance. Strong evidence suggests that membrane damage is the main mechanism of toxicity [42, 43], and previous metabolic engineering efforts successfully increased membrane integrity by modulating membrane composition to alter free fatty acid tolerance in E. coli [44–46].
Free radical scavenging ability and membrane integrity test
Since the cellular membrane is a vital factor that allows cells to acclimate to external stresses and is also one of the components that is strongly affected by organic solvents , many studies have proposed that the plasma membrane is the most affected target of organic solvents and plays a significant role in adapting to stress. Additionally, the length of the carbon backbone of the organic solvent can alter the toxicity mechanism; increasing the hydrophobicity of the solvent increases the level of toxicity . Long- and short-chain alcohols cause stress during biofuel production by altering membrane fluidity (also known as Overton’s Rule). Ethanol and n-butanol, respectively, decrease and increase membrane fluidity [47, 49, 50].
Toxicity that causes membrane damage is likely a key limitation for butanologenic microbes to produce butanol, suggesting that our engineering strategy of expressing OmpC-TMT is effective in protecting membrane integrity and scavenging free radicals to improve n-butanol tolerance and production in butanologenic E. coli strains.
To test whether OmpC-TMT increased membrane integrity during fermentation, bacteria were stained with DAPI and SYTOX Green nucleic acid stains after 48 h of incubation (Fig. 4). Viable bacteria or cells with membrane intact cells were stained and appeared blue only (a, d), whereas nonviable bacteria or cells with damaged membranes would appeared blue and green (b, e and c, f). The ratio of cell membrane damage was determined as the number of SYTOX stain-labeled bacteria cells/the number of DAPI stain-labeled bacteria cells and was 1.57% ± 0.33%, 48.59% ± 5.98%, and 19.69% ± 4.16%, for PGETS118-DE, BUT1-DE and BUT3-DE, respectively. The engineered E. coli strain BUT3-DE, expressing the outer-membrane fusion protein OmpC-TMT, exhibited better membrane integrity than that of E. coli strain BUT1-DE. These results clarified the relationship between membrane integrity and butanol resistance.
MTs are antioxidants that scavenge radicals, and alcohols are known to cause oxidative stress in E. coli . In the present study, we examined the ability of OmpC-TMT to scavenge free radicals when the host cells were producing n-butanol. We then analyzed the ROS content of the cells using 5(6)-carboxy-2′, 7′- dichlorodihydro-fluorescein diacetate (carboxy-H2DCFDA). During the first 60 h of incubation, the levels of free radicals in all strains increased with increasing concentration of n-butanol, as shown by the increasing levels of fluorescence (Fig. 5a). After fermentation for 20 h, BUT3-DE displayed 2-fold higher free radical content and 2.8-fold greater n-butanol production than BUT1-DE (Figs. 3b and 5a). In contrast, both PGETS118-DE and BUT1-DE had lower levels of radicals than the OmpC-TMT-expressing strains, because no or less n-butanol was produced. Moreover, we examined the capacity for OmpC-TMT to scavenge free radicals when the host cells produced n-butanol. As shown in Fig. 3b, BUT3-DE (320.24 mg/L) produced more n-butanol than BUT2-DE (151.03 mg/L). However, both BUT2-DE and BUT3-DE produced similarly large amounts of ROS after 60 h of incubation (Fig. 5a). The effect of n-butanol on each strain was defined as [specific fluorescence of each strain (A535/A600)/(n-butanol production by each strain)]. BUT3-DE showed approximately 2-fold higher free radical scavenging capacity (25.27) than BUT1-DE (54.07) and BUT2-DE (51.25) (Fig. 5b). Interestingly, this result is consistent with the result from the membrane integrity test. Therefore, OmpC-TMT overexpression enhanced n-butanol production by lowering the levels of free radicals and increasing membrane integrity.
Expression profiles of n-butanol and tolerance genes in E. coli
As shown in Table 3, the relative expression of the four genetic regions (thil, crt-bcd-etfB-etfA-hbd, adhe, ompC-tmt) with respect to the rrsA gene in the two sample types (control and experimental) for three cases (1. BUT1-DE and BUT2-DE; 2. BUT1-DE and BUT3-DE; and 3. BUT2-DE and BUT3-DE) was determined using qPCR assays. In this study, differential expression was observed between the OmpC-tMT strains and the control. In case 1, expression of the four genetic regions was lower in the control (BUT1-DE). BUT1-DE expressed ompC-tmt, but we did not introduce this gene into the plasmid. This phenomenon was observed in another study in which an E. coli strain was observed to tolerate solvents through the upregulation of OmpC [38, 52]. Therefore, gene expression may be derived from the host. Nevertheless, BUT2-DE exhibited high free radical scavenging capacity and high n-butanol yield compared with BUT1-DE (Fig. 5b), indicating that lower oxidative stress and higher n-butanol production may result from TMT protein expression. In case 2 with BUT3-DE, the expression of n-butanol production genes was still lower than in the control (BUT1-DE), but the ompC-tmt gene showed high expression due to the T7 promoter. Hence, we examined the ability of BUT3-DE to scavenge free radicals, and the result showed BUT3-DE had the greatest ability to reduce oxidative stress and the best n-butanol production (Figs. 3b and 5a). Finally, in case 3, BUT3-DE expressed more n-butanol production genes than BUT2-DE. These results again suggest that OmpC-TMT overexpression improved the robustness of the host to achieve better solvent production by decreasing the production of ROS.
Transcriptomic analysis in engineered E. coli
Understanding the membrane stress response to solvents and alcohols may facilitate the engineering of microorganisms for improved toxin tolerance. Therefore, the stress responses of organisms such as E. coli to ethanol exposure have been widely studied [53–55], and information from these studies has been successfully adapted to engineer improved ethanologenic hosts [56–59]. Although both strains appeared to grow into stationary phase after 36 h of incubation, meaningful increases (P-value < 0.01) in their butanol titers were still observed after 48 h (Fig. 3b, BUT1-DE: 74.49 ± 1.09 mg/l (36 h), 81.55 ± 1.15 mg/l (48 h) and 95.19 ± 1.62 mg/l (60 h); and BUT3-DE: 285.94 ± 4.10 mg/l (36 h), 299.81 ± 2.92 mg/l (48 h) and 320.24 ± 2.56 mg/l (60 h)). Meanwhile, as shown by the levels of fluorescence after 36 h and 60 h, ROS levels increased in BUT1-DE but decreased in BUT3-DE (Fig. 5a). Because BUT3-DE and BUT1-DE exhibited opposite changes in oxidative stress after 60 h (Fig. 5b), the effects of OmpC-TMT expression on butanologenic E. coli strains were explored via transcriptomic analysis using next-generation sequencing technology. Among the differentially expressed genes, a total of 147 genes (Additional file 1: Table S3), 73 upregulated and 74 downregulated genes in BUT3-DE vs. Control (BUT1-DE) met the selection criteria of (1) genes that changed by fragments per kilobase of exon per million fragments mapped reads (FPKM) > 0.3  and (2) a minimum of a twofold change in normalized read counts between groups. The P-value (P < 0.05) was estimated for each gene and was corrected for multiple testing using the Benjamini-Hochberg correction. The fold-change was used to partition the genes into up- and downregulated groups. Then, all the significantly expressed genes were subjected to KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway analysis (database from KEGG E. coli str. K12 substr. W3110). We identified five significant KEGG pathways in the BUT3-DE OmpC-TMT-expressing strain, including oxidative phosphorylation (4 genes), fructose and mannose metabolism (3 genes), ribosome (4 genes), glycolysis/gluconeogenesis (2 genes) and nicotinate and nicotinamide metabolism (2 genes). All 15 significantly expressed genes were compared between BUT3-DE and BUT1-DE (Table 4). We observed 4 downregulated genes (nuoI, nuoG, nuoF and nuoC) in oxidative phosphorylation (Table 4). The nuo operon encodes the NADH ubiquinone oxidoreductase (complex I) that couples the transfer of electrons from NADH to ubiquinone with the translocation of protons across the cytoplasmic membrane and forms an integral part of oxidative phosphorylation . Quinones are electron carriers with an isoprenoid side chain that anchors them to the membrane; they function as the primary electron carriers for respiration and are thought to regulate the ArcA–ArcB two-component system in response to the cellular redox state [62–64].
However, the past results showed that isobutanol disrupts the cell membrane, leading to quinone malfunction, which results in the release of quinone inhibition on ArcB and the subsequent autophosphorylation of ArcB and activation of ArcA, which adapts cellular metabolism for growth with decreased respiratory efficiency, such as downregulation of NADH dehydrogenase I (nuo operon) or upregulation of cytochrome d oxidase (cydAB operon) . This hypothesis coincides with our results. In addition, microbial fermentation pathways involve many redox reactions, which usually require NADH and NAD+ as cofactors. However, in the clostridial synthetic pathway, the synthesis of n-butanol from glucose can cause an NADH/NAD+ redox imbalance because more NADH is required in the synthetic pathway than is generated in the glycolytic pathway . Therefore, compared with the control strain, in the modified strain, a number of upregulated genes, such as fucA, srlE and srlA, may help supplement metabolic intermediates for fructose and mannose metabolism to generate NADH and assist in regulating NADH/NAD+ redox [66, 67]. These results suggested that maintaining the redox balance of NADH and NAD+ is a key to ensure the continued operation of cellular metabolism under n-butanol fermentation. The same strategy was used in previous research [8, 65]. Moreover, isobutanol disrupts quinone/quinol function, resulting in the malfunction of NADH ubiquinone oxidoreductase (complex I), which would diminish the ability to convert NADH to NAD+ and generate the proton motive force (PMF) . This process may affect the biosynthesis of ATP in the cell. We also found that two genes (pfkB, pgm) involved in glycolysis/gluconeogenesis were downregulated. One possible explanation is that the downregulation of pfkB, which encodes 6-phosphofructokinase II, results in the decreased conversion of fructose-6-phosphate (F6P) to fructose-1, 6-biphosphate, thereby reducing the consumption of ATP. Additionally, pgm encodes phosphoglucomutase and facilitates the interconversion of glucose-1-phosphate and glucose-6-phosphate (G6P). One possible explanation for the downregulation of pgm is that glucose-1-phosphate is further converted to ADP-glucose by ATP depletion. This result indicated that BUT3-DE tried to reduce ATP consumption following NADH dehydrogenase I malfunction. Moreover, both the pfkB and pgm genes were downregulated, which may result in the accumulation of G6P and F6P. Both G6P and F6P, as starting metabolites for the pentose phosphate pathway, produce precursors and provide a major source of NADPH for biosynthesis, which is involved in most of reductive pathways [68, 69]. However, previous studies show the energy-dependent reduction of NADP+ with NADH by PntAB, but increased NADPH formation will downregulate pntA for genetic and environmental manipulation .
In general, cells exposed to n-butanol show elevated levels of oxidative stress, and many oxidative-related genes have been identified as being significantly upregulated . However, in our case, although BUT3-DE showed elevated butanol production and higher radical scavenging capacity than BUT1-DE (Figs. 3b and 5b), we did not detect significant expression of oxidative stress-related genes. In fact, our analysis found that genes related to oxidative stress (sodA) and cell envelope stress (cpxP) were expressed (Additional file 1: Table S3). Meanwhile, KEGG pathway analysis revealed no significant differences between BUT1-DE and BUT3-DE with respect to these genes, and thus was not included in Table 4. Given that BUT1-DE experienced higher oxidative stress than BUT3-DE (Fig. 5b), these results suggest that OmpC-TMT overexpression decreased ROS production, as well as the expression of oxidative stress-related genes, in BUT3-DE.
This study uses a novel approach to develop E. coli strains that express both the n-butanol synthesis pathway and membrane-targeted MTs to improve cell fermentation via ROS scavenging. OmpC-TMT was able to decrease the production of free radicals and to improve host robustness during fermentation. Furthermore, a balanced intracellular redox state in microbes is recognized as essential to ensure the efficient production of fermentation products. Unfortunately, our studies suggest that although OmpC-TMT expression increased n-butanol production compared to that of the control strain, the accumulation of fermented products (e.g., n-butanol) in the cell or external environment cause malfunctions in the NADH and ATP synthesis system. Therefore, this fermentation system requires further improvement. In principle, the expression of membrane-targeted TMT has the potential to improve the alcohol fermentation bioprocess and may also serve as a practical strategy for the construction of platform E. coli strains for biofuel production.
All chemicals and reagents were purchased from Sigma-Aldrich Co., USA, unless otherwise noted. The reagents, when available, were molecular biology grade. All solutions were prepared using these reagents and sterile distilled water.
Construction of plasmids via the OGAB method
The OGAB method was described previously . The primer sequences used in this study are listed in Additional file 2: Table S1. To construct an expression plasmid containing the clostridial n-butanol synthesis pathway, thil (1.2 kb) and the crt-bcd-etfB-etfA-hbd cluster (4.8 kb) were first amplified from the genomic DNA of C. acetobutylicum ATCC824, and the adhel (2.7 kb) gene was amplified from C. acetobutylicum ATCC824 megaplasmid pSOL1. ompC-tmt (1.3 kb) was amplified from pET-OmTmt . The PCR products were subcloned into the yT&A cloning vector (YEASTERN) and confirmed by DNA sequencing. All PCR products were digested and then assembled as one transcriptional unit into pGETS118 in a polycistronic manner using the OGAB method, resulting in pBUT, pBUT-rO-tMT and pBUT-T7-rO-tMT (Fig. 1). All genes and an assembly vector were prepared as DNA fragments with arbitrary 3’-nucleotide protrusions that allowed the genes to be linked together in a specific order and direction. The protrusions were generated by type II restriction endonucleases that recognize two separate recognition sites, such as SfiI, which recognizes 5’-GGCCANNNAGGCC-3’, where N is any nucleotide in a random sequence. SfiI and BglI were chosen in this study since SfiI and BglI sites are absent in the relevant genes and vectors. The constructed plasmids were transferred into E. coli JM109 and JM109 (DE3) to obtain the recombinant E. coli strains BUT1, BUT2, BUT3, BUT1-DE, BUT2-DE and BUT3-DE (Table 1) and to assay the production of n-butanol.
Bacterial strains, culture media and growth conditions
MT-expressing engineered strains, protein expression and their locations in recombinants E. coli hosts were confirmed in our previous study . All batch cultures were grown at 37 °C in a rotary shaker at 200 rpm. Batch cultures were grown in 150 mL PYG medium under anaerobic conditions (per liter: 5 g of peptone, 10 g of yeast extract, 60 g of glucose, 5 g of tryptone, 40 mg of K2HPO4, 19.2 mg of MgSO4.7 (H2O), 8 mg of CaCl2, 40 mg of KH2PO4, 0.4 g of NaHCO3, 80 mg of NaCl and 1.1 mg of FeSO4.7 H2O) or in M9 (AMRESCO-J863) media. Each engineered E. coli strain, including PGETS118, BUT1, BUT2, BUT3, PGETS118-DE, BUT1-DE, BUT2-DE, and BUT3-DE (Table 1), was grown in medium supplemented with 50 μg/mL chloramphenicol. When the culture density reached OD 0.6, isopropyl-β-D-thiogalactopyranoside (IPTG) was added to a final culture concentration of 0.6 mM. After 12 h of incubation, cells were harvested for ROS and n-butanol production experiments. All solvent concentrations in media are reported as mg/L.
Cell density was measured at 600 nm using a UV–vis spectrophotometer (GENESYS 10S, Thermo Scientific, USA). Fermentation samples were analyzed by gas chromatography (GC) on an Agilent 7890A equipped with a split/splitless injector, flame ionization detector, and a DP-FFAP capillary column (30 m, 0.32 mm i.d., 0.25 μm film thickness). The temperatures of the injector and detector were both maintained at 225 °C. The temperature profile of the column oven was as follows: (i) initial hold at 50 °C for 4 min, increase to 100 °C at 20 °C/min and hold for 1 min; (ii) increase to 170 °C at 30 °C/min and hold for 2.5 min; and (iii) increase to 220 °C at 20 °C/min and hold for 4 min. Nitrogen gas was used as the carrier gas. Split injection mode was used, with a split ratio of 1:10.
Reactive oxygen species detected by carboxy-H2DCFDA under n-butanol stress
The engineered E. coli strains were cultured in PYG medium, and different concentrations of n-butanol were produced after 0, 20, 36 and 60 h of incubation. Aliquots of 100 μL of the cultured strains were re-suspended in 5 mL of M9 medium, and 140 μL of each diluted sample was transferred to a 96-well plate, followed by incubation at 37 °C. The assay method was adapted from a previous study . All samples were treated with 10 μL of 25 mM carboxy-H2DCFDA (Invitrogen, Co., Carlsbad, CA) and incubated at 37 °C for 15 min. The optical density at 600 nm and the fluorescence excitation/emission at 535/600 nm of each sample were measured on a plate reader. Tert-butyl hydroperoxide (TBHP) (Invitrogen, Carlsbad, CA) is a known stressor that produces intracellular H2O2; a set of positive controls for the ROS assay were prepared with strains cultured after inoculation and were treated as above, except with an initial dilution step with 10 μL of 7.78 M TBHP.
Staining of bacterial suspensions with DAPI and SYTOX Green
Bacterial suspensions containing either 1.0 × 106 or 1.2 × 106 organisms/mL were stained with 0.5 μM SYTOX Green and 1 μg/ml DAPI. The fluorescence emissions from E. coli were compared with the background fluorescence of each stain alone at the optimal excitation wavelength for each nucleic acid stain (DAPI, 358 nm; SYTOX Green, 504 nm). The mixture of cells and DAPI was incubated at 37 °C for 10 min, and the cells were washed twice with TBS buffer (pH 7.6). Next, the DAPI-stained cells were stained with SYTOX Green nucleic acid again, incubated at 37 °C for 10 min, and washed three times with TBS buffer (pH 7.6). Fluorescence was observed on a Nikon ECLIPSE 80i (Nikon, Tokyo, Japan).
Next-generation sequencing and analysis
After 60 h of incubation, total RNA was isolated from both butanologenic strains (BUT1-DE and BUT3-DE) for next-generation sequencing (RNA-Seq), which was performed by Welgene Biotech Co., Ltd. (Taipei, Taiwan). Total RNA was extracted using the Trizol® Reagent (Invitrogen, USA) according to the manufacturer’s instructions. Purified RNA was quantified at 260 nm using an ND-1000 spectrophotometer (Nanodrop Technology, USA) and was qualitatively analyzed using a Bioanalyzer 2100 (Agilent Technology, USA) with an RNA 6000 labchip kit (Agilent Technologies, USA). All procedures were performed according to the Illumina protocol. For all samples, library construction was performed using TruSeq RNA Sample Prep Kits v2 for 75 bp (Single-End) sequencing and the Solexa platform (Illumina Inc.). Sequences were directly determined using sequencing-by-synthesis technology via the TruSeq SBS Kit. Raw sequences were obtained from the Illumina Pipeline software bcl2fastq v2.0 and were expected to generate 10 M (million reads or Gb) per sample. The generated sequences were filtered to obtain qualified reads. Trimmomatics was implemented to trim or remove reads according to the quality score . Qualified reads after filtering low-quality data were analyzed using TopHat/Cufflinks for gene expression estimation . The gene expression level was calculated as FPKM (fragments per kilobase of transcript per million mapped reads). For differential expression analysis, CummeRbund was employed to perform statistical analysis of the gene expression profiles. The reference genome and gene annotations were retrieved from the Ensembl database for E. coli str. K-12 substr. W3110. The cuffdiff tool from the cufflinks package was run to calculate expression changes and the associated P-values for each gene between strain BUT1-DE (control) and strain BUT3-DE (Table 1) after 60 h of incubation. The output files of cuffdiff were annotated by adding gene functional descriptions. KEGG term enrichment analysis and fold-change enrichment for the gene lists of significantly upregulated and downregulated genes in two butanologenic strains were performed. KEGG pathways were analyzed using NIH DAVID Bioinformatics Resources 6.7  to identify the regulated biological themes.
Quantitative Real-Time PCR Analysis
Total RNA for quantitative reverse transcriptase PCR (qPCR) was extracted as described for next-generation sequencing and analysis. Two-step quantitative real-time PCR (qPCR) was conducted using LightCycler ® FastStart DNA Master PLUS SYBR Green I (Roche Life Science) on a Light Cycler 480 Real-Time PCR 3System (Roche Co., Germany) and was performed by Welgene Biotech Co., Ltd. (Taipei, Taiwan). The four genetic regions involved in n-butanol production pathway (thil, crt-bcd-etfB-etfA-hbd, adhe) and ROS-targeted scavenging (ompC-tmt) were analyzed. R16s (rrsA gene) was used as a reference in all analyses. The primers are provided in Additional file 3: Table S2 in the supplemental material. For the relative quantification of gene expression, the comparative CT method was employed. The averaged CT was subtracted from the corresponding averaged r16S value for each sample; resulting in ∆CT. ∆∆CT was obtained by subtracting the average control ∆CT value from the average experimental ∆CT. The fold increase was determined by calculating log2 (2−∆∆CT) for the experimental vs. control samples.
Medium containing peptone, yeast extract and glucose.
Reactive oxygen species
Ezeji TC, Qureshi N, Blaschek HP. Bioproduction of butanol from biomass: from genes to bioreactors. Curr Opin Biotechnol. 2007;18(3):220–7.
Green EM. Fermentative production of butanol--the industrial perspective. Curr Opin Biotechnol. 2011;22(3):337–43.
Procentese A, Raganati F, Olivieri G, Russo ME, De La Feld M, Marzocchella A. Renewable feedstocks for biobutanol production by fermentation. N Biotechnol. 2016. doi: 10.1016/j.nbt.2016.10.010.
Harvey BG, Meylemans HA. The role of butanol in the development of sustainable fuel technologies. J Chem Technol Biotechnol. 2011;86(1):2–9.
Atsumi S, Cann AF, Connor MR, Shen CR, Smith KM, Brynildsen MP, et al. Metabolic engineering of Escherichia coli for 1-butanol production. Metab Eng. 2008;10(6):305–11.
Inui M, Suda M, Kimura S, Yasuda K, Suzuki H, Toda H, et al. Expression of Clostridium acetobutylicum butanol synthetic genes in Escherichia coli. Appl Microbiol Biotechnol. 2008;77(6):1305–16.
Dellomonaco C, Clomburg JM, Miller EN, Gonzalez R. Engineered reversal of the beta-oxidation cycle for the synthesis of fuels and chemicals. Nature. 2011;476(7360):355–9.
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.
Ye Q, Bao J, Zhong J-J. Bioreactor Engineering Research and Industrial Applications I: Cell Factories. 1st ed. Springer: Berlin; 2016. https://link.springer.com/book/10.1007%2F978-3-662-49161-4.
Bowles LK, Ellefson WL. Effects of butanol on Clostridium acetobutylicum. Appl Environ Microbiol. 1985;50(5):1165–70.
Knoshaug EP, Zhang M. Butanol tolerance in a selection of microorganisms. Appl Biochem Biotechnol. 2009;153(1-3):13–20.
Mao S, Luo Y, Zhang T, Li J, Bao G, Zhu Y, et al. Proteome reference map and comparative proteomic analysis between a wild type Clostridium acetobutylicum DSM 1731 and its mutant with enhanced butanol tolerance and butanol yield. J Proteome Res. 2010;9(6):3046–61.
Liu XB, Gu QY, Yu XB. Repetitive domestication to enhance butanol tolerance and production in Clostridium acetobutylicum through artificial simulation of bio-evolution. Bioresour Technol. 2013;130:638–43.
Tomas CA, Beamish J, Papoutsakis ET. Transcriptional analysis of butanol stress and tolerance in Clostridium acetobutylicum. J Bacteriol. 2004;186(7):2006–18.
Gonzalez R, Tao H, Purvis JE, York SW, Shanmugam KT, Ingram LO. Gene array-based identification of changes that contribute to ethanol tolerance in ethanologenic Escherichia coli: comparison of KO11 (parent) to LY01 (resistant mutant). Biotechnol Prog. 2003;19(2):612–23.
Brynildsen MP, Liao JC. An integrated network approach identifies the isobutanol response network of Escherichia coli. Mol Syst Biol. 2009;5:277.
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(6):1935–45.
Sigler K, Chaloupka J, Brozmanova J, Stadler N, Hofer M. Oxidative stress in microorganisms--I. Microbial vs. higher cells--damage and defenses in relation to cell aging and death. Folia Microbiol (Praha). 1999;44(6):587–624.
Perrone GG, Tan SX, Dawes IW. Reactive oxygen species and yeast apoptosis. Biochim Biophys Acta. 2008;1783(7):1354–68.
Jones DP. Radical-free biology of oxidative stress. Am J Physiol Cell Physiol. 2008;295(4):C849–68.
Fridovich I. Superoxide radical and superoxide dismutases. Annu Rev Biochem. 1995;64:97–112.
Imlay JA. Pathways of oxidative damage. Annu Rev Microbiol. 2003;57:395–418.
Coyle P, Philcox JC, Carey LC, Rofe AM. Metallothionein: the multipurpose protein. Cell Mol Life Sci. 2002;59(4):627–47.
Vallee BL. The function of metallothionein. Neurochem Int. 1995;27(1):23–33.
Park JD, Liu Y, Klaassen CD. Protective effect of metallothionein against the toxicity of cadmium and other metals(1). Toxicology. 2001;163(2-3):93–100.
Sato M, Kondoh M. Recent studies on metallothionein: protection against toxicity of heavy metals and oxygen free radicals. Tohoku J Exp Med. 2002;196(1):9–22.
Lin KH, Chien MF, Hsieh JL, Huang CC. Mercury resistance and accumulation in Escherichia coli with cell surface expression of fish metallothionein. Appl Microbiol Biotechnol. 2010;87(2):561-69.
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.
Tsuge K, Matsui K, Itaya M. One step assembly of multiple DNA fragments with a designed order and orientation in Bacillus subtilis plasmid. Nucleic Acids Res. 2003;31(21):e133.
Zhu L, Dong H, Zhang Y, Li Y. Engineering the robustness of Clostridium acetobutylicum by introducing glutathione biosynthetic capability. Metab Eng. 2011;13(4):426–34.
Boyarskiy S, Davis Lopez 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.
Kempf B, Bremer E. Uptake and synthesis of compatible solutes as microbial stress responses to high-osmolality environments. Arch Microbiol. 1998;170(5):319–30.
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.
Zingaro KA, Terry PE. 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.
Yuan Y, Bi C, Nicolaou SA, Zingaro KA, Ralston M, Papoutsakis ET. Overexpression of the Lactobacillus plantarum peptidoglycan biosynthesis murA2 gene increases the tolerance of Escherichia coli to alcohols and enhances ethanol production. Appl Microbiol Biotechnol. 2014;98(19):8399–411.
Zingaro KA, Nicolaou SA, Yuan Y, Papoutsakis ET. Exploring the heterologous genomic space for building, stepwise, complex, multicomponent tolerance to toxic chemicals. ACS Synth Biol. 2014;3(7):476–86.
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.
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.
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.
Mialon L, Vanderhenst R, Pemba AG, Miller SA. Polyalkylenehydroxybenzoates (PAHBs): biorenewable aromatic/aliphatic polyesters from lignin. Macromol Rapid Commun. 2011;32(17):1386–92.
Polen T, Spelberg M, Bott M. Toward biotechnological production of adipic acid and precursors from biorenewables. J Biotechnol. 2013;167(2):75–84.
Steen EJ, Kang Y, Bokinsky G, Hu Z, Schirmer A, McClure A, et al. Microbial production of fatty-acid-derived fuels and chemicals from plant biomass. Nature. 2010;463(7280):559–62.
Lennen RM, Kruziki MA, Kumar K, Zinkel RA, Burnum KE, Lipton MS, et al. Membrane stresses induced by overproduction of free fatty acids in Escherichia coli. Appl Environ Microbiol. 2011;77(22):8114–28.
Jarboe LR, Royce LA, Liu P. Understanding biocatalyst inhibition by carboxylic acids. Front Microbiol. 2013;4:272.
Lennen RM, Pfleger BF. Modulating membrane composition alters free fatty acid tolerance in Escherichia coli. PLoS One. 2013;8(1):e54031.
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.
Tan Z, Yoon JM, Nielsen DR, Shanks JV, Jarboe LR. Membrane engineering via trans unsaturated fatty acids production improves Escherichia coli robustness and production of biorenewables. Metab Eng. 2016;35:105–13.
Ingram L. Microbial tolerance to alcohols: role of the cell membrane. Trends Biotechnol. 1986;4(2):40–4.
Aono R, Nakajima H. Organic solvent tolerance in Escherichia coli. Tanpakushitsu Kakusan Koso. 1997;42(15):2532–41.
Ingram LO, Buttke TM. Effects of alcohols on micro-organisms. Adv Microb Physiol. 1984;25:253–300.
Kabelitz N, Santos PM, Heipieper HJ. Effect of aliphatic alcohols on growth and degree of saturation of membrane lipids in Acinetobacter calcoaceticus. FEMS Microbiol Lett. 2003;220(2):223–7.
Trautwein K, Kuhner S, Wohlbrand L, Halder T, Kuchta K, Steinbuchel A, et al. Solvent stress response of the denitrifying bacterium “Aromatoleum aromaticum” strain EbN1. Appl Environ Microbiol. 2008;74(8):2267–74.
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.
Ingram LO. Ethanol tolerance in bacteria. Crit Rev Biotechnol. 1990;9(4):305–19.
Goodarzi H, Bennett BD, Amini S, Reaves ML, Hottes AK, Rabinowitz JD, et al. Regulatory and metabolic rewiring during laboratory evolution of ethanol tolerance in E. coli. Mol Syst Biol. 2010;6:378.
Woodruff LB, Pandhal J, Ow SY, Karimpour-Fard A, Weiss SJ, Wright PC, et al. Genome-scale identification and characterization of ethanol tolerance genes in Escherichia coli. Metab Eng. 2013;15:124–33.
Asako H, Nakajima H, Kobayashi K, Kobayashi M, Aono R. Organic solvent tolerance and antibiotic resistance increased by overexpression of marA in Escherichia coli. Appl Environ Microbiol. 1997;63(4):1428–33.
Dien BS, Cotta MA, Jeffries TW. Bacteria engineered for fuel ethanol production: current status. Appl Microbiol Biotechnol. 2003;63(3):258–66.
Luo LH, Seo PS, Seo JW, Heo SY, Kim DH, Kim CH. Improved ethanol tolerance in Escherichia coli by changing the cellular fatty acids composition through genetic manipulation. Biotechnol Lett. 2009;31(12):1867–71.
Jarboe LR, Liu P, Royce LA. Engineering inhibitor tolerance for the production of biorenewable fuels and chemicals. Curr Opin Chem Eng. 2011;1(1):38–42.
Ramskold D, Wang ET, Burge CB, Sandberg R. An abundance of ubiquitously expressed genes revealed by tissue transcriptome sequence data. PLoS Comput Biol. 2009;5(12):e1000598.
Ingledew WJ, Poole RK. The respiratory chains of Escherichia coli. Microbiol Rev. 1984;48(3):222–71.
Georgellis D, Kwon O, Lin EC. Quinones as the redox signal for the arc two-component system of bacteria. Science. 2001;292(5525):2314–6.
Malpica R, Franco B, Rodriguez C, Kwon O, Georgellis D. Identification of a quinone-sensitive redox switch in the ArcB sensor kinase. Proc Natl Acad Sci U S A. 2004;101(36):13318–23.
Malpica R, Sandoval GR, Rodriguez C, Franco B, Georgellis D. Signaling by the arc two-component system provides a link between the redox state of the quinone pool and gene expression. Antioxid Redox Signal. 2006;8(5-6):781–95.
Saini M, Li SY, Wang ZW, Chiang CJ, Chao YP. Systematic engineering of the central metabolism in Escherichia coli for effective production of n-butanol. Biotechnol Biofuels. 2016;9:69.
Du Toit PJ, Kotze JP. The isolation and characterization of sorbitol-6-phosphate dehydrogenase from Clostridium pasteurianum. Biochim Biophys Acta. 1970;206(3):333–42.
Hacking AJ, Lin EC. Disruption of the fucose pathway as a consequence of genetic adaptation to propanediol as a carbon source in Escherichia coli. J Bacteriol. 1976;126(3):1166–72.
White D. The physiology and biochemistry of prokaryotes. 3rd ed. New York: Oxford University Press; 2007.
Sundara Sekar B, Seol E, Mohan Raj S, Park S. Co-production of hydrogen and ethanol by pfkA-deficient Escherichia coli with activated pentose-phosphate pathway: reduction of pyruvate accumulation. Biotechnol Biofuels. 2016;9:95.
Sauer U, Canonaco F, Heri S, Perrenoud A, Fischer E. The soluble and membrane-bound transhydrogenases UdhA and PntAB have divergent functions in NADPH metabolism of Escherichia coli. J Biol Chem. 2004;279(8):6613–9.
Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30(15):2114–20.
Trapnell C, Roberts A, Goff L, Pertea G, Kim D, Kelley DR, et al. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat Protoc. 2012;7(3):562–78.
NIH DAVID Bioinformatics Resources 6.7. Available online: http://david.abcc.ncifcrf.gov/(accessed on 28 Dec 2014).
Nishizaki T, Tsuge K, Itaya M, Doi N, Yanagawa H. Metabolic engineering of carotenoid biosynthesis in Escherichia coli by ordered gene assembly in Bacillus subtilis. Appl Environ Microbiol. 2007;73(4):1355–61.
We thank Jian-Rong Li and Tzu-Hang Yuan for help with the next-generation sequencing data analysis and in gathering information.
This work was supported by grants (MOST 104-2621-M-005-003-MY3) from the National Science Council, Taiwan.
Availability of data and materials
The datasets supporting the conclusions of this article are included within the manuscript and its additional files. The data set supporting the results of this article are available in the NCBI Sequence Read Archive (http://www.ncbi.nlm.nih.gov/sra/) repositories, SRR5352318, and SRR5352319.
WCC performed the experiments, analyzed the data and drafted the manuscript. KHL performed some of the experimental work. KHL contributed to data interpretation. WCC and KHL wrote and revised the manuscript. WCC and CCH designed the study. CCH coordinated the study. All authors read and approved the final manuscript.
The authors declare no competing interests.
Consent for publication
Ethics approval and consent to participate
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
List of the 147 differentially expressed genes. In all 147 differentially expressed genes, including 132 down- and 4 up-regulated, were detected in the BUT3-DE, compared to the BUT1-DE. (DOCX 37 kb)
Primes used for OGAB method. The primer sequences used in this study are listed. (DOCX 23 kb)
Primes used for qPCR. The primer sequences used in this study are listed. (DOCX 22 kb)