Translational biomedical research is generating exponentially more data: thousands of whole-genome sequences (WGS) are now available; brain data are doubling every two years. Analyses of Big Data, including imaging, genomic, phenotypic, and clinical data, present qualitatively new challenges as well as opportunities. Among the challenges is a proliferation in ways analyses can fail, due largely to the increasing length and complexity of processing pipelines. Anomalies in input data, runtime resource exhaustion or node failure in a distributed computation can all cause pipeline hiccups that are not necessarily obvious in the output. Flaws that can taint results may persist undetected in complex pipelines, a danger amplified by the fact that research is often concurrent with the development of the software on which it depends. On the positive side, the huge sample sizes increase statistical power, which in turn can shed new insight and motivate innovative analytic approaches. We have developed a framework for Big Data Quality Control (BDQC) including an extensible set of heuristic and statistical analyses that identify deviations in data without regard to its meaning (domain-blind analyses). BDQC takes advantage of large sample sizes to classify the samples, estimate distributions and identify outliers. Such outliers may be symptoms of technology failure (e.g., truncated output of one step of a pipeline for a single genome) or may reveal unsuspected “ signal” in the data (e.g., evidence of aneuploidy in a genome). We have applied the framework to validate real-world WGS analysis pipelines. BDQC successfully identified data outliers representing various failure classes, including genome analyses missing a whole chromosome or part thereof, hidden among thousands of intermediary output files. These failures could then be resolved by reanalyzing the affected samples. BDQC both identified hidden flaws as well as yielded new insights into the data. BDQC is designed to complement quality software development practices. There are multiple benefits from the application of BDQC at all pipeline stages. By verifying input correctness, it can help avoid expensive computations on flawed data. Analysis of intermediary and final results facilitates recovery from aberrant termination of processes. All these computationally inexpensive verifications reduce cryptic analytical artifacts that could seriously preclude clinical-grade genome interpretation. BDQC is available at https://github.com/ini-bdds/bdqc.
[1]
Benjamin Haibe-Kains,et al.
BatchQC: interactive software for evaluating sample and batch effects in genomic data
,
2016,
Bioinform..
[2]
Jared C. Roach,et al.
Identification of copy number variants in whole-genome data using Reference Coverage Profiles
,
2015,
Front. Genet..
[3]
James A. Eddy,et al.
Human whole genome genotype and transcriptome data for Alzheimer’s and other neurodegenerative diseases
,
2016,
Scientific Data.
[4]
Jeffrey E. Barrick,et al.
Detecting rare structural variation in evolving microbial populations from new sequence junctions using breseq
,
2014,
Front. Genet..
[5]
Måns Magnusson,et al.
MultiQC: summarize analysis results for multiple tools and samples in a single report
,
2016,
Bioinform..
[6]
Christian Cole,et al.
AlmostSignificant: Simplifying quality control of high-throughput sequencing data
,
2016,
Bioinformatics.