MOTIVATION: Genomic instability in cancer leads to abnormal genome copy number alterations (CNA) that are associated with the development and behavior of tumors. Advances in microarray technology have allowed for greater resolution in detection of DNA copy number changes (amplifications or deletions) across the genome. However, the increase in number of measured signals and accompanying noise from the array probes present a challenge in accurate and fast identification of breakpoints that define CNA. This article proposes a novel detection technique that exploits the use of piece wise constant (PWC) vectors to represent genome copy number and sparse Bayesian learning (SBL) to detect CNA breakpoints. METHODS: First, a compact linear algebra representation for the genome copy number is developed from normalized probe intensities. Second, SBL is applied and optimized to infer locations where copy number changes occur. Third, a backward elimination (BE) procedure is used to rank the inferred breakpoints; and a cut-off point can be efficiently adjusted in this procedure to control for the false discovery rate (FDR). RESULTS: The performance of our algorithm is evaluated using simulated and real genome datasets and compared to other existing techniques. Our approach achieves the highest accuracy and lowest FDR while improving computational speed by several orders of magnitude. The proposed algorithm has been developed into a free standing software application (GADA, Genome Alteration Detection Algorithm). AVAILABILITY: http://biron.usc.edu/~piquereg/GADA

Alterations in the number of DNA copies are very common in tumor cells and may have a very important role in cancer development and progression. New array platforms provide means to analyze the copy number by comparing the hybridization intensities of thousands of DNA sections along the genome. However, detecting and locating the copy number changes from this data is a very challenging task due to the large amount of biological processes that affect hybridization and cannot be controlled. This paper proposes a new technique that exploits the key characteristic that the DNA copy number is piecewise-constant along the genome. First, wavelet footprints are used to obtain a basis for representing the DNA copy number that is maximally sparse in the number of copy number change points. Second, sparse Bayesian learning is applied to infer the copy number changes from noisy array probe intensities. Results demonstrate that sparse Bayesian learning has better performance than matching pursuits methods for this high coherence dictionary. Finally, our results are also shown to be very competitive in performance as compared to state-of-the-art methods for copy number detection.