A study of biases of DNA copy number estimation based on PICR model

Abstract

Affymetrix single-nucleotide polymorphism (SNP) arrays have been widely used for SNP genotype calling and copy number variation (CNV) studies, both of which are dependent on accurate DNA copy number estimation significantly. However, the methods for copy number estimation may suffer from kinds of difficulties: probe dependent binding affinity, crosshybridization of probes, and the whole genome amplification (WGA) of DNA sequences. The probe intensity composite representation (PICR) model, one former established approach, can cope with most complexities and achieve high accuracy in SNP genotyping. Nevertheless, the copy numbers estimated by PICR model still show array and site dependent biases for CNV studies. In this paper, we propose a procedure to adjust the biases and then make CNV inference based on both PICR model and our method. The comparison indicates that our correction of copy numbers is necessary for CNV studies.

Keywords

single-nucleotide polymorphism (SNP) array / copy number variation (CNV) / cross-hybridization

References

[1] Barnes C., Plagnol V., Fitzgerald T., Redon R., Marchini J., Clayton D., Hurles M. E. A robust statistical method for case-control association testing with copy number variation. Nat Genet, 2008, 40 10 1245-1252

[2] Bengtsson H., Irizarry R., Carvalho B., Speed T. P. Estimation and assessment of raw copy numbers at the single locus level. Bioinformatics, 2008, 24 6 759-767

[3] Bengtsson H., Wirapati P., Speed T. P. A single-array preprocessing method for estimating full-resolution raw copy numbers from all Affymetrix genotyping arrays including GenomeWideSNP 5 & 6. Bioinformatics, 2009, 25 17 2149-2156

[4] Bignell G. R., Huang J., Greshock J., Watt S., Butler A., West S., Grigorova M., Jones K. W., Wei W., Stratton M. R. High-resolution analysis of DNA copy number using oligonucleotide microarrays. Genome Res, 2004, 14 2 287-295

[5] Carter N. P. Methods and strategies for analyzing copy number variation using DNA microarrays. Nat Genet, 2007, 39 7Suppl S16-21

[6] Di X., Matsuzaki H., Webster T. A., Hubbell E., Liu G., Dong S., Bartell D., Huang J., Chiles R., Yang G. Dynamic model based algorithms for screening and genotyping over 100 K SNPs on oligonucleotide microarrays. Bioinformatics, 2005, 21 9 1958-1963

[7] Greenman C. D., Bignell G., Butler A., Edkins S., Hinton J., Beare D., Swamy S., Santarius T., Chen L., Widaa S., Futreal P. A., Stratton M. R. PICNIC: an algorithm to predict absolute allelic copy number variation with microarray cancer data. Biostatistics, 2010, 11 1 164-175

[8] Held G. A., Grinstein G., Tu Y. Modeling of DNA microarray data by using physical properties of hybridization. Proc Natl Acad Sci USA, 2003, 100 13 7575-7580

[9] Held G. A., Grinstein G., Tu Y. Relationship between gene expression and observed intensities in DNA microarrays-a modeling study. Nucleic Acids Res, 2006, 34 9 e70

[10] Huang J., Wei W., Chen J., Zhang J., Liu G., Di X., Mei R., Ishikawa S., Aburatani H., Jones K. W. CARAT: a novel method for allelic detection of DNA copy number changes using high density oligonucleotide arrays. BMC Bioinformatics, 2006, 7, 83

[11] Iafrate A. J., Feuk L., Rivera M.N., Listewnik M.L., Donahoe P. K., Qi Y., Scherer S.W., Lee C. Detection of large-scale variation in the human genome. Nat Genet, 2004, 36 9 949-951

[12] Johnson W. E., Li W., Meyer C. A., Gottardo R., Carroll J. S., Brown M., Liu X. S. Model-based analysis of tiling-arrays for ChIP-chip. Proc Natl Acad Sci USA, 2006, 103 33 12457-12462

[13] Kapur K., Jiang H., Xing Y., Wong W. H. Cross-hybridization modeling on Affymetrix exon arrays. Bioinformatics, 2008, 24 24 2887-2893

[14] Korn J. M., Kuruvilla F. G., McCarroll S. A., Wysoker A., Nemesh J., Cawley S., Hubbell E., Veitch J., Collins P. J., Darvishi K. Integrated genotype calling and association analysis of SNPs, common copy number polymorphisms and rare CNVs. Nat Genet, 2008, 40 10 1253-1260

