Display options
Share it on

Ann Appl Stat. 2010 Dec 01;4(4):1749-1773. doi: 10.1214/10-AOAS357.

RECONSTRUCTING DNA COPY NUMBER BY PENALIZED ESTIMATION AND IMPUTATION.

The annals of applied statistics

Zhongyang Zhang, Kenneth Lange, Roel Ophoff, Chiara Sabatti

Affiliations

  1. Department of Statistics University of California, Los Angeles Los Angeles, California 90095 USA [email protected].

PMID: 21572975 PMCID: PMC3092301 DOI: 10.1214/10-AOAS357

Abstract

Recent advances in genomics have underscored the surprising ubiquity of DNA copy number variation (CNV). Fortunately, modern genotyping platforms also detect CNVs with fairly high reliability. Hidden Markov models and algorithms have played a dominant role in the interpretation of CNV data. Here we explore CNV reconstruction via estimation with a fused-lasso penalty as suggested by Tibshirani and Wang [Biostatistics 9 (2008) 18-29]. We mount a fresh attack on this difficult optimization problem by the following: (a) changing the penalty terms slightly by substituting a smooth approximation to the absolute value function, (b) designing and implementing a new MM (majorization-minimization) algorithm, and (c) applying a fast version of Newton's method to jointly update all model parameters. Together these changes enable us to minimize the fused-lasso criterion in a highly effective way.We also reframe the reconstruction problem in terms of imputation via discrete optimization. This approach is easier and more accurate than parameter estimation because it relies on the fact that only a handful of possible copy number states exist at each SNP. The dynamic programming framework has the added bonus of exploiting information that the current fused-lasso approach ignores. The accuracy of our imputations is comparable to that of hidden Markov models at a substantially lower computational cost.

References

  1. Nucleic Acids Res. 2007;35(6):2013-25 - PubMed
  2. Ann Appl Stat. 2008 Jun 1;2(2):687-713 - PubMed
  3. Science. 2004 Jul 23;305(5683):525-8 - PubMed
  4. Genome Res. 2007 Nov;17(11):1665-74 - PubMed
  5. Nucleic Acids Res. 2008 Nov;36(19):e126 - PubMed
  6. Nat Genet. 2004 Sep;36(9):949-51 - PubMed
  7. Am J Hum Genet. 2008 Oct;83(4):504-10 - PubMed
  8. Nature. 2008 Sep 11;455(7210):232-6 - PubMed
  9. Nature. 2008 Feb 21;451(7181):998-1003 - PubMed
  10. Hum Hered. 2009;68(1):1-22 - PubMed
  11. Nat Genet. 2008 Oct;40(10):1253-60 - PubMed
  12. Nature. 2006 Nov 23;444(7118):444-54 - PubMed
  13. Biostatistics. 2008 Jan;9(1):18-29 - PubMed
  14. Bioinformatics. 2007 Sep 15;23(18):2470-6 - PubMed

Publication Types

Grant support