Volume 14 Supplement 2
Selected articles from the Eleventh Asia Pacific Bioinformatics Conference (APBC 2013): Bioinformatics
CoNVEX: copy number variation estimation in exome sequencing data using HMM
 Kaushalya C Amarasinghe^{1}Email author,
 Jason Li^{2} and
 Saman K Halgamuge^{1}
DOI: 10.1186/1471210514S2S2
© Amarasinghe et al.; licensee BioMed Central Ltd. 2013
Published: 21 January 2013
Abstract
Background
One of the main types of genetic variations in cancer is Copy Number Variations (CNV). Whole exome sequenicng (WES) is a popular alternative to whole genome sequencing (WGS) to study disease specific genomic variations. However, finding CNV in Cancer samples using WES data has not been fully explored.
Results
We present a new method, called CoNVEX, to estimate copy number variation in whole exome sequencing data. It uses ratio of tumour and matched normal average read depths at each exonic region, to predict the copy gain or loss. The useful signal produced by WES data will be hindered by the intrinsic noise present in the data itself. This limits its capacity to be used as a highly reliable CNV detection source. Here, we propose a method that consists of discrete wavelet transform (DWT) to reduce noise. The identification of copy number gains/losses of each targeted region is performed by a Hidden Markov Model (HMM).
Conclusion
HMM is frequently used to identify CNV in data produced by various technologies including Array Comparative Genomic Hybridization (aCGH) and WGS. Here, we propose an HMM to detect CNV in cancer exome data. We used modified data from 1000 Genomes project to evaluate the performance of the proposed method. Using these data we have shown that CoNVEX outperforms the existing methods significantly in terms of precision. Overall, CoNVEX achieved a sensitivity of more than 92% and a precision of more than 50%.
Keywords
CNV detection Cancer Genome Targeted resequencing Whole exome sequencing Hidden Markov Models Discrete Wavelet TransformBackground
Commercial products of Next Generation Sequencing (NGS) Technologies such as Roche/454 FLX, Illumina Genome Analyzer/HiSeq, Applied Biosystems SOLiD™System and Helicos Heliscope™have enabled the sequencing of DNA much faster and cheaper than before [1]. These have shifted the paradigm of biological sequence analysis to a new level. Currently these are not only being used for the sequencing of whole genome, but also for sequencing of known exons and transcriptomes as well. The main motivations behind the technology of targeted resequencing (TR) include the following among others. The actual coding regions or the exons of the human genome account only for ~1% of the total sequences, which consequently gives about 30 Mb data compared to 3 Gb data in WGS [2]. Currently, getting higher coverage of targeted regions using NGS technologies is about six times [3] cheaper and faster compared to achieving the same coverage of whole genome. On the other hand approximately 85% of disease causing mutations lie in the coding regions [4]. Targeted resequencing has been mainly used in medical sequencing to find disease causing genetic variations (a review can be found in [5]). Recent studies on TR and WES data have successfully detected cancer specific mutations (somatic mutations) in breast cancer [6–8], ovarian cancer [8] and prostate cancer [9]. Although, exome sequencing has been successfully used to find small variations in cancer genomes, its potential to find large structural variations such as CNV has not yet been fully explored.
Cancer arises due to the acquisition of many somatic variations by the DNA of cancer cells [10]. Copy Number Alterations (CNA) play a major part in the progression of this deadly disease [11]. Until recently, the most common method to detect CNV in cancer DNA was to use micro array based technologies. However, during last 2  3 years many algorithms have been developed to identify CNV in cancerous data generated by whole genome sequencing [11–15], making use of the vast amount of data produced by NGS technology. The higher resolution that can be achieved through NGS data has helped to detect new variations that were undetectable previously and include CNVs which are as small as 50 bp [16]. These methods use the number of reads mapped to a particular region in the genome, to find copy number varying regions in one genome compared to one or more other genomes. Some of these methods have been adapted from the methods used in aCGH. For example Circular Binary Segmentation (CBS) [17] and Hidden Markov Model (HMM) [18]. However, methods in whole genome sequencing cannot be directly applied to whole exome sequencing data due to the small size and sparseness of these data [19]. On the other hand, the useful signal will be hindered by the intrinsic noise present in the exome sequencing data itself due to various biases introduced in target capturing and sequencing phases. To address these issues and to utilize the advantages provided by targeted resequencing, new algorithms have to be developed. Since the end of 2011, very few bioinformatics methods for detecting copy number variations in targeted resequencing data have been published. The method in [20] describes the use of TR data to detect CNV in cancer samples. However, the targeted regions in this method are larger in size (~ 40 kb), where as exons are much shorter (200 bp 300 bp). Methods such as [21, 22] are developed to find CNV in non cancerous exome data, such as in population studies. CONTRA [19] is a recent method proposed to evaluate cancer TR data using a pooled or a matched normal sample. ExomeCNV [23] and Var Scan 2 [24] are specifically designed for whole exome sequencing of cancer samples. A limitation in these approaches is that they have a higher number of false positives which result in a very low precision (further discussed in Results and discussion section for ExomeCNV).
In this work, we present CoNVEX, a method that evaluates exon level depth of coverage ratios to assess variation in copy number of whole exome capture data produced from cancer samples. We propose to use Discrete Wavelet Transformation denoising to reduce the variability of coverage ratios and then use HMM to detect copy number variations. Our method reduces the number of false positives by efficient preprocessing of the data, which results in a mean precision of more than 50%.
Methods
Data preprocessing
Depth of coverage ratios at each targeted region
Where ${N}_{{T}_{i}}$ and ${N}_{{C}_{i}}$are the mean normalized DOC of tumour and control respectively.
DWT smoothing of the data
Here, ${\overline{R}}_{i}$ is the true signal of copy number variation with additive noise, ϵ_{ i }. This noise can be assumed to be iid with N(0, σ^{2})where σ is the standard deviation of the distribution. We have used DWT smoothing [26] on R_{ i }, to detect true signal ${\overline{R}}_{i}$ to increases the ability of actual copy number prediction. The DWT smoothing procedure starts by first taking discrete wavelet transformation of ratios using "HAAR" wavelet. The fundamental assumption behind discrete wavelet transform is that, there is a correlation between the two neighbouring samples or data points. This is very much true in predicting CNVs as they span multiple successive exons. The selection of HAAR wavelet family was based on the fact that it computes the wavelet coefficients as the difference between two near by blocks of data points. This feature helps to retain the information regarding copy number aberration points. The shrinking of the DWT coefficients were done using soft thresholding function and the threshold value was calculated by Stein's unbiased risk estimator (SURE) for each level of DWT. Finally, the modified coefficients were used to reconstruct the denoised signal at i^{ th }location of chromosome j, ${\overline{R}}_{ij}$, by taking the inverse transform.
CNV prediction using a Hidden Markov Model
 1.
