Skip to main content

A comparative multidimensional LC-MS proteomic analysis reveals mechanisms for furan aldehyde detoxification in Thermoanaerobacter pseudethanolicus 39E



Chemical and physical pretreatment of lignocellulosic biomass improves substrate reactivity for increased microbial biofuel production, but also restricts growth via the release of furan aldehydes, such as furfural and 5-hydroxymethylfurfural (5-HMF). The physiological effects of these inhibitors on thermophilic, fermentative bacteria are important to understand; especially as cellulolytic strains are being developed for consolidated bioprocessing (CBP) of lignocellulosic feedstocks. Identifying mechanisms for detoxification of aldehydes in naturally resistant strains, such as Thermoanaerobacter spp., may also enable improvements in candidate CBP microorganisms.


Thermoanaerobacter pseudethanolicus 39E, an anaerobic, saccharolytic thermophile, was found to grow readily in the presence of 30 mM furfural and 20 mM 5-HMF and reduce these aldehydes to their respective alcohols in situ. The proteomes of T. pseudethanolicus 39E grown in the presence or absence of 15 mM furfural were compared to identify upregulated enzymes potentially responsible for the observed reduction. A total of 225 proteins were differentially regulated in response to the 15 mM furfural treatment with 152 upregulated versus 73 downregulated. Only 87 proteins exhibited a twofold or greater change in abundance in either direction. Of these, 54 were upregulated in the presence of furfural and 33 were downregulated. Two oxidoreductases were upregulated at least twofold by furfural and were targeted for further investigation. Teth39_1597 encodes a predicted butanol dehydrogenase (BdhA) and Teth39_1598, a predicted aldo/keto reductase (AKR). Both genes were cloned from T. pseudethanolicus 39E, with the respective enzymes overexpressed in E. coli and specific activities determined against a variety of aldehydes. Overexpressed BdhA showed significant activity with all aldehydes tested, including furfural and 5-HMF, using NADPH as the cofactor. Cell extracts with AKR also showed activity with NADPH, but only with four-carbon butyraldehyde and isobutyraldehyde.


T. pseudethanolicus 39E displays intrinsic tolerance to the common pretreatment inhibitors furfural and 5-HMF. Multidimensional proteomic analysis was used as an effective tool to identify putative mechanisms for detoxification of furfural and 5-HMF. T. pseudethanolicus was found to upregulate an NADPH-dependent alcohol dehydrogenase 6.8-fold in response to furfural. In vitro enzyme assays confirmed the reduction of furfural and 5-HMF to their respective alcohols.


Thermophilic bacteria, such as Clostridium thermocellum and Caldicellulosiruptor species have gained interest for their possible use as biocatalysts for converting lignocellulosic biomass into renewable fuels and chemicals [1]-[4]. The potential advantages of thermal bioprocessing include improved kinetics, reduced viscosities of concentrated slurries, lower oxygen solubility, and reduced process cooling requirements [2]. In addition, several bacterial phyla include thermophiles that are able to utilize plant cell walls directly through the action of complex (hemi)cellulase systems expressed either as free enzymes, cellulosomes, or multifunctional enzymes [5]. Relying on these native enzymes in a bioprocessing scheme could substantially reduce or even eliminate the need for exogenous enzymes for cellulose solubilization with a resulting improvement in process economics [3],[6],[7].

To render plant material more reactive to microbial or enzymatic digestion, physical and chemical pretreatment methods are generally applied, and improvement in pretreatment technologies remains a highly active field of research [8],[9]. Pretreatment with dilute acid at high temperatures has the benefit of solubilizing the hemicellulose fraction of biomass, which produces fermentable C5 oligomers and monomers [9]. However, one disadvantage of dilute acid pretreatment is that the process is non-specific and, depending on its severity, generates a number of toxic by-products [10]. Inhibitory compounds generated by dilute acid pretreatment typically fall into four categories: organic acids (acetic acid, ferulic acid), ketones (acetovanillone), phenolics (coniferyl alcohol, catechol), and aldehydes (furfural, hydroxymethylfurfural, vanillin). Mixtures of inhibitors, especially those including the furan aldehyde furfural, often have a synergistic effect on inhibiting cell growth and fermentation. For example, furfural increases the toxicity of acetate in yeast [11] and phenols in Escherichia coli[12],[13]. Furfural is estimated to be responsible for 33% of the toxic effect of sugar cane hydrolysate on E. coli LYO1 [14].

In order to compete with more robust ethanologens such as Saccharomyces cerevisiae, several limitations inherent to fermentative thermophilic bacteria must be overcome. These limitations include relatively low ethanol titer and yield from mixed-acid fermentation pathways, although breakthroughs in metabolic engineering have improved the yield of ethanol from carbohydrates in some thermophiles to near theoretical limits [15]. Another major hurdle for thermophilic bioprocessing is growth inhibition by a wide range of compounds encountered in biomass fermentations. Insights into overcoming end-product inhibition in C. thermocellum have recently emerged [16],[17]; however, growth inhibition from other biomass-derived compounds remains underexplored in thermophilic microbes relative to S. cerevisiae[18]-[20] or engineered strains of E. coli[12],[21]-[23]. Interestingly, members of the genus Thermoanaerobacter have been shown to tolerate pretreated biomass hydrolysates [24],[25], and engineered strains give improved ethanol yields from both C5 and C6 sugars [26]. These properties have encouraged the development of several Thermoanaerobacter species for bioethanol production from hydrolysates (primarily xylose) and on cellulose when paired with a cellulolytic partner [27]. While surveying thermophilic bacteria for intrinsic tolerance to furfural (unpublished), we observed robust growth and rapid reduction of the compound by Thermoanaerobacter pseudethanolicus 39E (formally known as T. ethanolicus[28]). This study aims to identify and characterize traits that enable this organism to grow in the presence of and to simultaneously detoxify furan aldehydes through reduction to less toxic alcohols.

Results and Discussion

T. pseudethanolicus 39E furan aldehyde tolerance

We initially investigated the growth tolerance of T. pseudethanolicus 39E to the furan aldehydes furfural and 5-hydroxymethylfurfural (5-HMF). The addition of 10 mM and 15 mM furfural increased specific growth rates to 0.52 ± 0.03 h-1 and 0.49 ± 0.01 h-1, respectively, versus the control at 0.38 ± 0.01 h-1 (Figure 1A). 5-HMF also stimulated growth at 10 mM compared to no addition (0.51 ± 0.03 versus 0.41 ± 0.02 h-1), while growth rates were similar to the control at 15 mM 5-HMF (0.45 ± 0.02 versus 0.41 ± 0.02 h-1; Figure 1B). Both 10 mM furfural and 5-HMF slightly increased cell yield at 12 h by approximately 11% and 12%, respectively (Figure 1). Higher growth rates and increased cell yield from the addition of subinhibitory concentrations of furfural and 5-HMF suggest that 39E metabolism is constrained by electron flow, which is relieved by the furan aldehydes serving as an alternative dissimilatory electron acceptor. The concentration resulting in 50% inhibition of growth (IC50) with furfural was 30 mM after 12 h and 30 to 40 mM after 24 h. The IC50 for 5-HMF was between 20 and 30 mM after both 12 and 24 h. As shown in Table 1, the determined values are comparable to or slightly higher than those of other thermophilic bacteria, while they are higher than reported values for E. coli, S. cerevisiae, and Zymomonas mobilis. Though direct comparisons are difficult due to differences in the growth conditions used in the various studies, these results suggest that T. pseudethanolicus 39E has a comparable if not higher tolerance to the furan aldehydes furfural and 5-HMF than other studied organisms.

