Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2022 Dec 28:12:giac125.
doi: 10.1093/gigascience/giac125. Epub 2023 Jan 16.

xAtlas: scalable small variant calling across heterogeneous next-generation sequencing experiments

Jesse Farek  1 , Daniel Hughes  1   2 , William Salerno  1   3 , Yiming Zhu  1 , Aishwarya Pisupati  1 , Adam Mansfield  1   3 , Olga Krasheninina  1   3 , Adam C English  1 , Ginger Metcalf  1 , Eric Boerwinkle  1   4 , Donna M Muzny  1 , Richard Gibbs  1 , Ziad Khan  1 , Fritz J Sedlazeck  1
Affiliations

xAtlas: scalable small variant calling across heterogeneous next-generation sequencing experiments

Jesse Farek et al. Gigascience. .

Abstract

Background: The growing volume and heterogeneity of next-generation sequencing (NGS) data complicate the further optimization of identifying DNA variation, especially considering that curated high-confidence variant call sets frequently used to validate these methods are generally developed from the analysis of comparatively small and homogeneous sample sets.

Findings: We have developed xAtlas, a single-sample variant caller for single-nucleotide variants (SNVs) and small insertions and deletions (indels) in NGS data. xAtlas features rapid runtimes, support for CRAM and gVCF file formats, and retraining capabilities. xAtlas reports SNVs with 99.11% recall and 98.43% precision across a reference HG002 sample at 60× whole-genome coverage in less than 2 CPU hours. Applying xAtlas to 3,202 samples at 30× whole-genome coverage from the 1000 Genomes Project achieves an average runtime of 1.7 hours per sample and a clear separation of the individual populations in principal component analysis across called SNVs.

Conclusions: xAtlas is a fast, lightweight, and accurate SNV and small indel calling method. Source code for xAtlas is available under a BSD 3-clause license at https://github.com/jfarek/xatlas.

PubMed Disclaimer

Conflict of interest statement

The authors declare no conflicts of interest.

Figures

Figure 1:
Figure 1:
The major stages of xAtlas variant calling process.
Figure 2:
Figure 2:
(A) Variant caller runtimes of xAtlas and 4 other small variant callers (single-threaded operation) on 3 Illumina HiSeq X NA12878 whole-genome samples at 30× average coverage. (B, C) Precision vs. sensitivity of passing variants called for these samples concordant with GIAB v4.2.1 NA12878 high-confidence variants. Marker shape indicates values per variant caller for SNVs only (“X”), indels only (“+”), or SNVs and indels together (filled circles).
Figure 3:
Figure 3:
1000 Genome comparison: Sensitivity and precision of xAtlas variant calls compared against SNVs and indels from a high-coverage Illumina phased panel generated on 3,202 samples from the 1000 Genomes Project, using (A) all passing variants on chromosomes 1–22 and X and (B) passing variants within GIAB HG002 high-confidence variant regions on chromosomes 1–22. PCA plots of chromosome 22 SNVs with MAF ≥ 0.05 called by xAtlas on these 3,202 samples, showing (C) PC1 vs. PC2, (D) PC1 vs. PC3, and (E) PC2 vs. PC3.
See this image and copyright information in PMC

Similar articles

See all similar articles

Cited by

See all "Cited by" articles

References

    1. Wagner J, Olson ND, Harris L, et al. Benchmarking challenging small variants with linked and long reads. Cell Genomics. 2022;2(5):100128. - PMC - PubMed
    1. Olson ND, Wagner J, McDaniel J, et al. PrecisionFDA Truth Challenge V2: calling variants from short and long reads in difficult-to-map regions. Cell Genomics. 2022;2(5):100129. - PMC - PubMed
    1. Nielsen R, Paul JS, Albrechtsen A, et al. Genotype and SNP calling from next-generation sequencing data. Nat Rev Genet. 2011;12:443–51. - PMC - PubMed
    1. Altman RB, Prabhu S, Sidow A, et al. A research roadmap for next-generation sequencing informatics. Sci Transl Med. 2016;8:335ps. - PMC - PubMed
    1. Koboldt DC. Best practices for variant calling in clinical sequencing. Genome Med. 2020;12:91. - PMC - PubMed

Publication types