SEARCH · Search NASA

Results for “ecological genomics”

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.

SourceResearch connectionFrom yearThrough year

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

At least 19 records

All sources NASA NTRS DOE OSTI

NASA NTRS2016

Permanent Draft Genome of Strain ESFC-1: Ecological Genomics of a Newly Discovered Lineage of Filamentous Diazotrophic Cyanobacteria

The nonheterocystous filamentous cyanobacterium, strain ESFC-1, is a recently described member of the order Oscillatoriales within the Cyanobacteria. ESFC-1 has been shown to be a major diazotroph in the intertidal microbial mat system at Elkhorn Slough, CA, USA. Based on phylogenetic analyses of the 16S RNA gene, ESFC-1 appears to belong to a unique, genus-level divergence; the draft genome sequence of this strain has now been determined. Here we report features of this genome as they relate to the ecological functions and capabilities of strain ESFC-1. The 5,632,035 bp genome sequence encodes 4914 protein-coding genes and 92 RNA genes. One striking feature of this cyanobacterium is the apparent lack of either uptake or bi-directional hydrogenases typically expected within a diazotroph. Additionally, a large genomic island is found that contains numerous low GC-content genes and genes related to extracellular polysaccharide production and cell wall synthesis and maintenance.

R Craig Everroad↗

DOE OSTI2024

The mycorrhizal symbiosis: research frontiers in genomics, ecology, and agricultural application

Summary Mycorrhizal symbioses between plants and fungi are vital for the soil structure, nutrient cycling, plant diversity, and ecosystem sustainability. More than 250 000 plant species are associated with mycorrhizal fungi. Recent advances in genomics and related approaches have revolutionized our understanding of the biology and ecology of mycorrhizal associations. The genomes of 250+ mycorrhizal fungi have been released and hundreds of genes that play pivotal roles in regulating symbiosis development and metabolism have been characterized. rDNA metabarcoding and metatranscriptomics provide novel insights into the ecological cues driving mycorrhizal communities and functions expressed by these associations, linking genes to ecological traits such as nutrient acquisition and soil organic matter decomposition. Here, we review genomic studies that have revealed genes involved in nutrient uptake and symbiosis development, and discuss adaptations that are fundamental to the evolution of mycorrhizal lifestyles. We also evaluated the ecosystem services provided by mycorrhizal networks and discuss how mycorrhizal symbioses hold promise for sustainable agriculture and forestry by enhancing nutrient acquisition and stress tolerance. Overall, unraveling the intricate dynamics of mycorrhizal symbioses is paramount for promoting ecological sustainability and addressing current pressing environmental concerns. This review ends with major frontiers for further research.

Plant Sciences↗

DOE OSTI2021

Plant ecological genomics at the limits of life in the Atacama Desert

The Atacama Desert in Chile—hyperarid and with high–ultraviolet irradiance levels—is one of the harshest environments on Earth. Yet, dozens of species grow there, including Atacama-endemic plants. Herein, we establish the Talabre–Lejía transect (TLT) in the Atacama as an unparalleled natural laboratory to study plant adaptation to extreme environmental conditions. We characterized climate, soil, plant, and soil–microbe diversity at 22 sites (every 100 m of altitude) along the TLT over a 10-y period. We quantified drought, nutrient deficiencies, large diurnal temperature oscillations, and pH gradients that define three distinct vegetational belts along the altitudinal cline. Further, we deep-sequenced transcriptomes of 32 dominant plant species spanning the major plant clades, and assessed soil microbes by metabarcoding sequencing. The top-expressed genes in the 32 Atacama species are enriched in stress responses, metabolism, and energy production. Moreover, their root-associated soils are enriched in growth-promoting bacteria, including nitrogen fixers. To identify genes associated with plant adaptation to harsh environments, we compared 32 Atacama species with the 32 closest sequenced species, comprising 70 taxa and 1,686,950 proteins. To perform phylogenomic reconstruction, we concatenated 15,972 ortholog groups into a supermatrix of 8,599,764 amino acids. Using two codon-based methods, we identified 265 candidate positively selected genes (PSGs) in the Atacama plants, 64% of which are located in Pfam domains, supporting their functional relevance. For 59/184 PSGs with an Arabidopsis ortholog, we uncovered functional evidence linking them to plant resilience. As some Atacama plants are closely related to staple crops, these candidate PSGs are a “genetic goldmine” to engineer crop resilience to face climate change.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Understanding local plant extinctions before it is too late: bridging evolutionary genomics with global ecology

