Dataset Information

Optimized application of penalized regression methods to diverse genomic data.

ABSTRACT:

Motivation

Penalized regression methods have been adopted widely for high-dimensional feature selection and prediction in many bioinformatic and biostatistical contexts. While their theoretical properties are well-understood, specific methodology for their optimal application to genomic data has not been determined.

Results

Through simulation of contrasting scenarios of correlated high-dimensional survival data, we compared the LASSO, Ridge and Elastic Net penalties for prediction and variable selection. We found that a 2D tuning of the Elastic Net penalties was necessary to avoid mimicking the performance of LASSO or Ridge regression. Furthermore, we found that in a simulated scenario favoring the LASSO penalty, a univariate pre-filter made the Elastic Net behave more like Ridge regression, which was detrimental to prediction performance. We demonstrate the real-life application of these methods to predicting the survival of cancer patients from microarray data, and to classification of obese and lean individuals from metagenomic data. Based on these results, we provide an optimized set of guidelines for the application of penalized regression for reproducible class comparison and prediction with genomic data.

Availability and implementation

A parallelized implementation of the methods presented for regression and for simulation of synthetic data is provided as the pensim R package, available at http://cran.r-project.org/web/packages/pensim/index.html.

Contact

chuttenh@hsph.harvard.edu; juris@ai.utoronto.ca

Supplementary information

Supplementary data are available at Bioinformatics online.

SUBMITTER: Waldron L

PROVIDER: S-EPMC3232376 | biostudies-literature | 2011 Dec

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Optimized application of penalized regression methods to diverse genomic data.

Waldron Levi L Pintilie Melania M Tsao Ming-Sound MS Shepherd Frances A FA Huttenhower Curtis C Jurisica Igor I

Bioinformatics (Oxford, England) 20111201 24

<h4>Motivation</h4>Penalized regression methods have been adopted widely for high-dimensional feature selection and prediction in many bioinformatic and biostatistical contexts. While their theoretical properties are well-understood, specific methodology for their optimal application to genomic data has not been determined.<h4>Results</h4>Through simulation of contrasting scenarios of correlated high-dimensional survival data, we compared the LASSO, Ridge and Elastic Net penalties for prediction ...[more]

PMID: 22156367

Similar Datasets

Project description:Omics data integration is becoming necessary to investigate the still unknown genomic mechanisms of complex diseases. During the integration process, many challenges arise such as data heterogeneity, the smaller number of individuals in comparison to the number of parameters, multicollinearity, and interpretation and validation of results due to their complexity and lack of knowledge about biological mechanisms. To overcome some of these issues, innovative statistical approaches are being developed. In this work, we applied penalized regression methods (LASSO and ENET) to explore relationships between common genetic variants, DNA methylation and gene expression measured in bladder tumor samples and have proposed a permutation-based method to concomitantly assess significance and correct by multiple testing with the MaxT algorithm. The overall analysis flow consisted of three steps: (1) SNPs/CpGs were selected per each gene probe within 1Mb window upstream and downstream the gene; (2) LASSO and ENET were applied to assess the association between each expression probe and the selected SNPs/CpGs in three multivariable models (SNP, CPG, and Global models, the latter integrating SNPs and CPGs); and (3) the significance of each model was assessed using the permutation-based MaxT method. We identified 48 genes whom expression levels were associated with both SNPs and GPGs. Importantly, we replicated results for 36 (75%) of them in an independent data set (TCGA). We checked the performance of the proposed method with a simulation study and further supported our results with a biological interpretation based on an enrichment analysis. The approach we propose allows reducing computational time and is flexibly and easy to implement when analyzing several omics data. Our results highlight the importance of integrating omics data by applying appropriate statistical strategies to discover new insights into the complexity of disease genetic mechanisms. Gene expression data were obtained from 44 tumor samples using the Affymetrix DNA Microarray Human Gene 1.0 ST Array with 32,321 probes.

Project description:Omics data integration is becoming necessary to investigate the still unknown genomic mechanisms of complex diseases. During the integration process, many challenges arise such as data heterogeneity, the smaller number of individuals in comparison to the number of parameters, multicollinearity, and interpretation and validation of results due to their complexity and lack of knowledge about biological mechanisms. To overcome some of these issues, innovative statistical approaches are being developed. In this work, we applied penalized regression methods (LASSO and ENET) to explore relationships between common genetic variants, DNA methylation and gene expression measured in bladder tumor samples and have proposed a permutation-based method to concomitantly assess significance and correct by multiple testing with the MaxT algorithm. The overall analysis flow consisted of three steps: (1) SNPs/CpGs were selected per each gene probe within 1Mb window upstream and downstream the gene; (2) LASSO and ENET were applied to assess the association between each expression probe and the selected SNPs/CpGs in three multivariable models (SNP, CPG, and Global models, the latter integrating SNPs and CPGs); and (3) the significance of each model was assessed using the permutation-based MaxT method. We identified 48 genes whom expression levels were associated with both SNPs and GPGs. Importantly, we replicated results for 36 (75%) of them in an independent data set (TCGA). We checked the performance of the proposed method with a simulation study and further supported our results with a biological interpretation based on an enrichment analysis. The approach we propose allows reducing computational time and is flexibly and easy to implement when analyzing several omics data. Our results highlight the importance of integrating omics data by applying appropriate statistical strategies to discover new insights into the complexity of disease genetic mechanisms. Bisulphite modification of 46 tumor DNA samples using EZ-96 DNA METHYLATIONGOLD KIT (Zymo Research, Irvin, CA, USA), CpG methylation data was generated using the Infinum Human Methylation 27 BeadChip Kit that detected the CpG sites with two probes, one designed against the unmethylated site (signal U) and the other against the methylated site (signal M).

Dataset Information

Optimized application of penalized regression methods to diverse genomic data.

Motivation

Results

Availability and implementation

Contact

Supplementary information

Publications

Optimized application of penalized regression methods to diverse genomic data.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure