Hybridization and amplification rate correction for affymetrix SNP arrays.
BACKGROUND
Copy number variation (CNV) is essential to understand the pathology of many complex diseases at the DNA level. Affymetrix SNP arrays, which are widely used for CNV studies, significantly depend on accurate copy number (CN) estimation. Nevertheless, CN estimation may be biased by several factors, including cross-hybridization and training sample batch, as well as genomic waves of intensities induced by sequence-dependent hybridization rate and amplification efficiency. Since many available algorithms only address one or two of the three factors, a high false discovery rate (FDR) often results when identifying CNV. Therefore, we have developed a new CNV detection pipeline which is based on hybridization and amplification rate correction (CNVhac).


METHODS
CNVhac first estimates the allelic concentrations (ACs) of target sequences by using the sample independent parameters trained through physicochemical hybridization law. Then the raw CN is estimated by taking the ratio of AC to the corresponding average AC from a reference sample set for one specific site. Finally, a hidden Markov model (HMM) segmentation process is implemented to detect CNV regions.


RESULTS
Based on public HapMap data, the results show that CNVhac effectively smoothes the genomic waves and facilitates more accurate raw CN estimates compared to other methods. Moreover, CNVhac alleviates, to a certain extent, the sample dependence of inference and makes CNV calling with appreciable low FDRs.


CONCLUSIONS
CNVhac is an effective approach to address the common difficulties in SNP array analysis, and the working principles of CNVhac can be easily extended to other platforms.
Center for Theoretical Biology, Peking University, Beijing 100871, People's Republic of China.
Quan
Q
Wang
Peichao
P
Peng
Minping
M
Qian
Lin
L
Wan
Minghua
M
Deng
D016477:Artifacts; D002138:Calibration; D056915:DNA Copy Number Variations / Q000235:genetics; D005260:Female; D060148:HapMap Project; D006801:Humans; D008297:Male; D021141:Nucleic Acid Amplification Techniques / Q000379:methods*; D009693:Nucleic Acid Hybridization / Q000379:methods*; D020411:Oligonucleotide Array Sequence Analysis / Q000379:methods*; D020641:Polymorphism, Single Nucleotide / Q000235:genetics*