Functional equivalence of genome sequencing analysis pipelines enables harmonized variant calling across human genetics projects
Authors:
- Allison A RegierARAllison A Regier
- Yossi FarjounYFYossi Farjoun
- David LarsonDLDavid Larson
- Olga KrashenininaOKOlga Krasheninina
- Hyun Min KangHKHyun Min Kang
- Daniel P. HowriganDHDaniel P. Howrigan
- Bo-Juen ChenBCBo-Juen Chen
- Manisha KherMKManisha Kher
- Eric BanksEBEric Banks
- Darren C. AmesDADarren C. Ames
- Adam C. EnglishAEAdam C. English
- Heng LiHLHeng Li
Preprint
Abstract
Hundreds of thousands of human whole genome sequencing (WGS) datasets will be generated over the next few years to interrogate a broad range of traits, across diverse populations. These data are more valuable in aggregate: joint analysis of genomes from many sources increases sample size and statistical power for trait mapping, and will enable studies of genome biology, population genetics and genome function at unprecedented scale. A central challenge for joint analysis is that different WGS data processing and analysis pipelines cause substantial batch effects in combined datasets, necessitating computationally expensive reprocessing and harmonization prior to variant calling. This approach is no longer tenable given the scale of current studies and data volumes. Here, in a collaboration across multiple genome centers and NIH programs, we define WGS data processing standards that allow different groups to produce “functionally equivalent” (FE) results suitable for joint variant calling with minimal batch effects. Our approach promotes broad harmonization of upstream data processing steps, while allowing for diverse variant callers. Importantly, it allows each group to continue innovating on data processing pipelines, as long as results remain compatible. We present initial FE pipelines developed at five genome centers and show that they yield similar variant calling results – including single nucleotide (SNV), insertion/deletion (indel) and structural variation (SV) – and produce significantly less variability than sequencing replicates. Residual inter-pipeline variability is concentrated at low quality sites and repetitive genomic regions prone to stochastic effects. This work alleviates a key technical bottleneck for genome aggregation and helps lay the foundation for broad data sharing and community-wide “big-data” human genetics studies.Keywords
All Related Versions
- Published version: Nature Communications, 9 (1), 1.
This publication has 14 references indexed in Scilit:
- A reference panel of 64,976 haplotypes for genotype imputationNature Genetics, 2016
- Analysis of protein-coding genetic variation in 60,706 humansNature, 2016
- A global reference for human genetic variationNature, 2015
- SpeedSeq: ultra-fast personal genome analysis and interpretationNature Methods, 2015
- An efficient and scalable analysis framework for variant extraction and refinement from population-scale DNA sequence dataGenome Research, 2015
- Extending reference assembly modelsGenome Biology, 2015
- SAMBLASTER: fast duplicate marking and structural variant read extractionBioinformatics, 2014
- An integrated map of genetic variation from 1,092 human genomesNature, 2012
- Efficient storage of high throughput DNA sequencing data using reference-based compressionGenome Research, 2011
- The Sequence Alignment/Map format and SAMtoolsBioinformatics, 2009