Tobin J. Verbeke

Linking genome content to biofuel production yields: a meta-analysis of major catabolic pathways among select H2 and ethanol-producing bacteria

BackgroundFermentative bacteria offer the potential to convert lignocellulosic waste-streams into biofuels such as hydrogen (H2) and ethanol. Current fermentative H2 and ethanol yields, however, are below theoretical maxima, vary greatly among organisms, and depend on the extent of metabolic pathways utilized. For fermentative H2 and/or ethanol production to become practical, biofuel yields must be increased. We performed a comparative meta-analysis of (i) reported end-product yields, and (ii) genes encoding pyruvate metabolism and end-product synthesis pathways to identify suitable biomarkers for screening a microorganism’s potential of H2 and/or ethanol production, and to identify targets for metabolic engineering to improve biofuel yields. Our interest in H2 and/or ethanol optimization restricted our meta-analysis to organisms with sequenced genomes and limited branched end-product pathways. These included members of the Firmicutes, Euryarchaeota, and Thermotogae.ResultsBioinformatic analysis revealed that the absence of genes encoding acetaldehyde dehydrogenase and bifunctional acetaldehyde/alcohol dehydrogenase (AdhE) in Caldicellulosiruptor, Thermococcus, Pyrococcus, and Thermotoga species coincide with high H2 yields and low ethanol production. Organisms containing genes (or activities) for both ethanol and H2 synthesis pathways (i.e. Caldanaerobacter subterraneus subsp. tengcongensis, Ethanoligenens harbinense, and Clostridium species) had relatively uniform mixed product patterns. The absence of hydrogenases in Geobacillus and Bacillus species did not confer high ethanol production, but rather high lactate production. Only Thermoanaerobacter pseudethanolicus produced relatively high ethanol and low H2 yields. This may be attributed to the presence of genes encoding proteins that promote NADH production. Lactate dehydrogenase and pyruvate:formate lyase are not conducive for ethanol and/or H2 production. While the type(s) of encoded hydrogenases appear to have little impact on H2 production in organisms that do not encode ethanol producing pathways, they do influence reduced end-product yields in those that do.ConclusionsHere we show that composition of genes encoding pathways involved in pyruvate catabolism and end-product synthesis pathways can be used to approximate potential end-product distribution patterns. We have identified a number of genetic biomarkers for streamlining ethanol and H2 producing capabilities. By linking genome content, reaction thermodynamics, and end-product yields, we offer potential targets for optimization of either ethanol or H2 yields through metabolic engineering.

Predicting relatedness of bacterial genomes using the chaperonin-60 universal target (cpn60 UT): Application to Thermoanaerobacter species

D.R. Zeigler determined that the sequence identity of bacterial genomes can be predicted accurately using the sequence identities of a corresponding set of genes that meet certain criteria [32]. This three-gene model for comparing bacterial genome pairs requires the determination of the sequence identities for recN, thdF, and rpoA. This involves the generation of approximately 4.2kb of genomic DNA sequence from each organism to be compared, and also normally requires that oligonucleotide primers be designed for amplification and sequencing based on the sequences of closely related organisms. However, we have developed an analogous mathematical model for predicting the sequence identity of whole genomes based on the sequence identity of the 542-567 base pair chaperonin-60 universal target (cpn60 UT). The cpn60 UT is accessible in nearly all bacterial genomes with a single set of universal primers, and its length is such that it can be completely sequenced in one pair of overlapping sequencing reads via di-deoxy sequencing. These mathematical models were applied to a set of Thermoanaerobacter isolates from a wood chip compost pile and it was shown that both the one-gene cpn60 UT-based model and the three-gene model based on recN, rpoA, and thdF predicted that these isolates could be classified as Thermoanaerobacter thermohydrosulfuricus. Furthermore, it was found that the genomic prediction model using cpn60 UT gave similar results to whole-genome sequence alignments over a broad range of taxa, suggesting that this method may have general utility for screening isolates and predicting their taxonomic affiliations.

