- Open Access
Simultaneous achievement of high ethanol yield and titer in Clostridium thermocellum
Biotechnology for Biofuels volume 9, Article number: 116 (2016)
Biofuel production from plant cell walls offers the potential for sustainable and economically attractive alternatives to petroleum-based products. Fuels from cellulosic biomass are particularly promising, but would benefit from lower processing costs. Clostridium thermocellum can rapidly solubilize and ferment cellulosic biomass, making it a promising candidate microorganism for consolidated bioprocessing for biofuel production, but increases in product yield and titer are still needed.
Here, we started with an engineered C. thermocellum strain where the central metabolic pathways to products other than ethanol had been deleted. After two stages of adaptive evolution, an evolved strain was selected with improved yield and titer. On chemically defined medium with crystalline cellulose as substrate, the evolved strain produced 22.4 ± 1.4 g/L ethanol from 60 g/L cellulose. The resulting yield was about 0.39 gETOH/gGluc eq, which is 75 % of the maximum theoretical yield. Genome resequencing, proteomics, and biochemical analysis were used to examine differences between the original and evolved strains.
A two step selection method successfully improved the ethanol yield and the titer. This evolved strain has the highest ethanol yield and titer reported to date for C. thermocellum, and is an important step in the development of this microbe for industrial applications.
Low-carbon liquid fuels and organic chemicals will likely be derived from plant biomass, of which lignocellulose is the most prominent component. An economic process to overcome the recalcitrance of cellulosic biomass will lead to more widespread utilization of this resource . Consolidated bioprocessing (CBP), in which biomass solubilization and fermentation are accomplished in one step without added enzymes, is a promising configuration for low-cost biological conversion of plant cell walls [2, 3]. Key factors that determine the economic viability of CBP, as with any chemical transformation, are yield, titer, and rate. Target performance metrics for cost-effective production of ethanol from lignocellulose are a yield of >90 % of theoretical, titer of >40 g/L, and rate of 1 g/L/h .
Clostridium thermocellum is a good candidate organism for CBP due to its ability to rapidly ferment cellulose and produce ethanol. Wild-type strains typically produce ethanol with a yield of about 10–35 % of the theoretical maximum . Under controlled fermentation conditions, 100 g/L cellulose can be converted to 15 g/L ethanol in 75 h . Several approaches have been pursued aimed at engineering C. thermocellum to produce ethanol at higher yield , including mutations that block acetate, lactate , H2 , and formate  production. The most successful of these to date involved simultaneous deletion of hpt, hydG, ldh, pfl, and pta-ack . The resulting pathway for conversion of cellobiose to ethanol is redox balanced (Fig. 1). This strain, AG553, produced ethanol at 63.5 % of the theoretical maximum yield from 5 g/L cellulose, with lower yields observed at higher initial cellulose concentrations. Decreasing yield with increasing titer has also been observed for the wild-type strain , as well as other engineered strains of C. thermocellum [7, 11]. The highest reported titer produced by C. thermocellum is 23.6 g/L, although this was for strain I-1-B, an environmental isolate from a hot spring in Japan which is not widely available . Growth of wild-type C. thermocellum is inhibited by ethanol at concentrations of 4–16 g/L ; however, there are several reports showing that this organism can be adapted to grow in the presence of 50 g/L ethanol by serial transfer [13, 14].
To increase yield and titer, we started with strain AG553, using an adaptive evolution strategy to enrich for strains with faster growth in the hope that, due to the constraints on metabolism from gene deletions, ethanol production would improve. After adaptation, the resulting strains were characterized by whole genome sequencing as well as physiological, proteomic, and biochemical approaches to better understand the processes behind improved ethanol production.
Results and discussion
First round of adaptive evolution
Clostridium thermocellum strain AG553 grows much more slowly than wild type , presumably due to metabolic bottlenecks or imbalances. To improve the growth rate, the strain was serially transferred approximately daily in rich medium (CTFUD) with 5 g/L cellobiose, such that faster growing cells would begin to dominate the culture. Each transfer consisted of a ~1000-fold dilution, or approximately 10 generations. After 150 transfers (~1500 generations), the increase in ethanol production had slowed (Fig. 2a), and growth rate had significantly improved (Table 1).
The nature of this type of laboratory evolution results in the cultures being a mixture of genetic backgrounds, as different cells within the population acquire different mutations. This diversity of genotypes makes the understanding of the causative mutation difficult and further genetic manipulation impossible. Therefore, the isolation of an individual genetic background strain was necessary. An individual member of the population was isolated via single colony purification after transfer #150 and named AG601.
To determine the effect on ethanol titer, AG601 was grown in minimal medium with increasing substrate concentration from 5 to 50 g/L. However, preliminary experiments showed a decrease in yield as substrate loading increased (Additional file 1). Therefore, we performed a second round of adaptive evolution with a higher initial substrate concentration.
Second round of adaptive evolution
In this second round of adaptive evolution, we increased the substrate concentration to 50 g/L and switched from a rich medium (CTFUD) to a defined medium (MTC5) to avoid potential problems with auxotrophy and more closely mimic industrial conditions. The strain was transferred once per week, and ethanol titer was measured every 3rd transfer (Fig. 2b). After 13 serial transfers (~1000 generations total), the ethanol titer had increased from 9.5 ± 0.8 to 22.1 ± 1.2 g/L. This culture was used for single colony purification to isolate a pure genetic background, and one isolate was named LL1210.
Comparison of evolved strains on cellobiose
To examine strain improvements via evolution, the wild type and resulting three strains (AG553, AG601, and LL1210) were cultivated in serum bottles in defined medium. Maximum growth rate was determined on 5 g/L cellobiose (Table 1). The wild-type strain had the fastest growth rate, while the growth rate of the unevolved strain AG553 was the slowest. Each round of selection resulted in a strain with a faster growth rate than the previous culture.
Fermentation products were then measured after growth in batch culture in serum bottles on 50 g/L cellobiose (Fig. 3). Wild-type C. thermocellum had low ethanol yield and titer and produced large quantities of glucose. This is evidence of metabolic inhibition, where enzymatic hydrolysis of cellobiose continues even after C. thermocellum sugar catabolism stops. The unevolved strain AG553 consumed very little cellobiose, so even though the yield was higher than the wild-type strain, the titer was similar to that of wild type. The first round of selection (from AG553 to AG601) resulted in dramatically improved cellobiose consumption, and improved ethanol production, but about half of the cellobiose was still being converted to glucose. The second round of selection (from AG601 to LL1210) showed further improvements in ethanol production as well as reduced accumulation of glucose. Compared to the unevolved strain (AG553), strain LL1210 showed a 5.8-fold higher ethanol titer.
Batch culture of final evolved strain (LL1210) grown on cellulose
Clostridium thermocellum is noted for its ability to rapidly consume cellulose , and we wanted to confirm that we had not impacted this ability during selection. In a bioreactor batch fermentation of 60 g/L cellulose (Avicel PH105), 95 % of the substrate was consumed (Fig. 4), with a final ethanol titer of 22.4 ± 1.4 g/L and yield of 0.39 gEtOH/gGlu eq, which is 75 % of the maximum theoretical yield. The production of by-products including pyruvate, lactate, formate, and acetate were all less than 0.5 g/L. The cellulose consumption rate was 1.4 ± 0.2 g/L/h which is close to the cellobiose consumption rate of 1.6 ± 0.3 g/L/h.
The carbon recovery based on the T = 122 h sample was calculated, and the total carbon recovery was 99.8 % (Table 2). Besides biomass (7.2 %), the main by-product of the fermentation was extracellular amino acids which accounted for 10.1 % of the carbon. Valine (1.63 ± 0.18 g/L), alanine (0.81 ± 0.08 g/L), glutamate (0.70 ± 0.07 g/L), and threonine (0.69 ± 0.10 g/L) were the most abundant amino acids. Together, this demonstrates that LL1210 has retained the ability to efficiently utilize cellulose as a carbon and energy source.
To try to reach even higher titers, strain LL1210 was further grown on 95 g/L Avicel in a bioreactor, achieving 98 % solubilization. The ethanol titer increased from 22.4 ± 1.4 to 26.7 ± 0.9 g/L, but the yield decreased to 0.29 gETOH/gGluc eq (Additional file 2). The fact that the titer from 95 g/L cellulose was similar to the titer from 60 g/L cellulose suggested that the decrease in yield might be due to inhibition by ethanol. To test this hypothesis, strain LL1210 was grown in serum bottle batch cultures with different initial ethanol concentrations from 0 to 20 g/L (Fig. 5). Regardless of the initial ethanol concentration, the final concentration was around 22 g/L, suggesting that ethanol tolerance is currently limiting the ethanol titer in this strain.
Identification of the genes responsible for LL1210 strain adaptive evolution
The key phenotypes we observed in the evolved strains were improved growth rate, reduced glucose accumulation, improved cellobiose consumption, and improved ethanol production. To better understand the genetic basis for these changes, we sequenced the genomes of strains AG553, AG601, and LL1210 and compared them with the wild-type strain (Additional file 3).
All of the targeted genes for deletions (ldh, pta-ack, pfl, hydG, and hpt) were absent in AG553 and its descendants, as expected. During normal strain construction, mutations that spontaneously occur within the strain become fixed during single colony purification. Strain AG553, which was not intentionally evolved, contained a single nucleotide variation (SNV) in the bifunctional alcohol (ADH) and aldehyde dehydrogenase (ALDH), adhE (Clo1313_1798). This SNV originated in strain LL350 (∆hpt ∆hydG) , which is part of the AG553 lineage, and it increases NADPH-dependent ADH activity . A mutation also occurred 114 bp upstream of the ech hydrogenase gene (Clo1313_0575) (Additional file 3). We subsequently saw an increase in abundance of some of the Ech protein subunits (Clo1313_0571-2) (Additional file 4) in proteomics studies described below, suggesting that this SNV may increase ech expression. Clo1313_1831 is annotated as a “ROK domain containing protein” and contains a helix-turn-helix domain and a sugar kinase domain. The presence of a ROK domain often indicates function as a transcriptional repressor . These mutations may be responsible for the observed change in the abundance of proteins Clo1313_1831-3 (Additional file 4; Fig. 6a).
During the evolution from strain AG553 to AG601, a number of mutations occurred in or near genes related to the stress response and could be related to the improvement in strain performance. One such mutation was in the region upstream of DnaK (Clo1313_0933), a chaperone involved in protein folding and the stress response . One might expect this mutation to alter expression of dnaK; however, there was no significant change in DnaK levels (see proteomics data below), suggesting that this mutation might be neutral. Another mutation introduces an amino acid change in one of the two homologs of Spo0A (Clo1313_0637), typically the master regulator of sporulation. While most spore-forming bacteria have a single homolog of spo0A, C. thermocellum has two, Clo1313_0637 and Clo1313_1409, which share 56 % amino acid identity with each other. The Clo1313_1409 protein is 58.3 % identical to the well-characterized Bacillus subtilis Spo0A, whereas Clo1313_0637 is only 49 % identical. Previous deletion of Clo1313_1409 completely eliminated sporulation in C. thermocellum , suggesting that it might be the primary controller of sporulation. Interestingly, a mutant strain of C. thermocellum ATCC27405 that was evolved to be tolerant to poplar hydrolysate was found to have a mutation in the second spo0A homolog Cthe_3087  (99 % amino acid identity to Clo1313_0637), and we also found a mutation in this gene, suggesting that the second Spo0A may also play a role in coordinating a stress response, such as sporulation. This mutation could allow the cell to continue to metabolize sugar and grow even under stressful conditions. Further investigation into the role of Clo1313_0637 may help elucidate its function. Another potentially important mutation in AG601 is in the termination factor rho (Clo1313_2832). While this gene is typically essential for growth, SNVs in rho in Escherichia coli helped to confer resistance to ethanol by enhancing transcriptional read-through and altering gene expression . A similar mechanism may be at play here. Interestingly, two mutations were identified in CRISPR/Cas genes Cas2 and Cas4, which are both predicted to be involved in the acquisition of new CRISPR spacers to confer resistance to new foreign DNA, though the relation of these mutations to improved fermentation performance, if any, is unknown.
Additional mutations occurred during the evolution from strain AG601 to LL1210. One such mutation is in spoIIR (Clo1313_2594). During sporulation in Bacillus subtilis, SpoIIR is a regulatory protein that is required for activation of SigE in the mother cell, eventually leading to mother cell lysis . In strain LL1210, this mutation could help prevent metabolic shutdown, allowing continued fermentation of sugars even under stressful conditions. Additionally, homologs of the Bacillus subtilis mRNA-processing genes ribonuclease Y (Clo1313_1122-1123; ymdAB) could alter mRNA stability, leading to altered protein production and perhaps a more robust stress response or better metabolic flux.
Proteomic and enzymatic profile comparisons between original and evolved strains
To understand how adaptation affected protein expression, proteomic analysis was performed on all samples (Additional file 4). Protein abundance was compared between the wild type and final evolved strain LL1210 (Fig. 6a). Near the bottom horizontal axis, a serial of points represents proteins present in wild type, but absent or expressed at low level in strain LL1210. The most abundant of these are, in fact, the targeted gene deletions. Near the left vertical axis, another serial of points represents proteins present in LL1210, but absent or expressed at low level in wild type. Two of the most interesting proteins in this group are Clo1313_1831 and Clo1313_1832. Clo1313_1831 is the ROK-domain protein described above, and Clo1313_1832 is annotated as a phosphofructokinase (Pfk) B-domain protein, suggesting that it may be involved in conversion of fructose 6-phosphate to fructose 1, 6-bisphosphate during glycolysis. Indeed, PFK activity is increased in the evolved strains (Table 3), which might explain their improved growth.
Interestingly, several cellulosomal proteins also show increased expression (Fig. 6a), including CipA (Clo1313_0627), OlpB (Clo1313_0628), Orf2 (Clo1313_0629), and OlpA (also called SdbA, Clo1313_0630). These proteins encode the primary and secondary scaffoldin components of the cellulosome and have been shown to play a key role in cellulose consumption [23–25]. It is possible that the increased abundance of these proteins plays a role in the ability of the evolved strains to consume high concentrations of cellulose.
To understand how changes in protein abundance may have affected metabolism, we selected a subset of proteins thought to participate in glycolysis, electron transfer, and ethanol production, and compared their abundance in the two evolved strains (AG601 and LL1210) against their abundance in the parent strain (AG553) (Fig. 6b). In general, these metabolic enzymes exhibited less variation compared with the complete set of proteins (Fig. 6a, b). Two proteins, in particular, AdhE (Clo1313_1798) and GapDH (Clo1313_2095), are notable. Both exhibit a twofold increase despite the fact that they are among the most highly abundant proteins to begin with. AdhE is a bifunctional alcohol and aldehyde dehydrogenase . In accordance with the proteomic measurements, alcohol dehydrogenase (ADH) activity was assessed and also found to have increased about twofold compared with the parent strain (Table 3). High levels of ALDH and ADH activity may be directly linked to the improvements in ethanol titer in these strains. High level of GAPDH activity may have a similar effect to the increased levels of PFK, allowing greater flux through glycolysis and therefore better ethanol production.
Another interesting observation is that many proteins whose abundance increased after the first round of selection (strain AG601) subsequently decreased after the second round of selection (strain LL1210). These cases are indicated with red colored circles. Expression of metabolic enzymes needs to be balanced for optimal flux and redox balancing , so the decrease in protein abundance may indicate optimization of these central metabolic pathways.
In this work, an engineered strain of C. thermocellum strain was selected for improved growth. This selection resulted in strains with improvements in both ethanol yield and titer. The final evolved strain produced ethanol at 75 % of the maximum theoretical yield and titer of 22.4 g/L. The main products produced other than ethanol and cells were amino acids. If the carbon and electron flux currently used for amino acid production was instead directed toward ethanol production, we would expect the yield to reach 85 % of the theoretical maximum. The most apparent changes in the adapted strain are in genes Clo1313_1831-2, AdhE, and GapDH, and mutations related to sporulation and transcription. Future experiments may help elucidate the relative importance of these changes. Based on ethanol addition experiments, ethanol titer seems to be limited currently by low ethanol tolerance. Thus, more significant gains in ethanol production might benefit from adaptive selection for improved ethanol tolerance.
Bacterial strains, media and cultivation
Strains used in this study are listed in Table 1. All chemicals were reagent grade and obtained from Sigma-Aldrich (St. Louis, MO) or Fisher Scientific (Pittsburgh, PA) unless indicated otherwise. CTFUD rich medium  and MTC5 defined medium  were used for routine strain maintenance, strain evolution, and fermentation as indicated.
Serum bottle batch cultures were incubated at 55 °C and shaken at 180 rpm. Serum bottles were purged with N2 and sealed with butyl rubber stoppers. There is no gas exit or pH regulation for serum bottle batch culture. Bioreactor fermentations were carried out in 1.5-L (1-L working volume) Sartorius Biostat A-plus Sartorius Stedim (Sartorius Stedim, Bohemia, NY) bioreactors in modified MTC5 medium without MOPS buffer and with 2 g/L urea as the nitrogen source, with the temperature maintained at 55 °C and stirred at 150 rpm. The pH was controlled at 7.0 with a Mettler-Toledo pH probe (Columbus, OH) by the addition of 8N KOH. For vitamin supplementation, aliquots of a 50-fold concentrated vitamin-solution were added freshly. The bioreactor was inoculated with 5 mL −80 °C freezer stock culture grown on 5 g/L Avicel PH105 in MTC (0.5 % v/v). The headspace of the bioreactor was flushed with N2 gas prior to inoculation. The 95 g/L Avicel bioreactor fermentation was performed as previously described . 16S rRNA gene sequences of cell pellets from each fermentation were used to verify culture purity.
For the first round of selection, strain AG553 was inoculated into 5 mL of CTFUD medium with 5 g/L cellobiose as the carbon source in a COY anaerobic chamber (COY Labs, Grass Lake, MI) with an atmosphere of (85 % N2, 10 % CO2, 5 % H2) at 51 °C. The culture was diluted 1000-fold (0.1 % inoculum) by transferring to fresh CTFUD medium, which was performed approximately daily for 150 transfers. Frozen stocks were stored at −80 °C after every 10 transfers. After 150 transfers, the culture was streaked on an agar plate to isolate a single colony from the evolved population. Colonies were visible and distinct after 2 days, and single colonies were selected into 5-mL CTUFD medium. Strain AG601 was one of these isolated colonies.
MTC5 chemically defined medium was used for the second round of selection. For second round of selection, all the experiments were performed in a COY anaerobic chamber (85 % N2, 10 % CO2, and 5 % H2) at 55 °C. Strain AG601 was firstly plated in MTC agar. 50 colonies were selected from the plate and then grown in liquid MTC medium. Their growth curves were measured by BioTek PowerWave XS plate reader (BioTek Instruments Inc., Winooski, VT). The top 20 fastest growing colonies were selected and grown in MTC medium with 50 g/L cellobiose, respectively. The cultures were transferred every week, and the inoculum was diluted ~1000-fold for a total of 13 transfers. Frozen stocks were stored at −80 °C after 12 transfers, and the culture was streaked on an agar plate to isolate a single colony from the evolved population.
Acetate, formate, ethanol, glucose, and residual cellobiose were determined by high pressure liquid chromatography (HPLC, Waters, Milford, MA) with refractive index detection using an Aminex HPX-87H column (Bio-Rad, Hercules, CA) with a 5-mM sulfuric acid solution eluent. Pellet nitrogen was determined using a Shimadzu TOC-VCPH total organic carbon analyzer with added total nitrogen unit (Shimadzu Scientific Instruments, Columbia, MD), calibrated using an acidified glycine standard . 1 mL samples were centrifuged at 15,000×g for 10 min, the supernatant was discarded, and the pellet was rinsed twice using equal volumes of MilliQ water. Residual Avicel PH105 concentration was quantified by quantitative saccharification as previously described . Based on that, Avicel concentration was converted to equivalent glucose concentration, and the ethanol yield was calculated based on the equivalents of glucose. Supernatant protein was determined with the Bradford assay (Thermo Scientific, Rockford, IL) with bovine serum albumin (BSA) (Thermo Scientific, Rockford, IL) as a standard. Secreted amino acids were measured using post-column derivatization with ninhydrin followed by separation and quantification with an Aracus Amino Acid Analyzer (membraPure, Berlin, Germany) using a T111 Li-cation exchange column as previously described .
Clostridium thermocellum cells were harvested by centrifugation at 6000×g for 10 min, the supernatant was decanted, and the pellet was anaerobically suspended in a suitable buffer (dependent on assay). The cells were lysed with Ready-Lyse Lysozyme (Epicentre), and DNase I (New England Biolabs) was added to reduce viscosity. The resulting solution was centrifuged at 10,000×g for 5 min at room temperature, and the supernatant was used as cell-free-extract for enzyme assays.
Enzymes were assayed in an anaerobic chamber (COY Labs, Grass Lake, MI) at 55 °C with a headspace of 85 % N2, 10 % CO2, and 5 % H2 using an Agilent 8453 spectrophotometer. The reaction volume was 1 ml, in reduced-volume quartz cuvettes (part number 29MES10; Precision Cells Inc., NY) with a 1.0-cm path length. All enzyme activities are expressed as μmol of product min−1 (mg of cell extract protein)−1. For each enzyme assay, at least three concentrations of cell extracts were used to confirm that specific activity was proportional to the amount of extract added. Protein concentration was determined using the Bradford method with BSA as the standard.
ADH (EC 220.127.116.11 and EC 18.104.22.168) and ALDH (EC 22.214.171.124 and EC 126.96.36.199) activities were measured based on previously described methods , and the oxidation of NAD(P)H was observed at 340 nm (ε = 6.2 mM-1 cm-1). Phosphofructokinase (EC 188.8.131.52 or EC 184.108.40.206) was assayed by the oxidation of NADH, and the method was described in previous work . GAPDH (EC 220.127.116.11) was measured by BioVision (BioVision Inc., Milpitas, CA) GAPDH Activity Assay Kit.
Carbon balance calculations
Carbon balances were calculated as described previously . The molar concentration of Avicel PH105 was calculated based on glucose monomers with a formula weight of 163 g/mole and 2 C3 equivalents. For glucose, extracellular sugar (non-glucose), and isobutanol, each mole was assumed to be equivalent to 2 C3 units. Ethanol, acetate, malate, and lactate were assumed to be equivalent to 1 C3 unit. For amino acids, the number of C3 units was based on calculations from Stephanopoulos et al. . Alanine (1 C3 equivalent) and valine (2 C3 equivalents) account for the majority of the amino acid carbon. For biomass, 1 mol of pellet carbon was assumed to be equivalent to 1/3 mol of C3 units and based upon measured pellet nitrogen values . Extracellular protein was converted to C3 equivalents by assuming that extracellular protein is 45 % carbon by mass.
Genome resequencing was performed as previously described . Briefly, genomic DNA was submitted to the Joint Genome Institute (JGI) for sequencing with an Illumina MiSeq instrument.
Unamplified libraries were generated using a modified version of Illumina’s standard protocol. 100 ng of DNA was sheared to 500 bp using a focused ultrasonicator (Covaris). The sheared DNA fragments were size selected using SPRI beads (Beckman Coulter). The selected fragments were then end repaired, A tailed, and ligated to Illumina compatible adapters (IDT, Inc) using KAPA-Illumina library creation kit (KAPA Biosystems). Libraries were quantified using KAPA Biosystem’s next-generation sequencing library qPCR kit and run on a Roche LightCycler 480 real-time PCR instrument. The quantified libraries were then multiplexed into pools for sequencing. The pools were loaded and sequenced on the Illumina MiSeq sequencing platform utilizing a MiSeq Reagent Kit v2 (300 cycle) following a 2 × 150 indexed run recipe.
Paired-end reads were generated, with an average read length of 150 bp and paired distance of 500 bp. Raw data were analyzed using CLC Genomics Workbench, version 8.5 (Qiagen, USA). Reads were mapped to the reference genome (NC_017992). Mapping was improved by two rounds of local realignment. The CLC probabilistic variant detection algorithm was used to determine small mutations (single and multiple nucleotide polymorphisms, short insertions, and short deletions). Variants occurring in less than 90 % of the reads and variants that were identical to those of the wild-type strain (i.e., due to errors in the reference sequence) were filtered out. The fraction of the reads containing the mutation is shown in Additional file 3. To determine larger mutations, the CLC indel and structural variant algorithm was run. This tool analyzes unaligned ends of reads and annotates regions where a structural variation may have occurred, which are called breakpoints. Since the read length averaged 150 bp and the minimum mapping fraction was 0.5, a breakpoint can have up to 75 bp of sequence data. The resulting breakpoints were filtered to eliminate those with fewer than ten reads or less than 20 % “not perfectly matched.” The breakpoint sequence was searched with the basic local alignment search tool (BLAST) algorithm  for similarity to known sequences. Pairs of matching left and right breakpoints were considered evidence for structural variations, such as transposon insertions and gene deletions. The fraction of the reads supporting the mutation (left and right breakpoints averaged) is presented in Additional file 3. Raw data are available from the JGI Sequence Read Archive. For strain AG553, the SRA accession number is SRA188084; for strain AG601, the SRA accession number is SRA188071; for strain LL1210, the SRA accession number is SRR3503855.
Cell pellets for the four strains, wild-type C. thermocellum, AG553, AG601, and LL1210 were harvested at log phase and prepared for LC–MS/MS-based proteomic analysis. Briefly, proteins extracted via SDS, boiling, and sonic disruption were precipitated with trichloroacetic acid (TCA), cold acetone-washed, and pelleted as previously described . The acetone-washed protein pellet was resolubilized in urea and treated with dithiothreitol and iodoacetamide to reduce and block disulfide bonds prior to digestion with sequencing-grade trypsin (Sigma-Aldrich). Following proteolysis, tryptic peptides were salted, acidified, and filtered through a 10 kDa MWCO spin column (Vivaspin 2; GE Healthcare) and quantified by BCA assay (Pierce).
For each LC–MS/MS run, 25 µg of peptides were loaded via pressure cell onto a biphasic MudPIT column  for online 2D HPLC separation (strong-cation exchange and reversed-phase) and concurrent analysis via nanospray MS/MS using a hybrid LTQ-Orbitrap XL mass spectrometer (Thermo Scientific) operating in data-dependent acquisition (one full scan at 15 k resolution followed by 10 MS/MS scans in the LTQ, all one µscan; monoisotopic precursor selection; rejection of analytes with an undecipherable charge; dynamic exclusion = 30 s). Eleven salt cuts (25, 30, 35, 40, 45, 50, 65, 80, 100, 175, and 500 mM ammonium acetate) were performed per sample run with each followed by 120 min organic gradient to separate peptides .
Resultant peptide fragmentation spectra (MS/MS) were searched against the C. thermocellum 1313 proteome database concatenated with common contaminants and reversed sequences to control false-discovery rates using MyriMatch v.2.1 . Peptide spectrum matches (PSM) were filtered by IDPicker v.3  using a peptide-level FDR of <1 % per sample run and assigned matched-ion intensities (MIT) based on observed peptide fragment peaks. PSM MITs were summed on a per-peptide basis, and only those uniquely and specifically matching a particular protein were moved onto subsequent analysis with InfernoRDN . Peptide intensity distributions were log2-transformed, normalized across biological replicates by LOESS, and standardized by median absolute deviation and median centering across samples as suggested. Peptide abundance data were then assembled to proteins, scaled appropriately, and outliers removed (RRollup). Protein abundances were then filtered to maintain at least two values in at least one replicate set and missing values imputed using a random distribution of low-level values using Perseus (http://www.perseus-framework.org). Protein abundances were then compared across strains to identify proteins with differential abundance patterns.
bovine serum albumin
glyceraldehyde 3-phosphate dehydrogenase
single nucleotide variation
Lynd LR, Weimer PJ, Van Zyl WH, Isak S. Microbial cellulose utilization: fundamentals and biotechnology. Microbiol Mol Biol Rev. 2002;66:506–77.
Lynd L, van Zyl W, McBride J, Laser M. Consolidated bioprocessing of cellulosic biomass: an update. Curr Opin Biotechnol. 2005;16:577–83.
Olson DG, McBride JE, Shaw AJ, Lynd LR. Recent progress in consolidated bioprocessing. Curr Opin Biotechnol. 2012;23:396–405.
Dien BS, Cotta MA, Jeffries TW. Bacteria engineered for fuel ethanol production: current status. Appl Microbiol Biotechnol. 2003;63:258–66.
Olson DG, Sparling R, Lynd LR. Ethanol production by engineered thermophiles. Curr Opin Biotechnol. 2015;33:130–41.
Holwerda EK, Thorne PG, Olson DG, Amador-Noguez D, Engle NL, Tschaplinski TJ, van Dijken JP, Lynd LR. The exometabolome of Clostridium thermocellum reveals overflow metabolism at high cellulose loading. Biotechnol Biofuels. 2014;7:155.
Argyros DA, Tripathi SA, Barrett TF, Rogers SR, Feinberg LF, Olson DG, Foden JM, Miller BB, Lynd LR, Hogsett DA, Caiazza NC. High ethanol titers from cellulose by using metabolically engineered thermophilic, anaerobic microbes. Appl Environ Microbiol. 2011;77:8288–94.
Biswas R, Zheng T, Olson DG, Lynd LR, Guss AM. Elimination of hydrogenase active site assembly blocks H2 production and increases ethanol yield in Clostridium thermocellum. Biotechnol Biofuels. 2015;8:20.
Rydzak T, Lynd LR, Guss AM. Elimination of formate production in Clostridium thermocellum. J Ind Microbiol Biotechnol. 2015;42:1263–72.
Papanek B, Biswas R, Rydzak T, Guss AM. Characterization of Clostridium thermocellum strains with disrupted fermentation end-product pathways. Metab Eng. 2015;32:49–54.
Van Der Veen D, Lo J, Brown SD, Johnson CM, Tschaplinski TJ, Martin M, Engle NL, Van Den Berg RA, Argyros AD, Caiazza NC, Guss AM, Lynd LR. Characterization of Clostridium thermocellum strains with disrupted fermentation end-product pathways. J Ind Microbiol Biotechnol. 2013;40:725–34.
Sato K, Tomita M, Yonemura S, Goto S, Sekine K, Okuma E, Takagi Y, Hon-Nami K, Saiki T. Characterization of and ethanol hyper-production by Clostridium thermocellum I-1-B. Biosci Biotechnol Biochem. 1993;57:2116–21.
Brown SD, Guss AM, Karpinets TV, Parks JM, Smolin N, Yang S, Land ML, Klingeman DM, Bhandiwad A, Rodriguez M, Raman B, Shao X, 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–7.
Shao X, Raman B, Zhu M, Mielenz JR, Brown SD, Guss AM, Lynd LR. Mutant selection and phenotypic and genetic characterization of ethanol-tolerant strains of Clostridium thermocellum. Appl Microbiol Biotechnol. 2011;92:641–52.
Holwerda EK, Lynd LR. Testing alternative kinetic models for utilization of crystalline cellulose (Avicel) by batch cultures of Clostridium thermocellum. Biotechnol Bioeng. 2013;110:2389–94.
Zheng T, Olson DG, Tian L, Bomble YJ, Himmel ME, Lo J, Hon S, Shaw AJ, van Dijken JP, Lynd LR. Cofactor specificity of the bifunctional alcohol and aldehyde dehydrogenase (AdhE) in wild-type and mutants of Clostridium thermocellum and Thermoanaerobacterium saccharolyticum. J Bacteriol. 2015;197(15):2610–9.
Kazanov MD, Li X, Gelfand MS, Osterman AL, Rodionov DA. Functional diversification of ROK-family transcriptional regulators of sugar catabolism in the Thermotogae phylum. Nucleic Acids Res. 2013;41:790–803.
Mayer MP, Bukau B. Hsp70 chaperones: cellular functions and molecular mechanism. Cell Mol Life Sci. 2005;62:670–84.
Mearls EB, Lynd LR. The identification of four histidine kinases that influence sporulation in Clostridium thermocellum. Anaerobe. 2014;28:109–19.
Linville JL, Rodriguez M, Land M, Syed MH, Engle NL, Tschaplinski TJ, Mielenz JR, Cox CD. Industrial robustness: understanding the mechanism of tolerance for the populus hydrolysate-tolerant mutant strain of Clostridium thermocellum. PLoS One. 2013;8:1–16.
Haft RJF, Keating DH, Schwaegler T, Schwalbach MS, Vinokur J, Tremaine M, Peters JM, Kotlajich MV, Pohlmann EL, Ong IM, Grass JA, Kiley PJ, Landick R. Correcting direct effects of ethanol on translation and transcription machinery confers ethanol tolerance in bacteria. Proc Natl Acad Sci USA. 2014;111:E2576–85.
Karow LM, Glaser P, Piggot PJ. Identification of a gene, spoIIR, which links the activation of sigma E to the transcriptional activity of sigma F during sporulation in Bacillus subtilis. Proc Natl Acad Sci USA. 1995;92:2012–6.
Demain A, Newcomb M, Wu J. Cellulase, clostridia, and ethanol. Microbiol Mol Biol Rev. 2005;69:124–54.
Olson DG, Giannone RJ, Hettich RL, Lynd LR. Role of the CipA scaffoldin protein in cellulose solubilization, as determined by targeted gene deletion and complementation in Clostridium thermocellum. J Bacteriol. 2013;195:733–9.
Xu Q, Resch MG, Podkaminer K, Yang S, Baker JO, Donohoe BS, Wilson C, Klingeman DM, Olson DG, Decker SR, Giannone RJ, Hettich RL, Brown SD, Lynd LR, Bayer EA, Himmel ME, Bomble YJ. Dramatic performance of Clostridium thermocellum explained by its wide range of cellulase modalities. Sci Adv. 2016;2(2):e1501254.
Cheng KK, Lee BS, Masuda T, Ito T, Ikeda K, Hirayama A, Deng L, Dong J, Shimizu K, Soga T, Tomita M, Palsson BO, Robert M. Global metabolic network reorganization by adaptive mutations allows fast growth of Escherichia coli on glycerol. Nat Commun. 2014;5:3233.
Olson DG, Lynd LR. Transformation of Clostridium thermocellum by electroporation, vol. 510. 1st ed. Amsterdam: Elsevier Inc; 2012.
Hogsett DA. Cellulose hydrolysis and fermentation by Clostridium thermocellum for the production of ethanol. Hanover: Dartmouth College; 1995.
Holwerda EK, Ellis LD, Lynd LR. Development and evaluation of methods to infer biosynthesis and substrate consumption in cultures of cellulolytic microorganisms. Biotechnol Bioeng. 2013;110:2380–8.
Zhou J, Olson DG, Argyros DA, Deng Y, van Gulik WM, van Dijken JP, Lynd LR. Atypical glycolysis in Clostridium thermocellum. Appl Environ Microbiol. 2013;79:3000–8.
Ellis LD, Holwerda EK, Hogsett D, Rogers S, Shao X, Tschaplinski T, Thorne P, Lynd LR. Closing the carbon balance for fermentation by Clostridium thermocellum (ATCC 27405). Bioresour Technol. 2012;103:293–9.
Stephanopoulos G, Aristidou AA, Nielsen J. Metabolic engineering: principles and methodologies. Cambridge: Academic press; 1998.
Zhou J, Olson DG, Lanahan AA, Tian L, Murphy SJL, Lo J, Lynd LR. Physiological roles of pyruvate ferredoxin oxidoreductase and pyruvate formate-lyase in Thermoanaerobacterium saccharolyticum JW/SL-YS485. Biotechnol Biofuels. 2015;8:138.
Altschul SF, Gish W, Pennsylvania T, Park U. Basic local alignment search tool. J Mol Biol. 1990;215(3):403–10.
Giannone RJ, Wurch LL, Podar M, Hettich RL. Rescuing those left behind: recovering and characterizing underdigested membrane and hydrophobic proteins to enhance proteome measurement depth. Anal Chem. 2015;87:7720–8.
Delahunty CM, Yates JR. MudPIT: multidimensional protein identification technology. Biotechniques. 2007;43:563–7 (p. 565 passim).
Giannone RJ, Wurch LL, Heimerl T, Martin S, Yang Z, Huber H, Rachel R, Hettich RL, Podar M. Life on the edge: functional genomic response of Ignicoccus hospitalis to the presence of Nanoarchaeum equitans. ISME J. 2014;9:101–14.
Tabb DL, Fernando CG, Chambers MC. MyriMatch:highly accurate tandem mass spectral peptide identificaiton by multivariate hypergeometric analysis. J Proteome Res. 2008;6:654–61.
Ma ZQ, Dasari S, Chambers MC, Litton MD, Sobecki SM, Zimmerman LJ, Halvey PJ, Schilling B, Drake PM, Gibson BW, Tabb DL. IDPicker 2.0: improved protein assembly with high discrimination peptide identification filtering. J Proteome Res. 2009;8:3872–81.
Taverner T, Karpievitch YV, Polpitiya AD, Brown JN, Dabney AR, Anderson GA, Smith RD. DanteR: an extensible R-based tool for quantitative analysis of -omics data. Bioinformatics. 2012;28:2404–6.
LT, BP and DGO designed experiments, performed experiments, analyzed data and wrote the manuscript. TR measured the secreted amino acids. EH advised on fermentation experiments and carried out the 95 g/L Avicel fermentation. TZ and JZ carried out the ADH, ALDH and PFOR enzyme assay. NJ assisted with the first strain evolution and initial characterization of AG601. MM prepared the genome samples for sequencing. RJG and RLH performed protein abundance measurements and assisted with interpretation of the data. AMG and LRL designed experiments, and wrote the manuscript. All authors read and approved the final manuscript.
The genomic resequencing work conducted by the U.S. Department of Energy Joint Genome Institute, a DOE Office of Science User Facility, is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231.
Notice: This manuscript has been authored by Dartmouth College under Subcontract No. 4000115284 and Contract No. DE-AC05-00OR22725 with 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 world-wide license to publish or reproduce the published form of this manuscript or allow others to do so, for United States Government purposes.
Lee R. Lynd is a founder of the Enchi Corporation, which has a financial interest in C. thermocellum. No non-financial competing interests exist for any of the authors.
The BioEnergy Science Center 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.
Liang Tian and Beth Papanek contributed equally to this work