Figure 1
figure 1

T. pseudethanolicus furan aldehyde tolerance. T. pseudethanolicus was grown at 65°C with 40 mM glucose and increasing concentrations of (A) furfural or (B) 5-HMF. Growth was measured by optical density at 600 nm. Error bars are the standard deviation of three replicate cultures.

Table 1 IC 50 (g/L)* for furfural and 5-hydroxymethylfurfural in thermophilic bacteria and ethanologenic microorganisms

Furan aldehyde reduction and glucose fermentation

In order to establish the mechanism of increased furan aldehyde tolerance, T. pseudethanolicus 39E was grown in the presence and absence of 15 mM furfural or 5-HMF and the furan aldehyde and respective furan alcohol concentrations were measured. As shown in Figure 2, concomitant with growth, furfural and 5-HMF concentrations decreased while 2,5-furandimethanol concentration increased, indicating that T. pseudethanolicus 39E reduced 5-HMF to 2,5-furandimethanol. Furfural was most likely also reduced to furfuryl alcohol; however, quantitation of this compound is complicated by its polymerization at the growth temperature of T. pseudethanolicus 39E.

Figure 2
figure 2

T. pseudethanolicus furan aldehyde reduction. T. pseudethanolicus was grown at 65°C with 40 mM glucose and 15 mM (A) furfural or (B) 5-HMF. Growth was measured by optical density at 600 nm. Furan aldehyde concentration was measured spectrophotometrically, while furan alcohol concentration was measured by gas chromatography-mass spectrometry. Error bars are the standard deviation of three replicate cultures.

The addition of furan aldehydes affected end-product formation by T. pseudethanolicus 39E, as determined by HPLC analysis (Figure 3). At 10 mM concentrations, more acetate and lactate are produced, while ethanol production remains constant. At furan aldehyde concentrations above 10 mM where growth is observed (15, 20, 30 mM furfural and 15, 20 mM 5-HMF), ethanol decrease and acetate increase are directly proportional, suggesting that acetyl coenzyme A is converted to acetate through phosphotransacetylase (PTA) and acetate kinase (AK) rather than serving as an electron acceptor for ethanol production via NAD(P)H-dependent bifunctional alcohol dehydrogenase (ADH) activity. T. pseudethanolicus 39E possesses seven ADHs, but ethanol is primarily produced from NADPH-dependent AdhB [34]-[36]. The oxidative branch of the pentose phosphate pathway is also present in 39E, which could supply NADPH [34]. Reduction of furfural/5-HMF to their corresponding alcohols during growth competes with ethanol production for electrons delivered by NADPH. This is also suggested by stoichiometric shifts in end products, where added aldehydes resulted in about a 0.5 times decrease in corresponding molar ethanol concentrations (that is, the 30 mM furfural addition resulted in a decrease of 15 mM ethanol versus the control). This shift in ethanol/acetate concentrations is consistent with an electron balance of one NAD(P)H per furan aldehyde and two NAD(P)H per ethanol. The reason for increased lactate with added furan aldehyde, especially 5-HMF, is less clear. Further redox imbalances from the presence of furfural or 5-HMF may direct more NADH generated from glycolysis to be oxidized via lactate dehydrogenase (LdhA). Furfural addition to a growing culture of the related thermophilic bacterium C. thermocellum also resulted in increased lactate production and cessation of ethanol production [37], although the reason for these changes is unknown. The mesophilic ethanologen S. cerevisiae has also been shown to remove furan aldehydes by reduction to their respective alcohols at the expense of ethanol production [38]. This has been shown to involve an upregulation of central carbon metabolism, especially the NADPH-generating pentose phosphate pathway [39], and downregulation of enzymes involved in ethanol formation, thereby increasing the availability of reducing equivalents for aldehyde detoxification [40].

Figure 3
figure 3

24-h substrate and end product analysis. T. pseudethanolicus was grown at 65°C with 40 mM glucose and increasing concentrations of (A) furfural or (B) 5-HMF for 24 h. Culture supernatants were analyzed for glucose and end product formation by HPLC. Glucose consumption was determined by comparison with an uninoculated control. Error bars are the standard deviation of three replicate cultures.

Cell-wide proteomic response to furfural

In order to identify potential enzymes involved in reducing furfural and 5-HMF and also to examine cell-wide physiological effects of furan aldehyde exposure, proteomic analysis was performed comparing T. pseudethanolicus 39E grown with and without furfural. A concentration of 15 mM furfural was used to challenge the cells since this was the maximum amount that displayed a stimulatory effect during the growth experiments (Figure 1A). Triplicate cultures were grown in parallel either in the presence or absence (control) of furfural and harvested at 8 h after inoculation, which corresponded to the time interval necessary for complete reduction of furfural to furfuryl alcohol (Figure 2A). Peptide samples were prepared and analyzed as described in the Methods section. The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium [41] via the PRIDE partner repository with the dataset identifier PXD001446. The complete proteomics dataset and statistical analysis are also provided in Additional file 1: Table S1, which includes an in-table, color-coded heat map corresponding to Figure 4. In total, 1,294 proteins were identified across both conditions with roughly 300,000 spectra (SpC) assigned to constituent peptides. Of the 1,294 proteins, only 918 passed the 99% SpC cutoff and were moved to the ANOVA analysis. Using a P-value cutoff of 0.05, 225 proteins were found to be differentially expressed with 152 upregulated in furfural treated cells versus 73 downregulated. Culling this list even further, only 87 proteins exhibited a twofold or greater change in abundance in either direction. Of these, 54 were upregulated in the presence of furfural (Table 2) and 33 were downregulated (Table 3). Significantly regulated proteins were grouped into 11 clusters based on abundance pattern across all replicates (Figure 4). Many cellular functions were affected by furfural, with the most highly downregulated proteins involved in cell wall biosynthesis or sporulation. Hydrogenase-related proteins were also downregulated, along with several redox proteins predicted to use NAD(P)H. Upregulated proteins fell into 12 general cellular functions, with those regulated fivefold or higher falling into three categories: polar amino acid biosynthesis (arginine, cysteine), nucleotide metabolism, and redox proteins.

Figure 4
figure 4

Heat map of proteomic analysis with and without 15 mM furfural. Proteins exhibiting a statistically significant (ANOVA, P ≤0.05) difference in abundance are included. Each protein (row) was independently normalized to recast spectral count values as standard deviations from the row mean. Protein abundance differences were clustered according to trends measured across all biological replicates. Red = increased; green = decreased abundance.

Table 2 Proteins upregulated twofold or more by 15 mM furfural ( P ≤0.05)
Table 3 Proteins downregulated twofold or more by 15 mM furfural ( P ≤0.05)

Energy production and carbohydrate metabolism

