Skip to main content

Transcriptome analysis of Zymomonas mobilis ZM4 reveals mechanisms of tolerance and detoxification of phenolic aldehyde inhibitors from lignocellulose pretreatment



Phenolic aldehydes generated from lignocellulose pretreatment exhibited severe toxic inhibitions on microbial growth and fermentation. Numerous tolerance studies against furfural, 5-hydroxymethyl-2-furaldehyde (HMF), acetate, and ethanol were reported, but studies on inhibition of phenolic aldehyde inhibitors are rare. For ethanologenic strains, Zymomonas mobilis ZM4 is high in ethanol productivity and genetic manipulation feasibility, but sensitive to phenolic aldehyde inhibitors. Molecular mechanisms of tolerance for Z. mobilis toward phenolic aldehydes are not known.


We took the first insight into genomic response of Z. mobilis ZM4 to the phenolic aldehyde inhibitors derived from lignocellulose pretreatment. The results suggest that the toxicity to cells is caused by the functional group of phenolic aldehyde, similar to furfural and HMF, rather than aromatic groups or phenolic hydroxyl groups. Transcriptome response against 4-hydroxybenzaldehyde, syringaldehyde, and vanillin, representing phenolic groups H, S, and G, respectively, was investigated. The atlas of the important genes responsible for significantly enhanced and repressed genes at the genomic level was illustrated. 272 genes with twofold greater expressions than non-treated controls and 36 gene clusters in response to challenges of these phenolic aldehydes were identified. Several reductases encoded by ZMO1116, ZMO1696, and ZMO1885 were found to play the key roles in reducing phenolic aldehydes into the corresponding phenolic alcohols. Reduction of phenolic aldehydes by overexpression of ZMO1116, ZMO1696, and ZMO1885 in Z. mobilis ZM4 resulted in the increased inhibitor conversion and ethanol productivity, especially for 4-hydroxybenzaldehyde and vanillin. Several transporter genes such as ZMO0282, ZMO0283, ZMO0798, ZMO0799, and ZMO0800 was also displayed significantly increased expressions against the phenolic aldehydes.


The genes encoding reductases are with potentials on phenolic aldehydes-tolerant genes contributing to the reduction of phenolic aldehydes into the corresponding phenolic alcohols forms for Z. mobilis ZM4. Overexpression of the key genes improved the conversion ratio and ethanol productivity of 4-hydroxybenzaldehyde and vanillin with high toxicity. New knowledge obtained from this research aids understanding the mechanisms of bacterial tolerance and the development of the next-generation biocatalysts for advanced biofuels production.


Pretreatment is the central step for releasing fermentable sugars from lignocellulose biomass. During harsh pretreatment processes, various small molecules are generated from partial over-degradation of lignocellulose and inhibit consequent microbial fermentations. Commonly observed inhibitory compounds include furan aldehydes such as 2-furylaldehyde (furfural) and 5-hydroxymethyl-2-furaldehyde (HMF) from dehydration of pentoses and hexoses, weak organic acids such as acetic acid, formic acid, and levulinic acid from carboxylate group hydrolysis or furans oxidations, as well as phenolic compounds degraded from partial breakdown of lignin components [1, 2]. Phenolic compounds are classified into three major groups according to their methoxyl and functional groups: (1) p-hydroxyphenyl group (H) represented by 4-hydroxybenzaldehyde, 4-hydroxybenzoate, and 4-hydroxybenzyl alcohol; (2) syringyl group (S) including syringaldehyde (4-hydroxy-3,5-dimethoxybenzaldehyde), syringate, and syringic alcohol; and (3) guaiacyl group (G) including guaiacol, vanillin (4-hydroxy-3-methoxybenzaldehyde), vanillate, vanillyl alcohol, ferulic acid, and coniferyl aldehyde. Most of these phenolic compounds had been recognized in pretreated lignocellulose materials, including 4-hydroxybenzaldehyde, 4-hydroxybenzoate, syringaldehyde, syringate, vanillin, vanillate, and coniferyl aldehyde [1]. Phenolic compounds showed severe toxic inhibitions on fermenting strains [3] and cellulase activities [4]. Yet, identification and accurate measurement of all phenolic compounds in lignocellulose hydrolysate remain challenging because of poor water solubility and large number of derivatives.

Ethanologenic strain Zymomonas mobilis ZM4 has demonstrated its potential in lignocellulose biorefinery applications for high ethanol productivity, ethanol tolerance, and genetic manipulation feasibility [5, 6]. It tolerates high levels of phenolic acids [7], but is sensitive to phenolic aldehydes. Existence of minimum vanillin and syringaldehyde will lead to 1/3 decrease of ethanol yield for Z. mobilis [8] due to the disruptions of cell membrane and enzyme hydrophobic sites by phenolic compounds [1, 2]. While numerous studies on tolerance of varied strains were reported against ethanol [9], acetate [10], furfural [11], studies on inhibition of phenolic aldehyde inhibitors are rare. Molecular mechanisms of tolerance for Z. mobilis toward phenolic aldehydes are not known.

In this study, we investigated the profiles of genome expression of Z. mobilis ZM4 using DNA microarray in response to three typical phenolic aldehydes, 4-hydroxybenzaldehyde representing phenolic group H, syringaldehyde for group S, and vanillin for group G. Candidate genes against varied phenolic aldehydes were identified from Z. mobilis ZM4 and were selectively expressed in Z. mobilis ZM4 for the confirmation of gene functions. This study provides the first insight into the genome response of Z. mobilis against phenolic aldehyde inhibitors. Tolerant genes identified in this study will serve as valuable resources for robust strain development for future biorefinery applications.


Cell growth and fermentation response of Z. mobilis ZM4 to phenolic aldehydes

