Metagenome-assembled genome extraction and analysis from microbiomes using KBase.

Uncultivated Bacteria and Archaea account for the vast majority of species on Earth, but obtaining their genomes directly from the environment, using shotgun sequencing, has only become possible recently. To realize the hope of capturing Earth's microbial genetic complement and to facilitate the investigation of the functional roles of specific lineages in a given ecosystem, technologies that accelerate the recovery of high-quality genomes are necessary. We present a series of analysis steps and data products for the extraction of high-quality metagenome-assembled genomes (MAGs) from microbiomes using the U.S. Department of Energy Systems Biology Knowledgebase (KBase) platform ( http://www.kbase.us/ ). Overall, these steps take about a day to obtain extracted genomes when starting from smaller environmental shotgun read libraries, or up to about a week from larger libraries. In KBase, the process is end-to-end, allowing a user to go from the initial sequencing reads all the way through to MAGs, which can then be analyzed with other KBase capabilities such as phylogenetic placement, functional assignment, metabolic modeling, pangenome functional profiling, RNA-Seq and others. While portions of such capabilities are available individually from other resources, the combination of the intuitive usability, data interoperability and integration of tools in a freely available computational resource makes KBase a powerful platform for obtaining MAGs from microbiomes. While this workflow offers tools for each of the key steps in the genome extraction process, it also provides a scaffold that can be easily extended with additional MAG recovery and analysis tools, via the KBase software development kit (SDK).

A genomic catalog of Earth's microbiomes.

The reconstruction of bacterial and archaeal genomes from shotgun metagenomes has enabled insights into the ecology and evolution of environmental and host-associated microbiomes. Here we applied this approach to >10,000 metagenomes collected from diverse habitats covering all of Earth's continents and oceans, including metagenomes from human and animal hosts, engineered environments, and natural and agricultural soils, to capture extant microbial, metabolic and functional potential. This comprehensive catalog includes 52,515 metagenome-assembled genomes representing 12,556 novel candidate species-level operational taxonomic units spanning 135 phyla. The catalog expands the known phylogenetic diversity of bacteria and archaea by 44% and is broadly available for streamlined comparative analyses, interactive exploration, metabolic modeling and bulk download. We demonstrate the utility of this collection for understanding secondary-metabolite biosynthetic potential and for resolving thousands of new host linkages to uncultivated viruses. This resource underscores the value of genome-centric approaches for revealing genomic properties of uncultivated microorganisms that affect ecosystem processes.

Single cell genomics indicates horizontal gene transfer and viral infections in a deep subsurface Firmicutes population.

A major fraction of Earth's prokaryotic biomass dwells in the deep subsurface, where cellular abundances per volume of sample are lower, metabolism is slower, and generation times are longer than those in surface terrestrial and marine environments. How these conditions impact biotic interactions and evolutionary processes is largely unknown. Here we employed single cell genomics to analyze cell-to-cell genome content variability and signatures of horizontal gene transfer (HGT) and viral infections in five cells of Candidatus Desulforudis audaxviator, which were collected from a 3 km-deep fracture water in the 2.9 Ga-old Witwatersrand Basin of South Africa. Between 0 and 32% of genes recovered from single cells were not present in the original, metagenomic assembly of Desulforudis, which was obtained from a neighboring subsurface fracture. We found a transposable prophage, a retron, multiple clustered regularly interspaced short palindromic repeats (CRISPRs) and restriction-modification systems, and an unusually high frequency of transposases in the analyzed single cell genomes. This indicates that recombination, HGT and viral infections are prevalent evolutionary events in the studied population of microorganisms inhabiting a highly stable deep subsurface environment.

High-Quality Draft Genome Sequence of Desulfovibrio carbinoliphilus FW-101-2B, an Organic Acid-Oxidizing Sulfate-Reducing Bacterium Isolated from Uranium(VI)-Contaminated Groundwater.

Desulfovibrio carbinoliphilus subsp. oakridgensis FW-101-2B is an anaerobic, organic acid/alcohol-oxidizing, sulfate-reducing δ-proteobacterium. FW-101-2B was isolated from contaminated groundwater at The Field Research Center at Oak Ridge National Lab after in situ stimulation for heavy metal-reducing conditions. The genome will help elucidate the metabolic potential of sulfate-reducing bacteria during uranium reduction.

Phylogenomically guided identification of industrially relevant GH1 ß-glucosidases through DNA synthesis and nanostructure-initiator mass spectrometry.