Genomic Evaluation of Thermoanaerobacter spp. for the Construction of Designer Co-Cultures to Improve Lignocellulosic Biofuel Production

The microbial production of ethanol from lignocellulosic biomass is a multi-component process that involves biomass hydrolysis, carbohydrate transport and utilization, and finally, the production of ethanol. Strains of the genus Thermoanaerobacter have been studied for decades due to their innate abilities to produce comparatively high ethanol yields from hemicellulose constituent sugars. However, their inability to hydrolyze cellulose, limits their usefulness in lignocellulosic biofuel production. As such, co-culturing Thermoanaerobacter spp. with cellulolytic organisms is a plausible approach to improving lignocellulose conversion efficiencies and yields of biofuels. To evaluate native lignocellulosic ethanol production capacities relative to competing fermentative end-products, comparative genomic analysis of 11 sequenced Thermoanaerobacter strains, including a de novo genome, Thermoanaerobacter thermohydrosulfuricus WC1, was conducted. Analysis was specifically focused on the genomic potential for each strain to address all aspects of ethanol production mentioned through a consolidated bioprocessing approach. Whole genome functional annotation analysis identified three distinct clades within the genus. The genomes of Clade 1 strains encode the fewest extracellular carbohydrate active enzymes and also show the least diversity in terms of lignocellulose relevant carbohydrate utilization pathways. However, these same strains reportedly are capable of directing a higher proportion of their total carbon flux towards ethanol, rather than non-biofuel end-products, than other Thermoanaerobacter strains. Strains in Clade 2 show the greatest diversity in terms of lignocellulose hydrolysis and utilization, but proportionately produce more non-ethanol end-products than Clade 1 strains. Strains in Clade 3, in which T. thermohydrosulfuricus WC1 is included, show mid-range potential for lignocellulose hydrolysis and utilization, but also exhibit extensive divergence from both Clade 1 and Clade 2 strains in terms of cellular energetics. The potential implications regarding strain selection and suitability for industrial ethanol production through a consolidated bioprocessing co-culturing approach are examined throughout the manuscript.

Thermoanaerobacter thermohydrosulfuricus WC1 Shows Protein Complement Stability during Fermentation of Key Lignocellulose-Derived Substrates

ABSTRACT Thermoanaerobacter spp. have long been considered suitable Clostridium thermocellum coculture partners for improving lignocellulosic biofuel production through consolidated bioprocessing. However, studies using “omic”-based profiling to better understand carbon utilization and biofuel producing pathways have been limited to only a few strains thus far. To better characterize carbon and electron flux pathways in the recently isolated, xylanolytic strain, Thermoanaerobacter thermohydrosulfuricus WC1, label-free quantitative proteomic analyses were combined with metabolic profiling. SWATH-MS proteomic analysis quantified 832 proteins in each of six proteomes isolated from mid-exponential-phase cells grown on xylose, cellobiose, or a mixture of both. Despite encoding genes consistent with a carbon catabolite repression network observed in other Gram-positive organisms, simultaneous consumption of both substrates was observed. Lactate was the major end product of fermentation under all conditions despite the high expression of gene products involved with ethanol and/or acetate synthesis, suggesting that carbon flux in this strain may be controlled via metabolite-based (allosteric) regulation or is constrained by metabolic bottlenecks. Cross-species “omic” comparative analyses confirmed similar expression patterns for end-product-forming gene products across diverse Thermoanaerobacter spp. It also identified differences in cofactor metabolism, which potentially contribute to differences in end-product distribution patterns between the strains analyzed. The analyses presented here improve our understanding of T. thermohydrosulfuricus WC1 metabolism and identify important physiological limitations to be addressed in its development as a biotechnologically relevant strain in ethanologenic designer cocultures through consolidated bioprocessing.

Comparative analysis of carbohydrate active enzymes in Clostridium termitidis CT1112 reveals complex carbohydrate degradation ability.