The maximum cell growth of Z. mobilis ZM4 was depressed approximately 5.0, 1.4, and 3.1-folds when 4-hydroxybenzaldehyde, syringaldehyde, and vanillin were separately added, comparing to the non-phenolic aldehydes-treated culture (Fig. 1a). Z. mobilis ZM4 completely consumed glucose and reached the maximum ethanol titer at 8 h, but delayed about 12, 16, and 24 h and led to decrease of ethanol productivity from 0.73 to 0.24, 0.70, and 0.61 g/L/h when 4-hydroxybenzaldehyde, syringaldehyde, and vanillin were separately added (Fig. 1b, c). Apparently, 4-hydroxybenzaldehyde is the most potent inhibitory compound to Z. mobilis ZM4 among the three evaluated phenolic aldehydes. Conversion of phenolic aldehydes into the corresponding phenolic alcohols by Z. mobilis ZM4 was identified by GC–MS (Fig. 1d, e, Additional file 1). No further degradation of phenolic alcohols to other metabolites was detected. Approximately, 80 % of 4-hydroxybenzaldehyde was reduced into 4-hydroxybenzyl alcohol and 70 % of vanillin degraded into vanillyl alcohol after incubated for 36 h. On the other hand, syringaldehyde was relatively stable with only a minimum degradation in 36 h. The results were partially in agreement with the previous studies of Z. mobilis ZM4 converting 4-hydroxybenzaldehyde, syringaldehyde, and vanillin into their alcoholic forms [12].

Fig. 1

Inhibition of performance for Z. mobilis strain ZM4 by selective phenolic aldehydes. a Cell growth as measured by optical density at 600 nm (OD600nm); b glucose consumption; c ethanol fermentation; d degradation of phenolic aldehydes; e formation of phenolic alcohols in RM medium separately amended with 5 mM of 4-hydroxybenzaldehyde, syringaldehyde, and vanillin

Transcriptome response of Z. mobilis ZM4 against phenolic aldehydes

The microarray results revealed that a total of 2020 gene fragments were amplified. To confirm microarray expression data, 20 genes involving in central carbon metabolic pathways (Additional file 2) were selected for quantitative expression using qRT-PCR under the stress of the three phenolic aldehydes. Results were consistent in detection of the trend of up- or down-regulated expressions for both assays, except for only one gene showing a marginal variation (Fig. 2). As anticipated, quantitative results of qRT-PCR were more sensitive in general with greater fold changes than that of microarray.

Fig. 2

Comparison of gene expression levels of Z. mobilis ZM4 between the DNA microarray and qRT-PCR. The gene expression ratios of both microarray data and qRT-PCR data for 20 genes were log transformed in base 2 (log2, treatment/control), and the microarray log2 ratio values were plotted against the qRT-PCR log2 values. The 20 selected genes included ZMO0152 (pyruvate kinase), ZMO0177 (16S rRNA), ZMO0179 (fructose-bisphosphate aldolase), ZMO0367 (glucose-6-phosphate 1-dehydrogenase), ZMO0368 (6-phosphogluconate dehydrogenase), ZMO0369 (glucokinase), ZMO0387 (hpcH/hpaI aldolase), ZMO0543 (aconitate hydratase), ZMO0544 (isocitrate dehydrogenase), ZMO0567 (succinyl-CoA synthetase), ZMO0569 (succinate dehydrogenase), ZMO0997 (KDPG aldolase), ZMO1237 (d-isomer specific 2-hydroxyacid dehydrogenase), ZMO1307 (fumarase), ZMO1360 (pyruvate decarboxylase), ZMO1478 (6-phosphogluconolactonase), ZMO1496 (phosphoenolpyruvate carboxylase), ZMO1596 (Fe-containing alcohol dehydrogenase), ZMO1608 (phosphopyruvate hydratase), and ZMO1963 (citrate synthase). The primers for qRT-PCR of the genes are listed in Additional file 1

Our hierarchical cluster analysis showed two clear distinguished expression patterns for 1908 genes in the genome of Z. mobilis ZM4 in response to challenges of 4-hydroxybenzaldehyde, syringaldehyde, and vanillin representing three groups of phenolic aldehyde inhibitors (Fig. 3). The genome response to syringaldehyde challenge was similar and very closely related to the non-treated control. In contrast, the expression patterns against 4-hydroxybenzaldehyde and vanillin were closely related and distinct from the control as well as the responses treated with syringaldehyde. Results from the replicated microarray experiments showed high levels of consistence with minimum variations within each treatment. Based on their expression responses, two major groups of genes were clustered separately. The control and the syringaldehyde-treated cultures shared a group of approximate 16 genes showing up-regulated response. In contrast, most genes in this group displayed a repressed or down-regulated expression in response to the 4-hydroxybenzaldehyde or vanillin treatment. On the other hand, many genes showing up-regulated expressions against 4-hydroxybenzaldehyde and vanillin displayed normal or down-regulated expressions for the syringaldehyde treated and the control. In response to challenge of a single phenolic aldehyde, 144 genes showed up-regulated expression and 299 genes repressed with a greater than twofold changes against 4-hydroxybenzaldehyde, a representative of p-hydroxyphenyl group (H), compared with the control (Additional file 3). For vanillin, a representative of guaiacyl group, 111 genes were up-regulated expressed and 198 genes down-regulated. For syringaldehyde representing syringyl group (S), only 17 genes were up-regulated and 63 down-regulated. We found 272 genes including 36 gene clusters showing up-regulated expressions against these phenolic aldehydes distributed throughout the entire genome of Z. mobilis ZM4 (Fig. 4). Another 560 genes including 52 gene clusters were repressed significantly by three phenolic aldehydes. We identified 72 candidate tolerance genes showing significantly up-regulated expressions against challenges of at least two phenolic aldehydes examined in this study (Table 1).

Fig. 3