Comparatively, the proteomic response of T. pseudethanolicus 39E to furfural showed similarities at the functional level to the responses of C. thermocellum, as well as S. cerevisiae and E. coli. While central carbon metabolism did not appear to be significantly impacted by furfural in T. pseudethanolicus 39E, upregulated carbohydrate-related proteins included beta-galactosidase, mannose-6-phosphate isomerase, and kojibiose phosphorylase, while dihydroxyacetone kinase and two fructose-specific, phosphoenolpyruvate-dependent sugar phosphotransferase system transport proteins were downregulated. Genes involved in energy production and conversion, as well as carbohydrate transport and metabolism, were also regulated in C. thermocellum ATCC 27405 [37]. Acetate kinase and phosphoacetyltransferase were both downregulated, though this is likely a general stress response rather than furfural-specific. As in T. pseudethanolicus 39E, beta-galactosidase was upregulated in C. thermocellum ATCC 27405 upon furfural exposure, as were several glycosyl transferase family proteins, though the reason for this regulation is unclear. On the other hand, central carbon metabolism is significantly upregulated in both the S. cerevisiae and E. coli response to furfural. In anaerobic S. cerevisiae fermentations, an 8 g/L furfural treatment repressed the synthesis of enzymes involved in glucose catabolism and the tricarboxylic acid (TCA) cycle. Conversely, addition of 17 g/L furfural to an aerobic S. cerevisiae culture increased expression of proteins involved in glycolysis and the TCA cycle, while repressing expression of proteins involved in glycerol and ethanol production [42]. Analysis of a single-gene disruption library of S. cerevisiae BY4741 against growth with furfural identified several genes in the pentose phosphate pathway as important in furfural tolerance [39], especially ZWF1, whose overexpression allowed for growth with 50 mM furfural. In an ethanologenic strain of E. coli (LY180), a moderate furfural challenge (0.5 g/L) perturbed the expression of about 400 genes at least twofold, 15 min after exposure, with central carbon and energy metabolism being among the pathways regulated [22]. It is interesting to note that central carbon metabolism is significantly regulated in the mesophilic S. cerevisiae and E. coli and is much less affected in the thermophilic C. thermocellum and T. pseudethanolicus, though the reason for this difference remains unclear.

Stress response

With 15 mM furfural, T. pseudethanolicus 39E did not display a typical stress phenotype indicated by a reduced growth rate. Nevertheless, some functions associated with stress were differentially regulated. The expression of eight predicted transporters was affected by furfural, including upregulation of three metal transporters and one efflux transporter. A variety of nucleotide-related genes were also upregulated, including two de novo purine biosynthesis genes (Teth39_1713 and Teth39_1803) and two genes involved in nucleoside degradation (Teth39_1828 and Teth39_1832). Additionally, Teth39_1216, predicted to be involved in isoprenoid biosynthesis, and Teth39_0175, predicted to function in cell wall turnover, are downregulated in response to furfural. Similarly, a number of stress responses are upregulated by furfural in C. thermocellum ATCC 27405, including many genes that are homologous to class I and class IV heat shock response genes in Bacillus subtilis, though these genes were also upregulated by heat [37] and ethanol treatment [43]. A number of uncharacterized transporters were also regulated, as well as genes involved in transcription, RNA processing and modification, chromatin structure and dynamics, and DNA replication, recombination, and repair. In S. cerevisiae, stress responses upregulated by furfural include osmotic and salt stress, DNA damage, and pH stress [40]. S. cerevisiae also responds to furan aldehydes by regulating cell adaptation and survival processes, especially with respect to drug resistance, transport, and cell membrane composition [38]. In E. coli LY180, transport functions, as well as cell structure, DNA, and lipid synthesis functions are also regulated by furfural [22]. A general stress response to furfural thus appears to include upregulation of transport functions and nucleotide metabolism.

Amino acid metabolism

In T. pseudethanolicus 39E grown with 15 mM furfural, both cysteine and arginine biosynthetic genes were upregulated (cysteine synthase A, Teth39_0279; cysteine desulfurase, Teth39_0280; arginine biosynthesis, Teth39_0223-0229). Amino acid metabolism is also affected by furfural stress in C. thermocellum, S. cerevisiae, and E. coli. In C. thermocellum ATCC 27405, arginine biosynthetic genes are upregulated upon furfural addition [37]. While sulfur amino acid biosynthesis is not directly regulated, genes involved in sulfate transport and sulfur assimilation are upregulated by furfural. In S. cerevisiae, proteins involved in sulfur amino acid biosynthesis are downregulated upon exposure to 8 g/L furfural under anaerobic conditions [40]. In E. coli LY180, 0.5 g/L furfural also repressed genes involved in arginine biosynthesis, but induced expression of sulfur-containing amino acid biosynthetic genes [22]. This is due to a decrease in NADPH availability, which is instead used by the aldehyde reductases YqhD and DkgA for furfural reduction [22]. Upregulation of cysteine biosynthetic pathways in T. pseudethanolicus 39E suggests that, as in E. coli, increased expression of NAD(P)H-dependent aldehyde reductases (described below) may decrease NADPH availability for sulfur amino acid biosynthesis.

Redox metabolism

Six alcohol dehydrogenases (ADHs) identified in T. pseudethanolicus 39E were differentially regulated. The three functionally characterized ADHs [36],[44], AdhA (Teth39_0220), AdhB (Teth39_0218), and AdhE (Teth39_0206), were all downregulated (P <0.05); however, none more than twofold (AdhA 1.27-fold, AdhB 1.81-fold, and AdhE 1.59-fold). Of the other three identified alcohol dehydrogenases in T. pseudethanolicus 39E (Teth39_0878, Teth39_1597, Teth39_1979), only Teth39_1597 was significantly upregulated (6.8-fold; P <0.001). Another oxidoreductase upregulated by T. pseudethanolicus 39E in response to furfural is Teth39_1598 (6.3-fold; P <0.001). These genes potentially encode enzymes involved in reducing furfural and 5-HMF and will be discussed in more detail below. In C. thermocellum ATCC 27405, the bifunctional alcohol/aldehyde dehydrogenase Cthe_0423 and the redox regulator Rex (Cthe_0422) were downregulated upon furfural addition [37]. No other alcohol dehydrogenase or aldehyde reductase was differentially regulated; however, a putative carbon monoxide dehydrogenase (Cthe_0281) was upregulated by furfural and may play a role in redox balance in C. thermocellum ATCC 27405 [37]. In S. cerevisiae, NAD(P)H-dependent aldehyde reductases or alcohol dehydrogenases (ADH) have also been shown to affect furan aldehyde tolerance [45] and be regulated by furan aldehydes [40]. Transcriptomic expression analysis of known reductase and dehydrogenase genes showed that ADH2 was highly expressed in hydrolysate-tolerant S. cerevisiae strain TMB3000 compared to the wild-type CBS8066 and was also induced by 5-HMF [45]. In a proteomic analysis of the response of S. cerevisiae to 17 g/L furfural, six ADHs showed differential regulation, with Adh1p, Adh5p, and Adh6p upregulated, Adh2p and Sfa1p downregulated, and Adh4p unregulated by furfural [42]. Adh6p and Adh7p have furfural and 5-HMF reductase activity, with the former using both NADH and NADPH and the latter only NADH [46]. In E. coli LY180, energy functional groups are also highly regulated [22]. As in S. cerevisiae, NADPH-dependent aldehyde reductases with furfural reducing capacity are upregulated, namely yqhD and dkgA.

Enzyme cloning and activity measurements