Understanding evolutionary genomic and population processes within a species range is key to anticipating the extinction of plant species before it is too late. However, most models of biodiversity risk under global change do not account for the genetic variation and local adaptation of different populations. Population diversity is critical to understanding extinction because different populations may be more or less susceptible to global change and, if lost, would reduce the total diversity within a species. Here, two new modeling frameworks advance our understanding of extinction from a population and evolutionary angle: Rapid climate change-driven disruptions in population adaptation are predicted from associations between genomes and local climates. Furthermore, losses of population diversity from global land-use transformations are estimated by scaling relationships of species' genomic diversity with habitat area. Overall, these global eco-evolutionary methods advance the predictability – and possibly the preventability – of the ongoing extinction of plant species.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2022

Phylogenomic structure and speciation in an emerging model: the Sphagnum magellanicum complex (Bryophyta)

Sphagnum magellanicum is one of two Sphagnum species for which a reference-quality genome exists to facilitate research in ecological genomics. Phylogenetic and comparative genomic analyses were conducted based on resequencing data from 48 samples and RADseq analyses based on 187 samples. Here, we report that there are four clades/species within the S. magellanicum complex in eastern North America and that the reference genome belongs to Sphagnum divinum . The species exhibit tens of thousands (RADseq) to millions (resequencing) of fixed nucleotide differences. Two species, however, referred to informally as S. diabolicum and S. magni because they have not been formally described, are differentiated by only 100 (RADseq) to 1000 (resequencing) of differences. Introgression among species in the complex is demonstrated using D-statistics and f 4 ratios. One ecologically important functional trait, tissue decomposability, which underlies peat (carbon) accumulation, does not differ between segregates in the S. magellanicum complex, although previous research showed that many closely related Sphagnum species have evolved differences in decomposability/carbon sequestration. Phylogenetic resolution and more accurate species delimitation in the S. magellanicum complex substantially increase the value of this group for studying the early evolutionary stages of climate adaptation and ecological evolution more broadly.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2020

Ecological and genomic analyses of candidate phylum WPS-2 bacteria in an unvegetated soil

Members of the bacterial candidate phylum WPS-2 (or Eremiobacterota) are abundant in several dry, bare soil environments. In a bare soil deposited by an extinct iron-sulfur spring, we found that WPS-2 comprised up to 24% of the bacterial community and up to 108 cells per g of soil based on 16S rRNA gene sequencing and quantification. A single genus-level cluster (Ca. Rubrimentiphilum) predominated in bare soils but was less abundant in adjacent forest. Nearly complete genomes of Ca. Rubrimentiphilum were recovered as single amplified genomes (SAGs) and metagenome-assembled genomes (MAGs). Surprisingly, given the abundance of WPS-2 in bare soils, the genomes did not indicate any capacity for autotrophy, phototrophy, or trace gas metabolism. Instead, they suggest a predominantly aerobic organoheterotrophic lifestyle, perhaps based on scavenging amino acids, nucleotides, and complex oligopeptides, along with lithotrophic capacity on thiosulfate. Network analyses of the entire community showed that some species of Chloroflexi, Actinobacteria, and candidate phylum AD3 (or Dormibacterota) co-occurred with Ca. Rubrimentiphilum and may represent ecological or metabolic partners. We propose that Ca. Rubrimentiphilum act as efficient heterotrophic scavengers. Combined with previous studies, these data suggest that the phylum WPS-2 includes bacteria with diverse metabolic capabilities.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2022

Ecological and genomic responses of soil microbiomes to high-severity wildfire: linking community assembly to functional potential

Increasing wildfire severity, which is common throughout the western United States, can have deleterious effects on plant regeneration and large impacts on carbon (C) and nitrogen (N) cycling rates. Soil microbes are pivotal in facilitating these elemental cycles, so understanding the impact of increasing fire severity on soil microbial communities is critical. Here, we assess the long-term impact of high-severity fires on the soil microbiome. We find that high-severity wildfires result in a multi-decadal (>25 y) recovery of the soil microbiome mediated by concomitant differences in aboveground vegetation, soil chemistry, and microbial assembly processes. Our results depict a distinct taxonomic and functional successional pattern of increasing selection in post-fire soil microbial communities. Changes in microbiome composition corresponded with changes in microbial functional potential, specifically altered C metabolism and enhanced N cycling potential, which related to rates of potential decomposition and inorganic N availability, respectively. Based on metagenome-assembled genomes, we show that bacterial genomes enriched in our earliest site (4 y since fire) harbor distinct traits such as a robust stress response and a high potential to degrade pyrogenic, polyaromatic C that allow them to thrive in post-fire environments. Taken together, these results provide a biological basis for previously reported process rate measurements and explain the temporal dynamics of post-fire biogeochemistry, which ultimately constrains ecosystem recovery.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2026