Clostridium termitidis strain CT1112 is an anaerobic, gram positive, mesophilic, cellulolytic bacillus isolated from the gut of the wood-feeding termite, Nasutitermes lujae. It produces biofuels such as hydrogen and ethanol from cellulose, cellobiose, xylan, xylose, glucose, and other sugars, and therefore could be used for biofuel production from biomass through consolidated bioprocessing. The first step in the production of biofuel from biomass by microorganisms is the hydrolysis of complex carbohydrates present in biomass. This is achieved through the presence of a repertoire of secreted or complexed carbohydrate active enzymes (CAZymes), sometimes organized in an extracellular organelle called cellulosome. To assess the ability and understand the mechanism of polysaccharide hydrolysis in C. termitidis, the recently sequenced strain CT1112 of C. termitidis was analyzed for both CAZymes and cellulosomal components, and compared to other cellulolytic bacteria. A total of 355 CAZyme sequences were identified in C. termitidis, significantly higher than other Clostridial species. Of these, high numbers of glycoside hydrolases (199) and carbohydrate binding modules (95) were identified. The presence of a variety of CAZymes involved with polysaccharide utilization/degradation ability suggests hydrolysis potential for a wide range of polysaccharides. In addition, dockerin-bearing enzymes, cohesion domains and a cellulosomal gene cluster were identified, indicating the presence of potential cellulosome assembly.

Enhanced whole genome sequence and annotation of Clostridium stercorarium DSM8532T using RNA-seq transcriptomics and high-throughput proteomics

BackgroundGrowing interest in cellulolytic clostridia with potential for consolidated biofuels production is mitigated by low conversion of raw substrates to desired end products. Strategies to improve conversion are likely to benefit from emerging techniques to define molecular systems biology of these organisms. Clostridium stercorarium DSM8532T is an anaerobic thermophile with demonstrated high ethanol production on cellulose and hemicellulose. Although several lignocellulolytic enzymes in this organism have been well-characterized, details concerning carbohydrate transporters and central metabolism have not been described. Therefore, the goal of this study is to define an improved whole genome sequence (WGS) for this organism using in-depth molecular profiling by RNA-seq transcriptomics and tandem mass spectrometry-based proteomics.ResultsA paired-end Roche/454 WGS assembly was closed through application of an in silico algorithm designed to resolve repetitive sequence regions, resulting in a circular replicon with one gap and a region of 2 kilobases with 10 ambiguous bases. RNA-seq transcriptomics resulted in nearly complete coverage of the genome, identifying errors in homopolymer length attributable to 454 sequencing. Peptide sequences resulting from high-throughput tandem mass spectrometry of trypsin-digested protein extracts were mapped to 1,755 annotated proteins (68% of all protein-coding regions). Proteogenomic analysis confirmed the quality of annotation and improvement pipelines, identifying a missing gene and an alternative reading frame. Peptide coverage of genes hypothetically involved in substrate hydrolysis, transport and utilization confirmed multiple pathways for glycolysis, pyruvate conversion and recycling of intermediates. No sequences homologous to transaldolase, a central enzyme in the pentose phosphate pathway, were observed by any method, despite demonstrated growth of this organism on xylose and xylan hemicellulose.ConclusionsComplementary omics techniques confirm the quality of genome sequence assembly, annotation and error-reporting. Nearly complete genome coverage by RNA-seq likely indicates background DNA in RNA extracts, however these preps resulted in WGS enhancement and transcriptome profiling in a single Illumina run. No detection of transaldolase by any method despite xylose utilization by this organism indicates an alternative pathway for sedoheptulose-7-phosphate degradation. This report combines next-generation omics techniques to elucidate previously undefined features of substrate transport and central metabolism for this organism and its potential for consolidated biofuels production from lignocellulose.

Pentose sugars inhibit metabolism and increase expression of an AgrD-type cyclic pentapeptide in Clostridium thermocellum