Teth39_1597 [GenBank GeneID:5874751] shares 36% identity/54% similarity with yqhD from E. coli LY180 [22]. Teth39_1597 belongs to the Fe-dependent alcohol dehydrogenase superfamily (pfam00465) with predicted butanol dehydrogenase activity (BDH, cd08187). The gene product appears to be a close homolog of BdhA in Thermoanaerobacter mathranii (88% identity/94% similarity), which has been experimentally verified to have BDH activity [47]. Thus, Teth39_1597 is considered to be a butanol dehydrogenase and will be referred to as Teth39 bdhA. Teth39_1598 [GenBank GeneID:5874752] has 27% identity/41% similarity to dkgA from E. coli LY180, another enzyme shown to have NADPH-dependent furfural reductase activity [22]. Teth39_1598 is a predicted aldo/keto oxidoreductase and will be referred to as Teth39 akr.

Since both Teth39 bdhA and Teth39 akr were significantly upregulated in response to furfural and are homologs to similarly upregulated E. coli genes yqhD and dkgA, further biochemical characterization was performed to determine their cofactor and substrate specificities. The coding regions for Teth39 bdhA and Teth39 akr were PCR amplified from T. pseudethanolicus 39E genomic DNA and cloned into pET-30a behind a T7-lac promoter and N-terminal 6xHis- and S-tags. Overexpression plasmids, as well as the pET-30a plasmid alone, were transformed into E. coli BL21 (DE3) and Teth39 BdhA and Teth39 AKR were overexpressed (Figure 5). Whole cell lysates were prepared and assayed for aldehyde reductase activity aerobically at 60°C with acetaldehyde, furfural, and 5-HMF. Furfural was also assayed under anaerobic conditions. Teth39 BdhA and Teth39 AKR activities were compared to the vector-only control (Table 4). Neither enzyme showed any activity with NADH as cofactor. Teth39 AKR had minimal activity with both furan aldehydes using NADPH as the cofactor, but not above the vector control. In contrast, Teth39 BdhA showed activity above the vector control using NADPH as the cofactor with both furan aldehydes. The specific activity was 4.97 ± 0.17 U with furfural and 10.06 ± 0.80 U with 5-HMF.

Figure 5
figure 5

Overexpression of Teth39_1597 and Teth39_1598 in Escherichia coli . Teth39_1597 (48.4 kDa) and Teth39_1598 (42.1 kDa) were cloned from T. pseudethanolicus into pET-30a and overexpressed from a T7 promoter after induction with IPTG. Overexpression was determined in whole cell extract by Coomassie-stained SDS-PAGE.

Table 4 Specific activity (μmol/min/mg protein) of putative T. pseudethanolicus 39E aldehyde reductases

Based on these results, overexpression of Teth39 bdhA may increase microbial tolerance to the furan aldehydes furfural and 5-HMF produced during biomass pretreatment, especially in a thermophilic, cellulolytic microbe, such as a Caldicellulosiruptor sp. or Clostridium thermocellum. This approach has proven successful in S. cerevisiae, where overexpression of ADH6 or ADH7 allowed growth with 40 mM 5-HMF, where none was seen in a control [46]. However, the increased expression of Teth39 bdhA in T. pseudethanolicus 39E may have influenced other aspects of its furfural proteomic response. For example, Teth39 BdhA is predicted to contain iron-sulfur clusters, and iron-sulfur cluster biosynthetic genes (Teth39_0116-0120) are also upregulated by furfural. Additionally, other oxidoreductases predicted to use NADPH as a cofactor are downregulated (Teth39_0216, Teth39_0217, Teth39_1916), perhaps to increase NADPH availability for furfural reduction. Thus, while overexpression of Teth39 BdhA alone may increase furan aldehyde tolerance, coexpression of other genes may be required to achieve the phenotype overall.

Teth39 BdhA and Teth39 AKR also exhibited butyraldehyde and isobutyraldehyde reduction activity (Table 4), presumably forming butanol and isobutanol, with the Teth39 BdhA activity 2.7-fold higher than that of Teth39 AKR with isobutyraldehyde and 44.2-fold higher with butyraldehyde. Expression of Teth39 bdhA might therefore alternatively be used to biologically produce these higher energy fuel compounds in thermophilic microbes.


Inhibitors, such as furfural and 5-HMF, are generated from common pretreatment methods used for improving the reactivity of lignocellulosic biomass toward enzymatic solubilization. The physiological response and tolerance to these inhibitors must be understood in order to develop improved microorganisms such as C. thermocellum or a Caldicellulosiruptor sp. for thermophilic consolidated bioprocessing of biomass. We measured rapid growth in the presence of 10 to 15 mM furan aldehydes and their detoxification in situ by a thermophilic anaerobe, T. pseudethanolicus 39E. Physiological effects resulting from higher furfural and 5-HMF concentrations included decreased ethanol yield with increases in acetate and lactate production from glucose. A bottom-up proteomics approach was applied to screen for potential enzymes or pathways directly involved in furan aldehyde detoxification. While a number of cellular functions were impacted, including a decrease in expression of ADHs involved in ethanol production, an ADH encoded by Teth39_1597 was upregulated nearly sevenfold in the presence of 15 mM furfural. The enzyme, a putative butanol dehydrogenase, was cloned and overexpressed in E. coli and displayed NADPH-dependent activity against furfural and 5-HMF, suggesting a direct role in detoxifying furan aldehyde inhibitors in situ.


Growth experiments

Thermoanaerobacter pseudethanolicus 39E (DSMZ 2355) was purchased from Deutsche Sammlung von Mikroorganismen und Zellkulturen GmbH (DSMZ, Braunschweig, Germany). All chemicals were purchased from Sigma-Aldrich (St. Louis, MO) unless otherwise indicated. The anaerobic growth medium was prepared using a modified Hungate technique [48] and consisted of 4.5 mM KCl, 4.7 mM NH4Cl, 2.5 mM MgSO4 ∙7H2O, 1.0 mM NaCl, 0.7 mM CaCl2 ∙2H2O, 0.25 mg/mL resazurin, 2.8 mM cysteine-HCl, 6.0 mM NaHCO3, 1 mM potassium phosphate buffer (pH 6.8), 10 mM 3-(N-morpholino)propanesulfonic acid (pH 6.8), 1x Wolfe’s trace minerals [49], 1× Wolfe’s vitamin supplement [49], 0.1% (w/v) yeast extract (Fisher Scientific, Pittsburgh, PA), and 40 mM glucose. Furan aldehydes were added from degassed concentrated stock solutions. Cultures were grown at 65°C from a 1% inoculum in Balch tubes (10 mL) or 125-mL serum bottles (50 mL). Cell growth was monitored by optical density at 600 nm, either directly in the Balch tube using a Spectronic 200 spectrophotometer (Thermo Fisher Scientific, Waltham, MA) or as 200-μL samples transferred to a 96-well plate and read on a Synergy Mx plate reader (BioTek, Winooski, VT). All growth experiments were performed in triplicate.

Small molecule measurement

Glucose, lactate, acetate, and ethanol were measured via HPLC as previously described using an HPX-87H column (Bio-Rad Laboratories, Hercules, CA) at 60°C with detection via refractive index and 5 mM H2SO4 as the mobile phase [50]. Furfural and 5-HMF were measured spectrophotometrically (DU 800, Beckman Coulter, Brea, CA) at 304 and 323 nm, respectively, and concentrations were determined using standard curves generated in growth medium. 2,5-furandimethanol was measured using gas chromatography-mass spectrometry (GC-MS) following trimethylsilylation, with an Agilent 5975C standard quadrupole GC-MS using electron impact ionization (970 eV), as described previously [51].

