Abstract:Microbes live in complex communities that are of major importance for environmental ecology, public health, and animal physiology and pathology. Short-read metagenomic shotgun sequencing is currently the state-of-the-art technique for exploring these communities. With the aid of metagenomics, our understanding of the microbiome is moving from composition toward functionality, even down to the genetic variant level. While the exploration of single-nucleotide variation in a genome is a standard procedure in genomics, and many sophisticated tools exist to perform this task, identification of genetic variation in metagenomes remains challenging. Major factors that hamper the widespread application of variant-calling analysis include low-depth sequencing of individual genomes (which is especially significant for the microorganisms present in low abundance), the existence of large genomic variation even within the same species, the absence of comprehensive reference genomes, and the noise introduced by next-generation sequencing errors. Some bioinformatics tools, such as metaSNV or InStrain, have been created to identify genetic variants in metagenomes, but the performance of these tools has not been systematically assessed or compared with the variant callers commonly used on single or pooled genomes. In this study, we benchmark seven bioinformatic tools for genetic variant calling in metagenomics data and assess their performance. To do so, we simulated metagenomic reads to mimic human microbial composition, sequencing errors, and genetic variability. We also simulated different conditions, including low and high depth of coverage and unique or multiple strains per species. Our analysis of the simulated data shows that probabilistic method-based tools such as HaplotypeCaller and Mutect2 from the GATK toolset show the best performance. By applying these tools to longitudinal gut microbiome data from the Human Microbiome Project, we show that the genetic similarity between longitudinal samples from the same individuals is significantly greater than the similarity between samples from different individuals. Our benchmark shows that probabilistic tools can be used to call metagenomes, and we recommend the use of GATK's tools as reliable variant callers for metagenomic samples.

GenMPI: Cluster Scalable Variant Calling for Short/Long Reads Sequencing Data

ADS-HCSpark: A Scalable HaplotypeCaller Leveraging Adaptive Data Segmentation to Accelerate Variant Calling on Spark.

VC@Scale: Scalable and high-performance variant calling on cluster environments

MinimapR: A Parallel Alignment Tool for the Analysis of Large-Scale Third-Generation Sequencing Data

Integer programming framework for pangenome-based genome inference

Sentieon DNASeq Variant Calling Workflow Demonstrates Strong Computational Performance and Accuracy

Accelerating Minimap2 for Accurate Long Read Alignment on GPUs

Cloud-native distributed genomic pileup operations

Computational performance and accuracy of Sentieon DNASeq variant calling workflow

A Benchmark of Genetic Variant Calling Pipelines Using Metagenomic Short-Read Sequencing

Scalable and efficient DNA sequencing analysis on different compute infrastructures aiding variant discovery

PyroHMMsnp: an SNP Caller for Ion Torrent and 454 Sequencing Data.

A statistical framework for SNP calling, mutation discovery, association mapping and population genetical parameter estimation from sequencing data

Joint Haplotype Phasing and Genotype Calling of Multiple Individuals Using Haplotype Informative Reads

Minimap and miniasm: fast mapping and de novo assembly for noisy long sequences

Minimap2: pairwise alignment for nucleotide sequences

GMAP and GSNAP for Genomic Sequence Alignment: Enhancements to Speed, Accuracy, and Functionality

Psi-Caller: A Lightweight Short Read-Based Variant Caller with High Speed and Accuracy

High-Performance Genomic Analysis Heterogeneous System Using OpenCL

SMaSH: A Benchmarking Toolkit for Human Genome Variant Calling

MaxSSmap: A GPU program for mapping divergent short reads to genomes with the maximum scoring subsequence