Hierarchical cluster analysis of genome expression of Z. mobilis strain ZM4 in response to selective phenolic aldehydes. Gene expression levels were clustered based on log2 transformation. Color code red, green, or black each represents up-, down-, or normally regulated expression, respectively. Letter H, G, and S, each stands for, phenolic group H represented by 4-hydroxybenzaldehyde, phenolic group G represented by vanillin, and phenolic group S represented by syringaldehyde, respectively. CK stands for a non-phenolic aldehyde-treated control. Replications of each treatment for the microarray experiment were marked in Roman numeral

Fig. 4

Genome map of Z. mobilis strain ZM4 showing expression profiling in response to phenolic aldehyde challenges. Genes and gene clusters shown significantly up-regulated expressions were boxed in red, and those significantly down-regulated expressions boxed in blue. The genome map representing 1998 genes was adapted from KEGG database. A colored code indicating annotated gene functional categories is also attached

Table 1 Tolerance genes of Z. mobilis ZM4 against at least two phenolic aldehydes identified by gene expression analysis

The results show that the transcriptome response to syringaldehyde stress for Z. mobilis was significantly different from that to 4-hydroxybenzaldehyde and vanillin stress. The possible reasons may come from the following two aspects: (1) low molecular weight. The molecular weights of 4-hydroxybenzaldehyde, syringaldehyde, and vanillin were 122.1, 182.2, and 152.2, respectively. Similar to S. cerevisiae to low molecular weights phenolic compounds [1], the low molecular weight syringaldehyde gave the least toxicity for Z. mobilis ZM4. (2) Substituent position. The methoxyl substituent at ortho position increases the toxicity of vanillins [13], while the methoxyl and hydroxyl substituents at meta and para positions or vice versa do not change the toxicity [1]. Vanillin has one more methoxyl group and syringaldehyde has two more methoxyl groups on 4-hydroxybenzaldehyde. The configuration of the two methoxyl groups may generate space shelling for the toxic methoxyl substituent at ortho position, and then reduced the toxicity of syringaldehyde.

Important reductase genes

Among many highly induced expressions, two reductase genes, ZMO1116 encoding oxidoreductase and ZMO1885 encoding NADH: flavin oxidoreductase/NADH oxidase, were outstanding in response to the challenge of all three phenolic aldehydes. ZMO1116 was up-regulated expression with 2.53-, 2.47-, and 3.60-fold changes against 4-hydroxybenzaldehyde, vanillin, and syringaldehyde, respectively (Table 1), and ZMO1885 was up-regulated by 11.37-, 2.67-, and 3.33-fold changes. Other genes, ZMO0157 encoded d-isomer specific 2-hydroxyacid dehydrogenase NAD binding, ZMO0788 encoded a gluconate 2-dehydrogenase, ZMO0833 encoded UDP-N-acetylenolpyruvoylglucosamine reductase, ZMO1222 encoded 3-oxoacyl-(acyl-carrier-protein) reductase, ZMO1254 encoded redoxin domain protein, ZMO1303 encoded a pyrroline-5-carboxylate reductase, ZMO1399 encoded fatty acid hydroxylase, ZMO1576 encoded short-chain dehydrogenase/reductase SDR, ZMO1696 encoding zinc-binding alcohol dehydrogenase, ZMO1984 encoded aldo/keto reductase, were also up-regulated over twofold change under at least two phenolic aldehyde inhibitors. ZMO0020, ZMO0021, ZMO0074, ZMO0268, ZMO0270, ZMO0384, ZMO0388, ZMO0391, ZMO0440, ZMO0486, ZMO1334, ZMO1346, ZMO1386, ZMO1602, ZMO1821, ZMO1880, ZZM4_0133, and pzmob1_p38 encoded hypothetical proteins, also showed high levels of enhanced expressions in response to at least two phenolic aldehyde inhibitors. Figure 5 illustrates the degradation pathways of the three phenolic aldehydes in Z. mobilis ZM4 with the information of significantly differentially expressed genes.

Fig. 5

Conversion pathways of phenolic aldehydes. a Reduction of 4-hydroxybenzaldehyde, b 4-hydroxy-3-methoxybenzaldehyde (vanillin), and c syringaldehyde into the corresponding phenolic alcohols by candidate genes using NADH as a cofactor. Genes colored red indicate the most significant shared by all three phenolic aldehydes

Significant transporter genes for detoxification

Efflux system of living cells is an efficient mechanism for detoxification of external toxic compounds and internal damaging intermediates. We found at least four transporter genes were significantly up-regulated in response to the challenge of phenolic aldehydes. Two genes ZMO0799 and ZMO0800 belonging to adenosine triphosphate (ATP)-binding cassette (ABC) superfamily showed consistently greater expressions to all the three phenolic aldehydes (Table 1). Other two transporter genes, ZMO0282 and ZMO0798 belonging to resistance-nodulation-cell division (RND) superfamily, showed up to sevenfold increased expression to the phenolic challenge. One gene ZMO1288 encoding a major facilitator superfamily (MFS) transporter protein was down-regulated more than twofold change. To limit the intracellular concentration of toxic inhibitors, Gram-negative bacteria can diminish its entry by developing a barrier of low permeability [14]. This non-specific phenomenon is also existing in Z. mobilis ZM4, such as the down-regulation of ZMO1288 encoding a major facilitator superfamily (MFS) transporter protein.

Overexpression and tolerance evaluation of the important genes in Z. mobilis ZM4

