Integrative analysis of gene expression and copy number alterations using canonical correlation analysis.
BACKGROUND
With the rapid development of new genetic measurement methods, several types of genetic alterations can be quantified in a high-throughput manner. While the initial focus has been on investigating each data set separately, there is an increasing interest in studying the correlation structure between two or more data sets. Multivariate methods based on Canonical Correlation Analysis (CCA) have been proposed for integrating paired genetic data sets. The high dimensionality of microarray data imposes computational difficulties, which have been addressed for instance by studying the covariance structure of the data, or by reducing the number of variables prior to applying the CCA. In this work, we propose a new method for analyzing high-dimensional paired genetic data sets, which mainly emphasizes the correlation structure and still permits efficient application to very large data sets. The method is implemented by translating a regularized CCA to its dual form, where the computational complexity depends mainly on the number of samples instead of the number of variables. The optimal regularization parameters are chosen by cross-validation. We apply the regularized dual CCA, as well as a classical CCA preceded by a dimension-reducing Principal Components Analysis (PCA), to a paired data set of gene expression changes and copy number alterations in leukemia.


RESULTS
Using the correlation-maximizing methods, regularized dual CCA and PCA+CCA, we show that without pre-selection of known disease-relevant genes, and without using information about clinical class membership, an exploratory analysis singles out two patient groups, corresponding to well-known leukemia subtypes. Furthermore, the variables showing the highest relevance to the extracted features agree with previous biological knowledge concerning copy number alterations and gene expression changes in these subtypes. Finally, the correlation-maximizing methods are shown to yield results which are more biologically interpretable than those resulting from a covariance-maximizing method, and provide different insight compared to when each variable set is studied separately using PCA.


CONCLUSIONS
We conclude that regularized dual CCA as well as PCA+CCA are useful methods for exploratory analysis of paired genetic data sets, and can be efficiently implemented also when the number of variables is very large.
Soneson
Charlotte
C
Centre for Mathematical Sciences, Lund University, Box 118, SE-221 00 Lund, Sweden. lottas@maths.lth.se
Lilljebjörn
Henrik
H
Fioretos
Thoas
T
Fontes
Magnus
M
D000465:Algorithms; D030541:Databases, Genetic; D018628:Gene Dosage; D015870:Gene Expression; D020869:Gene Expression Profiling / Q000379:methods; D023281:Genomics / Q000379:methods*; D006801:Humans; D007938:Leukemia / Q000145:classification / Q000235:genetics; D025341:Principal Component Analysis
Cancer Res. 2002 Nov 1;62(21):6240-5
12414653
Genes Chromosomes Cancer. 2009 Jan;48(1):22-38
18803328
Nat Rev Cancer. 2006 Feb;6(2):107-16
16491070
Blood. 2003 Oct 15;102(8):2951-9
12730115
BMC Bioinformatics. 2009 Jan 26;10:34
19171069
Genes Chromosomes Cancer. 2009 Aug;48(8):637-60
19415723
Bioinformatics. 2004 May 22;20(8):1233-40
14871870
Nature. 2007 Apr 12;446(7137):758-64
17344859
Stat Appl Genet Mol Biol. 2009;8:Article 1
19222376
Cancer Res. 2002 Feb 15;62(4):1134-8
11861394
Proc Natl Acad Sci U S A. 2005 Dec 27;102(52):19069-74
16354839
Curr Opin Biotechnol. 2008 Feb;19(1):36-40
18162393
Nat Rev Cancer. 2007 Apr;7(4):233-45
17361217
Stat Appl Genet Mol Biol. 2008;7(1):Article3
18241193
Ann Appl Stat. 2010 Mar;4(1):53-77
24489618
Blood. 2006 Sep 1;108(5):1733-43
16705090
Nature. 2009 Apr 9;458(7239):719-24
19360079
Genome Biol. 2004;5(10):R80
15461798
Science. 2007 Feb 9;315(5813):848-53
17289997
Stat Appl Genet Mol Biol. 2008;7(1):Article 35
19049491
Mol Cancer Ther. 2006 Apr;5(4):853-67
16648555
Proc Natl Acad Sci U S A. 2002 Oct 1;99(20):12963-8
12297621
Leukemia. 2005 Jul;19(7):1224-8
15902281
Biostatistics. 2009 Jul;10(3):515-34
19377034
Neural Comput. 2004 Dec;16(12):2639-64
15516276