Clostridium thermocellum LL1210 pH homeostasis mechanisms informed by transcriptomics and metabolomics
© The Author(s) 2018
Received: 11 November 2017
Accepted: 24 March 2018
Published: 5 April 2018
Clostridium (Ruminiclostridium) thermocellum is a model fermentative anaerobic thermophile being studied and engineered for consolidated bioprocessing of lignocellulosic feedstocks into fuels and chemicals. Engineering efforts have resulted in significant improvements in ethanol yields and titers although further advances are required to make the bacterium industry-ready. For instance, fermentations at lower pH could enable co-culturing with microbes that have lower pH optima, augment productivity, and reduce buffering cost. C. thermocellum is typically grown at neutral pH, and little is known about its pH limits or pH homeostasis mechanisms. To better understand C. thermocellum pH homeostasis we grew strain LL1210 (C. thermocellum DSM1313 Δhpt ΔhydG Δldh Δpfl Δpta-ack), currently the highest ethanol producing strain of C. thermocellum, at different pH values in chemostat culture and applied systems biology tools.
Clostridium thermocellum LL1210 was found to be growth-limited below pH 6.24 at a dilution rate of 0.1 h−1. F1F0-ATPase gene expression was upregulated while many ATP-utilizing enzymes and pathways were downregulated at pH 6.24. These included most flagella biosynthesis genes, genes for chemotaxis, and other motility-related genes (> 50) as well as sulfate transport and reduction, nitrate transport and nitrogen fixation, and fatty acid biosynthesis genes. Clustering and enrichment of differentially expressed genes at pH values 6.48, pH 6.24 and pH 6.12 (washout conditions) compared to pH 6.98 showed inverse differential expression patterns between the F1F0-ATPase and genes for other ATP-utilizing enzymes. At and below pH 6.24, amino acids including glutamate and valine; long-chain fatty acids, their iso-counterparts and glycerol conjugates; glycolysis intermediates 3-phosphoglycerate, glucose 6-phosphate, and glucose accumulated intracellularly. Glutamate was 267 times more abundant in cells at pH 6.24 compared to pH 6.98, and intercellular concentration reached 1.8 μmol/g pellet at pH 5.80 (stopped flow).
Clostridium thermocellum LL1210 can grow under slightly acidic conditions, similar to limits reported for other strains. This foundational study provides a detailed characterization of a relatively acid-intolerant bacterium and provides genetic targets for strain improvement. Future studies should examine adding gene functions used by more acid-tolerant bacteria for improved pH homeostasis at acidic pH values.
Clostridium (Ruminiclostridium) thermocellum is a model fermentative anaerobic thermophile being studied and engineered for consolidated bioprocessing of lignocellulosic feedstocks into fuels and chemicals [1–3]. C. thermocellum forms specialized biofilms for growth on cellulose  and produces mobile  and cell-bound enzyme complexes termed cellulosomes  that mediate deconstruction of lignocellulosic substrates to short, β1,4-linked glucose oligosaccharides for fermentation to ethanol and organic acids. Wild type C. thermocellum fermentation products include acetic acid, lactic acid, formic acid, H2, ethanol, and amino acids such as valine , with additional products being made at high substrate loadings as a result of overflow metabolism .
Organic acid production decreases biofuel yields and acidifies culture medium, each of which are undesirable for an industrial process. Therefore, metabolic engineering efforts have focused on modifying carbon metabolism to decrease production of acidic co-products and increase ethanol production [9–17]. In the highest yielding C. thermocellum strain published to date, strain LL1210, pathways for the synthesis of lactic acid, acetic acid, formic acid, and most H2 production were eliminated , followed by adaptive laboratory evolution to improve growth rate and ethanol titer . This allowed the strain to produce 22 g/L ethanol from 60 g/L cellulose, with a maximum theoretical yield for ethanol of 75% although further advances are required for the bacterium to be industry-ready. Ethanol yields and titers of up to 80% and 38 g/L, respectively, have also been achieved with an engineered C. thermocellum strain in co-culture with an engineered Thermoanaerobacterium saccharolyticum strain . While these strains synthesize essentially no organic acids as end products, fermentation of sugars to reduced compounds such as ethanol or butanol results in the production of a more oxidized compound such as bicarbonate, which will also result in acidification of the medium.
Organisms must be tolerant to the fermentation conditions and the products formed during growth, especially at the high substrate levels required for industrial fermentations. The response to ethanol stress and ethanol-tolerant strains of C. thermocellum have been characterized [19, 20]. C. thermocellum inhibition by pentose sugars  and compounds generated within switchgrass fermentations  has also been characterized. However, deep understanding of the C. thermocellum response to low pH remains underexplored.
While the anaerobic cellulolytic bacterium Fibrobacter succinogenes has been evolved for steady-state growth at pH 5.75 , known anaerobic cellulolytic bacteria typically do not grow at pH values lower than pH 6.0 . High-yielding- and highly productive ethanologens such as Saccharomyces cerevisiae and Zymomonas mobilis, which are not natively cellulolytic, have broad pH ranges (pH 2–6.5  and pH 3.5–7.5 , respectively). C. thermocellum is typically cultured in medium maintained at neutral pH and it maintains its intercellular pH between ~ 7.3–8.5 [24, 25]. For different C. thermocellum strains, growth has been reported between pH 5.9 and 8.1 , as low as pH 6.0 , and between pH 6.2 and 7.7 . The lower pH limit of the only C. thermocellum strain that is currently genetically tractable, strain DSM1313, has not been reported.
A major mechanism of low pH toxicity is dissipation of the proton gradient across the cytoplasmic membrane. Weak acids such as acetic acid become protonated at the low pH in the supernatant, resulting in an uncharged molecule that can cross the cytoplasmic membrane. The cytoplasm is more alkaline, so the weak acid becomes deprotonated, effectively transporting a proton from outside the cell to the inside, collapsing the ∆pH and acidifying the cytoplasm. Cytoplasmic pH can be maintained by a variety of mechanisms that include proton pumps, antiporters, production of ammonia via urease and arginine deiminase pathways, and organic acid decarboxylation reactions that cleave a carboxylic acid side chain (pKa ~ 4.7) to bicarbonate (pKa = 6.4 and more volatile); tolerance to lower pH can be further enhanced by macromolecule repair or protection, lipid changes and biofilm formation [29–31]. Understanding and enhancing the mechanisms by which C. thermocellum responds to decreased pH will be critical for the industrial application of C. thermocellum for lignocellulosic biofuel production, but those mechanisms are not currently known. Therefore, to better understand the physiological response of C. thermocellum to acidic conditions, we cultured strain LL1210 at different pH values in chemostats and applied systems biology tools.
C. thermocellum LL1210 chemostat growth and lower pH limits in MTC medium
Samples at dilution rate of 0.1 h−1 were withdrawn from each bioreactor at pH 7 (actual pH = 6.98), pH 6.5 (actual pH = 6.48), pH 6.25 (actual pH = 6.24), and at apparent washout pH 6.15 (actual pH = 6.12) for analysis by transcriptomics and metabolomics (Fig. 1). Additional metabolomics samples were taken at pH 5.80 at a lower dilution rate (0.05 h−1) to help elucidate the physiological changes.
pH homeostasis mechanisms inferred by differential gene expression
Previous studies have identified genes commonly expressed during acid stress by Gram-positive neutrophilic bacteria , which includes genes for amino acid decarboxylases/antiporter systems, DnaK, GroEL, ClpP protease, urease, F1F0-ATPase, arginine deiminase (ADI), agmatine deiminase (AgDi), and cyclopropane fatty acid (CFA) synthase. Most of these gene systems were not differentially expressed in C. thermocellum and genes for ADI, AgDi, or CFA are not annotated in the genome. Rather than increasing in expression at low pH, several decarboxylases were downregulated at pH 6.24 and 6.12, including aspartate 1-decarboxylase (Clo1313_1318), diaminopimelate decarboxylase (Clo1313_1540), orotidine 5′-phosphate decarboxylase (Clo1313_1266), putative oxaloacetate decarboxylase (Clo1313_1523), and a putative sodium pump decarboxylase γ-subunit (Clo1313_1525). The only decarboxylase that was upregulated at any pH was the S-adenosylmethionine decarboxylase proenzyme (Clo1313_1509), which produces the aminopropyl substrate need for synthesis of the polyamines spermidine and spermine. However, expression of the agmatinase and spermidine synthase (Clo1313_1529-1530) that synthesize the spermidine and spermine were either not differentially expressed or showed decreased expression at lower pH. Then again, genes for a putative spermidine/putrescine ABC transporter (Clo1313_1471-5) were upregulated at pH 6.24 and 6.12 (Additional file 2). Because a polyamine transporter and biosynthetic gene were upregulated, and polyamines can decrease membrane permeability and may protect cells in acid environments , their potential to augment C. thermocellum LL1210 growth at lower pH values was evaluated. Medium was supplemented with 100 μM of the polyamines putrescine, spermidine, spermine, or metabolic precursor arginine and the terminal pH was determined after 144 h. C. thermocellum LL1210 fermentations with and without amendments initiated at pH 7 and 6.75 had similar growth profiles (Additional file 3).
Downregulation of the GLDH (Clo1313_1847) and GOGAT (Clo1313_2032-3, Clo1313_2035-6) genes were observed at pH 6.12, but no significant expression changes were observed for the Type I glutamine synthetase (GS, Clo1313_2031), two Type III glutamine synthetases (Clo1313_1357 and Clo1313_2038), or the glutamate synthetase (Clo1313_1849) (Additional file 2). The type III glutamine synthase (Clo1313_2303) was significantly upregulated at pH 6.24 (Additional file 2). To examine the role of nitrogen metabolism in relation to pH and physiology, deletion strains for GLDH, GOGAT, GS, GS–GOGAT, and an annotated NifH (Clo1313_2339) were characterized using unbuffered (MOPS free) carbon-replete medium (10 g/L) (Additional file 4). Consistent with an earlier study , the parental strain had a higher specific growth rate compared to strain LL1210 (Additional file 4a). The final pH for the GS-mutant- and LL1210 cultures were about pH 6.3 and pH 5.9, respectively, whereas the range for the other strains was ~ pH 5.6–5.8 (Additional file 4b). In the unbuffered system, the GS mutant achieved the highest culture turbidity and produced the second highest titer of ethanol after strain LL1210 (2.0 and 2.7 g/L, respectively) (Additional file 4c). The higher final ethanol titers and pH values for strains GS and LL1210 coincided lower levels of residual substrate.
Genes for two orphan histidine kinases previously predicted to be part of the sporulation cascade in C. thermocellum , a pro-σE processing protease and BofA, inhibitor of the stage IV pro-σK processing protease SpoIVFB, were upregulated in C. thermocellum LL1210 at pH 6.24 compared to 6.98 (Additional file 5). To better understand the role of sporulation in the pH response of C. thermcellum, growth and sporulation of strains LL1210, asporogenous strains M1726 (DSM1313 ∆hpt ∆spo0A1) and M1725 (DSM1313 ∆hpt ∆spo0A1 ∆ldh ∆pta), and the DSM1313 ∆hpt parental strain were examined. C. thermocellum sporulation was inefficient, consistent with an earlier study  and sporulation-deficient mutant strains had neither a growth advantage or disadvantage at lower pH values for the conditions used in this study (Additional file 5).
The molecular chaperone GrpE (Clo1313_0932) and heat shock protein Hsp 20 (Clo1313_0678) were upregulated (2.0- and 2.2-fold) at pH 6.12, which were previously found to be upregulated after furfural shock and heat shock in C. thermocellum ATCC 27405 . Hsp33 (Clo1313_2544) and MutS (Clo1313_1982) were upregulated at pH 6.24 (2.2- and 2.0-fold, respectively). Hsp33 was previously found to be more highly expressed in cellulose-adhered C. thermocellum ATCC 27405 compared to planktonic cells . At pH 6.12, downregulated genes included one of the annotated Clp protease genes (Clo1313_1116) (0.48-fold), an annotated DNA repair proteins RadC (Clo1313_1465) (0.36-fold), and UvrD (Clo1313_2004) (0.31-fold) and a small subunit of an exodeoxyribonuclease VII (Clo1313_1389) (0.50-fold).
Metabolite profiles for low pH conditions
Intercellular metabolites that were significantly higher or lower in concentration at below-standard pH values
Metabolite [retention time (min); key m/z]
pH 6.48 vs 6.98
pH 6.24 vs 6.98α
pH 6.12 vs 6.98
pH 5.80 vs 6.98
10.34 256 156 358 373 N-metabolite
11.93 288 198 172 390 M+ N-metabolite
9.76 174 276 186 248 N-metabolite
10.81 218 202 320 100 N-metabolite
13.30 299 315 357 328 211 P-metabolite
15.23 218 203 244 N-metabolite
8.92 259 274 184
15.2 347 glycoside
6.96 245 260 102 organic acid
10.89 218 191 100 362
11.33 299 211 328 415 175 P-metabolite
7.50 216 231 188 172 N-metabolite
7.23 117 74
15.90 292 uronic acid conjugate
Extracellular amino acid accumulation at growth-limiting pH
At pH 6.24, only valine, alanine, and threonine extracellular concentrations were significantly higher than when pH was 6.98 (α = 0.05) (Additional file 7). Asparagine was significantly lower at pH 6.24 compared to 6.98 (α = 0.05). Extracellular glutamate was not significantly different at pH 6.24 compared to 6.98, although it was significantly lower at pH 6.48 compared to 6.24 (α = 0.05). When cells were shocked with pH 6.12 compared to when they were growing at pH 6.98, there were significantly lower concentrations of all detectable extracellular amino acids except methionine (α = 0.05). The amino acids that had the highest average extracellular concentrations in chemostats were valine (1.5 mM at pH 6.24), asparagine (147 μM at pH 6.98), and alanine (104 μM at pH 6.24). The highest average concentration of glutamate was 32 μM (near the limit of detection) at pH 6.24 (Additional file 7). An ABC-type branch-chain amino acid transporter gene (Clo1313_0822) was upregulated at pH 6.24 and 6.12 and ABC-type polar amino acid transport genes (Clo1313_0794-5) were upregulated at pH 6.24 (Additional file 2).
Metabolic inhibition at a cellular level
Clostridium thermocellum strain LL1210 has deletions in key genes for metabolic pathways to produce acetate, lactate, formate, and H2, and it has among the highest reported ethanol yields and titers for this bacterium. LL1210 grew at pH 6.24 at a dilution rate of 0.1 h−1. At pH values at or below 6.24, we observed gene expression patterns and metabolite levels consistent with five pH homeostasis mechanisms used by C. thermocellum LL1210, namely (i) proton pumping by F1F0-ATPase, (ii) macromolecule protection by chaperones, and (iii) prevention of proton influx through MotAB were mechanisms supported by changes in gene expression (Fig. 2, Additional file 2); (iv) proton pumping by PPi-ases was supported by metabolomics (Table 1); and (v) changes in membrane fatty acid composition was supported by both techniques (Fig. 2, Additional file 2, Table 1).
Differential expression analysis at pH 6.24 revealed that the F1F0-ATPase genes were upregulated (Fig. 2), consistent with a model that the bacterium uses F1F0-ATPases to pump-protons at the expense of ATP (Fig. 3). Increased expression and activity of F1F0-ATPases are associated with pH tolerance in several other bacteria challenged with sublethal pH conditions [29, 30, 38]. In addition to the F1F0-ATPase system, proton-pumping may be possible via a constitutively and highly expressed putative proton-pumping PPi-ase (Clo1313_0823) (top 9% of genes) [39, 40]. Rhodospirillum rubrum and plant proton-pumping PPi-ases are reversible  and the C. thermocellum proton-pumping PPi-ase has been proposed to be important in energy metabolism in this organism ; further studies are required to experimentally test the function in C. thermocellum. Many ATP-utilizing enzymes and pathways were downregulated at lower pH, which included most of the flagella biosynthesis genes, genes for chemotaxis, and other motility-related genes (> 50) as well as sulfate transport and reduction, nitrate transport and nitrogen fixation, and fatty acid biosynthesis genes (Fig. 2, Additional file 2). Genes for sulfate transport/reduction and for nitrogen metabolism have shown differential expression in response to stresses [36, 52] and their roles in maintaining homeostasis and for channeling electrons to augment biofuels productively could better understood through modeling, biochemistry and genetic studies. In this study, clustered and enriched gene ontologies (GO) analyses indicate genes for enzymes using ATP had inverse differential expression patterns compared to F1F0-ATPase genes, which suggests metabolism is redirected to conserve energy (Fig. 2). Indeed, maintaining ∆pH and protecting macromolecules takes precedence over many reactions since critical enzymes for growth cannot function if they are damaged. In addition to a direct burden on ATP, flagella motility may have an indirect burden since protons channeled into the cell through MotAB would need to be pumped back out by to maintain pH homeostasis. Therefore, motAB downregulation at pH 6.24 would support pH homeostasis (Additional file 2, Fig. 3).
Genes involved in polyamine transport and biosynthesis were differentially expressed at lower pH, suggesting a potential functional role for these compounds in pH tolerance. Specifically, the transporter (Clo1313_1471-5) and SAM decarboxylase gene (Clo1313_1509) were upregulated (Additional file 2). Medium supplementation with polyamines or arginine (a polyamine precursor) did not significantly improve growth or enable C. thermocellum to ferment to a lower terminal pH value compared to unsupplemented fermentations (Additional file 3). Instead, a functional role for stabilization of tRNAs during translation  matches well with polyamine biosynthesis and transport gene expression profiles clustering with expression profiles of genes for tRNAs, tRNA modification proteins, rRNA, ribosomal proteins, ribosomal maturation factors, and translational initiation factors (Fig. 2, Additional file 2). Upregulation of translational machinery, tRNAs, and the polyamines that stabilize them, is likely a regulatory response to the abundance of available carbon at pH-limited growth.
Genes involved in sporulation were upregulated at pH 6.24 in chemostat culture. Interestingly, strain LL1210 has a mutation in one of the two copies of the spo0A gene that results in an amino acid change ; Spo0A is generally the master regulator of sporulation, and the role of this mutation has not been established. However, examination of sporulation in different genetic backgrounds at acidic pH fermentation demonstrated that C. thermocellum is inefficient at sporulation, and that LL1210 is similar to the wild type strain. Further, a growth advantage for sporulation mutants in acidic, MOPS-free, carbon-replete medium was not observed (Additional file 5), suggesting that a programmed metabolic shutdown via sporulation is not a major driver of sensitivity to low pH.
Nitrogen metabolism was strongly effected at lower pH values, resulting in increased secretion of amino acids including valine and alanine, which has been seen in other studies [8, 18], and also intracellular accumulation of glutamate with minimal secretion (Table 1, Additional files 5, 7). Some neutrophilic bacteria express glutamate decarboxylase in response to pH stress to convert glutamate into γ-aminobutyric acid; the decarboxylation reaction consumes protons from the cytoplasm which reduces intercellular pH . Heterologous expression of a glutamate decarboxylase in C. thermocellum LL1210 could make use of the accumulating intracellular glutamate, possibly enabling steady-state growth at pH values less than 6.25 with a dilution rate of 0.1 h−1. Glutamate decarboxylases with pH and temperature ranges suitable for C. thermocellum have already been characterized [44, 45], and heterologous expression of glutamate decarboxylase genes with enhanced production of γ-aminobutyric acid in distantly related organisms has been demonstrated . Urease is another target for improvement of pH tolerance in C. thermocellum. Urease activity in neutrophilic bacteria produces ammonia which is protonated intercellularly or exported and protonated to form ammonium thereby raising intracellular or extracellular pH . The urease genes had low expression in LL1210 at all pH values sampled (Additional file 2) in urea-containing growth medium. Previous work has showed that the parental strain (DSM1313 Δhpt) also had low urease expression in urea-containing fermentation medium, but that deletion of the GS gene caused significant upregulation of urease genes . Consistent with higher expression of the urease genes, higher terminal pH was also previously observed for the GS-mutant fermentations compared to the parental strain fermentations . It remains to be determined if higher urease expression increases C. thermocellum’s intracellular pH and enables steady-state growth at pH values below 6.25 with a dilution rate of 0.1 h−1. The GS mutant had less residual cellobiose and glucose in unbuffered medium than LL1210 and other nitrogen assimilation mutant strains tested (Additional file 4). Therefore, deletion of the GS in LL1210 might also improve substrate utilization and reduce or eliminate the need for buffering fermentation media.
Clostridium thermocellum LL1210 can grow under slightly acidic conditions, similar to limits reported for related C. thermocellum strains. In nature, C. thermocellum may be exposed to acidic pH values, and if it is, metabolism will slow or the bacterium could sporulate, which are not aligned with industrial processes. This foundational study provides a detailed characterization of a relatively acid-intolerant bacterium and provides genetic targets for strain improvement. Future studies should examine adding gene functions used by more acid-tolerant bacteria for improved pH homeostasis at acidic pH values.
Bacterial strains, media and growth
Clostridium thermocellum strains DSM1313 ∆hpt , LL1210 , M1726 (LL376) , M1725 (LL375) , and AG1329  were used in this study. The GLDH (Clo1313_1847), GOGAT (Clo1313_2032-6), GS–GOGAT (Clo1313_2031-6), and a nifH (Clo1313_2339) strain were deleted in a C. thermocellum DSM1313 ∆hpt background to create AG1327, AG1328, AG2068, and AG2078, respectively. Mutants were constructed using the procedures described in detail [49, 50]. Strains are available upon request. Primers used to construct strain AG2078 and other strains used in this study as well as SNPs identified from genome sequencing are provided in Additional file 8. Strains were stored frozen at − 80 °C in a modified version of DSM 122 medium, referred to as CTFUD , that contains 50 mM of MOPS (morpholinepropanesulfonic acid) and 10 mM of sodium citrate. Actively growing cultures were transferred twice from freezer stocks prior to batch fermentation experiments and inocula were 2% of the final volume. Cells were routinely cultured at 55 °C in either 50 ml (160 ml serum bottles) or 10 ml (26 ml Balch tubes) of MTC medium as described earlier , except that MOPS was omitted in this study. The LL1210 strain was also cultured in 1 ml volumes in a 48-well plate reader (Biotek Eon, Winooski, VT, USA) inside a Coy anaerobic chamber kept at 55 °C with growth data collected at OD600 nm every 15 min. Growth in MTC-containing 100 μM of spermine (AC132750010, Fisher Scientific, Hampton, NH, USA), spermidine (S0266, Sigma-Aldrich, St. Louis, MO, USA), putrescine (AC112120250, Fisher Scientific) or arginine (A5006, Sigma-Aldrich) was assessed in plate reader assays.
Strain LL1210 was grown at 55 °C in duplicate 750 ml (total vessel capacity 1.3 L) chemostat cultures using water-jacketed BioFlo110 bioreactors (New Brunswick Scientific, Edison, NJ, USA), essentially as described previously in another C. thermocellum chemostat study , except that the medium contained 3 g/L of cellobiose as the carbon source. Briefly, the MTC medium was MOPS free and fed at a dilution rate of 0.1 h−1, except where noted, for 1058 h. Temperature, pH, and agitation speed were monitored and controlled during fermentation using the BioCommand software (version 2.62). Gel-reference pH electrodes Mettler-Toledo (Woburn, MA) calibration, accuracy checks and pH control were maintained as described earlier , except that 1 N HCl and 3 N KOH were used for acid and base additions, respectively. Recalibration of the bioreactor probes was performed at 546 h after pH drifted by 0.08 and 0.12. Samples for omics were taken in 50 ml volumes; then harvested and stored as described in . For transcriptomics, two samples were taken from different chemostats at pH 6.48 and 6.12, and for pH values 6.98 and 6.24 biological duplicates and two technical replicates separated by at least one change of vessel medium were collected. For intracellular metabolomics, samples were taken immediately after transcriptomic samples, and additional biological duplicate samples were taken at pH 5.8 for metabolomics only. Culture turbidity was measured taking optical density readings at 600 nm using a Genesys 20 spectrophotometer (Thermo Fisher Scientific Inc., Waltham, MA, USA).
Cells were enumerated using a modified version of the direct cell counting procedure described in . Briefly, 1 ml of culture samples were preserved for up to 1 week in 0.025 mg/ml (final concentration) paraformaldehyde (Sigma-Aldrich,). Cells were diluted 1:1000 with 0.9% (m/v) saline solution and stained with 25 nM 4′,6-diamidino-2-phenylindole (DAPI) (final concentration), immobilized to a 1.5 cm diameter section of a black MilliQ filter, and then imaged with a Zeiss Axio Imager.M1 with a Plan Apochromat 63×/1.4 Oil Immersion objective and X-Cite Series 120 Exfo fluorescence power source (Carl Zeiss, Germany). 8-bit images were imported into ImageJ (http://imagej.nih.gov/ij/) for manual or automated counting. Manual counting was performed with the cell counter plugin, and automated counting was performed with JAVA scripts (Additional file 1) which assumed cells were visible at a % Area (saturation threshold) of 1% (± 0.1%) with a particle size greater than 5 pixels after binary images had been processed with Watershed. Spherical morphologies were assumed to have a particle size range of 20–200 pixels and circularity of 0.85–1.00 for conservative estimates, and a particle size range of 50–100 pixels and circularity of 0.95–1.00 for liberal estimates. Automated parameters were validated with manual counts of 24 images. Stepwise calculations of standard deviation of progressively more sample counts were performed to ensure the count variance plateaued. Highly refractive mature spores were also observed by light microscopy after heat fixing to 20 μl of stationary cells to glass slides.
Substrate and metabolite analyses
Extracellular fermentation product samples for high-performance liquid chromatography (HPLC) were collected and processed as described in . HPLC data were generated using a Shimadzu Prominence LC-20A Series system (Shimadzu Scientific Instruments, Columbia, MD) fitted with a refractive index detector (model RID-20A) and an Aminex HPX 87H HPLC column (300 × 7.8 mm) (Bio-Rad, Hercules, Dallas, TX, USA) against known standards for cellobiose, glucose, ethanol, acetate, lactate, and formate. Intracellular metabolite fold changes and select metabolite concentrations were determined using a 5975C inert XL gas chromatography–mass spectrometer (Agilent, Santa Clara, CA) . Extracellular amino acid concentrations were quantified by an Aracus Amino Acid Analyzer (membraPure, Berlin, Germany) . Two-tailed paired T tests were performed to determine if amino acid concentrations were significantly different. Homoscedasticity was confirmed by performing Breusch–Pagan and Abridged White’s tests prior to performing T tests.
RNA isolation, cDNA library preparation, sequencing and RNA-seq analysis
Total RNA was isolated from 50 ml of chemostat sample using the procedures described in , and quantity and quality were determined using a Nanodrop instrument (ThermoFisher Scientific Waltham, MA). Depletion, library preparation and sequencing were completed at the Joint Genome Institute (JGI, Walnut Creek, CA). In short, 100 ng of total RNA was depleted of ribosomal material using a Ribo-Zero (TM) rRNA Removal Kit (Epicentre). The rRNA depleted RNA was fragmented and reversed transcribed using random hexamers and SSII (Invitrogen) followed by second strand synthesis. Stranded cDNA libraries were generated using the Illumina Truseq Stranded RNA LT kit, following the manufactures protocols, using ten cycles of PCR. qPCR was used to determine the concentration of the libraries. Libraries were pooled and sequenced on an Illumina HiSeq 2500 platform with version 2 chemistry in a 2 × 151 bp configuration (JGI, Walnut Creek, CA; Project ID 503106). Fastq files were verified for integrity by checksum and read quality checked using FastQC . RNAseq reads were trimmed using Trimmomatic MAXINFO method with target length and strictness parameters set to 40 and 0.8, respectively . Reads were mapped to C. thermocellum DSM1313 (NC_017304, last modified 2015/07/30) using Bowtie2 with the same parameters as the very-sensitive preset option expect the number of mismatches was set to 1 [57, 58], and reads were counted with HTSeq . Principle component analysis (PCA) was used to check for sample variation using JMP Genomics 8 (SAS Institute, Cary, NC, USA). Differential expression analysis was performed with DESeq 2 . Genes were considered significantly differential expressed when their adjusted p values were less than or equal to 0.05 and log2-fold change values were ≥ 1 or ≤ 1. Raw RNA-seq data have been deposited in NCBI SRA under project number PRJNA395926. PCA plots, raw counts, RPKM normalized gene counts, and differential expression are provided in Additional file 2. Genes found to be differentially expressed at pH 6.48, 6.23, and pH 6.12 compared to pH 6.98 were clustered by their expression patterns using the K-means cluster procedure in JMP Genomics 8. Fisher’s exact test was then performed on 7 of the 9 K-means clusters to identify enriched GO terms. Prior to this analysis, GO terms were identified for the entire C. thermocellum DSM1313 genome using the standard workflow in Blast2GO 4.1.5 .
JMW designed and conducted experiments, analyzed data and wrote the manuscript. J-WM contributed to chemostat cultivation and sampling. MRJr contributed to chemostat cultivation and sampling, and HPLC analysis. NLE generated metabolomics data. DMK extracted RNA. TR generated and analyzed extracellular amino acid data and constructed the nifH mutant. MMA contributed nifH mutant construction. TJT analyzed metabolomics data. AMG oversaw nifH mutant construction, assisted with data interpretation and manuscript preparation. SDB designed experiments, analyzed data and wrote the manuscript. All authors read and approved the final manuscript.
This work is supported by the BioEnergy Science Center (BESC), which is a US Department of Energy Bioenergy Research Center supported by the Office of Biological and Environmental Research in the DOE Office of Science. ORNL is managed by UT-Battelle, LLC, Oak Ridge, TN, USA, for the DOE under contract DE-AC05-00OR22725. RNA-Seq data were generated by the US Department of Energy Joint Genome Institute, a DOE Office of Science User Facility, is supported by the Office of Science of the US Department of Energy under Contract No. DE-AC02-05CH11231.
The manuscript has been authored by UT-Battelle, LLC, under Contract No. DE-AC05-00OR22725 with the US Department of Energy. The funders had no role in study design, data collection and interpretation, preparation of the manuscript, or the decision to submit the work for publication. 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. The Department of Energy will provide public access to these results of federally sponsored research in accordance with the DOE Public Access Plan (http://energy.gov/downloads/doe-public-access-plan).
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
This work is supported by the BioEnergy Science Center (BESC), which is a US Department of Energy Bioenergy Research Center supported by the Office of Biological and Environmental Research in the DOE Office of Science. ORNL is managed by UT-Battelle, LLC, Oak Ridge, TN, USA, for the DOE under contract DE-AC05-00OR22725.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Blumer-Schuette SE, Brown SD, Sander KB, Bayer EA, Kataeva I, Zurawski JV, et al. Thermophilic lignocellulose deconstruction. Rev: FEMS Microbiol; 2014.Google Scholar
- Lynd LR, Weimer PJ, van Zyl WH, Pretorius IS. Microbial cellulose utilization: fundamentals and biotechnology. Microbiol Mol Biol Rev. 2002;66:506–77.View ArticleGoogle Scholar
- Lin PP, Mi L, Morioka AH, Yoshino KM, Konishi S, Xu SC, et al. Consolidated bioprocessing of cellulose to isobutanol using Clostridium thermocellum. Metab Eng. 2015;31:44–52.View ArticleGoogle Scholar
- Dumitrache A, Klingeman DM, Natzke J, Rodriguez M Jr, Giannone RJ, Hettich RL, et al. Specialized activities and expression differences for Clostridium thermocellum biofilm and planktonic cells. Sci Rep. 2017;7:43583.View ArticleGoogle Scholar
- Xu Q, Resch MG, Podkaminer K, Yang S, Baker JO, Donohoe BS, et al. Dramatic performance of Clostridium thermocellum explained by its wide range of cellulase modalities. Sci. Adv. 2016;2:e1501254.View ArticleGoogle Scholar
- Bayer EA, Lamed R, White BA, Flints HJ. From cellulosomes to cellulosomics. Chem Rec. 2008;8:364–77.View ArticleGoogle Scholar
- Ellis LD, Holwerda EK, Hogsett D, Rogers S, Shao X, Tschaplinski T, et al. Closing the carbon balance for fermentation by Clostridium thermocellum (ATCC 27405). Bioresour Technol. 2012;103:293–9.View ArticleGoogle Scholar
- Holwerda EK, Thorne PG, Olson DG, Amador-Noguez D, Engle NL, Tschaplinski TJ, et al. The exometabolome of Clostridium thermocellum reveals overflow metabolism at high cellulose loading. Biotechnol Biofuels. 2014;7:155.View ArticleGoogle Scholar
- Tripathi SA, Olson DG, Argyros DA, Miller BB, Barrett TF, Murphy DM, et al. Development of pyrF-based genetic system for targeted gene deletion in Clostridium thermocellum and creation of a pta mutant. Appl Environ Microbiol. 2010;76:6591–9.View ArticleGoogle Scholar
- Argyros DA, Tripathi SA, Barrett TF, Rogers SR, Feinberg LF, Olson DG, et al. High ethanol titers from cellulose by using metabolically engineered thermophilic, anaerobic microbes. Appl Environ Microbiol. 2011;77:8288–94.View ArticleGoogle Scholar
- Rydzak T, Lynd LR, Guss AM. Elimination of formate production in Clostridium thermocellum. J Ind Microbiol Biotechnol. 2015;42:1263–72.View ArticleGoogle Scholar
- Papanek B, Biswas R, Rydzak T, Guss AM. Elimination of metabolic pathways to all traditional fermentation products increases ethanol yields in Clostridium thermocellum. Metab Eng. 2015;32:49–54.View ArticleGoogle Scholar
- Biswas R, Prabhu S, Lynd LR, Guss AM. Increase in ethanol yield via elimination of lactate production in an ethanol-tolerant mutant of Clostridium thermocellum. PLoS ONE. 2014;9:1–7.Google Scholar
- Olson DG, Hörl M, Fuhrer T, Cui J, Zhou J, Maloney MI, et al. Glycolysis without pyruvate kinase in Clostridium thermocellum. Metab Eng. 2017;39:169–80.View ArticleGoogle Scholar
- Biswas R, Wilson CM, Zheng T, Giannone RJ, Dawn M, Olson DG, et al. Elimination of hydrogenase active site assembly blocks H2 production and increases ethanol yield in Clostridium thermocellum. Biotechnol Biofuels. 2015;8:20.View ArticleGoogle Scholar
- Hon S, Olson DG, Holwerda EK, Lanahan AA, Murphy SJL, Maloney MI, et al. The ethanol pathway from Thermoanaerobacterium saccharolyticum improves ethanol production in Clostridium thermocellum. Metab Eng. 2017;42:175–84.View ArticleGoogle Scholar
- Lo J, Olson DG, Murphy SJL, Tian L, Hon S, Lanahan A, et al. Engineering electron metabolism to increase ethanol production in Clostridium thermocellum. Metab Eng. 2017;39:71–9.View ArticleGoogle Scholar
- Tian L, Papanek B, Olson DG, Rydzak T, Holwerda EK, Zheng T, et al. Simultaneous achievement of high ethanol yield and titer in Clostridium thermocellum. Biotechnol Biofuels. 2016;9:116.View ArticleGoogle Scholar
- Brown SD, Guss AM, Karpinets TV, Parks JM, Smolin N, Yang S, et al. Mutant alcohol dehydrogenase leads to improved ethanol tolerance in Clostridium thermocellum. Proc Natl Acad Sci. 2011;108:13752–7.View ArticleGoogle Scholar
- Shao X, Raman B, Zhu M, Mielenz JR, Brown SD, Guss AM, et al. Mutant selection and phenotypic and genetic characterization of ethanol-tolerant strains of Clostridium thermocellum. Appl Microbiol Biotechnol. 2011;92:641–52.View ArticleGoogle Scholar
- Verbeke TJ, Giannone RJ, Klingeman DM, Engle NL, Rydzak T, Guss AM, et al. Pentose sugars inhibit metabolism and increase expression of an AgrD-type cyclic pentapeptide in Clostridium thermocellum. Sci Rep. 2017;7:43355.View ArticleGoogle Scholar
- Poudel S, Giannone RJ, Rodriguez M, Raman B, Martin MZ, Engle NL, et al. Integrated omics analyses reveal the details of metabolic adaptation of Clostridium thermocellum to lignocellulose-derived growth inhibitors released during the deconstruction of switchgrass. Biotechnol Biofuels. 2017;10:14.View ArticleGoogle Scholar
- Wu C-W, Spike T, Klingeman DM, Rodriguez M, Bremer VR, Brown SD. Generation and characterization of acid tolerant Fibrobacter succinogenes S85. Sci Rep. 2017;7:2277.View ArticleGoogle Scholar
- Nochur SV, Demain AL, Roberts MF. Carbohydrate utilization by Clostridium thermocellum: importance of internal pH in regulating growth. Enzyme Microb Technol. 1992;14:338–49.View ArticleGoogle Scholar
- Herrero AA, Gomez RF, Snedecor B, Tolman CJ, Roberts MF. Growth inhibition of Clostridium thermocellum by carboxylic acids: a mechanism based on uncoupling by weak acids. Appl Microbiol Biotechnol. 1985;22:53–62.View ArticleGoogle Scholar
- Mori Y. Characterization of symbiotic coculture of Clostridium thermohydrosulfuricum YM3 and Clostridium thermocellum YM4. Appl Environ Microbiol. 1990;56:37–42.Google Scholar
- Koeck DE, Zverlov VV, Liebl W, Schwarz WH. Comparative genotyping of Clostridium thermocellum strains isolated from biogas plants: genetic markers and characterization of cellulolytic potential. Microbiol: Syst Appl; 2014.Google Scholar
- Freier D, Mothershed CP, Wiegel J. Characterization of Clostridium thermocellum JW20. Appl Environ Microbiol. 1988;54:204–11.Google Scholar
- Cotter PD, Hill C. Surviving the acid test: responses of Gram-positive bacteria to low pH. Microbiol Mol Biol Rev. 2003;67:429–53.View ArticleGoogle Scholar
- Lund P, Tramonti A, De Biase D. Coping with low pH: molecular strategies in neutralophilic bacteria. FEMS Microbiol Rev. 2014;38:1091–125.View ArticleGoogle Scholar
- Kobayashis H, Murakami N, Unemoto T. Regulation of the Cytoplasmic pH in. 1982;257:13246–52.Google Scholar
- Terahara N, Noguchi Y, Nakamura S, Kami-ike N, Ito M, Namba K, et al. Load- and polysaccharide-dependent activation of the Na+-type MotPS stator in the Bacillus subtilis flagellar motor. Sci Rep. 2017;7:46081.View ArticleGoogle Scholar
- Rhee HJ, Kim EJ, Lee JK. Physiological polyamines: simple primordial stress molecules. J Cell Mol Med. 2007;11:685–703.View ArticleGoogle Scholar
- Mearls EB, Lynd LR. The identification of four histidine kinases that influence sporulation in Clostridium thermocellum. Anaerobe. 2014;28:109–19.View ArticleGoogle Scholar
- Mearls EB, Izquierdo JA, Lynd LR. Formation and characterization of non-growth states in Clostridium thermocellum: spores and L-forms. BMC Microbiol. 2012;12:180.View ArticleGoogle Scholar
- Wilson CM, Yang S, Rodriguez M Jr, Ma Q, Johnson CM, Dice L, et al. Clostridium thermocellum transcriptomic profiles after exposure to furfural or heat stress. Biotechnol Biofuels. 2013;6:131.View ArticleGoogle Scholar
- Sohlenkamp C. Membrane homeostasis in bacteria upon pH challenge. Biogenesis of fatty acids, lipids and membranes. Berlin: Springer; 2017. p. 1–13.View ArticleGoogle Scholar
- Kobayashi H, Suzuki T, Unemoto T. Streptococcal cytoplasmic pH is regulated by changes in amount and activity of a proton-translocating ATPase. J Biol Chem. 1986;261:627–30.Google Scholar
- Zhou J, Olson DG, Argyros DA, Deng Y, van Gulik WM, van Dijken JP, et al. Atypical glycolysis in Clostridium thermocellum. Appl Environ Microbiol. 2013;79:3000–8.View ArticleGoogle Scholar
- Rydzak T, McQueen PD, Krokhin OV, Spicer V, Ezzati P, Dwivedi RC, et al. Proteomic analysis of Clostridium thermocellum core metabolism: relative protein expression profiles and growth phase-dependent changes in protein expression. BMC Microbiol. 2012;12:214.View ArticleGoogle Scholar
- Baltsche M, Schultz A, Baltsche H. H+-PPases : a tightly membrane-bound family. FEBS Lett. 1999;457:527–33.View ArticleGoogle Scholar
- Lynd LR, Guss AM, Himmel ME, Beri D, Herring C, Holwerda EK, et al. Advances in consolidated bioprocessing using Clostridium thermocellum and Thermoanaerobacter saccharolyticum. In: Industrial biotechnology. Weinheim: Wiley-VCH; 2016. p. 365–94.Google Scholar
- Ouameur AA, Bourassa P, Tajmir-Riahi H-A. Probing tRNA interaction with biogenic polyamines. RNA. 2010;16:1968–79.View ArticleGoogle Scholar
- Tsuchiya K, Nishimura K, Iwahara M. Purification and characterization of glutamate decarboxylase from Aspergillus oryzae. Food Sci Technol Res. 2003;9:283–7.View ArticleGoogle Scholar
- Seo M-J, Nam Y-D, Lee S-Y, Park S-L, Yi S-H, Lim S-I. Expression and characterization of a glutamate decarboxylase from Lactobacillus brevis 877G producing γ-aminobutyric acid. Biosci Biotechnol Biochem. 2013;77:853–6.View ArticleGoogle Scholar
- Park K-B, Ji G-E, Park M-S, Oh S-H. Expression of rice glutamate decarboxylase in Bifidobacterium longum enhances gamma-aminobutyric acid production. Biotechnol Lett. 2005;27:1681–4.View ArticleGoogle Scholar
- Rydzak T, Garcia D, Stevenson DM, Sladek M, Klingeman DM, Holwerda EK, et al. Deletion of type I glutamine synthetase deregulates nitrogen metabolism and increases ethanol production in Clostridium thermocellum. Metab Eng. 2017;41:182–91.View ArticleGoogle Scholar
- van der Veen D, Lo J, Brown SD, Johnson CM, Tschaplinski TJ, Martin M, et al. Characterization of Clostridium thermocellum strains with disrupted fermentation end-product pathways. J Ind Microbiol Biotechnol. 2013;40:725–34.View ArticleGoogle Scholar
- Olson DG, Lynd LR. Transformation of Clostridium thermocellum by electroporation. Methods enzymol. 1st ed. Waltham: Academic Press; 2012. p. 317–30.Google Scholar
- Guss AM, Olson DG, Caiazza NC, Lynd LR. Dcm methylation is detrimental to plasmid transformation in Clostridium thermocellum. Biotechnol Biofuels. 2012;5:30.View ArticleGoogle Scholar
- Dumitrache A, Akinosho H, Rodriguez M, Meng X, Yoo CG, Natzke J, et al. Consolidated bioprocessing of Populus using Clostridium (Ruminiclostridium) thermocellum: a case study on the impact of lignin composition and structure. Biotechnol Biofuels. 2016;9:31.View ArticleGoogle Scholar
- Sander K, Wilson CM, Rodriguez M, Klingeman DM, Rydzak T, Davison BH, et al. Clostridium thermocellum DSM 1313 transcriptional responses to redox perturbation. Biotechnol Biofuels. 2015;8:211.View ArticleGoogle Scholar
- Dumitrache A, Wolfaardt GM, Allen DG, Liss SN, Lynd LR. Tracking the cellulolytic activity of Clostridium thermocellum biofilms. Biotechnol Biofuels. 2013;6:175.View ArticleGoogle Scholar
- Wilson CM, Klingeman DM, Schlachter C, Syed MH, Wu C, Guss AM, et al. LacI transcriptional regulatory networks in Clostridium thermocellum DSM1313. Appl Environ Microbiol. 2017;83:e02751–816.View ArticleGoogle Scholar
- Andrews S. FastQC: a quality control tool for high throughput sequence data. 2010. http://www.bioinformatics.babraham.ac.uk/projects/fastqc. Accessed 5 Feb 2016.
- Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30:2114–20.View ArticleGoogle Scholar
- Langmead B, Trapnell C, Pop M, Salzberg SL. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009;10:R25.View ArticleGoogle Scholar
- Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9:357–9.View ArticleGoogle Scholar
- Anders S, Pyl PT, Huber W. HTSeq-A python framework to work with high-throughput sequencing data. Bioinformatics. 2015;31:166–9.View ArticleGoogle Scholar
- Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15:550.View ArticleGoogle Scholar
- Götz S, García-Gómez JM, Terol J, Williams TD, Nagaraj SH, Nueda MJ, et al. High-throughput functional annotation and data mining with the Blast2GO suite. Nucleic Acids Res. 2008;36:3420–35.View ArticleGoogle Scholar