Genes with strong response to phenolic aldehydes on the transcriptional level were screened and overexpressed in Z. mobilis ZM4 to confirm their functions. Four genes, ZMO1116, ZMO1288, ZMO1696, and ZMO1885, were selected, including (1) significantly up-regulated genes: ZMO1116 encoding an oxidoreductase with up-regulation by 2.53-, 2.47-, and 3.60-fold change, and ZMO1885 encoding a NADH:flavin oxidoreductase/NADH oxidase with 11.37-, 2.67-, and 3.33-fold change under the stress of 4-hydroxybenzaldehyde, syringaldehyde, and vanillin, respectively (data shown in Table 1); (2) important genes on the degradation pathway with relatively high up-regulation: ZMO1696 encoding zinc-binding alcohol dehydrogenase with up-regulation by 2.46-, 1.52-, and 2.94-fold change under the same stress (data shown in Table 1); (3) significantly down-regulated gene: ZMO1288 encoding a major facilitator superfamily (MFS) transporter protein with down-regulation by 4.07-, 4.18-, and 3.44-fold change under the above same stress (see Additional file 4). The selection of ZMO1288 for overexpression is to validate whether the down-regulation of ZMO1288 contributed to the resistance of phenolic aldehyde inhibitors.

Four pHW20a plasmids harboring ZMO1116, ZMO1288, ZMO1696, or ZMO1885 (Fig. 6) were introduced into Z. mobilis ZM4 to generate the recombinants according to the protocol in Dong et al. [15]. The shuttle expression vector pHW20a, as well as pHW20a-gfp with gfp and a short linker were used as the two controls. The successful expression of ZMO1116, ZMO1696, and ZMO1885 were confirmed by vanillin degradation by the recombinants at the rate of 0.112, 0.125, and 0.023 μmol per minute, respectively. The expression of the transporter gene ZMO1288 was confirmed by the fluorescence intensity of 10,340 at 507 nm for the recombinant with the fused gfp, comparing to 9497 of the control containing only gfp reporter gene. The recombinants were cultured in RM medium separately amended with 5 mM of 4-hydroxybenzaldehyde, syringaldehyde, or vanillin. The cell growth, glucose consumption rate, ethanol productivity, and phenolic aldehyde inhibitors conversion rate were measured as shown in Table 2. Under the stress of 4-hydroxybenzaldehyde, the overexpression of ZMO1116 increased the cell growth, glucose consumption, ethanol productivity, and 4-hydroxybenzaldehyde conversion by 3.92, 23.39, 13.64, and 10.92 %, respectively. Overexpression of ZMO1288 led to the increases of the above parameters by 10.42, 60.84, 50.00, and 86.84 %, respectively. Overexpression of ZMO1696 led to the increases of the above parameters by 3.92, 20.00, 13.64, and 133.48 %, respectively. Overexpression of ZMO1885 led to the increases of the above parameters by 5.88, 31.61, 18.18, and 52.11 %, respectively. On the other hand, the expression of all the genes on cell growth increased a little. Under the stress of syringaldehyde, the overexpression of ZMO1116, ZMO1288, ZMO1696, and ZMO1885 increased syringaldehyde conversion by 18.40, 51.87, 20.43, and 30.47 %, respectively, but the cell growth, glucose consumption, and ethanol productivity only showed limited change. Under the stress of vanillin, the overexpression of ZMO1116 increased the cell growth, glucose consumption, ethanol productivity, and vanillin conversion by 40.68, 16.32, 11.76, and 19.49 %, respectively. Overexpression of ZMO1288 led to the increases of 45.90, 19.10, 14.29, and 210.23 %, respectively. Overexpression of ZMO1696 led to the increases of 8.47, 20.52, 7.84, and 193.93 %, respectively. Overexpression of ZMO1885 led to the increases of 16.95, 8.16, 13.73, and 97.43 %, respectively.

Fig. 6

Construction of the plasmids for expression of functional genes with strong response to phenolic aldehydes. a pHW20a. b pHW20a-gfp. c pHW20a-ZMO1116. d pHW20a-ZMO1288-gfp. e pHW20a-ZMO1696. f pHW20a-ZMO1885

Table 2 Evaluation of selective Z. mobilis recombinants with genes showing differential expressions in genome expression analysis

The low toxicity of syringaldehyde to Z. mobilis led to the limited inhibitory influence on cell growth and ethanol productivity, so the reduced level of syringaldehyde by the overexpression of the key genes also did not increase the cell growth and ethanol productivity. But the obvious improvements of cell growth and ethanol fermentation were obtained for the more toxic 4-hydroxybenzaldehyde and vanillin for inhibitor conversion.

Overexpression of the reductase genes ZMO1116, ZMO1696, and ZMO1885 in Z. mobilis ZM4 significantly improved the tolerance to phenolic aldehydes, indicating the importance of reductases for phenolic aldehydes conversion. Müller et al. [16] found the reduction activity of NAD(P)H-dependent reductase encoded by ZMO1885 from Z. mobilis on the carbon double bond, and this study also verified the reduction activity of phenolic aldehydes into phenolic alcohols. Overexpression of ZMO1288 encoding a major facilitator superfamily (MFS) protein contributed to the highest inhibitor conversion, suggesting that the increase of toxin transport into intracellular space by enhancing ZMO1288 activity did not harm the intracellular activity. Instead, the enhanced transport of inhibitors to intracellular space facilitated the fast degradation and reduction of the toxic phenolic aldehydes.


In this study, our transcriptome analysis using microarray revealed the first insight into genomic response of Z. mobilis ZM4 to the three groups of phenolic aldehyde inhibitors derived from lignocellulose pretreatment. We identified 272 up-regulated genes, including 36 gene clusters, against phenolic aldehydes as the potential candidate tolerance genes. We also identified significantly repressed genes and gene clusters to the phenolic aldehydes and provided their physical locations in the genome map. Furthermore, we illustrated mechanisms of bacterial tolerance and detoxification of these phenolic aldehydes through a group of genes encoding varied reductases and transporters. New knowledge obtained from this research aids understanding of the bacterial tolerance and the development of strong inhibitor tolerant fermenting strains.

