Search NASA⌕ Search

SEARCH · Search NASA

Results for “genome binning”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

ATLAS (Automatic Tool for Local Assembly Structures) - A Comprehensive Infrastructure for Assembly, Annotation, and Genomic Binning of Metagenomic and Metaranscripomic Data

ATLAS (Automatic Tool for Local Assembly Structures) is a comprehensive multiomics data analysis pipeline that is massively parallel and scalable. ATLAS contains a modular analysis pipeline for assembly, annotation, quantification and genome binning of metagenomics and metatranscriptomics data and a framework for reference metaproteomic database construction. ATLAS transforms raw sequence data into functional and taxonomic data at the microbial population level and provides genome-centric resolution through genome binning. ATLAS provides robust taxonomy based on majority voting of protein coding open reading frames rolled-up at the contig level using modified lowest common ancestor (LCA) analysis. ATLAS provides robust taxonomy based on majority voting of protein coding open reading frames rolled-up at the contig level using modified lowest common ancestor (LCA) analysis. ATLAS is user-friendly, easy install through bioconda maintained as open-source on GitHub, and is implemented in Snakemake for modular customizable workflows.

genome binning↗
DOE OSTI2016

Moleculo Long-Read Sequencing Facilitates Assembly and Genomic Binning from Complex Soil Metagenomes

ABSTRACT Soil metagenomics has been touted as the “grand challenge” for metagenomics, as the high microbial diversity and spatial heterogeneity of soils make them unamenable to current assembly platforms. Here, we aimed to improve soil metagenomic sequence assembly by applying the Moleculo synthetic long-read sequencing technology. In total, we obtained 267 Gbp of raw sequence data from a native prairie soil; these data included 109.7 Gbp of short-read data (~100 bp) from the Joint Genome Institute (JGI), an additional 87.7 Gbp of rapid-mode read data (~250 bp), plus 69.6 Gbp (>1.5 kbp) from Moleculo sequencing. The Moleculo data alone yielded over 5,600 reads of >10 kbp in length, and over 95% of the unassembled reads mapped to contigs of >1.5 kbp. Hybrid assembly of all data resulted in more than 10,000 contigs over 10 kbp in length. We mapped three replicate metatranscriptomes derived from the same parent soil to the Moleculo subassembly and found that 95% of the predicted genes, based on their assignments to Enzyme Commission (EC) numbers, were expressed. The Moleculo subassembly also enabled binning of >100 microbial genome bins. We obtained via direct binning the first complete genome, that of “ Candidatus Pseudomonas sp. strain JKJ-1” from a native soil metagenome. By mapping metatranscriptome sequence reads back to the bins, we found that several bins corresponding to low-relative-abundance Acidobacteria were highly transcriptionally active, whereas bins corresponding to high-relative-abundance Verrucomicrobia were not. These results demonstrate that Moleculo sequencing provides a significant advance for resolving complex soil microbial communities. IMPORTANCE Soil microorganisms carry out key processes for life on our planet, including cycling of carbon and other nutrients and supporting growth of plants. However, there is poor molecular-level understanding of their functional roles in ecosystem stability and responses to environmental perturbations. This knowledge gap is largely due to the difficulty in culturing the majority of soil microbes. Thus, use of culture-independent approaches, such as metagenomics, promises the direct assessment of the functional potential of soil microbiomes. Soil is, however, a challenge for metagenomic assembly due to its high microbial diversity and variable evenness, resulting in low coverage and uneven sampling of microbial genomes. Despite increasingly large soil metagenome data volumes (>200 Gbp), the majority of the data do not assemble. Here, we used the cutting-edge approach of synthetic long-read sequencing technology (Moleculo) to assemble soil metagenome sequence data into long contigs and used the assemblies for binning of genomes. Author Video : An author video summary of this article is available.

De Novo assembly↗
DOE OSTI2020

ATLAS: a Snakemake workflow for assembly, annotation, and genomic binning of metagenome sequence data