The total number of hidden states in the model is given by K and those are denoted by S = S_{1}, S_{2},..., S_{ K }. If there are L exons in the sample of consideration, the state of l^{ th }exon (e_{ l }) equals to S_{ k }where 1 ≤ l ≤ L and 1 ≤ k ≤ K.
 2.The initial state distribution π = {π_{ k }} where${\pi}_{k}=P\left({e}_{1}={S}_{k}\right),\phantom{\rule{1em}{0ex}}1\le k\le K$
 3.The state transition probability distribution A = a_{ mp }where${a}_{mp}=P\left({e}_{l+1}={S}_{p}{e}_{l}={S}_{m}\right),\phantom{\rule{1em}{0ex}}1\le m,p\le K$
 4.The emission probability distribution is given by B = {b_{ k }(O)} where$\left\{{b}_{k}\left(\mathbf{O}\right)\right\}=\mathcal{N}\left({\mathbf{O}}_{l},{\mu}_{k},{\sigma}^{2}\right),\phantom{\rule{1em}{0ex}}1\le l\le L\phantom{\rule{0.3em}{0ex}}\phantom{\rule{0.3em}{0ex}}and\phantom{\rule{0.3em}{0ex}}\phantom{\rule{0.3em}{0ex}}1\le k\le K$
Here, $\mathcal{N}$ represents the Gaussian distribution. Mean (μ_{ k }) of that distribution vary with different states and the provided normal cell contamination percentage and ploidy. We used a common standard deviation, σ, to all states.
The above HMM can be represented compactly as λ = (A, B, π) where A, B and π represent transition probability matrix, emission probability distribution and initial state distribution. When fitting the above HMM, the K states must be fixed at first and normal contamination and tumour ploidy must be given as inputs.
The optimal λ is selected by optimizing the negative loglikelihood [27, 28]. The initial state distribution π is chosen such that higher probability is attached to the most abundantly expected state or the normal state (i.e. copy 2 in normal humans). Similarly, the transition probability matrix A, is chosen such that, a higher probability is assigned to remain in the same state and lower probability is assigned to transition to another state. Also the transition to normal state has higher probability than transition to a CNV state. Then we used Viterbi Algorithm to assign the most appropriate copy number state for each exon.
Relationship between DOC ratio and copy number
where P_{ T } is the most abundant ploidy in the tumour sample.
Data from 1000 Genome Project
We randomly selected six samples, NA18536, NA18543, NA18544, NA18548, NA18557, NA18558, from 1000 Genome project, which share some common attributes, to evaluate the performance of the proposed method. These selected individuals have been studied by the HapMap project http://www.hapmap.org. The common features in these individuals are (i) exome sequencing was performed by the Beijing Genome Institute, hence a common exomecapture (NimbleGen V2) has been performed, (ii) male individuals and (iii) from CHB population.
Simulated data with known copy number variations
We used depth of coverage data at each exon of 1000 Genome samples to simulate CNV. This ensures that we retain as much intrinsic noise present in non copy number varying regions. The simulation procedure is as follows,
1. First, we retain only the copy number neutral regions in each sample. The CNV information were downloaded from the HapMap project genotype file.
2. We selected one sample (NA18536) as the Control and others as the tumour with known CNV.
3. To do the simulation of gains and losses we randomly selected a region in the Chr 1 and reduce (e.g: multiplied by 0.05) or amplified (e.g: multiplied by 2) the number of reads in that particular region. For each variation type, we perform 100 simulations.
4. When we evaluated the performance using only one sample (NA18543), we used 100 simulations for each variation type. When we used 5 samples in simulations, 20 variations were simulated in each individual sample.
5. To incorporate contamination in the simulation, we mix the control sample and simulated sample as per the relationship (α∗ Control + (1  α) ∗ Tumor) where α is the contamination proportion.
Results and discussion
Exon level depth of coverage ratios to detect CNV in whole exome data
Different methods have been proposed to reduce the experimental biases present in TR data. These include GC content bias reduction using regression methods [20, 21], taking base level ratios between normal and control samples [19, 23] and bait capture bias reduction using log transformation [20]. Those methods, adapted from aCGH or whole genome sequencing based approaches, try to reduce different experimental biases separately. Hsu et. al. [29] proposed DWT smoothing as an effective method to extract true copy number variations from aCGH data.
In this work, we propose to combine the strengths of both DWT and HMM to robustly predict copy number variations in cancer samples. The main novelty of our approach is the use of DWT smoothing to reduce experimental biases present in whole exome sequencing data prior to applying a Hidden Markov Model. These experimental biases are modelled here as additive noise to the true signal. The wavelet coefficients, which are the differences between two nearby data blocks, can be used to reduce noise. This is achieved through approximating some coefficients that do not by pass a certain threshold to zero. After thresholding step when the inverse transform is performed on these wavelet coefficients, we can generate a smoother version of the input signal. Exon level ratios, before and after DWT smoothing, for data downloaded from 1000 Genome project (http://www.1000genomes.org) are given in Figure 1.
After smoothing, we applied an HMM described in Methods section to detect copy gains and losses. Hidden Markov Models have been previously used to detect CNV in exome data (an R package called ExomeCopy) [21], but not used in this manner to detect CNV in tumour samples. The differences between ExomeCopy and CoNVEX are,

ExomeCopy uses HMM to identify CNVs in male patients with Xlinked Intellectual Disabilities (XLID)

They have used depth of coverage of exons as observations or emissions of hidden states

The robustness in copy identification is achieved by pooling coverage data from all patients
Therefore, it fails to identify relative copy number in cancer samples against a matched normal.
Comparison of the performance of CoNVEX against other methods
Comparison against ExomeCNV using simulated data
We carried out a comparison between the proposed method and the existing method, ExomeCNV [23]. Using simulated data, we were able to assess the performance of CoNVEX and ExomeCNV for different size ranges.
A true positive (TP) is identified when the gain or loss of an exon is correctly identified by the algorithm and a false positive (FP) identification is defined in the same manner. When using ExomeCNV, we used their primary CNV detection method (here after referred to as ExomeCNV1) and the extension which combines DNACopy [17] (here after referred to as ExomeCNV2) separately on our simulated data sets. The DNACopy version of ExomeCNV is applied to make sure that we get results for all exons that pass the default cutoff level of the coverage (this is a direction given by the authors of the paper). We used default parameters given in ExomeCNV R package for CNV prediction, except for read length and admixture rate, which we set to 90 and 0.0 in our evaluation.
Performance of proposed method for 100 simulations.
Type  Proposed Method  

Sensitivity  Specificity  Precision  Accuracy  
Deletions (1 k 1 M bp)  97.82 ± 12.37%  99.94 ± 0.081%  79.25 ± 23.23%  99.94 ± 0.081% 
Duplications (1 k 1 M bp)  95.25 ± 19.64%  99.93 ± 0.082%  77.04 ± 26.43%  99.93 ± 0.085% 
Performance of ExomeCNV1 for 100 simulations.
Type  ExomeCNV1  

Sensitivity  Specificity  Precision  Accuracy  
Deletions (1 k  1 M bp)  97.91 ± 2.81%  86.20 ± 1.57%  8.76 ± 6.54%  86.24 ± 1.56% 
Duplications (1 k  1 M bp)  90.68 ± 9.02%  86.26 ± 1.55%  8.96 ± 8.57%  86.28 ± 1.54% 
Performance of ExomeCNV2 for 100 simulations.
Type  ExomeCNV2  

Sensitivity  Specificity  Precision  Accuracy  
Deletions (1 k  1 M bp)  99.26 ± 2.11%  96.00 ± 1.67%  8.69 ± 6.50%  96.01 ± 1.66% 
Duplications (1 k  1 M bp)  99.98 ± 0.16%  96.06 ± 1.65%  9.62 ± 9.25%  96.08 ± 1.64% 
When compared with ExomeCNV2, our method showed superior performance in terms of specificity, precision and accuracy. Slight decrease in sensitivity was observed in CoNVEX, this is mainly due to the detecting short variations involving 1 or 2 exons. This can be attributed to the smoothing step we performed using DWT. Because of this we separately tested the performance of CoNVEX for shorter variations sizes as described below. Both versions of ExomeCNV, showed very poor performance when it comes to precision, as it tries to detect as many as possible variations to maintain a higher sensitivity rate.
Performance assessment of other methods against CoNVEX
Performance of CoNVEX against other methods.
Method  True positives  False positives 

CoNVEX  9/10  10/15850 
Var Scan2  6/7  4983/15283 
ExomeCopy  0/10  9/15850 
CONTRA  0/10  0/15847 
ExomeCopy and CONTRA did not identify any of the variations present in the test sample. This can be attributed to the fact that these are specifically designed for using a background sample [21] or a robust baseline [19]. VarScan 2 was able to identify the hemizygous duplication in the region with 60% sensitivity, however the number of false positives reported by the method was very high (false positive rate of 32%). CoNVEX performed well with 90% sensitivity and 0.06% false positive rate.
Performance of proposed method at different duplication and deletion sizes
Median sensitivity of CoNVEX for small variation detection is 100%. Every deletion of size, more than 200 bp was detected by our method. Hence, giving a mean sensitivity of 100% for detecting deletions. Mean sensitivity of detecting each duplication size was more than 85%. As seen in the graph, almost every variation of size of more than 800 bases can be detected by the proposed method. Also, a median precision of more than 30% can be achieved.
Performance assessment at different levels of contamination
Conclusions
Exome sequencing data can be used to detect copy number variations as an initial screening procedure. It is a cheap and time efficient method. We have successfully applied the proposed method on exome data to identify CNVs spanning one to thousands of exons. However, actual breakpoint of the CNV would not necessarily lie in the coding region. This limits the use of WES in identifying actual breakpoints of the CNV.
As discussed in the Results and Discussion section, we have achieved a higher precision than existing methods in detecting variations due to the data smoothing step. However, detection of some of the small variations may be missed by this smoothing step, as these can be recognised as noise. Further analysis is needed in order to better detect these variations among higher level of noise.
Although, we have used a matched normal sample to detect CNVs, the CNV identification can be done based on a pooled normal sample as described in [19]. This might give an advantage in finding CNVs in familial studies assuming all members have a median copy number of two.
Notes
Declarations
Acknowledgements
We thank Dr. Isaam Saeed and Dr. Suhinthan Maheswararajah for initial discussions on HMM. We used resources from both University of Melbourne and Peter MacCallum Cancer Centre for data processing and analysis. This work is partially funded by Australian Research Council (grant DP1096296).
This article has been published as part of BMC Bioinformatics Volume 14 Supplement 2, 2013: Selected articles from the Eleventh Asia Pacific Bioinformatics Conference (APBC 2013): Bioinformatics. The full contents of the supplement are available online at http://www.biomedcentral.com/bmcbioinformatics/supplements/14/S2.
Declarations
The funding for open access charges were provided by The University of Melbourne.
Authors’ Affiliations
References
 Holt RA, Jones SJM: The new paradigm of flow cell sequencing. Genome Research. 2008, 18 (6): 839846. 10.1101/gr.073262.107.View ArticlePubMedGoogle Scholar
 Ng SB, Turner EH, Robertson PD, Flygare SD, Bigham AW, Lee C, Shaffer T, Wong M, Bhattacharjee A, Eichler EE, Bamshad M, Nickerson DA, Shendure J: Targeted capture and massively parallel sequencing of 12 human exomes. NATURE. 2009, 461 (7261): 272U153. 10.1038/nature08250.PubMed CentralView ArticlePubMedGoogle Scholar
 Biesecker LG, Shianna KV, Mullikin JC: Exome sequencing: the expert view. GENOME BIOLOGY. 2011, 12 (9, SI):
 Choi M, Scholl UI, Ji W, Liu T, Tikhonova IR, Zumbo P, Nayir A, Bakkaloglu A, Ozen S, Sanjad S, NelsonWilliams C, Farhi A, Mane S, Lifton RP: Genetic diagnosis by whole exome capture and massively parallel DNA sequencing. Proceedings of the National Academy of Sciences of the United States of America. 2009, 106 (45): 1909619101. 10.1073/pnas.0910672106.PubMed CentralView ArticlePubMedGoogle Scholar
 Teer JK, Mullikin JC: Exome sequencing: the sweet spot before whole genomes. Human Molecular Genetics. 2010, 19 (R2): R145R151. 10.1093/hmg/ddq333. [http://hmg.oxfordjournals.org/content/19/R2/R145.abstract]PubMed CentralView ArticlePubMedGoogle Scholar
 Stephens PJ, Tarpey PS, Davies H, Van Loo P, Greenman C, Wedge DC, Zainal SN, Martin S, Varela I, Bignell GR, Yates LR, Papaemmanuil E, Beare D, Butler A, Cheverton A, Gamble J, Hinton J, Jia M, Jayakumar A, Jones D, Latimer C, Lau KW, McLaren S, McBride DJ, Menzies A, Mudie L, Raine K, Rad R, Spencer Chapman M, Teague J, Easton D, Langerod A, Lee MTM, Shen CY, Tee BTK, Huimin BW, Broeks A, Vargas AC, Turashvili G, Martens J, Fatima A, Miron P, Chin SF, Thomas G, Boyault S, Mariani O, Lakhani SR, van de Vijver M, van/'t Veer L, Foekens J, Desmedt C, Sotiriou C, Tutt A, Caldas C, ReisFilho JS, Aparicio SAJR, Salomon AV, BorresenDale AL, Richardson A, Campbell PJ, Futreal PA, Stratton MR: The landscape of cancer genes and mutational processes in breast cancer. Nature. 2012, advance online publication:, http://dx.doi.org/10.1038/nature11017Google Scholar
 Banerji S, Cibulskis K, RangelEscareno C, Brown KK, Carter SL, Frederick AM, Lawrence MS, Sivachenko AY, Sougnez C, Zou L, Cortes ML, FernandezLopez JC, Peng S, Ardlie KG, Auclair D, BautistaPina V, Duke F, Francis J, Jung J, MaffuzAziz A, Onofrio RC, Parkin M, Pho NH, QuintanarJurado V, Ramos AH, RebollarVega R, RodriguezCuevas S, RomeroCordoba SL, Schumacher SE, Stransky N, Thompson KM, UribeFigueroa L, Baselga J, Beroukhim R, Polyak K, Sgroi DC, Richardson AL, JimenezSanchez G, Lander ES, Gabriel SB, Garraway LA, Golub TR, MelendezZajgla J, Toker A, Getz G, HidalgoMiranda A, Meyerson M: Sequence analysis of mutations and translocations across breast cancer subtypes. Nature. 2012, 486 (7403): 405409. 10.1038/nature11154. [http://dx.doi.org/10.1038/nature11154]PubMed CentralView ArticlePubMedGoogle Scholar
 Walsh T, Lee MK, Casadei S, Thornton AM, Stray SM, Pennil C, Nord AS, Mandell JB, Swisher EM, King MC: Detection of inherited mutations for breast and ovarian cancer using genomic capture and massively parallel sequencing. PROCEEDINGS OF THE NATIONAL ACADEMY OF SCIENCES OF THE UNITED STATES OF AMERICA. 2010, 107 (28): 1262912633. 10.1073/pnas.1007983107.PubMed CentralView ArticlePubMedGoogle Scholar
 Kumar A, White TA, MacKenzie AP, Clegg N, Lee C, Dumpit RF, Coleman I, Ng SB, Salipante SJ, Rieder MJ, Nickerson DA, Corey E, Lange PH, Morrissey C, Vessella RL, Nelson PS, Shendure J: Exome sequencing identifies a spectrum of mutation frequencies in advanced and lethal prostate cancers. Proceedings of the National Academy of Sciences. 2011, 108 (41): 1708717092. 10.1073/pnas.1108745108. [http://www.pnas.org/content/108/41/17087.abstract]View ArticleGoogle Scholar
 Stratton MR, Campbell PJ, Futreal PA: The cancer genome. NATURE. 2009, 458 (7239): 719724. 10.1038/nature07943.PubMed CentralView ArticlePubMedGoogle Scholar
 Ivakhno S, Royce T, Cox AJ, Evers DJ, Cheetham RK, Tavare S: CNAsega novel framework for identification of copy number changes in cancer from secondgeneration sequencing data. Bioinformatics. 2010, 26 (24): 30513058. 10.1093/bioinformatics/btq587. [http://bioinformatics.oxfordjournals.org/content/26/24/3051.abstract]View ArticlePubMedGoogle Scholar
 Campbell PJ, Stephens PJ, Pleasance ED, O'Meara S, Li H, Santarius T, Stebbings LA, Leroy C, Edkins S, Hardy C, Teague JW, Menzies A, Goodhead I, Turner DJ, Clee CM, Quail MA, Cox A, Brown C, Durbin R, Hurles ME, Edwards PAW, Bignell GR, Stratton MR, Futreal PA: Identification of somatically acquired rearrangements in cancer using genomewide massively parallel pairedend sequencing. Nat Genet. 2008, 40 (6): 722729. 10.1038/ng.128. [http://dx.doi.org/10.1038/ng.128]PubMed CentralView ArticlePubMedGoogle Scholar
 Yoon S, Xuan Z, Makarov V, Ye K, Sebat J: Sensitive and accurate detection of copy number variants using read depth of coverage. GENOME RESEARCH. 2009, 19 (9): 15861592. 10.1101/gr.092981.109.PubMed CentralView ArticlePubMedGoogle Scholar
 Xie C, Tammi M: CNVseq, a new method to detect copy number variation using highthroughput sequencing. BMC Bioinformatics. 2009, 10: 8010.1186/147121051080. [http://www.biomedcentral.com/14712105/10/80]PubMed CentralView ArticlePubMedGoogle Scholar
 Boeva V, Popova T, Bleakley K, Chiche P, Cappo J, Schleiermacher G, JanoueixLerosey I, Delattre O, Barillot E: ControlFREEC: a tool for assessing copy number and allelic content using next generation sequencing data. Bioinformatics. 2011, [http://bioinformatics.oxfordjournals.org/content/early/2011/12/05/bioinformatics.btr670.abstract]Google Scholar
 Alkan C, Coe BP, Eichler EE: APPLICATIONS OF NEXTGENERATION SEQUENCING Genome structural variation discovery and genotyping. NATURE REVIEWS GENETICS. 2011, 12 (5): 363375. 10.1038/nrg2958.PubMed CentralView ArticlePubMedGoogle Scholar
 Olshen AB, Venkatraman ES, Lucito R, Wigler M: Circular binary segmentation for the analysis of arraybased DNA copy number data. Biostatistics. 2004, 5 (4): 557572. 10.1093/biostatistics/kxh008. [http://biostatistics.oxfordjournals.org/content/5/4/557.abstract]View ArticlePubMedGoogle Scholar
 Fridlyand J, Snijders AM, Pinkel D, Albertson DG, Jain AN: Hidden Markov models approach to the analysis of array CGH data. Journal of Multivariate Analysis. 2004, 90: 132153. 10.1016/j.jmva.2004.02.008. [¡ce:title¿Special Issue on Multivariate Methods in Genomic Data Analysis¡/ce:title¿], [http://www.sciencedirect.com/science/article/pii/S0047259X04000260]View ArticleGoogle Scholar
 Li J, Lupat R, Amarasinghe KC, Thompson ER, Doyle MA, Ryland GL, Tothill RW, Halgamuge SK, Campbell IG, Gorringe KL: CONTRA: copy number analysis for targeted resequencing. Bioinformatics. 2012, 28 (10): 13071313. 10.1093/bioinformatics/bts146. [http://bioinformatics.oxfordjournals.org/content/28/10/1307.abstract]PubMed CentralView ArticlePubMedGoogle Scholar
 Nord AS, Lee M, King MC, Walsh T: Accurate and exact CNV identification from targeted highthroughput sequence data. BMC GENOMICS. 2011, 12:Google Scholar
 Sun R, Kalscheuer V, Vingron M, Haas SA: Modeling Read Counts for CNV Detection in Exome Sequencing Data. Statistical Applications in Genetics and Molecular Biology. 2011, 10 (52): Love A Michael Iand Mysickov Ã¡, [http://www.bepress.com/sagmb/vol10/iss1/art52]Google Scholar
 Krumm N, Sudmant PH, Ko A, O'Roak BJ, Malig M, Coe BP, NHLBI Exome Sequencing Project N, Quinlan AR, Nickerson DA, Eichler EE: Copy number variation detection and genotyping from exome sequence data. Genome Research. 2012, [http://genome.cshlp.org/content/early/2012/05/14/gr.138115.112.abstract]Google Scholar
 Sathirapongsasuti JF, Lee H, Horst BAJ, Brunner G, Cochran AJ, Binder S, Quackenbush J, Nelson SF: Exome sequencingbased copynumber variation and loss of heterozygosity detection: ExomeCNV. Bioinformatics. 2011, 27 (19): 26482654. 10.1093/bioinformatics/btr462. [http://bioinformatics.oxfordjournals.org/content/27/19/2648.abstract]PubMed CentralView ArticlePubMedGoogle Scholar
 Koboldt DC, Zhang Q, Larson DE, Shen D, McLellan MD, Lin L, Miller CA, Mardis ER, Ding L, Wilson RK: VarScan 2: Somatic mutation and copy number alteration discovery in cancer by exome sequencing. Genome Research. 2012, 22 (3): 568576. 10.1101/gr.129684.111. [http://genome.cshlp.org/content/22/3/568.abstract]PubMed CentralView ArticlePubMedGoogle Scholar
 Quinlan AR, Hall IM: BEDTools: a flexible suite of utilities for comparing genomic features. BIOINFORMATICS. 2010, 26 (6): 841842. 10.1093/bioinformatics/btq033.PubMed CentralView ArticlePubMedGoogle Scholar
 Percival DB, Walden AT: Wavelet Methods for Time Series Analysis (Cambridge Series in Statistical and Probabilistic Mathematics). 2006, Cambridge University Press, [http://www.worldcat.org/isbn/0521685087g]Google Scholar
 Rabiner L: A tutorial on hidden Markov models and selected applications in speech recognition. Proceedings of the IEEE. 1989, 77 (2): 257286. 10.1109/5.18626.View ArticleGoogle Scholar
 Zucchini W, MacDonald IL: Hidden Markov models for time series: an introduction using R. 2009View ArticleGoogle Scholar
 Hsu L, Self SG, Grove D, Randolph T, Wang K, Delrow JJ, Loo L, Porter P: Denoising arraybased comparative genomic hybridization data using wavelets. Biostatistics. 2005, 6 (2): 211226. 10.1093/biostatistics/kxi004. [http://biostatistics.oxfordjournals.org/content/6/2/211.abstract]View ArticlePubMedGoogle Scholar
Copyright
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.