Phenolic aldehydes are among the most toxic inhibitory compounds for microbial growth and fermentation. The toxicity is often thought caused by aromatic ring and the mechanism of inhibition is not clear. It is generally accepted that these phenolic compounds are classified into three groups of p-hydroxyphenyl group (H), guaiacyl group (G), and syringyl group (S) based on their methoxyl groups [1]. We used three phenolic aldehydes of 4-hydroxybenzaldehyde, vanillin, and syringaldehyde to represent each of these three different phenolic groups. We observed in the Z. mobilis ZM4 cultures, phenolic aldehydes concentrations were reduced and the corresponding alcohols were increased. As a result, when the phenolic aldehydes were reduced to certain lower levels, Z. mobilis ZM4 cells were able to accelerate the rate of glucose consumption and cell growth and eventually complete the fermentation. As demonstrated by our GC–MS analysis, 4-hydroxybenzaldehyde was converted into 4-hydroxybenzyl alcohol, and vanillin into vanillyl alcohol, and syringaldehyde into syringic alcohol. In these corresponding alcohols, a typical aromatic ring still existed but it apparently did not cause inhibition to the cell growth and metabolism of Z. mobilis ZM4. Functional group of aldehyde is apparently the toxic component of these phenolic aldehydes, regardless of their classification of phenolic group H, G, or S, based on the mass balance of phenolic aldehydes and its alcohol or acid derivatives, as well as the almost constant total phenolic concentrations. In response to the challenge of these phenolic aldehydes, the gene expression of ZMO0157, ZMO0788, ZMO0833, ZMO1116, ZMO1222, ZMO1254, ZMO1303, ZMO1399, ZMO1576, ZMO1696, ZMO1885, and ZMO1984 was significantly up-regulated under the stress of at least two phenolic aldehydes. These genes were found primarily encoding oxidoreductase, alcohol dehydrogenase, oxidoreductase/NADH oxidase, aldo/keto reductase, and a member of short-chain dehydrogenase/reductase (SDR), respectively. Our results suggest these genes are commonly involved in the conversion of phenolic aldehydes into the corresponding phenolic alcohols. However, since the complexity of the compound structure differs, the fine toned mechanism of the conversion pathway can be different. For example, syringaldehyde is more stable due to its complex structure and may require additional conversion pathways. It indicated that reduction and transport contributed to the conversion of the phenolic aldehydes into the corresponding phenolic alcohols.

Similarly, furfural and HMF were conventionally called furan inhibitors that hindered studies on microbial tolerance to these common biomass fermentation inhibitors. He et al. [9, 11] has found 127 genes and 433 genes differentially expressed in Z. mobilis in response to ethanol and furfural by DNA microarray, respectively. Metabolic conversion pathways of furfural and HMF into furanmethanol and furan-2,5-dimethanol by tolerant yeast Saccharomyces cerevisiae NRRL Y-50049 were presented and the toxic functional group of aldehyde, not furan group, was clarified [1719]. It appeared there might be an aldehyde reductase gene family exist [20, 21]. In this study, we observed multiple enzymes were involved in the conversion of phenolic aldehydes into corresponding alcohols that is consistent with previous reported. The overexpression of oxidoreductase (ZMO1116), zinc-binding alcohol dehydrogenase (ZMO1696), and NADH:flavin oxidoreductase/NADH oxidase (ZMO1885) with significantly enhanced conversion of phenolic aldehydes into phenolic alcohols gave further supports the multiple gene involvement in microbe tolerance to aldehyde inhibitors.

In our study, we also observed expression levels of at least eight transporter genes were significantly increased in response to varied phenolic aldehydes, including ZMO0143, ZMO0799, and ZMO0800 from ATP-binding cassette (ABC) superfamily; ZMO0282, ZMO0798, and ZMO1529 from RND (resistance-nodulation-cell division) superfamily; and ZMO0283 and ZMO0965 from efflux system proteins. In yeast, numerous transporter genes were identified to be involved in response to aldehyde including ABC transporter genes and other drug and toxin transporter genes [22]. Efflux pumps from RND family are considered one of the most efficient mechanisms for tolerance to toxin in Gram-negative bacteria [23]. Results of this study in Z. mobilis ZM4 in response to phenolic aldehyde inhibitors concur with the conclusion obtained from S. cerevisiae.

Our comprehensive analysis of the genome expression in Z. mobilis ZM4 provided the first important gene atlas for significantly enhanced and repressed genes at the genomic level. This illustrated distribution of potential candidate tolerance genes and gene clusters at the genomic landscape can serve as a valuable resource for genetic engineering efforts in genomic manipulation of Z. mobilis ZM4 for the continued tolerant strain development. In general, enhanced expressed genes under the inhibitor challenge pressure are considered as tolerance genes. However, this has to be in the context of consistent dynamics over a time course. In this study, we took merely a snapshot in the mid of the logarithm growth phase, thus, a more conserved approach should be taken in selection of candidate tolerance genes. On the other hand, it needs to point out that repressed genes are also important, particularly for those recovered from the early repression. Such genes often contribute to adaptation, redirecting gene interactions, and rewired networks for survival and tolerance development as demonstrated in S. cerevisiae [24, 25].


The ethanologenic bacterium Z. mobilis ZM4 degrades phenolic aldehydes into their corresponding phenolic alcohols forms, and the toxicity of phenolic aldehydes to cell growth and ethanol fermentation is caused by the functional group of phenolic aldehyde rather than the phenolic aromatic ring. Among 272 up-regulated genes, including 36 gene clusters, reductases are the potential candidate tolerance genes for the reduction of phenolic aldehydes. Tolerant genes identified in this study will serve as valuable resources for robust strain development for future biorefinery applications.


Bacterial strain, plasmid, and reagents

Zymomonas mobilis ZM4 (ATCC 31821) was purchased from American Type Culture Collection (Manassas, VA, USA). Shuttle vector pHW20a was used for construction of Z. mobilis recombinants [15]. Yeast extract was purchased from Oxoid (Hampshire, UK). Tetracycline (Tc) and nalidixic acid were from Sigma-Aldrich, St. Louis, MO, USA. 4-Hydroxybenzaldehyde, syringaldehyde, and vanillin were from Sangon Biotech Co., Shanghai, China. All other analytical grade chemicals were from Sinopharm Chemical Reagents, Shanghai, China.