Background: Metagenomics and metatranscriptomics studies provide valuable insight into the composition and function of microbial populations from diverse environments, however the data processing pipelines that rely on mapping reads to gene catalogs or genome databases for cultured strains yield results that underrepresent the genes and functional potential of uncultured microbes. Recent improvements in sequence assembly methods have eased the reliance on genome databases, thereby allowing the recovery of genomes from uncultured microbes. However, configuring these tools, linking them with advanced binning and annotation tools, and maintaining provenance of the processing continues to be challenging for researchers. Results: Here we present ATLAS, a software package for customizable data processing from raw sequence reads to functional and taxonomic annotations using state-of-the-art tools to assemble, annotate, quantify, and bin metagenome and metatranscriptome data. Genome-centric resolution and abundance estimates are provided for each sample in a dataset. ATLAS is written in Python and the workflow implemented in Snakemake; it operates in a Linux environment, and is compatible with Python 3.5+ and Anaconda 3+ versions. The source code for ATLAS is freely available, distributed under a BSD-3 license. Conclusions: ATLAS provides a user-friendly, modular and customizable Snakemake workflow for metagenome and metatranscriptome data processing; it is easily installable with conda and maintained as open-source on GitHub at https://github.com/metagenome-atlas/atlas.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2017

Metagenomic and Metatranscriptomic Analyses Reveal the Structure and Dynamics of a Dechlorinating Community Containing Dehalococcoides mccartyi and Corrinoid-Providing Microorganisms under Cobalamin-Limited Conditions

The aim of this paper is to obtain a systems-level understanding of the interactions between Dehalococcoides and corrinoid-supplying microorganisms by analyzing community structures and functional compositions, activities, and dynamics in trichloroethene (TCE)-dechlorinating enrichments. Metagenomes and metatranscriptomes of the dechlorinating enrichments with and without exogenous cobalamin were compared. Seven putative draft genomes were binned from the metagenomes. At an early stage (2 days), more transcripts of genes in the Veillonellaceae bin-genome were detected in the metatranscriptome of the enrichment without exogenous cobalamin than in the one with the addition of cobalamin. Among these genes, sporulation-related genes exhibited the highest differential expression when cobalamin was not added, suggesting a possible release route of corrinoids from corrinoid producers. Other differentially expressed genes include those involved in energy conservation and nutrient transport (including cobalt transport). The most highly expressed corrinoid de novo biosynthesis pathway was also assigned to the Veillonellaceae bin-genome. Targeted quantitative PCR (qPCR) analyses confirmed higher transcript abundances of those corrinoid biosynthesis genes in the enrichment without exogenous cobalamin than in the enrichment with cobalamin. Furthermore, the corrinoid salvaging and modification pathway of Dehalococcoides was upregulated in response to the cobalamin stress. Finally, this study provides important insights into the microbial interactions and roles played by members of dechlorinating communities under cobalamin-limited conditions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2016

Identification and genome reconstruction of abundant distinct taxa in microbiomes from one thermophilic and three mesophilic production-scale biogas plants

Biofuel production from conversion of biomass is indispensable in the portfolio of renewable energies. Complex microbial communities are involved in the anaerobic digestion process of plant material, agricultural residual products and food wastes. Analysis of the genetic potential and microbiology of communities degrading biomass to biofuels is considered to be the key to develop process optimisation strategies. Hence, due to the still incomplete taxonomic and functional characterisation of corresponding communities, new and unknown species are of special interest. Three mesophilic and one thermophilic production-scale biogas plants (BGPs) were taxonomically profiled using high-throughput 16S rRNA gene amplicon sequencing. All BGPs shared a core microbiome with the thermophilic BGP featuring the lowest diversity. However, the phyla Cloacimonetes and Spirochaetes were unique to BGPs 2 and 3, Fusobacteria were only found in BGP3 and members of the phylum Thermotogae were present only in the thermophilic BGP4. Taxonomic analyses revealed that these distinctive taxa mostly represent so far unknown species. The only exception is the dominant Thermotogae OTU featuring 16S rRNA gene sequence identity to Defluviitoga tunisiensis L3, a sequenced and characterised strain. To further investigate the genetic potential of the biogas communities, corresponding metagenomes were sequenced in a deepness of 347.5 Gbp in total. A combined assembly comprised 80.3 % of all reads and resulted in the prediction of 1.59 million genes on assembled contigs. Genome binning yielded genome bins comprising the prevalent distinctive phyla Cloacimonetes, Spirochaetes, Fusobacteria and Thermotogae. Comparative genome analyses between the most dominant Thermotogae bin and the very closely related Defluviitoga tunisiensis L3 genome originating from the same BGP revealed high genetic similarity. This finding confirmed applicability and reliability of the binning approach. The four highly covered genome bins of the other three distinct phyla showed low or very low genetic similarities to their closest phylogenetic relatives, and therefore indicated their novelty. In this study, the 16S rRNA gene sequencing approach and a combined metagenome assembly and binning approach were used for the first time on different production-scale biogas plants and revealed insights into the genetic potential and functional role of so far unknown species.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2017

