Main

Current estimates are that virus particles massively outnumber live cells in most habitats1,2, but only a tiny fraction of viruses have been cultivated in the laboratory. An unprecedented diversity of viruses are being discovered through culture-independent sequencing3. Progress has been made in reconstructing genomes of uncultivated viruses de novo, from biotic and abiotic environments, without laboratory isolation of the virus–host system. For example, in the past 2 years, more than 750,000 uncultivated virus genomes (UViGs) have been identified in metagenome and metatranscriptome datasets4,5,6,7,8,9, five times the total number of genomes sequenced from virus isolates (Fig. 1), and UViGs already represent ≥95% of the taxonomic diversity in publicly available virus sequences10,11. Although double-stranded DNA (dsDNA) genomes are over-represented in UViGs because most metagenomic protocols exclusively target dsDNA, UViGs nonetheless enable an assessment of global virus diversity and an evaluation of structure and drivers of viral communities. UViGs also contribute to improving our understanding of the evolutionary history of viruses and virus–host interactions.

Figure 1: Size of virus genome databases over time4,7,22,45,83,84,85,86,87,88,89.
figure 1

Genome sequences from isolates (blue and green) or from UViGs (yellow) are shown. For genomes from isolates, the total number of genomes (blue) and the number of 'reference' genomes (green) are shown. Data were downloaded using the queries “Viruses[Organism] AND srcdb_refseq[PROP] NOT wgs[PROP] NOT cellular organisms[ORGN] NOT AC_000001:AC_999999[PACC]” for reference genomes and “Viruses[Organism] NOT cellular organisms[ORGN] NOT wgs[PROP] NOT AC_000001:AC_999999[pacc] NOT gbdiv syn[prop] AND nuccore genome samespecies[Filter]” for total number of virus genomes, on the NCBI nucleotide database portal (https://www.ncbi.nlm.nih.gov/nuccore) in January 2018. Genomes from the influenza virus database (https://www.ncbi.nlm.nih.gov/genomes/FLU/Database/nph-select.cgi?go=genomeset) were also added to the total number of virus genomes. UViGs can be assembled from metagenomes, from proviruses identified in microbial genomes, or from single-virus genomes, and estimated total UViG numbers were obtained by compiling data from the literature and from the total number of sequences in the IMG/VR database in January 2017, January 2018 and July 2018 (https://img.jgi.doe.gov/vr/)11. UpViG, uncultivated provirus.

Analysis and interpretation of standalone genomes present substantial challenges, whether the genomes are eukaryotic, bacterial, archaeal or viral. To address these challenges, MISAG and MIMAG standards were drafted to improve the quality of reporting of microbial genomes derived from single cell or metagenome sequences, which are often incomplete12. Although some aspects of MISAG and MIMAG can be applied to UViGs, the extraordinary diversity of viral genome composition and content, replication strategies, and hosts means that the completeness, quality, taxonomy and ecology of UViGs need to be evaluated via virus-specific metrics.

The Genomic Standards Consortium (http://gensc.org) maintains metadata checklists for MIxS, encompassing genome and metagenome sequences13, marker gene sequences14 and single amplified and metagenome-assembled bacterial and archaeal genomes12. Here we present a set of standards that extend the MIxS checklists to include identification, quality assessment, analysis and reporting of UViGs (Table 1 and Supplementary Tables 1 and 2), together with recommendations on how to perform these analyses. We provide a metadata checklist for database submission and publication of UViGs designed to be flexible enough to accommodate technological and methodological changes over time (Table 1 and Supplementary Table 1). The information gathered through the MIUViG checklist can be directly submitted with new UViG sequences to International Nucleotide Sequence Database Collaboration (INSDC) member databases—the DNA Database of Japan (DDBJ), the European Molecular Biology Laboratory–European Bioinformatics Institute (EMBL-EBI) and US National Center for Biotechnology Information (NCBI)—which will host and display checklist metadata alongside the UViG sequence. These MIUViG standards should also be used along with existing guidelines for virus genome analysis, including those issued by the International Committee on Taxonomy of Viruses (ICTV), which recently endorsed the incorporation of UViGs into the official virus classification scheme15 (https://talk.ictvonline.org). Although MIUViG standards and best practices were designed for genomes of viruses infecting microorganisms, they can also be applied to viruses infecting animals, fungi and plants, and are compatible with standards that are already in place for epidemiological analysis of these viruses16 (Supplementary Table 3).

Table 1 List of mandatory metadata for UViGs

Recovery of UViGs after virus enrichment

UViGs can be retrieved from datasets enriched for virus genomes, namely viral metagenomes and single-virus genomes (Fig. 2). Viral metagenomes are usually obtained through a combination of filtration steps, DNase or RNase treatments, and RNA or DNA extraction depending on the targeted viruses, then reverse transcription (to find RNA viruses) and shotgun sequencing3,17,18,19. Targeted sequence capture methods can be applied to recover specific virus groups (Fig. 2), and these methods have proven especially useful when viruses are present in small amounts (for example, clinical samples)20. Single-virus methods use flow cytometry to sort individual viral particles before genome amplification and sequencing, to produce viral single amplified genomes (SAGs)9,21,22,23 (Fig. 2). Viral metagenomes and single-virus genomes are usually sequenced with short-read, high-throughput technologies, such as Illumina sequencing, and assembled by algorithms similar to those used for microbial genomes and metagenomes. However, owing to their relatively small genome size (92% of virus genomes in the NCBI Viral RefSeq database are <100 kb)10, short read-based genome assemblies could soon be superseded by long-read sequencing technologies24 (for example, PacBio zero-mode waveguide technology or Oxford Nanopore Technology nanopore sequencing; Fig. 2). Sequencing virus genomes from a single template would notably enable the identification of individual genotypes in mixed populations.

Figure 2: Identification of UViGs.
figure 2

Schematic of methods used to obtain UViGs. Steps that have been adapted from those used to assemble MAGs and SAGs12 or added for UViG are shown for sample preparation (orange) and bioinformatics analysis (blue). Steps specifically required for virus targeting and identification are highlighted in bold. *For viruses with short genomes, long-read technologies can provide complete genomes from shotgun sequencing in a single read, bypassing the assembly step24. **Targeted sequence capture can be used to recover viral genomes from a known virus group. These genomes can be recovered from samples in which they represent a small fraction of the templates (for example, clinical samples20).

The main advantages of datasets produced after enrichment for viruses are good de novo assembly of both abundant and rare viruses, increased confidence that the sequence is of viral origin, and the ability to sequence both active and 'inactive' or 'cryptic' viruses (i.e., viruses that are present in the sample but cannot infect). However, virus-enriched datasets can have over-representation of virulent viruses with high burst size (high number of virus particles released from each infected cell) and under-representation of larger viruses with capsids ≥0.2 μm, such as giant viruses, as a result of the selective filtration steps used25. Furthermore, in silico approaches are often the only option available to determine the host range of UViGs obtained from virus-enriched samples.

Recovery of UViGs without enrichment

Virus sequences are also present in non-virus-enriched datasets, including sorted cells, tissues, or environmental samples collected on 0.2 μm filters4,26,27,28. These sequences could originate from viruses that are replicating in cells, from temperate viruses (proviruses or prophages) that are either integrated into host genomes or present as episomal elements in the host cell, or from free virus particles present in samples.

Analyzing datasets without virus enrichment has several advantages. It can detect lytic, temperate and persistent infection, it overcomes some of the biases arising from the size-based selection of virus particles, and it can be applied to any metagenome. However, UViGs from non-virus-enriched datasets may be biased toward viruses that infect the dominant host cell in the sample, and rare viruses or those infecting rare hosts could be under-represented or absent. Finally, comparisons between virus-enriched and non-virus-enriched datasets suggest that analyzing UViGs across different size fractions and sample types is valuable for exploring the virus genome sequence space29 (Supplementary Fig. 1 and Supplementary Note 1).

Computational identification of viral sequences

Regardless of the type of dataset, the viral origin of UViGs must be validated because even samples enriched for virus particles still contain a substantial amount of cellular DNA30. Contamination can arise either from difficulty in separating virus particles from cellular fractions (for example, ultra-small bacteria31) or from the capture of extracellular DNA in the virus fraction. Cellular sequences can also derive from cell genome fragments that are encased in virus capsids or comparable particles (for example, via transduction), DNA-containing membrane vesicles, or gene transfer agents32,33,34.

Several bioinformatic tools and protocols have been developed to identify sequences from bacteriophages and archaeal viruses35,36,37,38; eukaryotic viruses39; or combinations of bacteriophages, archaeal viruses and large eukaryotic viruses40 (Supplementary Table 4). These approaches rely on a few characteristics, such that a sequence is considered viral if it is significantly similar to known viruses (in terms of gene content or nucleotide usage pattern) or if it is unrelated to any known virus and cellular genome but contains one or more hallmark virus genes. UViGs must therefore be accompanied by a list of virus detection tool(s) and protocol(s) used, together with any thresholds applied (Table 1 and Supplementary Table 1).

Identification of integrated proviruses and their precise boundaries in the host genome is problematic (Box 1). Notably, no high-throughput approach can accurately distinguish active proviruses (still able to replicate and produce virions) from inactive proviral remnants of a past infection28. Thus, although prediction methods are improving, UViGs identified as proviruses should be clearly marked as such, so that these caveats are clear (Table 1 and Supplementary Table 1).

Estimating quality of UViGs

We propose three categories of UViG sequences: genome fragment(s), high-quality draft genomes and finished genomes (Fig. 3 and Table 2). These categories mirror those in MISAG and MIMAG12, and they are matched to categories already proposed for complete-genome sequencing of small viruses in epidemiology and surveillance16 (Supplementary Table 3). UViG quality is more challenging to evaluate than metagenome-assembled genomes (MAGs) or SAGs because most viruses lack conserved sets of single-copy marker genes that can be used to estimate draft genome completeness. However, exceptions exist, such as large eukaryotic dsDNA viruses. To date, researchers have estimated UViG sequence completeness by identifying circular contigs or contigs with inverted terminal repeats as putative complete genomes. For linear contigs, completeness is estimated by comparison to reference genome sequences and typically requires a taxonomic assignment to a (candidate) (sub)family or genus because genome length is relatively homogeneous at these ranks (±10%; Supplementary Fig. 2 and Supplementary Table 5). This assignment can be based on the detection of specific marker genes, such as clade-specific viral orthologous groups (Supplementary Table 6), or based on genome-based classification tools (see “Taxonomy of UViGs”). Estimating completeness is more difficult for segmented genomes, which require either a closely related reference genome or additional in vitro experiments16. A detailed example of how this quality tier classification can be performed on the Global Ocean Virome dataset7 is presented in Supplementary Note 2 and Supplementary Table 7.

Figure 3: UViG classification and associated sequence analyses.
figure 3

“Functional potential” is functional annotation used in gene content analysis. “Host prediction” is the application of different in silico host prediction tools. “Taxonomic classification” is classification of the contig to established groups using marker genes or gene content comparison. “Diversity and distribution” includes vOTU clustering and relative abundance estimation through metagenome read mapping, at the geographical scale or across anatomical sites for host-associated datasets. “New taxonomic groups” concerns the delineation of new proposed groups (for example, families or genera) based exclusively on UViG sequences. “New reference species” refers to the proposal of a new entry in ICTV (https://talk.ictvonline.org/files/taxonomy-proposal-templates/). *Some of these approaches require a minimum contig size—for example, contigs ≥10 kb for taxonomic classification based on gene content59 or diversity estimation47—and will not be applicable to every genome fragment.

Table 2 Summary of required characteristics for each category

Contigs or genome bins representing <90% of the expected genome length, or for which no expected genome length can be determined, would be considered genome fragments. This category might include UViG fragments large enough to be assigned to known virus groups on the basis of gene content and average nucleotide identity. However, high-quality draft or finished genomes are required to establish new taxa (Fig. 3). Sequences from UViG fragments can be used in phylogenetic and diversity studies, either as references for virus operational taxonomic units (see Supplementary Note 4), or through the analysis of virus marker genes encoded in these genome fragments; for example, capsid proteins, terminases, ribonucleotide reductases and DNA- or RNA-dependent RNA polymerases41,42,43,44,45,46. Similarly, UViG fragments can be analyzed to assess the functional gene complement of unknown viruses or link them to potential hosts. Importantly, current methods for automatic virus sequence identification35,36,37,38,39,40 cannot reliably identify short (<10 kb) viral sequences, which should be interpreted with utmost caution.

Contigs or genome bins either predicted as complete or representing ≥90% of the expected genome sequence are high-quality drafts, consistent with standards for microbial genomes12. Repeat regions may lead to erroneous assembly of partial genomes as circular contigs47. Thus, the length of the assembled circular contig should be considered when assessing UViG completeness (Box 1). For UViGs not derived from a consensus assembly, such as single long reads, base calling quality >99% on average (phred score >20) is needed to assign a “high-quality draft” label. Genome sequences assembled into a single contig, or one per segment, with extensive manual review and annotation, can be labeled “finished genomes.” Annotation must include identification of putative gene functions; structural, replication or lysogeny modules; and transcriptional units. The “finished genomes” category is reserved for only the highest quality, manually curated UViGs and is required for the establishment of new virus species (Fig. 3 and Table 2).

Unlike that of SAGs and MAGs12, quality estimation of UViGs does not include a genome contamination threshold. Contamination issues are most prominent in the case of genome bins, whereas most UViGs are represented by a single contig for which in silico simulations have shown that chimeric sequences are rare and present at <2% (ref. 47). In addition, no tools exist to automatically estimate UViG contamination, and thus this information is not included in the current MIUViG checklist. A future updated version of the MIUViG checklist may, however. For include contamination thresholds if such a tool were to be developed. For example, such a tool might exploit single-copy marker genes (once these have been defined for a broader range of viruses) or it might use coverage by metagenome reads, which should in principle be evenly distributed along the genome with no major deviance, except for highly conserved genes.

Annotation of UViGs

Functional annotation of UViGs comprises the following tasks: predicting features in the genome sequence, such as protein-coding genes, tRNAs and integration sites; assigning functions to as many predicted features as possible; and assigning the remaining hypothetical proteins to uncharacterized protein families. Annotation pipelines have been established for different types of viruses48,49, and large differences between viral genome types likely preclude the development of a single tool able to annotate every virus50. Therefore, we recommend that software used to annotate UViGs be reported (Supplementary Table 1).

The choice of methods and reference databases used to annotate predicted proteins should be clearly stated. Homologs of novel virus genes may not be detected with standard methods for pairwise sequence similarity detection, such as BLAST, but instead require the use of more sensitive profile similarity approaches, such as HMMER51, PSI-BLAST52 or HHPred53 (Supplementary Table 8; reviewed in ref. 54). Although sequence profiles for many protein families have been collected, they frequently remain unassociated with any specific function. Therefore, UViG analyses should always report (i) feature prediction method(s), (ii) sequence similarity search method(s), and (iii) database(s) searched (Box 1 and Supplementary Table 1).

Taxonomy of UViGs

Taxonomic classification can provide information on the relationship of a UViG with known viruses. Although the information and criteria used for virus classification have changed over time, virus classification has now converged to genome-based analyses15 (Box 2). The ICTV established specific demarcation criteria for each virus group (Supplementary Table 9) owing to the vast range of viral genomes, mutation rates and evolution. Recently, a consensus has emerged on using whole-genome average nucleotide identity for classification at the species rank, which is used in downstream ecological, evolutionary and functional studies. This consensus was reached through analysis of published population genetics studies55,56 and gene content comparison of NCBI RefSeq10 virus genomes57,58,59 (Supplementary Note 3 and Supplementary Fig. 3). We propose to formalize the use of species-rank virus groups and to name these “virus operational taxonomic units” (vOTUs) to avoid confusion because species groups have been variously named “viral population,” “viral cluster” or “contig cluster” in the literature4,7,60. We suggest standard thresholds of 95% average nucleotide identity over 85% alignment fraction (relative to the shorter sequence) on the basis of a comparison of sequences currently available in NCBI RefSeq10 and IMG/VR11 (Supplementary Note 3 and Supplementary Figs. 3 and 4). Although partial genomes remain challenging to classify, these common thresholds will enable comparative analyses (Supplementary Fig. 5). In addition, vOTU reports should include the clustering method and cutoff, the reference database used (if any), and the genome alignment approach because small differences have been observed between different methods61 (Supplementary Table 1).

For higher taxonomic ranks than species, no consensus has been reached on which approach should be used, although several have been proposed58,59,62,63,64,65,66. Keeping this in mind, UViG reports including taxonomy must clearly indicate the methods and cutoffs applied, and any new taxon must be highlighted as preliminary (for example, “genus-rank cluster,” “putative genus” or “candidate genus,” but not simply “genus,” as this category is reserved for ICTV-recognized groups; Supplementary Table 1). Authors should submit formal taxonomic proposals to the ICTV for consideration (https://talk.ictvonline.org/files/taxonomy-proposal-templates/).

Finally, information about the nature of the genome and mode of expression (i.e., Baltimore classification67) should be included in the UViG description. Similarly, the predicted segmentation state of the genome (segmented or nonsegmented) should be reported, typically derived from taxonomic classification and comparison with the closest references (Supplementary Table 1).

In silico host prediction

Once a new virus genome has been assembled, an important step toward understanding the ecological role of the associated virus is to predict its host(s). In silico approaches are often the only option for UViGs (reviewed in ref. 68; Supplementary Table 10). These can be separated into four main types. First, hosts can be predicted with relatively high precision on the basis of sequence similarity between the UViG and a reference virus genome when a closely related virus is available69,70. Second, hosts can be predicted on the basis of sequence similarities between a UViG and a host genome. These sequence similarities can range from short exact matches (20–100 bp), which include CRISPR spacers4,7,68,71, to longer (>100 bp) nucleotide sequence matches, including proviruses integrated into a larger host contig26,68,72,73 (Supplementary Table 10). Host-range predictions based on sequence similarity are the most reliable but require that a closely related host genome has been sequenced68. Third, host taxonomy from domain down to genus rank can be predicted from nucleotide usage signatures reflecting coevolution between virus and host genomes in terms of G+C content, k-mer frequency and codon usage26,74,75. These approaches are usually less specific than sequence similarity–based ones and cannot reliably predict host range below the genus rank, but can provide a predicted host for a larger number of UViGs7 (Supplementary Table 10). Finally, host predictions can be computed from a comparison of abundance profiles of host and virus sequences across spatial or temporal scales, either through abundance correlation25,76,77,78 or through more sophisticated model-based interaction predictors79. Although few datasets are available for robust evaluation of host prediction based on comparison of abundance profiles, we expect this approach to become more powerful and relevant as high-resolution time-series metagenomics becomes more common.

As all these bioinformatic approaches remain predictive, it is crucial that robust false-discovery rate estimations are reported (Supplementary Table 1). Moreover, computational tools do not predict quantitative infection characteristics (for example, infection rate or burst size), which are important for understanding the impacts of viruses on host biology, and thus far only apply to viruses infecting bacteria or archaea. Nevertheless, these predictions are important guides for subsequent in silico, in vitro and in vivo studies, including experimental validation to unequivocally demonstrate a viral infection of a given microbial host. Host predictions should be reported along with details regarding the specific tool(s) used and, importantly, their estimated accuracy as derived either from published benchmarks or from tests conducted in the study (Supplementary Table 1). This information will allow virus–host databases69,80 to progressively incorporate UViGs while still controlling for the sensitivity and accuracy of the predictions provided to users.

Reporting UViGs

We recommend the following best practice for sharing and archiving UViGs and UViG-related data: data publication should center on the data resources of INSDC (http://www.insdc.org/) through one of the member databases, at DDBJ (https://www.ddbj.nig.ac.jp/index-e.html), EMBL-EBI's European Nucleotide Archive (ENA; https://www.ebi.ac.uk/ena) or NCBI (GenBank and the Sequence Read Archive; https://www.ncbi.nlm.nih.gov/nucleotide). If needed, INSDC database curators can be contacted directly for large-scale batch dataset submissions. Where new datasets are generated as part of a UViG study, sequenced samples should be described according to the environment-relevant MIxS checklists and raw read data should be submitted. High-quality and finished UViGs should be submitted as assemblies, the former reported as “draft” accompanied by the required metadata (Table 1). Incomplete assemblies may be submitted, but they must be accompanied by the required metadata (Table 1 and Supplementary Table 1).

Where available, annotation and taxonomic classification should be submitted to INSDC, and occurrence and abundance data reported as 'Analysis' records in the ENA. Reports of abundance data estimated by short-read metagenome mapping should include information about the nucleotide identity and coverage thresholds used, with corresponding estimates of false-positive and false-negative rates either computed de novo or extracted from the literature (for example, from refs. 47, 81; Supplementary Note 4). All INSDC accession codes must be cited in publications. For ICTV classification, only coding-complete genomes (complete high-quality and finished draft UViGs) are currently considered82.

Conclusions

MIUViG standards and best practices for UViG analysis are the virus-specific counterparts to MISAG and MIMAG12. Virus genomics and metagenomics are rapidly expanding and improving as sequencing technologies emerge and mature. At the same time, the development of genome-based virus taxonomy methods as well as unified, comprehensive, and annotated reference databases of virus genomes and/or proteins continues apace. Community adoption of these standards, including through ongoing collaborations with other virus committees (ICTV) and data centers (DDBJ, EMBL-EBI and NCBI), will provide a framework for a systematic exploration of viral genome sequence space and enable the research community to better utilize and report UViGs.

Additional information

Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.