Harnessing the biotechnological potential of the large number of proteins available in sequence databases requires scalable methods for functional characterization. Here we propose a workflow to address this challenge by combining phylogenomic guided DNA synthesis with high-throughput mass spectrometry and apply it to the systematic characterization of GH1 ß-glucosidases, a family of enzymes necessary for biomass hydrolysis, an important step in the conversion of lignocellulosic feedstocks to fuels and chemicals. We synthesized and expressed 175 GH1s, selected from over 2000 candidate sequences to cover maximum sequence diversity. These enzymes were functionally characterized over a range of temperatures and pHs using nanostructure-initiator mass spectrometry (NIMS), generating over 10,000 data points. When combined with HPLC-based sugar profiling, we observed GH1 enzymes active over a broad temperature range and toward many different ß-linked disaccharides. For some GH1s we also observed activity toward laminarin, a more complex oligosaccharide present as a major component of macroalgae. An area of particular interest was the identification of GH1 enzymes compatible with the ionic liquid 1-ethyl-3-methylimidazolium acetate ([C2mim][OAc]), a next-generation biomass pretreatment technology. We thus searched for GH1 enzymes active at 70 °C and 20% (v/v) [C2mim][OAc] over the course of a 24-h saccharification reaction. Using our unbiased approach, we identified multiple enzymes of different phylogentic origin with such activities. Our approach of characterizing sequence diversity through targeted gene synthesis coupled to high-throughput screening technologies is a broadly applicable paradigm for a wide range of biological problems.

Metagenomes of tropical soil-derived anaerobic switchgrass-adapted consortia with and without iron.

Tropical forest soils decompose litter rapidly with frequent episodes of anoxia, making it likely that bacteria using alternate terminal electron acceptors (TEAs) such as iron play a large role in supporting decomposition under these conditions. The prevalence of many types of metabolism in litter deconstruction makes these soils useful templates for improving biofuel production. To investigate how iron availability affects decomposition, we cultivated feedstock-adapted consortia (FACs) derived from iron-rich tropical forest soils accustomed to experiencing frequent episodes of anaerobic conditions and frequently fluctuating redox. One consortium was propagated under fermenting conditions, with switchgrass as the sole carbon source in minimal media (SG only FACs), and the other consortium was treated the same way but received poorly crystalline iron as an additional terminal electron acceptor (SG + Fe FACs). We sequenced the metagenomes of both consortia to a depth of about 150 Mb each, resulting in a coverage of 26× for the more diverse SG + Fe FACs, and 81× for the relatively less diverse SG only FACs. Both consortia were able to quickly grow on switchgrass, and the iron-amended consortium exhibited significantly higher microbial diversity than the unamended consortium. We found evidence of higher stress in the unamended FACs and increased sugar transport and utilization in the iron-amended FACs. This work provides metagenomic evidence that supplementation of alternative TEAs may improve feedstock deconstruction in biofuel production.

MetaMicrobesOnline: phylogenomic analysis of microbial communities.