Cell growth conditions

Cells of Z. mobilis ZM4 were cultured and maintained in RM medium containing 20 g/L glucose, 2 g/L KH2PO4, and 10 g/L yeast extract. The initial culture was incubated at 30 °C without agitation. RM medium was amended separately by adding 5 mM each of 4-hydroxybenzaldehyde, syringaldehyde, or vanillin for testing the inhibitor tolerance. A non-phenolic aldehyde-treated culture was served as control. Samples were taken at a 4 h interval till 36 h. For DNA microarray, 10 mL culture was transferred into 100 mL fresh RM medium when the culture was 2.0 at the optical density at 600 nm. Cells were harvested at 4 h for RNA extraction, DNA microarray, and quantitative real-time PCR (qRT-PCR).

DNA microarray

Total RNA of Z. mobilis ZM4 cells was extracted using a Trizol reagent kit from Invitrogen (Carlsbad, CA, USA). Purity and concentrations of the RNA samples were measured by a ratio of OD260/280 readings using a NanoDrop ND-1000 spectrophotometer. DNA microarray was performed by CapitalBio Co., Beijing, China. cDNA fragments were hybridized with the probes and labeled with Cy3-dCTP fluorescent dye [26]. Data were analyzed using GeneSpring V12 Agilent Technologies (Santa Clara, CA, USA) with proper normalization and quality control procedures. Raw data were transformed into algorism phase based on Log2 and statistical analysis was carried out using software CLUSTER 3.0 [27]. Significance of differential expression levels compared with the control was examined using False Discovery Rate (FDR) test (p < 0.05) [28] and analysis of variance (ANOVA). A threshold of fluorescence intensity with more or less than twofold changes was regarded as significantly differentially expressed genes. Finally, a hierarchical clustering analysis was performed using Java Treeview (Stanford University, Stanford, CA, USA).

qRT-PCR assays

To confirm expression data from microarray, the quantitative mRNA expressions were performed using qRT-PCR on a CFX96TM Real-Time System with C10000TM Thermal Cycler from Bio-Rad Laboratories (Hercules, CA, USA). Primers for qRT-PCR of the selected genes are shown in Additional file 2. The first strand of cDNA was synthesized using a cDNA synthesis kit (Torobo Co., Osaka, Japan). PCR amplification was prepared using an SYBR Green Real-time PCR Master Mix. PCR amplification profile used was follows: 94 °C for 5 min, then 35 cycles at 94 °C for 2 min and 62 °C for 30 s, and 72 °C for 30 s. 16S rRNA gene ZMOr003 was used as an internal control for data acquisition and normalization.

Construction of Z. mobilis recombinants

Genomic DNA of Z. mobilis ZM4 was extracted using Omega-Biotek Bacterial DNA kit (Norcross, GA, USA). Primers for amplification of ZMO1116, ZMO1288, ZMO1696, and ZMO1885 genes were listed in Additional file 3. Construction of the plasmids harboring ZMO1116, ZMO1288, ZMO1696, and ZMO1885 genes is illustrated in Fig. 6. The transporter gene ZMO1288 was fused with a reporter gene gfp which was amplified from the plasmid pMG36E [29], with a short linker encoding seven glycines at the C-terminal of gfp between gfp and ZMO1288 for confirmation of its expression. The plasmids were introduced into E. coli S17-1 λπ, and then into Z. mobilis ZM4 through biparental transconjugation method [15].

Evaluation of Z. mobilis ZM4 recombinants

The oxidoreductase activity of the recombinants was evaluated in 200 μL of reaction mixture containing 0.1 M potassium phosphate buffer (pH 7.5), 0.329 mM vanillin, 2.25 mM NADH, and 20 μL cell-free extract [30] at 30 °C for 30 min, and then stopped by addition of 10 μL of 3 M trichloroacetic acid solution and measured vanillin reduction by HPLC. One oxidoreductase unit (U) was defined as the amount of the crude enzyme required to reduce 1 μmol of vanillin per min at pH 7.5 and 25 °C. The expression of transporter gene ZMO1288 fused with gfp gene was detected by monitoring fluorescence using Synergy H1 Hybrid Microplate Reader at 507 nm (Biotek Instruments, Winooski, VT, USA). Fermentation performance of the newly generated Z. mobilis recombinants was evaluated as described above using RM medium separately amended with 5 mM each of 4-hydroxybenzaldehyde, syringaldehyde, or vanillin. Samples were periodically taken from the culture broth and supernatant obtained by centrifuging at 12,000g for 5 min and filtering through 0.22-μm filters.

HPLC and GC–MS analysis

Glucose and ethanol were analyzed using HPLC equipped with a refractive index detector RID-10A (Shimadzu, Kyoto, Japan) with an Aminex HPX-87H column (Bio-rad, Hercules, CA, USA) at 65 °C using a mobile phase of 5 mM H2SO4 at a rate of 0.6 mL/min. 4-Hydroxybenzaldehyde, syringaldehyde, and vanillin were analyzed using reverse-phase HPLC (SPD-20A, Shimadzu, Kyoto, Japan) equipped with a YMC-Pack ODS-A column (Tokyo, Japan) at 35 °C. 4-Hydroxybenzaldehyde and syringaldehyde were detected at 270 nm using a mobile phase of 30 % acetonitrile solution at a rate of 1.0 mL/min, and vanillin was detected at 320 nm using 50 % acetonitrile solution at a rate of 0.8 mL/min.

