Association analysis using next-generation sequence data from publicly available control groups: the robust variance score statistic

Andriy Derkach, Theodore Chiang, Jiafen Gong, Laura Addis, Sara Dobbins, Ian Tomlinson, Richard Houlston, Deb K Pal, Lisa J Strug

Research output: Contribution to journalArticlepeer-review

24 Citations (Scopus)


MOTIVATION: Sufficiently powered case-control studies with next-generation sequence (NGS) data remain prohibitively expensive for many investigators. If feasible, a more efficient strategy would be to include publicly available sequenced controls. However, these studies can be confounded by differences in sequencing platform; alignment, single nucleotide polymorphism and variant calling algorithms; read depth; and selection thresholds. Assuming one can match cases and controls on the basis of ethnicity and other potential confounding factors, and one has access to the aligned reads in both groups, we investigate the effect of systematic differences in read depth and selection threshold when comparing allele frequencies between cases and controls. We propose a novel likelihood-based method, the robust variance score (RVS), that substitutes genotype calls by their expected values given observed sequence data.

RESULTS: We show theoretically that the RVS eliminates read depth bias in the estimation of minor allele frequency. We also demonstrate that, using simulated and real NGS data, the RVS method controls Type I error and has comparable power to the 'gold standard' analysis with the true underlying genotypes for both common and rare variants.

AVAILABILITY AND IMPLEMENTATION: An RVS R script and instructions can be found at, and at

Original languageEnglish
Pages (from-to)2179-2188
Number of pages10
Issue number15
Early online date14 Apr 2014
Publication statusPublished - 1 Aug 2014


  • Algorithms
  • Analysis of Variance
  • Case-Control Studies
  • Child
  • Computational Biology
  • Control Groups
  • Data Interpretation, Statistical
  • Epilepsy, Rolandic
  • Gene Frequency
  • Genotype
  • High-Throughput Nucleotide Sequencing
  • Human Genome Project
  • Humans
  • Likelihood Functions
  • Polymorphism, Single Nucleotide


Dive into the research topics of 'Association analysis using next-generation sequence data from publicly available control groups: the robust variance score statistic'. Together they form a unique fingerprint.

Cite this