Search

26 pages, 5 figures, supporting information https://doi.org/10.1371/journal.pbio.3002083.-- Data Availability: All data directly relevant are within the paper and its Supporting Information files. Benchmarking analysis were based on the publicly available IMG/VR v3 database (doi 10.1093/nar/gkaa946 - https://genome.jgi.doe.gov/portal/IMG_VR/IMG_VR.home.html), The extraordinary diversity of viruses infecting bacteria and archaea is now primarily studied through metagenomics. While metagenomes enable high-throughput exploration of the viral sequence space, metagenome-derived sequences lack key information compared to isolated viruses, in particular host association. Different computational approaches are available to predict the host(s) of uncultivated viruses based on their genome sequences, but thus far individual approaches are limited either in precision or in recall, i.e., for a number of viruses they yield erroneous predictions or no prediction at all. Here, we describe iPHoP, a two-step framework that integrates multiple methods to reliably predict host taxonomy at the genus rank for a broad range of viruses infecting bacteria and archaea, while retaining a low false discovery rate. Based on a large dataset of metagenome-derived virus genomes from the IMG/VR database, we illustrate how iPHoP can provide extensive host prediction and guide further characterization of uncultivated viruses, BED was supported by the European Research Council (ERC) Consolidator grant 865694: DiversiPHI, the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under Germany’s Excellence Strategy – EXC 2051 – Project-ID 390713860, the Alexander von Humboldt Foundation in the context of an Alexander von Humboldt Professorship funded by the German Federal Ministry of Education and Research, and the European Union’s Horizon 2020 research and innovation program, under the Marie Skłodowska-Curie Actions Innovative Training Networks grant agreement no. 955974 (VIROINF). FHC was supported by a Juan de la Cierva - Incoporación fellowship (Grant IJC2019-039859-I), and had the institutional support of the “Severo Ochoa Centre of Excellence'' accreditation (CEX2019-000928-S). This work was supported by the U.S. Department of Energy, Office of Science, Biological and Environmental Research, Early Career Research Program (SR) awarded under UC-DOE Prime Contract DE-AC02-05CH11231. The work conducted by the U.S. Department of Energy Joint Genome Institute (https://ror.org/04xm1d337), a DOE Office of Science User Facility, is supported by the Office of Science of the U.S. Department of Energy operated under Contract No. DE-AC02-05CH11231 (SR, APC, SN)

Bacteriophages have important roles in the ecology of the human gut microbiome but are under-represented in reference databases. To address this problem, we assembled the Metagenomic Gut Virus catalogue that comprises 189,680 viral genomes from 11,810 publicly available human stool metagenomes. Over 75% of genomes represent double-stranded DNA phages that infect members of the Bacteroidia and Clostridia classes. Based on sequence clustering we identified 54,118 candidate viral species, 92% of which were not found in existing databases. The Metagenomic Gut Virus catalogue improves detection of viruses in stool metagenomes and accounts for nearly 40% of CRISPR spacers found in human gut Bacteria and Archaea. We also produced a catalogue of 459,375 viral protein clusters to explore the functional potential of the gut virome. This revealed tens of thousands of diversity-generating retroelements, which use error-prone reverse transcription to mutate target genes and may be involved in the molecular arms race between phages and their bacterial hosts.

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., Nature Biotechnology, 39 (4), ISSN:1546-1696, ISSN:1087-0156

This paper was originally published under standard Springer Nature copyright (© The Author(s), under exclusive licence to Springer Nature America, Inc.). It is now available as an open-access paper under a Creative Commons Attribution 4.0 International license. The error has been corrected in the print, HTML and PDF versions of the article.

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.

In the version of this article initially published, four people were missing from the alphabetical list of IMG/M Data Consortium members: Lauren V. Alteio of the Centre for Microbiology and Environmental Systems Science, University of Vienna, Vienna, Austria; Jeffrey L. Blanchard of the Biology Department, University of Massachusetts Amherst, Amherst, MA, USA; Kristen M. DeAngelis of the Department of Microbiology, University of Massachusetts Amherst, Amherst, MA, USA; and William Rodriguez-Reillo of the Research Computing Division, Harvard Medical School, Boston, MA, USA. The error has been corrected in the PDF and HTML versions of the article.

This is the open peer reviewers comments and recommendations regarding the submitted GigaScience article and/or dataset.

