Abstract:BackgroundIdentifying viral sequences in mixed metagenomes containing both viral and host contigs is a critical first step in analyzing the viral component of samples. Current tools for distinguishing prokaryotic virus and host contigs primarily use gene-based similarity approaches. Such approaches can significantly limit results especially for short contigs that have few predicted proteins or lack proteins with similarity to previously known viruses.MethodsWe have developed VirFinder, the first k-mer frequency based, machine learning method for virus contig identification that entirely avoids gene-based similarity searches. VirFinder instead identifies viral sequences based on our empirical observation that viruses and hosts have discernibly different k-mer signatures. VirFinder’s performance in correctly identifying viral sequences was tested by training its machine learning model on sequences from host and viral genomes sequenced before 1 January 2014 and evaluating on sequences obtained after 1 January 2014.ResultsVirFinder had significantly better rates of identifying true viral contigs (true positive rates (TPRs)) than VirSorter, the current state-of-the-art gene-based virus classification tool, when evaluated with either contigs subsampled from complete genomes or assembled from a simulated human gut metagenome. For example, for contigs subsampled from complete genomes, VirFinder had 78-, 2.4-, and 1.8-fold higher TPRs than VirSorter for 1, 3, and 5 kb contigs, respectively, at the same false positive rates as VirSorter (0, 0.003, and 0.006, respectively), thus VirFinder works considerably better for small contigs than VirSorter. VirFinder furthermore identified several recently sequenced virus genomes (after 1 January 2014) that VirSorter did not and that have no nucleotide similarity to previously sequenced viruses, demonstrating VirFinder’s potential advantage in identifying novel viral sequences. Application of VirFinder to a set of human gut metagenomes from healthy and liver cirrhosis patients reveals higher viral diversity in healthy individuals than cirrhosis patients. We also identified contig bins containing crAssphage-like contigs with higher abundance in healthy patients and a putative Veillonella genus prophage associated with cirrhosis patients.ConclusionsThis innovative k-mer based tool complements gene-based approaches and will significantly improve prokaryotic viral sequence identification, especially for metagenomic-based studies of viral ecology.

Eukfinder: a pipeline to retrieve microbial eukaryote genomes from metagenomic sequencing data

Accurate and sensitive detection of microbial eukaryotes from whole metagenome shotgun sequencing

Probing the eukaryotic microbes of ruminants with a deep-learning classifier and comprehensive protein databases

Metagenomic discovery of microbial eukaryotes in stool microbiomes

Microbiome Single Cell Atlases Generated with a Commercial Instrument

UnigeneFinder: An automated pipeline for gene calling from transcriptome assemblies without a reference genome

MitoFinder: Efficient automated large‐scale extraction of mitogenomic data in target enrichment phylogenomics

DeepMicrobeFinder Sorts Metagenomes into Prokaryotes, Eukaryotes and Viruses, with Marine Applications

MetaMap: an Atlas of Metatranscriptomic Reads in Human Disease-Related RNA-seq Data

EUKulele: Taxonomic annotation of the unsung eukaryotic microbes

Soil Nematode Community Profiling Using Reference-Free Mito-Metagenomics

Reference-free mito-metagenomics using deep-learning method

VirusSeeker, a computational pipeline for virus discovery and virome composition analysis

Bin Chicken: targeted metagenomic coassembly for the efficient recovery of novel genomes

VirFinder: a novel k-mer based tool for identifying viral sequences from assembled metagenomic data

Recovery of genomes from metagenomes via a dereplication, aggregation and scoring strategy

Metaxa: a software tool for automated detection and discrimination among ribosomal small subunit (12S/16S/18S) sequences of archaea, bacteria, eukaryotes, mitochondria, and chloroplasts in metagenomes and environmental sequencing datasets

Metabolome searcher: a high throughput tool for metabolite identification and metabolic pathway mapping directly from mass spectrometry and using genome restriction

GutMeta: online microbiome analysis and interactive visualization with build-in curated human gut microbiome database

Annotating microbial functions with ProkFunFind

: A Feature-Rich Metagenomic Analysis Tool for Accessible and Comprehensive Metagenomic Exploration