Abstract:Evaluation of survival models to predict cancer patient prognosis is one of the most important areas of emphasis in cancer research. A binary classification approach has difficulty directly predicting survival due to the characteristics of censored observations and the fact that the predictive power depends on the threshold used to set two classes. In contrast, the traditional Cox regression approach has some drawbacks in the sense that it does not allow for the identification of interactions between genomic features, which could have key roles associated with cancer prognosis. In addition, data integration is regarded as one of the important issues in improving the predictive power of survival models since cancer could be caused by multiple alterations through meta-dimensional genomic data including genome, epigenome, transcriptome, and proteome. Here we have proposed a new integrative framework designed to perform these three functions simultaneously: (1) predicting censored survival data; (2) integrating meta-dimensional omics data; (3) identifying interactions within/between meta-dimensional genomic features associated with survival. In order to predict censored survival time, martingale residuals were calculated as a new continuous outcome and a new fitness function used by the grammatical evolution neural network (GENN) based on mean absolute difference of martingale residuals was implemented. To test the utility of the proposed framework, a simulation study was conducted, followed by an analysis of meta-dimensional omics data including copy number, gene expression, DNA methylation, and protein expression data in breast cancer retrieved from The Cancer Genome Atlas (TCGA). On the basis of the results from breast cancer dataset, we were able to identify interactions not only within a single dimension of genomic data but also between meta-dimensional omics data that are associated with survival. Notably, the predictive power of our best meta-dimensional model was 73% which outperformed all of the other models conducted based on a single dimension of genomic data. Breast cancer is an extremely heterogeneous disease and the high levels of genomic diversity within/between breast tumors could affect the risk of therapeutic responses and disease progression. Thus, identifying interactions within/between meta-dimensional omics data associated with survival in breast cancer is expected to deliver direction for improved meta-dimensional prognostic biomarkers and therapeutic targets.

Two-level Bayesian interaction analysis for survival data incorporating pathway information

Gene-gene interaction analysis incorporating network information via a structured Bayesian approach

Incorporating biological information into linear models: A Bayesian approach to the selection of pathways and genes

BHAFT: Bayesian heredity-constrained accelerated failure time models for detecting gene-environment interactions in survival analysis

Bayesian Phase II Adaptive Randomization by Jointly Modeling Time-to-event Efficacy and Binary Toxicity

Link-based identification of survival time-related biological pathways

Integration of Pathway Structure Information into a Reweighted Partial Cox Regression Approach for Survival Analysis on High-Dimensional Gene Expression Data

Gene-Based Association Analysis for Censored Traits Via Fixed Effect Functional Regressions

Joint longitudinal and time-to-event models for multilevel hierarchical data

Overlapping group screening for detection of gene-gene interactions: application to gene expression profiles with survival trait

Identification of Gene-Environment Interactions by Non-Parametric Kendall's Partial Correlation with Application to TCGA Ultrahigh-Dimensional Survival Genomic Data

Predicting censored survival data based on the interactions between meta-dimensional omics data in breast cancer

Penalized Bayesian forward continuation ratio model with application to high-dimensional data with discrete survival outcomes

Multicategory Survival Outcomes Classification via Overlapping Group Screening Process Based on Multinomial Logistic Regression Model With Application to TCGA Transcriptomic Data

Gene-level pharmacogenetic analysis on survival outcomes using gene-trait similarity regression

Kernel-based hierarchical structural component models for pathway analysis on survival phenotype

Gene-based Association Analysis for Bivariate Time-to-event Data through Functional Regression with Copula Models

Bayesian Variable Selection for Parametric Survival Model with Applications to Cancer Omics Data

Dependent Dirichlet Processes for Analysis of a Generalized Shared Frailty Model

Bridging the gap between two‐stage and joint models: The case of tumor growth inhibition and overall survival models

Survival-supervised latent Dirichlet allocation models for genomic analysis of time-to-event outcomes