Proteomic analysis: sample preparation

Cell pellets (10 to 50 mg) from cultures grown for 8 h (early stationary phase) with and without 15 mM furfural were frozen at -80°C prior to preparation. The thawed pellets were resuspended in 1 mL lysis buffer (4% SDS, 100 mM Tris-HCl pH 8.0, 50 mM dithiothreitol) and boiled for 5 min. Samples were then pulse-sonicated (10 s on, 10 s off) for 2 min with an ultrasonic disruptor (Branson, Danbury, CT) at 20% amplitude. The samples were boiled again for 5 min, cleared by centrifugation (21,000 × g, 10 min, RT), and immediately precipitated with 20% trichloroacetic acid overnight at -20°C. Precipitated proteins were washed twice with ice-cold acetone, air dried, and resuspended in 8 M urea in 100 mM Tris-HCl, pH 8.0. The samples were sonicated as before and incubated for 30 min at room temperature (RT). Samples were adjusted to 10 mM dithiothreitol (10 min, RT) and then 10 mM iodoacetamide (10 min, RT, in the dark) to both reduce and block cysteine residues. Sample aliquots containing about 1.5 mg of crude protein were diluted 1:1 (v/v) with 100 mM Tris-HCl, pH 8.0 and 20 mM CaCl2 and digested with sequencing-grade trypsin (Promega, San Luis Obispo, CA) at a 1:75 (w/w) enzyme:protein ratio (16 h, RT). The samples were again diluted 1:1 (v/v) and digested with a second aliquot of trypsin (1:75; w/w) for an additional 4 h. Following digestion, each sample was adjusted to 200 mM NaCl and 0.1% formic acid and filtered through a 10-kDa cutoff spin column filter (Vivaspin 2, GE Healthcare, Pittsburgh, PA). The peptide-enriched flow-through was then quantified by the bicinchoninic acid assay.

Proteomic analysis: measurement of peptides by two-dimensional liquid chromatography-tandem mass spectrometry (LC-MS/MS)

For each sample, 100 μg of peptides were bomb-loaded onto a biphasic MudPIT back column [52] packed with about 5 cm of strong cation exchange (SCX) resin for charge-based separation of peptides followed by about 3 cm C18 reversed phase (RP) for online desalting (Luna and Aqua respectively, Phenomenex, Torrance, CA). Once loaded, the sample columns were washed offline with solvent A (5% acetonitrile, 95% HPLC-grade water, 0.1% formic acid) for 15 min, followed by a gradient to 100% solvent B (70% acetonitrile, 30% HPLC-grade water, 0.1% formic acid) over 30 min. The washed samples were then placed in-line with an in-house pulled nanospray emitter (100-μm inner diameter) packed with 15 cm of C18 RP material and analyzed via 24-h MudPIT two-dimensional LC-MS/MS (eleven salt pulses: 5, 7.5, 10, 12.5, 15, 17.5, 20, 25, 35, 50, 100% of 500 mM ammonium acetate followed by a 100-min gradient to 50% solvent B) with an LTQ XL mass spectrometer (Thermo Fisher Scientific) operating in data-dependent mode. A total of three biological replicate measurements were obtained for each sample.

Proteomic analysis: MS data analysis and evaluation

Acquired MS/MS spectra were assigned to specific peptide sequences using the SEQUEST search algorithm [53] with a FASTA proteome database specific to T. pseudethanolicus. The database contained common contaminant protein entries as well as reversed decoy sequences to assess protein-level false discovery rates. SEQUEST-scored peptide sequence data were filtered and assembled into protein loci using DTASelect [54] with the following conservative criteria: XCorr: +1 = 1.8, +2 = 2.5, +3 = 3.5, DeltCN 0.08, and two peptides per protein identification with at least one required to be unique.

Prior to the semiquantitative analysis, spectral counts were rebalanced to properly distribute non-unique/shared peptides between their potential parent proteins, as previously described [55]. To represent proteins that were sporadically identified across runs (that is, blank/zero values in a portion of the six sample runs), a fraction of a spectral count (0.33) was added to the entire dataset. This distributional shift maintains the originally measured spectral count differential but allows for blank/zero values to be considered in the ensuing statistical analysis [56]. These adjusted values were then converted to normalized spectral counts (nSpC), an extension of the widely recognized normalized spectral abundance factor (NSAF) [57] that is calculated by multiplying the NSAF values by an arbitrary number representative of the number of spectra collected for each run. In this case, the number 50,000 was used for facile data interpretation. Once calculated, an SpC cutoff was applied to all proteins identified in the dataset so that 99% of the total raw SpC assigned to each (summed across all replicates and conditions) remained.

These remaining proteins were log2 transformed, and statistically assessed by ANOVA with JMP Genomics ver. 4.1 (SAS Institute, Cary, NC) to identify proteins in the furfural treated samples that were significantly (P ≤0.05) up- or downregulated relative to the control. These differentially expressed proteins were then hierarchically clustered based on their abundance patterns across all replicates and conditions using the “Fast Ward” algorithm. To remove differences based on raw magnitude differences in nSpC, each protein’s abundance was standardized to represent the number of standard deviations away from the row mean.

Aldehyde reductase cloning and overexpression

T. pseudethanolicus genomic DNA was isolated using the Wizard Genomic DNA Purification kit (Promega, Madison, WI). Teth39_1597 and Teth39_1598 were PCR amplified using Phusion Polymerase (New England Biolabs, Ipswich, MA), cloned into pET-30a (EMD Millipore, Billerica, MA) behind 6xHis- and S-tags, and the final constructs were sequence verified (University of Tennessee, Knoxville, Molecular Biology Resource Facility). Expression plasmids were transformed in BL21 (DE3) Escherichia coli according to the manufacturer’s protocol (Invitrogen, Grand Island, NY). Cells were grown in 50 mL 2×YT medium at 37°C to OD600 of 0.8-1.0, then induced with 100 μM isopropyl-β-D-thiogalactopyranoside (IPTG) and switched to 30°C for 16 h. Cells were harvested at 4°C (3,000 × g, 30 min), washed in 50 mL 100 mM sodium phosphate buffer, pH 7 (buffer A), and resuspended in 5 mL buffer A. The cell suspension (450 μL) was added to 0.1 mm zirconia beads (300 μL) and vortexed 4 × 60 s with 30 s on ice in between. The samples were centrifuged (14,000 × g, 2 min), and the resulting supernatant was used for enzyme assays.

Enyzme assays

The in vitro aldehyde reductase activity was measured as previously described [23] in 100 mM sodium phosphate buffer (pH 7) with 0.2 mM NAD(P)H and 20 mM substrate. The assay mix was added to a sealed 2-mL quartz cuvette and equilibrated to 60°C. Assays were read at 340 nm (DU 800) for 150 s to establish a baseline slope before whole cell lysate was added (1 to 5 μL). The cuvettes were inverted once to mix and read an additional 450 s. The decrease in absorbance over time was calculated and the baseline slope was subtracted. The NAD(P)H concentration was determined using the extinction coefficient (NADH: 6,220 M-1 cm-1, NADPH: 6,270 M-1 cm-1), and the specific activity was calculated as the change in μmol NAD(P)H/min/mg of whole cell lysate protein. Protein concentration was determined using the Bradford assay (Bio-Rad, Hercules, CA) with bovine serum albumin as a standard. The specific activity was measured for the pET-30a vector (control) and overexpressed Teth39_1597 and Teth39_1598 with acetaldehyde, butyraldehyde, isobutyraldehyde, furfural, and 5-HMF.