[15] Laframboise T., Harrington D., Weir B. A. PLASQ: a generalized linear model-based procedure to determine allelic dosage in cancer cells from SNP array data. Biostatistics, 2007, 8 2 323-336

[16] McCarroll S. A., Kuruvilla F. G., Korn J. M., Cawley S., Nemesh J., Wysoker A., Shapero M. H., de Bakker P. I., Maller J. B., Kirby A. Integrated detection and populationgenetic analysis of SNPs and copy number variation. Nat Genet, 2008, 40 10 1166-1174

[17] Nannya Y., Sanada M., Nakazaki K., Hosoya N., Wang L., Hangaishi A., Kurokawa M., Chiba S., Bailey D. K., Kennedy G. C. A robust algorithm for copy number detection using high-density oligonucleotide single nucleotide polymorphism genotyping arrays. Cancer Res, 2005, 65 14 6071-6079

[18] Olshen A. B., Venkatraman E. S., Lucito R., Wigler M. Circular binary segmentation for the analysis of array-based DNA copy number data. Biostatistics, 2004, 5 4 557-572

[19] Ono N., Suzuki S., Furusawa C., Agata T., Kashiwagi A., Shimizu H., Yomo T. An improved physico-chemical model of hybridization on high-density oligonucleotide microarrays. Bioinformatics, 2008, 24 10 1278-1285

[20] Pugh T. J., Delaney A. D., Farnoud N., Flibotte S., Griffith M., Li H. I., Qian H., Farinha P., Gascoyne R. D., Marra M. A. Impact of whole genome amplification on analysis of copy number variants. Nucleic Acids Res, 2008, 36 13 e80

[21] Rabbee N., Speed T. P. A genotype calling algorithm for affymetrix SNP arrays. Bioinformatics, 2006, 22 1 7-12

[22] Redon R., Ishikawa S., Fitch K. R., Feuk L., Perry G. H., Andrews T. D., Fiegler H., Shapero M. H., Carson A. R., Chen W. Global variation in copy number in the human genome. Nature, 2006, 444 7118 444-454

[23] Scherer S. W., Lee C., Birney E., Altshuler D. M., Eichler E. E., Carter N. P., Hurles M. E., Feuk L. Challenges and standards in integrating surveys of structural variation. Nat Genet, 2007, 39 7Suppl S7-15

[24] Sebat J., Lakshmi B., Troge J., Alexander J., Young J., Lundin P., Maner S., Massa H., Walker M., Chi M. Y. Large-scale copy number polymorphism in the human genome. Science, 2004, 305 5683 525-528

[25] Slater H. R., Bailey D. K., Ren H., Cao M., Bell K., Nasioulas S., Henke R., Choo K. H., Kennedy G. C. High-resolution identification of chromosomal abnormalities using oligonucleotide arrays containing 116,204 SNPs. Am J Hum Genet, 2005, 77 5 709-726

[26] Wan L., Sun K., Ding Q., Cui Y., Li M., Wen Y., Elston R. C., Qian M., Fu W.J. Hybridization modeling of oligonucleotide SNP arrays for accurate DNA copy number estimation. Nucleic Acids Res, 2009, 37 17 e117

[27] Wan L., Xiao Y., Chen Q., Deng M., Qian M. The analysis of biases of copy numbers from Affymetrix SNP arrays. Communications in Information and Systems, 2010, 10 2 81-96

[28] Weir B. A., Woo M.S., Getz G., Perner S., Ding L., Beroukhim R., Lin W.M., Province M.A., Kraja A., Johnson L. A. Characterizing the cancer genome in lung adenocarcinoma. Nature, 2007, 450 7171 893-898

[29] Xiao Y., Segal M. R., Yang Y. H., Yeh R. F. A multi-array multi-SNP genotyping algorithm for Affymetrix SNP microarrays. Bioinformatics, 2007, 23 12 1459-1467

[30] Zhang L., Miles M. F., Aldape K. D. A model of molecular interactions on short oligonucleotide microarrays. Nat Biotechnol, 2003, 21 7 818-821

[31] Zhang L., Wu C., Carta R., Zhao H. Free energy of DNA duplex formation on short oligonucleotide microarrays. Nucleic Acids Res, 2007, 35 3 e18