Clostridium thermocellum could potentially be used as a microbial biocatalyst to produce renewable fuels directly from lignocellulosic biomass due to its ability to rapidly solubilize plant cell walls. While the organism readily ferments sugars derived from cellulose, pentose sugars from xylan are not metabolized. Here, we show that non-fermentable pentoses inhibit growth and end-product formation during fermentation of cellulose-derived sugars. Metabolomic experiments confirmed that xylose is transported intracellularly and reduced to the dead-end metabolite xylitol. Comparative RNA-seq analysis of xylose-inhibited cultures revealed several up-regulated genes potentially involved in pentose transport and metabolism, which were targeted for disruption. Deletion of the ATP-dependent transporter, CbpD partially alleviated xylose inhibition. A putative xylitol dehydrogenase, encoded by Clo1313_0076, was also deleted resulting in decreased total xylitol production and yield by 41% and 46%, respectively. Finally, xylose-induced inhibition corresponds with the up-regulation and biogenesis of a cyclical AgrD-type, pentapeptide. Medium supplementation with the mature cyclical pentapeptide also inhibits bacterial growth. Together, these findings provide new foundational insights needed for engineering improved pentose utilizing strains of C. thermocellum and reveal the first functional Agr-type cyclic peptide to be produced by a thermophilic member of the Firmicutes.

Chromatographic behavior of peptides containing oxidized methionine residues in proteomic LC-MS experiments: Complex tale of a simple modification.

On average, the oxidation of a single Met residue to Mso (methionine S-oxide, methionine sulfoxide) and Msn (methionine S,S-dioxide, methionine sulfone) decreases peptide retention in RP HPLC by 2.37 and 1.95 Hydrophobicity Index units (% acetonitrile), respectively. At the same time, the magnitude of the retention shift varies greatly (-9.1 to +0.4% acetonitrile for Mso) depending on peptide sequence. The latter effects are mostly associated with the stabilization of secondary structures upon peptide interaction with the hydrophobic stationary phase: when an oxidized residue is located in the hydrophobic face of an amphipathic helix, the decrease in retention is profound. The same amino acid positioning leads to complete or partial resolution of pairs of peptides containing diastereomeric Mso residues. Contrary to all previously reported observations, and the nature of this modification, we also demonstrate for the first time that methionine oxidation may increase peptide hydrophobicity. This behavior is characteristic for Met residues in the N3 position of an N-capping box stabilization motif prior to the amphipathic helix. All these findings indicate that the prediction of peptide secondary structures upon interaction with hydrophobic surfaces must become an integral part of peptide retention modeling in proteomic applications going forward.

Omics Approaches for Designing Biofuel Producing Cocultures for Enhanced Microbial Conversion of Lignocellulosic Substrates

Abstract Consolidated bioprocessing (CBP) is a system in which enzyme production, substrate hydrolysis, and fermentation are accomplished in a single process step by lignocellulolytic microorganisms. CBP offers the potential for lower biofuel production costs due to simpler feedstock processing, lower energy inputs, and higher conversion efficiencies than separate hydrolysis and fermentation processes, and it is an economically attractive near-term goal for “next-generation” cellulosic biofuel production. Although monocultures of cellulolytic bacteria have proven useful, they do have limitations, and cocultures of mesophilic or thermophilic bacteria that can directly ferment lignocellulosic biomass have been reported to display increased rates of cellulose hydrolysis and higher ethanol titers than observed in monocultures. Development of “designer consortia” for novel biorefining processes may be a highly productive area of research for H 2 and/or ethanol production coupled to synthesis of value-added co-products. A comprehensive understanding of lignocellulose fermentation may require, in future studies, the utilization of holistic approaches that involve genomic, proteomic, and transcriptomic methodologies. In this chapter, we discuss the application of “omics” (bioinformatics, genomics, proteomics, and transcriptomics) to identify synergistic cocultures of bacteria that display higher rates of substrate conversion and enhanced yields of desired end products.