Bioinformatics Advance Access first published online on December 3, 2008
This version published online on December 4, 2008
Bioinformatics, doi:10.1093/bioinformatics/btn585
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Aneuploidy Prediction and Tumor Classification with Heterogeneous Hidden Conditional Random Fields
1Department of Computer Science, Princeton University, 35 Olden Street, Princeton, NJ 08540
2Lewis-Sigler Institute for Integrative Genomics, Carl Icahn Laboratory, Princeton University, Princeton, NJ 08544
3Institute of Community Medicine, Tromso University, Tromso, Norway
*To whom correspondence should be addressed. Dr. Olga G. Troyanskaya, E-mail: ogt{at}cs.princeton.edu
| Abstract |
|---|
Motivation: The heterogeneity of cancer cannot always be recognized by tumor morphology, but may be reflected by the underlying genetic aberrations. Array-CGH methods provide highthroughput data on genetic copy numbers, but determining the clinically relevant copy number changes remains a challenge. Conventional classification methods for linking recurrent alterations to clinical outcome ignore sequential correlations in selecting relevant features. Conversely, existing sequence classification methods can only model overall copy number instability, without regard to any particular position in the genome.
Results: Here we present the Heterogeneous Hidden Conditional Random Field, a new integrated array-CGH analysis method for jointly classifying tumors, inferring copy numbers, and identifying clinically relevant positions in recurrent alteration regions. By capturing the sequentiality as well as the locality of changes, our integrated model provides better noise reduction, and achieves more relevant gene retrieval and more accurate classification than existing methods. We provide an efficient L1-regularized discriminative training algorithm, which notably selects a small set of candidate genes most likely to be clinically relevant and driving the recurrent amplicons of importance. Our method thus provides unbiased starting points in deciding which genomic regions and which genes in particular to pursue for further examination. Our experiments on synthetic data and real genomic cancer prediction data show that our method is superior, both in prediction accuracy and relevant feature discovery, to existing methods. We also demonstrate that it can be used to generate novel biological hypotheses for breast cancer.
Contact: Olga G. Troyanskaya (ogt{at}cs.princeton.edu)
Associate Editor: Dr. Jonathan Wren
Received on July 29, 2008; revised on November 9, 2008; accepted on November 9, 2008