Additional file





alcohol dehydrogenase


Aldo/keto reductase


analysis of variance


butanol dehydrogenase


consolidated bioprocessing


gas chromatography-mass spectrometry


high performance liquid chromatography




liquid chromatography-mass spectrometry


normalized spectral counts


optical density


reverse phase


spectral counts


tricarboxylic acid


  1. Wiegel J: Formation of ethanol by bacteria - a pledge for the use of extreme thermophilic anaerobic bacteria in industrial ethanol fermentation processes. Experientia 1980, 36: 1434-1446. 10.1007/BF01960144

    Article  Google Scholar 

  2. Taylor MP, Eley KL, Martin S, Tuffin MI, Burton SG, Cowan DA: Thermophilic ethanologenesis: future prospects for second-generation bioethanol production. Trends Biotechnol 2009, 27: 398-405. 10.1016/j.tibtech.2009.03.006

    Article  Google Scholar 

  3. Lynd LR, Weimer PJ, van Zyl WH, Pretorius IS: Microbial cellulose utilization: fundamentals and biotechnology. Microbiol Mol Biol Rev 2002, 66: 739-739. 10.1128/MMBR.66.4.739.2002

    Article  Google Scholar 

  4. Chung D, Cha M, Guss AM, Westpheling J: Direct conversion of plant biomass to ethanol by engineered Caldicellulosiruptor bescii . Proc Natl Acad Sci USA 2014, 111: 8931-8936. 10.1073/pnas.1402210111

    Article  Google Scholar 

  5. Bayer E, Shoham Y, Lamed R: Cellulose-decomposing bacteria and their enzyme systems. In The Prokaryotes. Volume 2. Part 1. Edited by: Dworkin M, Falkow S, Rosenberg E, Schleifer K-H, Stackebrandt E. Springer, New York; 2006:578-617. 10.1007/0-387-30742-7_19

    Chapter  Google Scholar 

  6. Lynd LR, Laser MS, Brandsby D, Dale BE, Davison B, Hamilton R, Himmel M, Keller M, McMillan JD, Sheehan , Wyman CE: How biotech can transform biofuels. Nat Biotechnol 2008, 26: 169-172. 10.1038/nbt0208-169

    Article  Google Scholar 

  7. Olson DG, McBride JE, Shaw AJ, Lynd LR: Recent progress in consolidated bioprocessing. Curr Opin Biotech 2012, 23: 396-405. 10.1016/j.copbio.2011.11.026

    Article  Google Scholar 

  8. Alvira P, Tomas-Pejo E, Ballesteros M, Negro MJ: Pretreatment technologies for an efficient bioethanol production process based on enzymatic hydrolysis: a review. Bioresour Technol 2010, 101: 4851-4861. 10.1016/j.biortech.2009.11.093

    Article  Google Scholar 

  9. Hendriks A, Zeeman G: Pretreatments to enhance the digestibility of lignocellulosic biomass. Bioresour Technol 2009, 100: 10-18. 10.1016/j.biortech.2008.05.027

    Article  Google Scholar 

  10. 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. 10.1007/s00253-004-1642-2

    Article  Google Scholar 

  11. Almarsdottir AR, Sigurbjornsdottir MA, Orlygsson J: Effect of various factors on ethanol yields from lignocellulosic biomass by Thermoanaerobacterium AK17. Biotechnol Bioeng 2012, 109: 686-694. 10.1002/bit.24346

    Article  Google Scholar 

  12. Zaldivar J, Martinez A, Ingram LO: Effect of selected aldehydes on the growth and fermentation of ethanologenic Escherichia coli . Biotechnol Bioeng 1999, 65: 24-33. 10.1002/(SICI)1097-0290(19991005)65:1<24::AID-BIT4>3.0.CO;2-2

    Article  Google Scholar 

  13. Zaldivar J, Martinez A, Ingram LO: Effect of alcohol compounds found in hemicellulose hydrolysate on the growth and fermentation of ethanologenic Escherichia coli . Biotechnol Bioeng 2000, 68: 524-530. 10.1002/(SICI)1097-0290(20000605)68:5<524::AID-BIT6>3.0.CO;2-T

    Article  Google Scholar 

  14. Martinez A, Rodriguez ME, York SW, Preston JF, Ingram LO: Effects of Ca(OH)2 treatments ("overliming") on the composition and toxicity of bagasse hemicellulose hydrolysates. Biotechnol Bioeng 2000, 69: 526-536. 10.1002/1097-0290(20000905)69:5<526::AID-BIT7>3.0.CO;2-E

    Article  Google Scholar 

  15. Shaw AJ, Podkaminer KK, Desai SG, Bardsley JS, Rogers SR, Thorne PG, Hogsett DA, Lynd LR: Metabolic engineering of a thermophilic bacterium to produce ethanol at high yield. Proc Natl Acad Sci USA 2008, 105: 13769-13774. 10.1073/pnas.0801266105

    Article  Google Scholar 

  16. Shao XJ, Raman B, Zhu MJ, Mielenz JR, Brown SD, Guss AM, Lynd LR: Mutant selection and phenotypic and genetic characterization of ethanol-tolerant strains of Clostridium thermocellum . Appl Microbiol Biot 2011, 92: 641-652. 10.1007/s00253-011-3492-z

    Article  Google Scholar 

  17. Brown SD, Guss AM, Karpinets TV, Parks JM, Smolin N, Yang SH, Land ML, Klingeman DM, Bhandiwad A, Rodriguez M, Raman B, Shao XJ, Mielenz JR, Smith JC, Keller M, Lynd LR: Mutant alcohol dehydrogenase leads to improved ethanol tolerance in Clostridium thermocellum . Proc Natl Acad Sci USA 2011, 108: 13752-13757. 10.1073/pnas.1102444108

    Article  Google Scholar 

  18. Liu ZL, Blaschek HP: Biomass conversion inhibitors and in situ detoxification. In Biomass to Biofuels: Strategies for Global Industries. Edited by Vertès AA, Oureshi A, Blaschek HP, Yukawa H. West Sussex, United Kingdom: Wiley; 2010:233–259.

  19. Liu ZL: Molecular mechanisms of yeast tolerance and in situ detoxification of lignocellulose hydrolysates. Appl Microbiol Biotechnol 2011, 90: 809-825. 10.1007/s00253-011-3167-9

    Article  Google Scholar 

  20. Jonsson LJ, Alriksson B, Nilvebrant N-O: Bioconversion of lignocellulose: inhibitors and detoxification. Biotechnol Biofuels 2013, 6: 16. 10.1186/1754-6834-6-16

    Article  Google Scholar 

  21. Mills TY, Sandoval NR, Gill RT: Cellulosic hydrolysate toxicity and tolerance mechanisms in Escherichia coli . Biotechnol Biofuels 2009, 2: 26. 10.1186/1754-6834-2-26

    Article  Google Scholar 

  22. Miller EN, Jarboe LR, Turner PC, Pharkya P, Yomano LP, York SW, Nunn D, Shanmugam KT, Ingram LO: Furfural inhibits growth by limiting sulfur assimilation in ethanologenic Escherichia coli strain LY180. Appl Environ Microbiol 2009, 75: 6132-6141. 10.1128/AEM.01187-09

    Article  Google Scholar 

  23. Miller EN, Jarboe LR, Yomano LP, York SW, Shanmugam KT, Ingram LO: Silencing of NADPH-dependent oxidoreductase genes ( yqhD and dkgA ) in furfural-resistant ethanologenic Escherichia coli . Appl Environ Microbiol 2009, 75: 4315-4323. 10.1128/AEM.00567-09

    Article  Google Scholar 

  24. Ahring BK, Licht D, Schmidt AS, Sommer P, Thomsen AB: Production of ethanol from wet oxidised wheat straw by Thermoanaerobacter mathranii . Bioresour Technol 1999, 68: 3-9. 10.1016/S0960-8524(98)00077-7

    Article  Google Scholar 

  25. Georgieva TI, Ahring BK: Evaluation of continuous ethanol fermentation of dilute-acid corn stover hydrolysate using thermophilic anaerobic bacterium Thermoanaerobacter BG1L1. Appl Microbiol Biotechnol 2007, 77: 61-68. 10.1007/s00253-007-1149-8

    Article  Google Scholar 

  26. Yao S, Mikkelsen MJ: Metabolic engineering to improve ethanol production in Thermoanaerobacter mathranii . Appl Microbiol Biotechnol 2010, 88: 199-208. 10.1007/s00253-010-2703-3

    Article  Google Scholar 

  27. He Q, Hemme CL, Jiang HL, He ZL, Zhou JZ: Mechanisms of enhanced cellulosic bioethanol fermentation by co-cultivation of Clostridium and Thermoanaerobacter spp. Bioresour Technol 2011, 102: 9586-9592. 10.1016/j.biortech.2011.07.098

    Article  Google Scholar 

  28. Onyenwoke RU, Kevbrin VV, Lysenko AM, Wiegel J: Thermoanaerobacter pseudethanolicus sp. nov., a thermophilic heterotrophic anaerobe from Yellowstone National Park. Int J Syst Evol Microbiol 2007, 57: 2191-2193. 10.1099/ijs.0.65051-0

    Article  Google Scholar 

  29. Bischoff KM, Liu SQ, Hughes SR, Rich JO: Fermentation of corn fiber hydrolysate to lactic acid by the moderate thermophile Bacillus coagulans . Biotechnol Lett 2010, 32: 823-828. 10.1007/s10529-010-0222-z

    Article  Google Scholar 

  30. de Vrije T, Bakker RR, Budde MAW, Lai MH, Mars AE, Claassen PAM: Efficient hydrogen production from the lignocellulosic energy crop Miscanthus by the extreme thermophilic bacteria Caldicellulosiruptor saccharolyticus and Thermotoga neapolitana . Biotechnol Biofuels 2009, 2: 12. 10.1186/1754-6834-2-12

    Article  Google Scholar 

  31. Cao GL, Ren NQ, Wang AJ, Guo WQ, Xu JF, Liu BF: Effect of lignocellulose-derived inhibitors on growth and hydrogen production by Thermoanaerobacterium thermosaccharolyticum W16. Int J Hydrogen Energy 2010, 35: 13475-13480. 10.1016/j.ijhydene.2009.11.127

    Article  Google Scholar 

  32. Delgenes JP, Moletta R, Navarro JM: Effects of lignocellulose degradation products on ethanol fermentations of glucose and xylose by Saccharomyces cerevisiae , Zymomonas mobilis , Pichia stipitis , and Candida shehatae . Enzyme Microb Technol 1996, 19: 220-225. 10.1016/0141-0229(95)00237-5

    Article  Google Scholar 

  33. Banerjee N, Bhatnagar R, Viswanathan L: Inhibition of glycolysis by furfural in Saccharomyces cerevisiae . Eur J Appl Microbiol Biotechnol 1981, 11: 226-228. 10.1007/BF00505872

    Article  Google Scholar 

  34. Hemme CL, Fields MW, He Q, Deng Y, Lin L, Tu QC, Mouttaki H, Zhou AF, Feng XY, Zuo Z, Ramsay BD, He Z, Wu L, Van Nostrand J, Xu J, Tang YJ, Wiegel J, Phelps TJ, Zhou J: Correlation of genomic and physiological traits of Thermoanaerobacter species with biofuel yields. Appl Environ Microbiol 2011, 77: 7998-8008. 10.1128/AEM.05677-11

    Article  Google Scholar 

  35. Burdette DS, Vieille C, Zeikus JG: Cloning and expression of the gene encoding the Thermoanaerobacter ethanolicus 39E secondary-alcohol dehydrogenase and biochemical characterization of the enzyme. Biochem J 1996, 316: 115-122.

    Article  Google Scholar 

  36. Burdette D, Zeikus JG: Purification of acetaldehyde dehydrogenase and alcohol dehydrogenases from Thermoanaerobacter ethanolicus 39E and characterization of the secondary-alcohol dehydrogenase (2° Adh) as a bifunctional alcohol-dehydrogenase acetyl-CoA reductive thioesterase. Biochem J 1994, 302: 163-170.

    Article  Google Scholar 

  37. Wilson CM, Yang S, Rodriguez M Jr, Ma Q, Johnson CM, Dice L, Xu Y, Brown SD: Clostridium thermocellum transcriptomic profiles after exposure to furfural or heat stress. Biotechnol Biofuels 2013, 6: 131. 10.1186/1754-6834-6-131

    Article  Google Scholar 

  38. Ma M, Liu ZL: Comparative transcriptome profiling analyses during the lag phase uncover YAP1, PDR1, PDR3, RPN4, and HSF1 as key regulatory genes in genomic adaptation to the lignocellulose derived inhibitor HMF for Saccharomyces cerevisiae . BMC Genomics 2010, 11: 660. 10.1186/1471-2164-11-660

    Article  Google Scholar 

  39. Gorsich SW, Dien BS, Nichols NN, Slininger PJ, Liu ZL, Skory CD: Tolerance to furfural-induced stress is associated with pentose phosphate pathway genes ZWF1, GND1, RPE1, and TKL1 in Saccharomyces cerevisiae . Appl Microbiol Biotechnol 2006, 71: 339-349. 10.1007/s00253-005-0142-3

    Article  Google Scholar 

  40. Lin F-M, Tan Y, Yuan YJ: Temporal quantitative proteomics of Saccharomyces cerevisiae in response to a nonlethal concentration of furfural. Proteomics 2009, 9: 5471-5483. 10.1002/pmic.200900100

    Article  Google Scholar 

  41. Vizcaíno JA, Deutsch EW, Wang R, Csordas A, Reisinger F, Ríos D, Dianes JA, Sun Z, Farrah T, Bandeira N, Binz PA, Xenarios I, Eisenacher M, Mayer G, Gatto L, Campos A, Chalkley RJ, Kraus HJ, Albar JP, Martinez-Bartolomé S, Apweiler R, Omenn GS, Martens L, Jones AR, Hermjakob H: ProteomeXchange provides globally coordinated proteomics data submission and dissemination. Nature Biotechnol 2014, 30: 223-226. 10.1038/nbt.2839

    Article  Google Scholar 

  42. Lin FM, Qiao B, Yuan YJ: Comparative proteomic analysis of tolerance and adaptation of ethanologenic Saccharomyces cerevisiae to furfural, a lignocellulosic inhibitory compound. Appl Environ Microbiol 2009, 75: 3765-3776. 10.1128/AEM.02594-08

    Article  Google Scholar 

  43. Yang S, Giannone RJ, Dice L, Yang ZK, Engle NL, Tschaplinski TJ, Hettich RL, Brown SD: Clostridium thermocellum ATCC27405 transcriptomic, metabolomic and proteomic profiles after ethanol stress. BMC Genomics 2012, 13: 336. 10.1186/1471-2164-13-336

    Article  Google Scholar 

  44. Peng H, Wu G, Shao W: The aldehyde/alcohol dehydrogenase (AdhE) in relation to the ethanol formation in Thermoanaerobacter ethanolicus JW200. Anaerobe 2008, 14: 125-127. 10.1016/j.anaerobe.2007.09.004

    Article  Google Scholar 

  45. Petersson A, Almeida JRM, Modig T, Karhumaa K, Hahn-Hagerdal B, Gorwa-Grauslund MF, Liden G: A 5-hydroxymethyl furfural reducing enzyme encoded by the Saccharomyces cerevisiae ADH6 gene conveys HMF tolerance. Yeast 2006, 23: 455-464. 10.1002/yea.1370

    Article  Google Scholar 

  46. Liu ZL, Moon J, Andersh BJ, Slininger PJ, Weber S: Multiple gene-mediated NAD(P)H-dependent aldehyde reduction is a mechanism of in situ detoxification of furfural and 5-hydroxymethylfurfural by Saccharomyces cerevisiae . Appl Microbiol Biotechnol 2008, 81: 743-753. 10.1007/s00253-008-1702-0

    Article  Google Scholar 

  47. Yao S, Mikkelsen MJ: Identification and overexpression of a bifunctional aldehyde/alcohol dehydrogenase responsible for ethanol production in Thermoanaerobacter mathranii . J Mol Microbiol Biotechnol 2010, 19: 123-133. 10.1159/000321498

    Article  Google Scholar 

  48. Miller TL, Wolin MJ: A serum bottle modification of the Hungate technique for cultivating obligate anaerobes. Appl Microbiol 1974, 27: 985-987.

    Google Scholar 

  49. Wolin EA, Wolin MJ, Wolfe RS: Formation of methane by bacterial extracts. J Biol Chem 1963, 238: 2882-2886.

    Google Scholar 

  50. Hamilton-Brehm SD, Mosher JJ, Vishnivetskaya T, Podar M, Carroll S, Allman S, Phelps TJ, Keller M, Elkins JG: Caldicellulosiruptor obsidiansis sp nov., an anaerobic, extremely thermophilic, cellulolytic bacterium isolated from Obsidian Pool, Yellowstone National Park. Appl Environ Microbiol 2010, 76: 1014-1020. 10.1128/AEM.01903-09

    Article  Google Scholar 

  51. Li YC, Tschaplinski TJ, Engle NL, Hamilton CY, Rodriguez M, Liao JC, Schadt CW, Guss AM, Yang YF, Graham DE: Combined inactivation of the Clostridium cellulolyticum lactate and malate dehydrogenase genes substantially increases ethanol yield from cellulose and switchgrass fermentations. Biotechnol Biofuels 2012, 5: 2. 10.1186/1754-6834-5-2

    Article  Google Scholar 

  52. McDonald WH, Ohi R, Miyamoto DT, Mitchison TJ, Yates JR: Comparison of three directly coupled HPLC MS/MS strategies for identification of proteins from complex mixtures: single-dimension LC-MS/MS, 2-phase MudPIT, and 3-phase MudPIT. Int J Mass Spectrom 2002, 219: 245-251. 10.1016/S1387-3806(02)00563-8

    Article  Google Scholar 

  53. Eng JK, McCormack AL, Yates JR: An approach to correlate tandem mass-spectral data of peptides with amino-acid-sequences in a protein database. J Am Soc Mass Spectrom 1994, 5: 976-989. 10.1016/1044-0305(94)80016-2

    Article  Google Scholar 

  54. Tabb DL, McDonald WH, Yates JR: DTASelect and Contrast: tools for assembling and comparing protein identifications from shotgun proteomics. J Proteome Res 2002, 1: 21-26. 10.1021/pr015504q

    Article  Google Scholar 

  55. Giannone RJ, Huber H, Karpinets T, Heimerl T, Kuper U, Rachel R, Keller M, Hettich RL, Podar M: Proteomic characterization of cellular and molecular processes that enable the Nanoarchaeum equitans - Ignicoccus hospitalis relationship. PLoS One 2011, 6: e22942. 10.1371/journal.pone.0022942

    Article  Google Scholar 

  56. Lochner A, Giannone RJ, Keller M, Antranikian G, Graham DE, Hettich RL: Label-free quantitative proteomics for the extremely thermophilic bacterium Caldicellulosiruptor obsidiansis reveal distinct abundance patterns upon growth on cellobiose, crystalline cellulose, and switchgrass. J Proteome Res 2011, 10: 5302-5314. 10.1021/pr200536j

    Article  Google Scholar 

  57. Zybailov B, Mosley AL, Sardiu ME, Coleman MK, Florens L, Washburn MP: Statistical analysis of membrane proteome expression changes in Saccharomyces cerevisiae . J Proteome Res 2006, 5: 2339-2347. 10.1021/pr060161n

    Article  Google Scholar 