Additional file of Influence of the polar light cycle on seasonal dynamics of an Antarctic lake microbial community

Over the last several years, metagenomics has enabled the assembly of millions of new viral sequences that have vastly expanded our knowledge of Earth’s viral diversity. However, these sequences range from small fragments to complete genomes and no tools currently exist for estimating their quality. To address this problem, we developed CheckV, which is an automated pipeline for estimating the completeness of viral genomes as well as the identification and removal of non-viral regions found on integrated proviruses. After validating the approach on mock datasets, CheckV was applied to large and diverse viral genome collections, including IMG/VR and the Global Ocean Virome, revealing that the majority of viral sequences were small fragments, with just 3.6% classified as high-quality (i.e. > 90% completeness) or complete genomes. Additionally, we found that removal of host contamination significantly improved identification of auxiliary metabolic genes and interpretation of viral-encoded functions. We expect CheckV will be broadly useful for all researchers studying and reporting viral genomes assembled from metagenomes. CheckV is freely available at: http://bitbucket.org/berkeleylab/CheckV.

Additional file 2: Analysis of ‘unassigned’ taxa. Supplementary text, figures and tables.

Additional file 1: Ace Lake carbon cycle, expedition information, and taxonomic and functional analyses. Supplementary text, figures, tables and dataset.

Additional file 4: Viral analyses. Supplementary tables and dataset.

Bacteriophages from the Inoviridae family (inoviruses) are characterized by their unique morphology, genome content, and infection cycle. To date, a relatively small number of inovirus isolates have been extensively studied, either for biotechnological applications such as phage display, or because of their impact on the toxicity of known bacterial pathogens including Vibrio cholerae and Neisseria meningitidis. Here we show that the current 56 members of the Inoviridae family represent a minute fraction of a highly diverse group of inoviruses. Using a new machine learning approach leveraging a combination of marker gene and genome features, we identified 10,295 inovirus-like genomes from microbial genomes and metagenomes. Collectively, these represent six distinct proposed inovirus families infecting nearly all bacterial phyla across virtually every ecosystem. Putative inoviruses were also detected in several archaeal genomes, suggesting that these viruses may have occasionally transferred from bacterial to archaeal hosts. Finally, we identified an expansive diversity of inovirus-encoded toxin-antitoxin and gene expression modulation systems, alongside evidence of both synergistic (CRISPR evasion) and antagonistic (superinfection exclusion) interactions with co-infecting viruses which we experimentally validated in a Pseudomonas model. Capturing this previously obscured component of the global virosphere sparks new avenues for microbial manipulation approaches and innovative biotechnological applications.

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 to10,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.

Figure S1. Map of sampling sites. Figure S2. TNF PCA plots for SAGs illustrating separation of Nanoarchaeota and putative host genome bins. Figure S3. Maximum likelihood phylogeny of phylum Nanoarchaeota based on 16S rRNA gene sequences at least 400Â nt in length. Figure S4. Identification of putative host genome bins based on ANI to reference genomes and metagenome bins. Figure S5. Scanning electron micrograph of multiple Nanoarchaeota cells attached to host cells. Figure S6. SNP type and density in individual clade 1 Nanoarchaeota SAGs. Figure S7. Variation in sSNP density in clade 1 Nanoarchaeota genes by functional category. (DOCX 24585 kb)

Deep metagenomic sequencing has the potential to illuminate the intra-species genomic variation of abundant microbial species. In this thesis, I develop a new tool MIDAS (Metagenomic Intra-species Diversity Analysis System) for rapidly and automatically quantifying species abundance, single nucleotide polymorphisms (SNPs), and gene copy number variants (CNVs) from metagenomes. To illustrate the utility of this approach, I reanalyze three public datasets with MIDAS. First, I re-analyze stool metagenomes from 98 mother-infant pairs and used rare SNPs to track strain transmission. I find that early colonizers are likely transmitted from the mother whereas late colonizers are likely transmitted from the environment. Second, I re-analyze >300 stool metagenomes from healthy adults and use SNPs to identify examples of both strain co-existence and strain coexclusion. Third, I re-analyze 198 globally distributed marine metagenomes and used gene copy number variants to show that many species have population structure that correlates with geographic location. Strain level genetic variants clearly reveal extensive structure and dynamics that are obscured when metagenomes are analyzed at coarser taxonomic resolution.