The metaMicrobesOnline database (freely available at http://meta.MicrobesOnline.org) offers phylogenetic analysis of genes from microbial genomes and metagenomes. Gene trees are constructed for canonical gene families such as COG and Pfam. Such gene trees allow for rapid homologue analysis and subfamily comparison of genes from multiple metagenomes and comparisons with genes from microbial isolates. Additionally, the genome browser permits genome context comparisons, which may be used to determine the closest sequenced genome or suggest functionally associated genes. Lastly, the domain browser permits rapid comparison of protein domain organization within genes of interest from metagenomes and complete microbial genomes.

Manipulation of the carbon storage regulator system for metabolite remodeling and biofuel production in Escherichia coli.

BACKGROUND: Microbial engineering strategies that elicit global metabolic perturbations have the capacity to increase organism robustness for targeted metabolite production. In particular, perturbations to regulators of cellular systems that impact glycolysis and amino acid production while simultaneously decreasing fermentation by-products such as acetate and CO(2) make ideal targets. Intriguingly, perturbation of the Carbon Storage Regulator (Csr) system has been previously implicated in large changes in central carbon metabolism in E. coli. Therefore, we hypothesized that perturbation of the Csr system through the CsrA-CsrB ribonucleoprotein complex might increase production of biofuels and their intermediates from heterologous pathways. RESULTS: We engaged the CsrA-CsrB ribonucleoprotein complex of E. coli via overexpression of CsrB. CsrB is a 350-nucleotide non-coding RNA that antagonizes CsrA, an RNA-binding protein that regulates translation of specific mRNA targets. By using shotgun proteomics and targeted metabolomics we established that elevation of CsrB levels leads to alterations in metabolite and protein levels in glycolysis, the TCA cycle and amino acid levels. Consequently, we show that such changes can be suitably applied to improve the production of desired compounds through the native fatty acid and heterologous n-butanol and isoprenoid pathways by up to two-fold. We also observed concomitant decreases in undesirable fermentation by-products such as acetate and CO(2). CONCLUSIONS: We have demonstrated that simple engineering of the RNA-based Csr global regulatory system constitutes a novel approach to obtaining pathway-independent improvements within engineered hosts. Additionally, since Csr is conserved across most prokaryotic species, this approach may also be amenable to a wide variety of production hosts.

Tracing determinants of dual substrate specificity in glycoside hydrolase family 5.

Enzymes are traditionally viewed as having exquisite substrate specificity; however, recent evidence supports the notion that many enzymes have evolved activities against a range of substrates. The diversity of activities across glycoside hydrolase family 5 (GH5) suggests that this family of enzymes may contain numerous members with activities on multiple substrates. In this study, we combined structure- and sequence-based phylogenetic analysis with biochemical characterization to survey the prevalence of dual specificity for glucan- and mannan-based substrates in the GH5 family. Examination of amino acid profile differences between the subfamilies led to the identification and subsequent experimental confirmation of an active site motif indicative of dual specificity. The motif enabled us to successfully discover several new dually specific members of GH5, and this pattern is present in over 70 other enzymes, strongly suggesting that dual endoglucanase-mannanase activity is widespread in this family. In addition, reinstatement of the conserved motif in a wild type member of GH5 enhanced its catalytic efficiency on glucan and mannan substrates by 175 and 1,600%, respectively. Phylogenetic examination of other GH families further indicates that the prevalence of enzyme multispecificity in GHs may be greater than has been experimentally characterized. Single domain multispecific GHs may be exploited for developing improved enzyme cocktails or facile engineering of microbial hosts for consolidated bioprocessing of lignocellulose.

Complete genome sequence of "Enterobacter lignolyticus" SCF1.

In an effort to discover anaerobic bacteria capable of lignin degradation, we isolated "Enterobacter lignolyticus" SCF1 on minimal media with alkali lignin as the sole source of carbon. This organism was isolated anaerobically from tropical forest soils collected from the Short Cloud Forest site in the El Yunque National Forest in Puerto Rico, USA, part of the Luquillo Long-Term Ecological Research Station. At this site, the soils experience strong fluctuations in redox potential and are net methane producers. Because of its ability to grow on lignin anaerobically, we sequenced the genome. The genome of "E. lignolyticus" SCF1 is 4.81 Mbp with no detected plasmids, and includes a relatively small arsenal of lignocellulolytic carbohydrate active enzymes. Lignin degradation was observed in culture, and the genome revealed two putative laccases, a putative peroxidase, and a complete 4-hydroxyphenylacetate degradation pathway encoded in a single gene cluster.

Optimization of a heterologous mevalonate pathway through the use of variant HMG-CoA reductases.

Expression of foreign pathways often results in suboptimal performance due to unintended factors such as introduction of toxic metabolites, cofactor imbalances or poor expression of pathway components. In this study we report a 120% improvement in the production of the isoprenoid-derived sesquiterpene, amorphadiene, produced by an engineered strain of Escherichia coli developed to express the native seven-gene mevalonate pathway from Saccharomyces cerevisiae (Martin et al. 2003). This substantial improvement was made by varying only a single component of the pathway (HMG-CoA reductase) and subsequent host optimization to improve cofactor availability. We characterized and tested five variant HMG-CoA reductases obtained from publicly available genome databases with differing kinetic properties and cofactor requirements. The results of our in vitro and in vivo analyses of these enzymes implicate substrate inhibition of mevalonate kinase as an important factor in optimization of the engineered mevalonate pathway. Consequently, the NADH-dependent HMG-CoA reductase from Delftia acidovorans, which appeared to have the optimal kinetic parameters to balance HMG-CoA levels below the cellular toxicity threshold of E. coli and those of mevalonate below inhibitory concentrations for mevalonate kinase, was identified as the best producer for amorphadiene (54% improvement over the native pathway enzyme, resulting in 2.5mM or 520 mg/L of amorphadiene after 48 h). We further enhanced performance of the strain bearing the D. acidovorans HMG-CoA reductase by increasing the intracellular levels of its preferred cofactor (NADH) using a NAD(+)-dependent formate dehydrogenase from Candida boidinii, along with formate supplementation. This resulted in an overall improvement of the system by 120% resulting in 3.5mM or 700 mg/L amorphadiene after 48 h of fermentation. This comprehensive study incorporated analysis of several key parameters for metabolic design such as in vitro and in vivo kinetic performance of variant enzymes, intracellular levels of protein expression, in-pathway substrate inhibition and cofactor management to enable the observed improvements. These metrics may be applied to a broad range of heterologous pathways for improving the production of biologically derived compounds.

GLAMM: Genome-Linked Application for Metabolic Maps.

The Genome-Linked Application for Metabolic Maps (GLAMM) is a unified web interface for visualizing metabolic networks, reconstructing metabolic networks from annotated genome data, visualizing experimental data in the context of metabolic networks and investigating the construction of novel, transgenic pathways. This simple, user-friendly interface is tightly integrated with the comparative genomics tools of MicrobesOnline [Dehal et al. (2010) Nucleic Acids Research, 38, D396-D400]. GLAMM is available for free to the scientific community at glamm.lbl.gov.

MicrobesOnline: an integrated portal for comparative and functional genomics.

Since 2003, MicrobesOnline (http://www.microbesonline.org) has been providing a community resource for comparative and functional genome analysis. The portal includes over 1000 complete genomes of bacteria, archaea and fungi and thousands of expression microarrays from diverse organisms ranging from model organisms such as Escherichia coli and Saccharomyces cerevisiae to environmental microbes such as Desulfovibrio vulgaris and Shewanella oneidensis. To assist in annotating genes and in reconstructing their evolutionary history, MicrobesOnline includes a comparative genome browser based on phylogenetic trees for every gene family as well as a species tree. To identify co-regulated genes, MicrobesOnline can search for genes based on their expression profile, and provides tools for identifying regulatory motifs and seeing if they are conserved. MicrobesOnline also includes fast phylogenetic profile searches, comparative views of metabolic pathways, operon predictions, a workbench for sequence analysis and integration with RegTransBase and other microbial genome resources. The next update of MicrobesOnline will contain significant new functionality, including comparative analysis of metagenomic sequence data. Programmatic access to the database, along with source code and documentation, is available at http://microbesonline.org/programmers.html.

Contribution of mobile genetic elements to Desulfovibrio vulgaris genome plasticity.

The genome of Desulfovibrio vulgaris strain DePue, a sulfate-reducing Deltaproteobacterium isolated from heavy metal-impacted lake sediment, was completely sequenced and compared with the type strain D. vulgaris Hildenborough. The two genomes share a high degree of relatedness and synteny, but harbour distinct prophage and signatures of past phage encounters. In addition to a highly variable phage contribution, the genome of strain DePue contains a cluster of open-reading frames not found in strain Hildenborough coding for the production and export of a capsule exopolysaccharide, possibly of relevance to heavy metal resistance. Comparative whole-genome microarray analysis on four additional D. vulgaris strains established greater interstrain variation within regions associated with phage insertion and exopolysaccharide biosynthesis.

Environmental genomics reveals a single-species ecosystem deep within Earth.

DNA from low-biodiversity fracture water collected at 2.8-kilometer depth in a South African gold mine was sequenced and assembled into a single, complete genome. This bacterium, Candidatus Desulforudis audaxviator, composes >99.9% of the microorganisms inhabiting the fluid phase of this particular fracture. Its genome indicates a motile, sporulating, sulfate-reducing, chemoautotrophic thermophile that can fix its own nitrogen and carbon by using machinery shared with archaea. Candidatus Desulforudis audaxviator is capable of an independent life-style well suited to long-term isolation from the photosphere deep within Earth's crust and offers an example of a natural ecosystem that appears to have its biological component entirely encoded within a single genome.

Structure prediction for CASP7 targets using extensive all-atom refinement with Rosetta@home.

We describe predictions made using the Rosetta structure prediction methodology for both template-based modeling and free modeling categories in the Seventh Critical Assessment of Techniques for Protein Structure Prediction. For the first time, aggressive sampling and all-atom refinement could be carried out for the majority of targets, an advance enabled by the Rosetta@home distributed computing network. Template-based modeling predictions using an iterative refinement algorithm improved over the best existing templates for the majority of proteins with less than 200 residues. Free modeling methods gave near-atomic accuracy predictions for several targets under 100 residues from all secondary structure classes. These results indicate that refinement with an all-atom energy function, although computationally expensive, is a powerful method for obtaining accurate structure predictions.