Ecological and genomic variation in ectomycorrhizal fungal exploration types

Ectomycorrhizal fungi (EMF) produce mycelia with variable extension and complexity, which can be classified according to soil ‘exploration types’ (ETs). ETs have received attention as one of the few mycorrhizal trait frameworks, but without an empirical classification of ET functional diversity and environmental preferences, understanding and interpreting EMF biogeographic patterns has been difficult. We conducted a synthesis combining: comparative EMF genomics to describe functional divergence in decomposition and nutrient cycling genes across ETs; and EMF trait distribution modeling across continental Europe, pairing soil and root EMF surveys to establish biogeographic ET niche profiles. We demonstrate a signature of ETs encoded in EMF genomes, which is independent from phylogeny and linked to biomass production strategies. EMF ET relative abundances were separated by soil, root, and dominant tree leaf type habitats and exhibited unique correlations with forest biotic (e.g. plant productivity and plant pathogen densities) and abiotic (e.g. nitrogen deposition and soil pH) conditions. These findings support a theory that EMF niche partitioning can be partially explained by extraradical mycelial traits, with underlying variation in ET biogeography likely arising from distinct decomposition and nutrient cycling potentials. We also identify important limitations to this trait framework and provide a guided outlook for future research.

biogeography↗

DOE OSTI2023

Methyl-Based Methanogenesis: an Ecological and Genomic Review

Methyl-based methanogenesis is one of three broad categories of archaeal anaerobic methanogenesis, including both the methyl dismutation (methylotrophic) pathway and the methyl-reducing (also known as hydrogen-dependent methylotrophic) pathway. Methyl-based methanogenesis is increasingly recognized as an important source of methane in a variety of environments. Here, we provide an overview of methyl-based methanogenesis research, including the conditions under which methyl-based methanogenesis can be a dominant source of methane emissions, experimental methods for distinguishing different pathways of methane production, molecular details of the biochemical pathways involved, and the genes and organisms involved in these processes. In this study, we also identify the current gaps in knowledge and present a genomic and metagenomic survey of methyl-based methanogenesis genes, highlighting the diversity of methyl-based methanogens at multiple taxonomic levels and the widespread distribution of known methyl-based methanogenesis genes and families across different environments.

methanogens↗

DOE OSTI2021

Metapangenomics reveals depth-dependent shifts in metabolic potential for the ubiquitous marine bacterial SAR324 lineage

Abstract Background Oceanic microbiomes play a pivotal role in the global carbon cycle and are central to the transformation and recycling of carbon and energy in the ocean’s interior. SAR324 is a ubiquitous but poorly understood uncultivated clade of Deltaproteobacteria that inhabits the entire water column, from ocean surface waters to its deep interior. Although some progress has been made in elucidating potential metabolic traits of SAR324 in the dark ocean, very little is known about the ecology and the metabolic capabilities of this group in the euphotic and twilight zones. To investigate the comparative genomics, ecology, and physiological potential of the SAR324 clade, we examined the distribution and variability of key genomic features and metabolic pathways in this group from surface waters to the abyss in the North Pacific Subtropical Gyre, one of the largest biomes on Earth. Results We leveraged a pangenomic ecological approach, combining spatio-temporally resolved single-amplified genome, metagenomic, and metatranscriptomic datasets. The data revealed substantial genomic diversity throughout the SAR324 clade, with distinct depth and temporal distributions that clearly differentiated ecotypes. Phylogenomic subclade delineation, environmental distributions, genomic feature similarities, and metabolic capacities revealed strong congruence. The four SAR324 ecotypes delineated in this study revealed striking divergence from one another with respect to their habitat-specific metabolic potentials. The ecotypes living in the dark or twilight oceans shared genomic features and metabolic capabilities consistent with a sulfur-based chemolithoautotrophic lifestyle. In contrast, those inhabiting the sunlit ocean displayed higher plasticity energy-related metabolic pathways, supporting a presumptive photoheterotrophic lifestyle. In epipelagic SAR324 ecotypes, we observed the presence of two types of proton-pumping rhodopsins, as well as genomic, transcriptomic, and ecological evidence for active photoheterotrophy, based on xanthorhodopsin-like light-harvesting proteins. Conclusions Combining pangenomic and both metagenomic and metatranscriptomic profiling revealed a striking divergence in the vertical distribution, genomic composition, metabolic potential, and predicted lifestyle strategies of geographically co-located members of the SAR324 bacterial clade. The results highlight the utility of metapangenomic approaches employed across environmental gradients, to decipher the properties and variation in function and ecological traits of specific phylogenetic clades within complex microbiomes.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2026