The biodegradation intermediates of phenolic aldehydes by Z. mobilis ZM4 was identified by GC–MS. Samples were taken at 4 h intervals after inoculation and concentrated by rotary evaporator with vacuum system, then dissolved in ethyl acetate and acetonitrile solution (2:1, v/v) and silylated with NO-bis-trimethylsilyl trifluoro-acetamide according to [31, 32]. The treated samples were analyzed using Agilent 6890 GC–MS fitted with an HP-5 MS column (30 m × 0.25 mm × 0.25 μm) (Agilent Technologies, Santa Clara, CA, USA) from 80 °C (held for 4 min) to 280 °C at 8 °C/min. 1 μL sample was injected and detected under splitless condition.



adenosine triphosphate


adenosine triphosphate (ATP)-binding cassette


complementary DNA


2′-deoxycytidine 5′-triphosphate


False Discovery Rate


guaiacyl group


gas chromatography–mass spectrometer


green fluorescent protein


p-hydroxyphenyl group




high-performance liquid chromatography


major facilitator superfamily


nicotinamide adenine dinucleotide


optical density


quantitative real-time polymerase chain reaction


rich medium


resistance-nodulation-cell division


syringyl group


Zymomonas mobilis open reading frame


short-chain dehydrogenase/reductase


  1. 1.

    Klinke HB, Thomsen AB, Ahring BK. Inhibition of ethanol-producing yeast and bacteria by degradation products produced during pre-treatment of biomass. Appl Microbiol Biotechnol. 2004;66:10–26.

    CAS  Article  Google Scholar 

  2. 2.

    Palmqvist E, Hahn-Hagerdal B. Fermentation of lignocellulosic hydrolysates. II: inhibitors and mechanisms of inhibition. Bioresour Technol. 2000;74:25–33.

    CAS  Article  Google Scholar 

  3. 3.

    Larsson S, Quintana-Sainz A, Reimann A, Nilvebrant NO, Jonsson LJ. Influence of lignocellulose-derived aromatic compounds on oxygen-limited growth and ethanolic fermentation by Saccharomyces cerevisiae. Appl Biochem Biotechnol. 2000;84–86:617–32.

    Article  Google Scholar 

  4. 4.

    Jing X, Zhang X, Bao J. Inhibition performance of lignocellulose degradation products on industrial cellulase enzymes during cellulose hydrolysis. Appl Biochem Biotechnol. 2009;159:696–707.

    CAS  Article  Google Scholar 

  5. 5.

    Rogers PL, Jeon YJ, Lee KJ, Lawford HG. Zymomonas mobilis for fuel ethanol and higher value products. Adv Biochem Eng Biotechnol. 2007;108:263–88.

    CAS  Google Scholar 

  6. 6.

    Zhang M, Eddy C, Deanda K, Finkelstein M, Picataggio S. Metabolic engineering of a pentose metabolism pathway in ethanologenic Zymomonas mobilis. Science. 1995;267:240–3.

    CAS  Article  Google Scholar 

  7. 7.

    Gu HQ, Zhang J, Bao J. High tolerance and physiological mechanism of Zymomonas mobilis to phenolic inhibitors in ethanol fermentation of corncob residue. Biotechnol Bioeng. 2015;112:1770–82.

    CAS  Article  Google Scholar 

  8. 8.

    Ranatunga TD, Jervis J, Helm RF, McMillan JD, Hatzis C. Identification of inhibitory components toxic toward Zymomonas mobilis CP4 (pZB5) xylose fermentation. Appl Biochem Biotechnol. 1997;67:185–98.

    CAS  Article  Google Scholar 

  9. 9.

    He MX, Wu B, Shui ZX, Hu QC, Wang WG, Tan FR, Tang XY, Zhu QL, Pan K, Li Q, Su XH. Transcriptome profiling of Zymomonas mobilis under ethanol stress. Biotechnol Biofuels. 2012;5:75–84.

    CAS  Article  Google Scholar 

  10. 10.

    Yang S, Pan C, Tschaplinski TJ, Hurst GB, Engle NL, Zhou W, Dam P, Xu Y, Rodriguez M Jr, Dice L, Johnson CM, Davison BH, Brown SD. Systems biology analysis of Zymomonas mobilis ZM4 ethanol stress responses. PLoS One. 2013;7:e68886–900.

    Article  Google Scholar 

  11. 11.

    He MX, Wu B, Shui ZX, Hu QC, Wang WG, Tan FR, Tang XY, Zhu QL, Pan K, Li Q, Su XH. Transcriptome profiling of Zymomonas mobilis under furfural stress. Appl Microbiol Biotechnol. 2012;95:189–99.

    CAS  Article  Google Scholar 

  12. 12.

    Franden MA, Pilath HM, Mohagheghi A, Pienkos PT, Zhang M. Inhibition of growth of Zymomonas mobilis by model compounds found in lignocellulosic hydrolysates. Biotechnol Biofuels. 2013;6:99–113.

    CAS  Article  Google Scholar 

  13. 13.

    Palmqvist E, Almeida JS, Hahn-Hagerdal B. Influence of furfural on anaerobic glycolytic kinetics of Saccharomyces cerevisiae in batch culture. Biotechnol Bioeng. 1999;62:447–54.

    CAS  Article  Google Scholar 

  14. 14.

    Nikaido H. Molecular basis of bacterial outer membrane permeability revisited. Microbiol Mol Biol Rev. 2003;67:593–656.

    CAS  Article  Google Scholar 

  15. 15.

    Dong HW, Bao J, Ryu DD, Zhong JJ. Design and construction of improved new vectors for Zymomonas mobilis recombinants. Biotechnol Bioeng. 2011;108:1616–27.

    CAS  Article  Google Scholar 

  16. 16.

    Müller A, Hauer B, Rosche B. Asymmetric alkene reduction by yeast old yellow enzymes and by a novel Zymomonas mobilis reductase. Biotechnol Bioeng. 2007;98:22–9.

    Article  Google Scholar 

  17. 17.

    Liu Z. Molecular mechanisms of yeast tolerance and in situ detoxification of lignocellulose hydrolysates. Appl Microbiol Biotechnol. 2011;90:809–25.

    CAS  Article  Google Scholar 

  18. 18.

    Liu Z. Genomics of yeast tolerance and in situ detoxification. In: Liu ZL, editor. Microbial stress tolerance for biofuels, vol. 22., Microbiology monographsBerlin: Springer; 2011. p. 1–28.

    Google Scholar 

  19. 19.

    Liu Z, Menggen M, Song MJ. Evolutionarily engineered ethanologenic yeast detoxifies lignocellulosic biomass conversion inhibitors by reprogrammed pathways. Mol Genet Genom. 2009;282:233–44.

    Article  Google Scholar 

  20. 20.

    Moon J, Liu ZL. Engineered NADH-dependent GRE2 from Saccharomyces cerevisiae by directed enzyme evolution enhances HMF reduction using additional cofactor NADPH. Enzym Microb Technol. 2012;50:115–20.

    CAS  Article  Google Scholar 

  21. 21.

    Moon J, Liu ZL. Direct enzyme assay evidence confirms aldehyde reductase function of Ydr541cp and Ygl039wp from Saccharomyces cerevisiae. Yeast. 2015;32:399–407.

    CAS  Article  Google Scholar 

  22. 22.

    Ma M, Liu Z. Mechanisms of ethanol tolerance in Saccharomyces cerevisiae. Appl Microbiol Biotechnol. 2010;87:829–45.

    CAS  Article  Google Scholar 

  23. 23.

    Segura A, Molina L, Fillet S, Krell T, Bernal P, Muñoz-Rojas J, Ramos JL. Solvent tolerance in Gram-negative bacteria. Curr Opin Biotechnol. 2012;23:415–21.

    CAS  Article  Google Scholar 

  24. 24.

    Liu ZL. Molecular mechanisms of yeast tolerance and in situ detoxification of lignocellulose hydrolysate. Appl Microbiol Biotechnol. 2011;90:809–25.

    CAS  Article  Google Scholar 

  25. 25.

    Zhang Y, Liu ZL, Song M. ChiNet uncovers rewired transcription subnetworks in tolerant yeast for advanced biofuels conversion. Nucl Acids Res. 2015;43:4393–407. doi:10.1093/nar/gkv358.

    Article  Google Scholar 

  26. 26.

    Patterson TA, Lobenhofer EK, Fulmer-Smentek SB, Collins PJ, Chu TM, Bao W, Fang H, Kawasaki ES, Hager J, Tikhonova IR, Walker SJ, Zhang L, Hurban P, de Longueville F, Fuscoe JC, Tong W, Shi L, Wolfinger RD. Performance comparison of one-color and two-color platforms within the MicroArray Quality Control (MAQC) project. Nat Biotechnol. 2006;24:1140–50.

    CAS  Article  Google Scholar 

  27. 27.

    Eisen MB, Spellman PT, Brown PO, Botstein D. Cluster analysis and display of genome-wide expression patterns. Proc Natl Acad Sci USA. 1998;95:14863–8.

    CAS  Article  Google Scholar 

  28. 28.

    Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc B. 1995;57:289–300.

    Google Scholar 

  29. 29.

    Van de Guchte M, van der Vossen JM, Kok J, Venema G. Construction of a lactococcal expression vector: expression of hen egg white lysozyme in Lactococcus lactis subsp. lactis. Appl Environ Microbiol. 1989;55:224–8.

    Google Scholar 

  30. 30.

    Abd E-Mawla AM, Beerhues L. Benzoic acid biosynthesis in cell cultures of Hypericum androsaemum. Planta. 2002;214:727–33.

    Article  Google Scholar 

  31. 31.

    Klinke HB, Ahring BK, Schmidt AS, Thomsen AB. Characterization of degradation products from alkaline wet oxidation of wheat straw. Bioresour Technol. 2002;82:15–26.

    CAS  Article  Google Scholar 

  32. 32.

    Varga E, Klinke HB, Réczey K, Thomsen AB. High solid simultaneous saccharification and fermentation of wet oxidized corn stover to ethanol. Biotechnol Bioeng. 2004;88:567–74.

    CAS  Article  Google Scholar 