Targeted in situ metatranscriptomics for selected taxa from mesophilic and thermophilic biogas plants

Summary Biogas production is performed anaerobically by complex microbial communities with key species driving the process. Hence, analyses of their in situ activities are crucial to understand the process. In a previous study, metagenome sequencing and subsequent genome binning for different production‐scale biogas plants ( BGP s) resulted in four genome bins of special interest, assigned to the phyla Thermotogae , Fusobacteria , Spirochaetes and Cloacimonetes , respectively, that were genetically analysed. In this study, metatranscriptome sequencing of the same BGP samples was conducted, enabling in situ transcriptional activity determination of these genome bins. For this, mapping of metatranscriptome reads on genome bin sequences was performed providing transcripts per million ( TPM ) values for each gene. This approach revealed an active sugar‐based metabolism of the Thermotogae and Spirochaetes bins and an active amino acid‐based metabolism of the Fusobacteria and Cloacimonetes bins. The data also hint at syntrophic associations of the four corresponding species with methanogenic Archaea .

Biotechnology & Applied Microbiology↗
DOE OSTI2015

Single-Cell-Genomics-Facilitated Read Binning of Candidate Phylum EM19 Genomes from Geothermal Spring Metagenomes

Here, the vast majority of microbial life remains uncatalogued due to the inability to cultivate these organisms in the laboratory. This “microbial dark matter” represents a substantial portion of the tree of life and of the populations that contribute to chemical cycling in many ecosystems. In this work, we leveraged an existing single-cell genomic data set representing the candidate bacterial phylum “ Calescamantes ” (EM19) to calibrate machine learning algorithms and define metagenomic bins directly from pyrosequencing reads derived from Great Boiling Spring in the U.S. Great Basin. Compared to other assembly-based methods, taxonomic binning with a read-based machine learning approach yielded final assemblies with the highest predicted genome completeness of any method tested. Read-first binning subsequently was used to extract Calescamantes bins from all metagenomes with abundant Calescamantes populations, including metagenomes from Octopus Spring and Bison Pool in Yellowstone National Park and Gongxiaoshe Spring in Yunnan Province, China. Metabolic reconstruction suggests that Calescamantes are heterotrophic, facultative anaerobes, which can utilize oxidized nitrogen sources as terminal electron acceptors for respiration in the absence of oxygen and use proteins as their primary carbon source. Despite their phylogenetic divergence, the geographically separate Calescamantes populations were highly similar in their predicted metabolic capabilities and core gene content, respiring O 2 , or oxidized nitrogen species for energy conservation in distant but chemically similar hot springs.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2018

Insights into microbial community structure and function from a shallow, simulated CO 2 -leakage aquifer demonstrate microbial selection and adaptation

Geological carbon storage is likely to be a part of a comprehensive strategy to minimize the atmospheric release of carbon dioxide (CO 2 ), raising concerns that injected CO 2 will leak into overlying freshwater aquifers. CO 2(aq) leakage may impact the dominant microbial community responsible for important ecosystem functions such as nutrient cycling, metal cycling and carbon conversion. Here, we examined the impact of an experimental in situ CO 2 leakage on a freshwater aquifer microbial community. High-throughput 16S rRNA gene sequencing demonstrated lower microbial diversity in freshwater wells with CO 2 concentrations above 1.15 g/L. Metagenomic sequencing and population genome binning were used to evaluate the metabolic potential of microbial populations across four CO 2 exposed samples and one control sample. Population genome binning resulted in the recovery and annotation of three metagenome assembled genomes (MAGs). Two of the MAGs, most closely related to Curvibacter and Sulfuricurvum, had the functional capacity for CO 2 utilization via carbon fixation coupled to sulfur and iron oxidation. The third draft genome was an Archaea, most closely related to Methanoregula, characterized by the metabolic potential for methanogenesis. Together, these findings show that CO 2 leakage in a freshwater aquifer poses a strong selection, driving both microbial community structure and metabolic function.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2018