We present the Metagenomic Intra-species Diversity Analysis System (MIDAS), which is an integrated computational pipeline for quantifying bacterial species abundance and strain-level genomic variation, including gene content and single-nucleotide polymorphisms (SNPs), from shotgun metagenomes. Our method leverages a database of more than 30,000 bacterial reference genomes that we clustered into species groups. These cover the majority of abundant species in the human microbiome but only a small proportion of microbes in other environments, including soil and seawater. We applied MIDAS to stool metagenomes from 98 Swedish mothers and their infants over one year and used rare SNPs to track strains between hosts. Using this approach, we found that although species compositions of mothers and infants converged over time, strain-level similarity diverged. Specifically, early colonizing bacteria were often transmitted from an infant's mother, while late colonizing bacteria were often transmitted from other sources in the environment and were enriched for spore-formation genes. We also applied MIDAS to 198 globally distributed marine metagenomes and used gene content to show that many prevalent bacterial species have population structure that correlates with geographic location. Strain-level genetic variants present in metagenomes clearly reveal extensive structure and dynamics that are obscured when data are analyzed at a coarser taxonomic resolution.

The phylum Actinobacteria includes important human pathogens like Mycobacterium tuberculosis and Corynebacterium diphtheriae and renowned producers of secondary metabolites of commercial interest, yet only a small part of its diversity is represented by sequenced genomes. Here, we present 824 actinobacterial isolate genomes in the context of a phylum-wide analysis of 6,700 genomes including public isolates and metagenome-assembled genomes (MAGs). We estimate that only 30%–50% of projected actinobacterial phylogenetic diversity possesses genomic representation via isolates and MAGs. A comparison of gene functions reveals novel determinants of host-microbe interaction as well as environment-specific adaptations such as potential antimicrobial peptides. We identify plasmids and prophages across isolates and uncover extensive prophage diversity structured mainly by host taxonomy. Analysis of >80,000 biosynthetic gene clusters reveals that horizontal gene transfer and gene loss shape secondary metabolite repertoire across taxa. Our observations illustrate the essential role of and need for high-quality isolate genome sequences. The work (proposal DOI[s]: https://doi.org/10.46936/10.25585/60001024; https://doi.org/10.46936/10.25585/60000886; https://doi.org/10.46936/10. 25585/60001401; https://doi.org/10.46936/10.25585/60001087; https:// doi.org/10.46936/10.25585/60001079; https://doi.org/10.46936/10.25585/ 60001044) conducted by the US Department of Energy Joint Genome Institute (https://ror.org/04xm1d337), a DOE Office of Science User Facility, is supported by the Office of Science of the US Department of Energy operated under contract no. DE-AC02-05CH11231. Design and synthesis of peptide constructs were supported by the Biosystems Design program (DOE Office of Science contract DE-SC0018260) and the US Department of Energy Joint Genome Institute (DOE Office of Science contract DEAC02-05CH1123), respectively. Characterization of the peptide was supported by the Secured Biosystem Design project entitled ‘‘Rapid Design and Engineering of Smart and Secure Microbiological Systems’’ (DOE Office of Science contract DE-AC02-05CH1123). We thank the following researchers for their support of this study by providing free use of their public genome data: Kristen De-Angelis, Grace Pold, Mallory Choudoir, Camila Carlos-Shanley, Paul Carini, H. Corby C. Kistler, James Elkins, Javier A. Izquierdo, Dimitris Hatzinikolaou, Daniel Schachtman, Paul R. Jensen, Aindrila Mukhopadhyay, John Vogel, Carolin Frank, Paul M. D’Agostino, Ann M. Hirsch, Satoshi Yuzawa, Regina Lamendella, Bernhard Fuchs, Dale Pelletier, Laila P. Partida-Martinez, Cameron Currie, Seth De-Bolt, Jeff Dangl, David Mead, Susannah Tringe, David A. Baltrus, Seung Bum Kim, Linda Kinkel, Kelly Wrighton, William Mohn, Ludmila Christoserdova, Sarah Lebeis, Janet Janssen, Sandra Baena Garzon, and Nicholas Coleman. Special thanks to Marie Louise Ballon at the JGI Communications and Outreach office for all her help with the graphical abstract and many other details.