Machine learning identifies novel signatures of antifungal drug resistance in Saccharomycotina yeasts

Antifungal drug resistance is a major challenge in fungal infection management. Numerous genomic changes are known to contribute to acquired drug resistance in clinical isolates of specific pathogens, but whether they broadly explain natural resistance across entire lineages is unknown. We leveraged genomic, ecological, and phenotypic trait data from naturally sampled strains from nearly all known species in subphylum Saccharomycotina to examine the evolution of resistance to eight antifungal drugs. The phylogenetic distribution of drug resistance varied by drug; fluconazole resistance was widespread, while 5-fluorocytosine resistance was rare, except in Lipomycetales. A random forest algorithm trained on genomic data predicted drug-resistant yeasts with 54–75% accuracy. Fluconazole resistance was consistently predicted with the highest accuracy (75.2%). Furthermore, fluconazole resistance prediction accuracy was similar between models trained on genome-wide variation in the presence and number of InterPro protein annotations across Saccharomycotina (75.2%) and those trained on amino acid sequence alignment data of Erg11, a protein known to be involved in fluconazole resistance (74.3-74.9%). Interestingly, the top Erg11 residues for predicting fluconazole resistance across Saccharomycotina do not overlap with, are not spatially close to, and are less conserved than those previously linked to resistance in clinical isolates of Candida albicans. In silico deep mutational scanning of the C. albicans Erg11 protein reveals that amino acid variants implicated in clinical cases of resistance are almost universally destabilizing while variants in our most informative residues are energetically more neutral, explaining why the latter are much more common than the former in natural populations. Importantly, previous experimental analyses of C. albicans Erg11 have shown that amino acid variation in our most informative residues, despite having never been directly implicated in clinical cases, can directly contribute to resistance. Our results suggest that studies of natural resistance in yeast species never encountered in the clinic will yield a fuller understanding of antifungal drug resistance.