“Candidatus Paraporphyromonas polyenzymogenes” encodes multi-modular cellulases linked to the type IX secretion system

Background In Nature, obligate herbivorous ruminants have a close symbiotic relationship with their gastrointestinal microbiome, which proficiently deconstructs plant biomass. Despite decades of research, lignocellulose degradation in the rumen has thus far been attributed to a limited number of culturable microorganisms. Here, we combine metaomics and enzymology to identify and describe a novel Bacteroidetes family (UMH11) composed entirely of uncultivated strains that are predominant in ruminants and only distantly related to previously characterized taxa. Results The first metabolic reconstruction of UMH11-affiliated genome bins, with a particular focus on the provisionally named UParaporphyromonas polyenzymogenes, illustrated their capacity to degrade various lignocellulosic substrates via comprehensive inventories of singular and multi-modular carbohydrate active enzymes (CAZymes). Closer examination revealed an absence of archetypical polysaccharide utilization loci found in human-gut microbiota. Instead, we identified many multi-modular CAZymes putatively secreted via the Bacteroidetes-specific Type 9 secretion system (T9SS). This included cellulases with two or more catalytic domains, which are modular arrangements that are unique to Bacteroidetes species studied to date. Core metabolic proteins from UP. polyenzymogenes were detected in metaproteomic data and were enriched in rumen-incubated plant biomass, indicating that active saccharification and fermentation of complex carbohydrates could be assigned to members of this novel family. Biochemical analysis of selected UP. polyenzymogenes CAZymes further iterated the cellulolytic activity of this hitherto uncultured bacterium towards linear polymers, such as amorphous and crystalline cellulose as well as mixed linkage β-glucans. Conclusion We propose that UP. olyenzymogenes genotypes and other UMH11 members actively degrade plant biomass in the rumen of cows, sheep, and most likely other ruminants, utilizing singular and multi-domain catalytic CAZymes secreted through the T9SS. The discovery of a prominent role of multi-modular cellulases in the Gramnegative Bacteroidetes, together with similar findings for Gram-positive cellulosomal bacteria (Ruminococcus flavefaciens) and anaerobic fungi (Orpinomyces sp.), suggests that complex enzymes are essential and have evolved within all major cellulolytic dominions inherent to the rumen.

Environmental Molecular Sciences Laboratory↗
DOE OSTI2019

Elucidating Syntrophic Butyrate-Degrading Populations in Anaerobic Digesters Using Stable-Isotope-Informed Genome-Resolved Metagenomics

Linking the genomic content of uncultivated microbes to their metabolic functions remains a critical challenge in microbial ecology. Resolving this challenge has implications for improving our management of key microbial interactions in biotechnologies such as anaerobic digestion, which relies on slow-growing syntrophic and methanogenic communities to produce renewable methane from organic waste. In this study, we combined DNA stable-isotope probing (SIP) with genome-centric metagenomics to recover the genomes of populations enriched in 13 C after growing on [ 13 C]butyrate. Differential abundance analysis of recovered genomic bins across the SIP metagenomes identified two metagenome-assembled genomes (MAGs) that were significantly enriched in heavy [ 13 C]DNA. Phylogenomic analysis assigned one MAG to the genus Syntrophomonas and the other MAG to the genus Methanothrix. Metabolic reconstruction of the annotated genomes showed that the Syntrophomonas genome encoded all the enzymes for beta-oxidizing butyrate, as well as several mechanisms for interspecies electron transfer via electron transfer flavoproteins, hydrogenases, and formate dehydrogenases. The Syntrophomonas genome shared low average nucleotide identity (<95%) with any cultured representative species, indicating that it is a novel species that plays a significant role in syntrophic butyrate degradation within anaerobic digesters. The Methanothrix genome contained the complete pathway for acetoclastic methanogenesis, indicating that it was enriched in 13 C from syntrophic acetate transfer. This study demonstrates the potential of stable-isotope-informed genome-resolved metagenomics to identify in situ interspecies metabolic cooperation within syntrophic consortia important to anaerobic waste treatment as well as global carbon cycling.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