Download references

Authors’ contributions

XY and JB designed the experiment. XY, ZLL, and JB wrote the manuscript. XY conducted the experiment. HQG analyzed part of the inhibitors. QQG analyzed the cloning experiment. ZLL developed the gene cluster analysis method. ZLL and JB conceived the study. All authors read and approved the final manuscript.


This research was supported by the National Basic Research Program of China (2011CB707406), the Natural Science Foundation of China (31300070), and the National High-Tech Program of China (2012AA022301/2014AA021901).

Compliance with ethical guidelines

Competing interests The authors declare that they have no competing interests.

Author information



Corresponding authors

Correspondence to Z. Lewis Liu or Jie Bao.

Additional files


Additional file 1. Degradation products of phenolic aldehydes for Z. mobilis ZM4 by GC–MS.


Additional file 2. Primers for qRT-PCR of the genes.


Additional file 3. Primers for the construction of the recombinant Z. mobilis ZM4.


Additional file 4. Differentially expressed genes (>2-fold change) under the stress of the phenolic aldehydes in Z. mobilis ZM4.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Yi, X., Gu, H., Gao, Q. et al. Transcriptome analysis of Zymomonas mobilis ZM4 reveals mechanisms of tolerance and detoxification of phenolic aldehyde inhibitors from lignocellulose pretreatment. Biotechnol Biofuels 8, 153 (2015).

Download citation


  • Zymomonas mobilis ZM4
  • Lignocellulose pretreatment
  • Phenolic aldehyde inhibitors
  • Stress tolerance
  • DNA microarray
  • Recombinants