Harrison, Marie-Claire [Vanderbilt Univ., Nashvill↗

DOE OSTI2021

Ontology-Enriched Specifications Enabling Findable, Accessible, Interoperable, and Reusable Marine Metagenomic Datasets in Cyberinfrastructure Systems

Marine microbial ecology requires the systematic comparison of biogeochemical and sequence data to analyze environmental influences on the distribution and variability of microbial communities. With ever-increasing quantities of metagenomic data, there is a growing need to make datasets Findable, Accessible, Interoperable, and Reusable (FAIR) across diverse ecosystems. FAIR data is essential to developing analytical frameworks that integrate microbiological, genomic, ecological, oceanographic, and computational methods. Although community standards defining the minimal metadata required to accompany sequence data exist, they haven’t been consistently used across projects, precluding interoperability. Moreover, these data are not machine-actionable or discoverable by cyberinfrastructure systems. By making ‘omic and physicochemical datasets FAIR to machine systems, we can enable sequence data discovery and reuse based on machine-readable descriptions of environments or physicochemical gradients. In this work, we developed a novel technical specification for dataset encapsulation for the FAIR reuse of marine metagenomic and physicochemical datasets within cyberinfrastructure systems. This includes using Frictionless Data Packages enriched with terminology from environmental and life-science ontologies to annotate measured variables, their units, and the measurement devices used. This approach was implemented in Planet Microbe, a cyberinfrastructure platform and marine metagenomic web-portal. Here, we discuss the data properties built into the specification to make global ocean datasets FAIR within the Planet Microbe portal. We additionally discuss the selection of, and contributions to marine-science ontologies used within the specification. Finally, we use the system to discover data by which to answer various biological questions about environments, physicochemical gradients, and microbial communities in meta-analyses. This work represents a future direction in marine metagenomic research by proposing a specification for FAIR dataset encapsulation that, if adopted within cyberinfrastructure systems, would automate the discovery, exchange, and re-use of data needed to answer broader reaching questions than originally intended.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2024

Multicellular magnetotactic bacteria are genetically heterogeneous consortia with metabolically differentiated cells

Consortia of multicellular magnetotactic bacteria (MMB) are currently the only known example of bacteria without a unicellular stage in their life cycle. Because of their recalcitrance to cultivation, most previous studies of MMB have been limited to microscopic observations. To study the biology of these unique organisms in more detail, we use multiple culture-independent approaches to analyze the genomics and physiology of MMB consortia at single-cell resolution. We separately sequenced the metagenomes of 22 individual MMB consortia, representing 8 new species, and quantified the genetic diversity within each MMB consortium. This revealed that, counter to conventional views, cells within MMB consortia are not clonal. Single consortia metagenomes were then used to reconstruct the species-specific metabolic potential and infer the physiological capabilities of MMB. To validate genomic predictions, we performed stable isotope probing (SIP) experiments and interrogated MMB consortia using fluorescence in situ hybridization (FISH) combined with nanoscale secondary ion mass spectrometry (NanoSIMS). By coupling FISH with bioorthogonal noncanonical amino acid tagging (BONCAT), we explored their in situ activity as well as variation of protein synthesis within cells. We demonstrate that MMB consortia are mixotrophic sulfate reducers and that they exhibit metabolic differentiation between individual cells, suggesting that MMB consortia are more complex than previously thought. These findings expand our understanding of MMB diversity, ecology, genomics, and physiology, as well as offer insights into the mechanisms underpinning the multicellular nature of their unique lifestyle.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2023

Advancing the central role of non-model biorepositories in predictive modeling of emerging pathogens

The COVID-19 pandemic demonstrated the insufficiency of a reactive approach to emerging zoonotic pathogens. With spillover increasing in frequency as environments change and the human footprint continues to grow, pandemic prevention will require predictive models that can identify (i) potential zoonoses with a high likelihood of emergence and (ii) environmental or other features that may trigger a shift in host, vector, or pathogen baselines associated with emergence and/or spillover. Artificial intelligence (AI), and particularly its machine learning and deep learning branches, holds enormous potential for detecting shifts in large-scale biodiversity and disease datasets (genomic, ecological, geospatial, etc.). Such algorithms can be trained to identify subtle patterns in large volumes of data to yield insights into complex phenomena for which we have limited knowledge of the true cause(s) or predictor(s), as is the case for emerging infectious diseases.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2026

Resilient Crop Production for the Bioeconomy: Frontier Science for the Bioeconomy Workshop Series

A bioeconomy based on energy, chemicals, and bioproducts derived from nonfood crops promotes U.S. prosperity, energy independence, and national security. Crop productivity (i.e., yield) affects profitability for farmers, whereas biomass and seed oil quality affects profitability for biorefineries through the conversion efficiency of feedstock crops into energy, chemicals, and bioproducts. Scientific breakthroughs in plant genetics and genomics, ecological processes, bioprocessing, microbial engineering and design, and conversion technologies have improved the potential market viability of products derived from nonfood crops specifically cultivated as biomass and seed oil feedstocks. Improving understanding of a plant’s complex responses to withstand or recover from stress is crucial to maintaining feedstock yield and quality, particularly in challenging production zones. Stress may be related to biotic or abiotic (i.e., nonliving) environmental conditions that negatively affect plant growth, development, and productivity. Research and development to enhance the environmental resilience of feedstock crops will enable improved production by minimizing losses during stress and improving production on currently economically nonviable marginal lands.

09 BIOMASS FUELS↗

DOE OSTI2023

Insights into the Ecological Diversification of the Hymenochaetales based on Comparative Genomics and Phylogenomics With an Emphasis on Coltricia

To elucidate the genomic traits of ecological diversification in the Hymenochaetales, we sequenced 15 new genomes, with attention to ectomycorrhizal (EcM) Coltricia species. Together with published data, 32 genomes, including 31 Hymenochaetales and one outgroup, were comparatively analyzed in total. Compared with those of parasitic and saprophytic members, EcM species have significantly reduced number of plant cell wall degrading enzyme genes, and expanded transposable elements, genome sizes, small secreted proteins, and secreted proteases. EcM species still retain some of secreted carbohydrate-active enzymes (CAZymes) and have lost the key secreted CAZymes to degrade lignin and cellulose, while possess a strong capacity to degrade a microbial cell wall containing chitin and peptidoglycan. There were no significant differences in secreted CAZymes between fungi growing on gymnosperms and angiosperms, suggesting that the secreted CAZymes in the Hymenochaetales evolved before differentiation of host trees into gymnosperms and angiosperms. Nevertheless, parasitic and saprophytic species of the Hymenochaetales are very similar in many genome features, which reflect their close phylogenetic relationships both being white rot fungi. Phylogenomic and molecular clock analyses showed that the EcM genus Coltricia formed a clade located at the base of the Hymenochaetaceae and divergence time later than saprophytic species. And Coltricia remains one to two genes of AA2 family. These indicate that the ancestors of Coltricia appear to have originated from saprophytic ancestor with the ability to cause a white rot. This study provides new genomic data for EcM species and insights into the ecological diversification within the Hymenochaetales based on comparative genomics and phylogenomics analyses.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2025

Hybridization breaks species barriers in long-term coevolution of a cyanobacterial population

Bacterial species often undergo rampant recombination yet maintain cohesive genomic identity. Ecological differences can generate recombination barriers between species and sustain genomic clusters in the short term. But can these forces prevent genomic mixing during long-term coevolution? Cyanobacteria in Yellowstone hot springs comprise several diverse species that have coevolved for hundreds of thousands of years, providing a rare natural experiment. By analyzing more than 300 single-cell genomes, we show that despite each species forming a distinct genomic cluster, much of the diversity within species is the result of hybridization driven by selection, which has mixed their ancestral genotypes. This widespread mixing is contrary to the prevailing view that ecological barriers can maintain cohesive bacterial species and highlights the importance of hybridization as a source of genomic diversity.

Evolutionary Biology↗

DOE OSTI2022

METABOLIC: high-throughput profiling of microbial genomes for functional traits, metabolism, biogeochemistry, and community-scale functional networks

Background Advances in microbiome science are being driven in large part due to our ability to study and infer microbial ecology from genomes reconstructed from mixed microbial communities using metagenomics and single-cell genomics. Such omics-based techniques allow us to read genomic blueprints of microorganisms, decipher their functional capacities and activities, and reconstruct their roles in biogeochemical processes. Currently available tools for analyses of genomic data can annotate and depict metabolic functions to some extent; however, no standardized approaches are currently available for the comprehensive characterization of metabolic predictions, metabolite exchanges, microbial interactions, and microbial contributions to biogeochemical cycling. Results We present METABOLIC (METabolic And BiogeOchemistry anaLyses In miCrobes), a scalable software to advance microbial ecology and biogeochemistry studies using genomes at the resolution of individual organisms and/or microbial communities. The genome-scale workflow includes annotation of microbial genomes, motif validation of biochemically validated conserved protein residues, metabolic pathway analyses, and calculation of contributions to individual biogeochemical transformations and cycles. The community-scale workflow supplements genome-scale analyses with determination of genome abundance in the microbiome, potential microbial metabolic handoffs and metabolite exchange, reconstruction of functional networks, and determination of microbial contributions to biogeochemical cycles. METABOLIC can take input genomes from isolates, metagenome-assembled genomes, or single-cell genomes. Results are presented in the form of tables for metabolism and a variety of visualizations including biogeochemical cycling potential, representation of sequential metabolic transformations, community-scale microbial functional networks using a newly defined metric “MW-score” (metabolic weight score), and metabolic Sankey diagrams. METABOLIC takes ~ 3 h with 40 CPU threads to process ~ 100 genomes and corresponding metagenomic reads within which the most compute-demanding part of hmmsearch takes ~ 45 min, while it takes ~ 5 h to complete hmmsearch for ~ 3600 genomes. Tests of accuracy, robustness, and consistency suggest METABOLIC provides better performance compared to other software and online servers. To highlight the utility and versatility of METABOLIC, we demonstrate its capabilities on diverse metagenomic datasets from the marine subsurface, terrestrial subsurface, meadow soil, deep sea, freshwater lakes, wastewater, and the human gut. Conclusion METABOLIC enables the consistent and reproducible study of microbial community ecology and biogeochemistry using a foundation of genome-informed microbial metabolism, and will advance the integration of uncultivated organisms into metabolic and biogeochemical models. METABOLIC is written in Perl and R and is freely available under GPLv3 at https://github.com/AnantharamanLab/METABOLIC.

59 BASIC BIOLOGICAL SCIENCES↗