MVP: a modular viromics pipeline to identify, filter, cluster, annotate, and bin viruses from metagenomes

While numerous computational frameworks and workflows are available for recovering prokaryote and eukaryote genomes from metagenome data, only a limited number of pipelines are designed specifically for viromics analysis. With many viromics tools developed in the last few years alone, it can be challenging for scientists with limited bioinformatics experience to easily recover, evaluate quality, annotate genes, dereplicate, assign taxonomy, and calculate relative abundance and coverage of viral genomes using state-of-the-art methods and standards. Here, we describe Modular Viromics Pipeline (MVP) v.1.0, a user-friendly pipeline written in Python and providing a simple framework to perform standard viromics analyses. MVP combines multiple tools to enable viral genome identification, characterization of genome quality, filtering, clustering, taxonomic and functional annotation, genome binning, and comprehensive summaries of results that can be used for downstream ecological analyses. Overall, MVP provides a standardized and reproducible pipeline for both extensive and robust characterization of viruses from large-scale sequencing data including metagenomes, metatranscriptomes, viromes, and isolate genomes. As a typical use case, we show how the entire MVP pipeline can be applied to a set of 20 metagenomes from wetland sediments using only 10 modules executed via command lines, leading to the identification of 11,656 viral contigs and 8,145 viral operational taxonomic units (vOTUs) displaying a clear beta-diversity pattern. Further, acting as a dynamic wrapper, MVP is designed to continuously incorporate updates and integrate new tools, ensuring its ongoing relevance in the rapidly evolving field of viromics. MVP is available at https://gitlab.com/ccoclet/mvp and as versioned packages in PyPi and Conda.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

GenomeFace v1.0

GenomeFace is meta-genome binning software. Metagenomic binning, the process of grouping DNA sequences into taxonomic units, is critical for understanding the functions, interactions, and evolutionary dynamics of microbial communities. We propose a deep learning approach to binning using two neural networks, one based on composition and another on environmental abundance, dynamically weighting the contribution of each based on characteristics of the input data. Trained on over 43,000 prokaryotic genomes, our network for composition-based binning is inspired by metric learning techniques used for facial recognition. Using a task-specific, multi-GPU accelerated algorithm to cluster the embeddings produced by our network, our binner leverages marker genes observed to be universally present in nearly all taxa to grade and select optimal clusters of sequences from a hierarchy of candidates. We evaluate our approach on four simulated datasets with known ground truth. Our linear time integration of marker genes recovers more near complete genomes than state of the art but computationally infeasible solutions using them, while being over an order of magnitude faster. Finally, we demonstrate the scalability and acuity of our approach by testing it on three of the largest metagenome assemblies ever performed. Compared to other binners, we produced 47%-183% more near complete genomes. From these datasets, we find over the genomes of over 3000 new candidate species which have never been previously cataloged, representing a potential 4% expansion of the known bacterial tree of life.