Download references


We thank Sarah Kaufmann for laboratory assistance and Adam M. Guss for providing helpful comments on the manuscript. This work was supported by the BioEnergy Science Center (BESC), which is a U.S. Department of Energy Bioenergy Research Center supported by the Office of Biological and Environmental Research in the DOE Office of Science, Oak Ridge National Laboratory. Oak Ridge National Laboratory is managed by UT-Battelle, LLC, for the U.S. Department of Energy under contract DE-AC05-00OR22725.

This manuscript has been authored by UT-Battelle, LLC, under Contract DE-AC05-00OR22725 with the U.S. Department of Energy. The United States Government retains, and the publisher by accepting the article for publication acknowledges, that the United States Government retains a non-exclusive, paid-up, irrevocable, worldwide license to publish or reproduce the published form of this manuscript, or allow others to do so, for United States Government purposes.

Author information

Authors and Affiliations


Corresponding author

Correspondence to James G Elkins.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

SMC, SDH-B, and JGE conceived and designed the study; SMC, SDH-B, and JGE performed growth experiments; SMC measured furan aldehydes; SMC and JGE performed HPLC analyses; SDH-B prepared samples for proteomic analysis; RJG and RLH performed proteomic experiments and statistical analyses; NLE and TJT performed GC-MS quantifications; SMC performed cloning and enzyme activity assays; SMC, RJG, and JGE analyzed the data and wrote the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material


Additional file 1: Table S1.: Description of data. Proteomics data file for 225 proteins that were differentially regulated in response to a 15 mM furfural treatment. Color-coded hierarchical cluster analysis, cluster number, cluster order, average nSpC per condition, fold change, and P-value are provided for all replicates. The raw output for all detected proteins and statistical analyses are also included. (XLSX 1 MB)

Authors’ original submitted files for images

Rights and permissions

Open Access  This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.

The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.

To view a copy of this licence, visit

The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Clarkson, S.M., Hamilton-Brehm, S.D., Giannone, R.J. et al. A comparative multidimensional LC-MS proteomic analysis reveals mechanisms for furan aldehyde detoxification in Thermoanaerobacter pseudethanolicus 39E. Biotechnol Biofuels 7, 165 (2014).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: