Detecting copy number variations from array CGH data based on a conditional random field model.
Array comparative genomic hybridization (aCGH) allows identification of copy number alterations across genomes. The key computational challenge in analyzing copy number variations (CNVs) using aCGH data or other similar data generated by a variety of array technologies is the detection of segment boundaries of copy number changes and inference of the copy number state for each segment. We have developed a novel statistical model based on the framework of conditional random fields (CRFs) that can effectively combine data smoothing, segmentation and copy number state decoding into one unified framework. Our approach (termed CRF-CNV) provides great flexibilities in defining meaningful feature functions. Therefore, it can effectively integrate local spatial information of arbitrary sizes into the model. For model parameter estimations, we have adopted the conjugate gradient (CG) method for likelihood optimization and developed efficient forward/backward algorithms within the CG framework. The method is evaluated using real data with known copy numbers as well as simulated data with realistic assumptions, and compared with two popular publicly available programs. Experimental results have demonstrated that CRF-CNV outperforms a Bayesian Hidden Markov Model-based approach on both datasets in terms of copy number assignments. Comparing to a non-parametric approach, CRF-CNV has achieved much greater precision while maintaining the same level of recall on the real data, and their performance on the simulated data is comparable.
Yin
Xiao-Lin
XL
Electrical Engineering and Computer Science Department, Case Western Reserve University, Cleveland, Ohio 44106, United States. xly@case.edu
Li
Jing
J
D000465:Algorithms; D001499:Bayes Theorem; D055028:Comparative Genomic Hybridization / Q000706:statistics & numerical data*; D019295:Computational Biology; D056915:DNA Copy Number Variations; D030561:Databases, Nucleic Acid; D006801:Humans; D016013:Likelihood Functions; D008390:Markov Chains; D008957:Models, Genetic; D015233:Models, Statistical
Genome Res. 2006 Aug;16(8):949-61
16809666
Nucleic Acids Res. 2007;35(6):2013-25
17341461
Bioinformatics. 2006 Jul 15;22(14):e431-9
16873504
J Am Stat Assoc. 2008 Jun 1;103(482):485-497
22375091
Nat Genet. 2001 Nov;29(3):263-4
11687795
Biostatistics. 2004 Oct;5(4):557-72
15475419
Cytogenet Genome Res. 2006;115(3-4):262-72
17124409
Bioinformatics. 2005 Nov 15;21(22):4084-91
16159913
BMC Genet. 2008 Mar 28;9:27
18373861
Nature. 2006 Nov 23;444(7118):444-54
17122850
Bioinformatics. 2004 Dec 12;20(18):3413-22
15381628
Nat Methods. 2009 Jan;6(1):99-103
19043412
Nat Genet. 2007 Jul;39(7 Suppl):S16-21
17597776
Bioinformatics. 2005 Oct 1;21(19):3763-70
16081473
Genome Res. 2007 Nov;17(11):1665-74
17921354