Lettich, Richard [Lawrence Berkeley National Labor↗
DOE OSTI2017

Critical Assessment of Metagenome Interpretation—a benchmark of metagenomics software

Methods for assembly, taxonomic profiling and binning are key to interpreting metagenome data, but a lack of consensus about benchmarking complicates performance assessment. The Critical Assessment of Metagenome Interpretation (CAMI) challenge has engaged the global developer community to benchmark their programs on highly complex and realistic data sets, generated from ~700 newly sequenced microorganisms and ~600 novel viruses and plasmids and representing common experimental setups. Assembly and genome binning programs performed well for species represented by individual genomes but were substantially affected by the presence of related strains. Taxonomic profiling and binning programs were proficient at high taxonomic ranks, with a notable performance decrease below family level. Parameter settings markedly affected performance, underscoring their importance for program reproducibility. The CAMI results highlight current challenges but also provide a roadmap for software selection to answer specific research questions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2016

Draft Genome Sequence of Methanohalophilus mahii Strain DAL1 Reconstructed from a Hydraulic Fracturing-Produced Water Metagenome

ABSTRACT We report here the 1,882,100-bp draft genome sequence of Methanohalophilus mahii strain DAL1, recovered from Marcellus Shale hydraulic fracturing-produced water using metagenomic contig binning. Genome annotation revealed several key methanogenesis genes and provides valuable information on archaeal activity associated with hydraulic fracturing-produced water environments.

Lipus, Daniel↗
DOE OSTI2020

Experimentally-validated correlation analysis reveals new anaerobic methane oxidation partnerships with consortium-level heterogeneity in diazotrophy

Abstract Archaeal anaerobic methanotrophs (“ANME”) and sulfate-reducing Deltaproteobacteria (“SRB”) form symbiotic multicellular consortia capable of anaerobic methane oxidation (AOM), and in so doing modulate methane flux from marine sediments. The specificity with which ANME associate with particular SRB partners in situ, however, is poorly understood. To characterize partnership specificity in ANME-SRB consortia, we applied the correlation inference technique SparCC to 310 16S rRNA amplicon libraries prepared from Costa Rica seep sediment samples, uncovering a strong positive correlation between ANME-2b and members of a clade of Deltaproteobacteria we termed SEEP-SRB1g. We confirmed this association by examining 16S rRNA diversity in individual ANME-SRB consortia sorted using flow cytometry and by imaging ANME-SRB consortia with fluorescence in situ hybridization (FISH) microscopy using newly-designed probes targeting the SEEP-SRB1g clade. Analysis of genome bins belonging to SEEP-SRB1g revealed the presence of a complete nifHDK operon required for diazotrophy, unusual in published genomes of ANME-associated SRB. Active expression of nifH in SEEP-SRB1g within ANME-2b—SEEP-SRB1g consortia was then demonstrated by microscopy using hybridization chain reaction (HCR-) FISH targeting nifH transcripts and diazotrophic activity was documented by FISH-nanoSIMS experiments. NanoSIMS analysis of ANME-2b—SEEP-SRB1g consortia incubated with a headspace containing CH4 and 15N2 revealed differences in cellular 15N-enrichment between the two partners that varied between individual consortia, with SEEP-SRB1g cells enriched in 15N relative to ANME-2b in one consortium and the opposite pattern observed in others, indicating both ANME-2b and SEEP-SRB1g are capable of nitrogen fixation, but with consortium-specific variation in whether the archaea or bacterial partner is the dominant diazotroph.

58 GEOSCIENCES↗
DOE OSTI2021

Competitive Exclusion and Metabolic Dependency among Microorganisms Structure the Cellulose Economy of an Agricultural Soil

Microorganisms that degrade cellulose utilize extracellular reactions that yield free by-products which can promote interactions with noncellulolytic organisms. We hypothesized that these interactions determine the ecological and physiological traits governing the fate of cellulosic carbon (C) in soil. We performed comparative genomics with genome bins from a shotgun metagenomic-stable isotope probing experiment to characterize the attributes of cellulolytic and noncellulolytic taxa accessing 13 C from cellulose. We hypothesized that cellulolytic taxa would exhibit competitive traits that limit access, while noncellulolytic taxa would display greater metabolic dependency, such as signatures of adaptive gene loss. We tested our hypotheses by evaluating genomic traits indicative of competitive exclusion or metabolic dependency, such as antibiotic production, growth rate, surface attachment, biomass degrading potential, and auxotrophy. The most 13 C-enriched taxa were cellulolytic Cellvibrio (Gammaproteobacteria) and Chaetomium (Ascomycota), which exhibited a strategy of self-sufficiency (prototrophy), rapid growth, and competitive exclusion via antibiotic production. Auxotrophy was more prevalent in cellulolytic Actinobacteria than in cellulolytic Proteobacteria, demonstrating differences in dependency among cellulose degraders. Noncellulolytic taxa that accessed 13 C from cellulose (Planctomycetales, Verrucomicrobia, and Vampirovibrionales) were also more dependent, as indicated by patterns of auxotrophy and 13 C labeling (i.e., partial labeling or labeling at later stages). Major 13 C-labeled cellulolytic microbes (e.g., Sorangium, Actinomycetales, Rhizobiales, and Caulobacteraceae) possessed adaptations for surface colonization (e.g., gliding motility, hyphae, attachment structures) signifying the importance of surface ecology in decomposing particulate organic matter. Our results demonstrated that access to cellulosic C was accompanied by ecological trade-offs characterized by differing degrees of metabolic dependency and competitive exclusion.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2019

Review, Evaluation, and Directions for Gene-Targeted Assembly for Ecological Analyses of Metagenomes

Shotgun metagenomics has greatly advanced our understanding of microbial communities over the last decade. Metagenomic analyses often include assembly and genome binning, computationally daunting tasks especially for big data from complex environments such as soil and sediments. In many studies, however, only a subset of genes and pathways involved in specific functions are of interest; thus, it is not necessary to attempt global assembly. In addition, methods that target genes can be computationally more efficient and produce more accurate assembly by leveraging rich databases, especially for those genes that are of broad interest such as those involved in biogeochemical cycles, biodegradation, and antibiotic resistance or used as phylogenetic markers. Here, we review six gene-targeted assemblers with unique algorithms for extracting and/or assembling targeted genes: Xander, MegaGTA, SAT-Assembler, HMM-GRASPx, GenSeed-HMM, and MEGAN. We tested these tools using two datasets with known genomes, a synthetic community of artificial reads derived from the genomes of 17 bacteria, shotgun sequence data from a mock community with 48 bacteria and 16 archaea genomes, and a large soil shotgun metagenomic dataset. We compared assemblies of a universal single copy gene (rplB) and two N cycle genes (nifH and nirK). We measured their computational efficiency, sensitivity, specificity, and chimera rate and found Xander and MegaGTA, which both use a probabilistic graph structure to model the genes, have the best overall performance with all three datasets, although MEGAN, a reference matching assembler, had better sensitivity with synthetic and mock community members chosen from its reference collection. Also, Xander and MegaGTA are the only tools that include post-assembly scripts tuned for common molecular ecology and diversity analyses. Additionally, we provide a mathematical model for estimating the probability of assembling targeted genes in a metagenome for estimating required sequencing depth.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Metagenome-assembled genomes measured at 3 depths during snowmelt period in East River, CO (March, May, and June, September 2017)

Snowmelt is a critical biogeochemical period that accounts for large nitrogen (N) export events from high-elevation watersheds. Soil microbial populations bloom and immobilize N during snowmelt, yet the population size crashes in spring, which releases a pulse of soil N. We sought to discover the N sources fueling this microbial bloom and determine the fate of N following microbial die-off. Here, focusing on the snowmelt period within a headwater catchment of the Upper Colorado River Basin (East River, CO), we deployed strain-resolved metagenomics to identify the metabolic pathways and processes that mobilize soil N during and after snowmelt. Soil metagenome samples were taken from 6 snowpits from 3 depths (0-5cm, 5-15cm, >15cm) at 4 time points during snowmelt period (March 2017, May 2017, and June 2017, September 2017) generating 48 metagenomes. We reconstructed 474 metagenome-assembled genomes (MAGs) across all metagenomes.All 48 metagenomes were sequenced at JGI and raw data can be found under JGI (Joint Genome Institute) GOLD Study Gs0135149. Metagenome assemblies from IMG under the same study were used for genome binning. This dataset (1) a zip file of 474 MAGs (as fasta files, Gs0135149_bins_tar.gz), (2) sample metadata file with sample IGSNs (International Generic Sample Numbers) (samples.csv), (3) bounding box coordinates for the sampled locations (Gs0135149.kml), (4) metagenome metadata file listing IMG/M (Integrated Microbial Genomes/Metagenomes) metagenome accessions linking samples to metagenomes (metagenomes.csv), (5) location metadata file (locations.csv), (6) file-level metadata file (flmd.csv) and (7) data dictionary (dd.csv) file.This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

54 ENVIRONMENTAL SCIENCES↗