An integrated omics approach highlights how epigenetic events can explain and predict response to neoadjuvant chemotherapy and bevacizumab in breast cancer.

Treatment with the anti-angiogenic drug bevacizumab in addition to chemotherapy has shown efficacy for breast cancer in some clinical trials, but better biomarkers are needed to optimally select patients for treatment. Here, we present an omics approach where DNA methylation profiles are integrated with gene expression and results from proteomic data in breast cancer patients to predict response to therapy and pinpoint response-related epigenetic events. Fresh-frozen tumor biopsies taken before, during, and after treatment from human epidermal growth factor receptor 2 negative non-metastatic patients receiving neoadjuvant chemotherapy with or without bevacizumab were subjected to molecular profiling. Here, we report that DNA methylation at enhancer CpGs related to cell cycle regulation can predict response to chemotherapy and bevacizumab for the estrogen receptor positive subset of patients (AUC = 0.874), and we validated this observation in an independent patient cohort with a similar treatment regimen (AUC = 0.762). Combining the DNA methylation scores with the scores from a previously published protein signature resulted in a slight increase in the prediction performance (AUC = 0.784). We also show that tumors receiving the combination treatment underwent more extensive epigenetic alterations. Finally, we performed an integrative expression-methylation quantitative trait loci analysis on alterations in DNA methylation and gene expression levels, showing that the epigenetic alterations that occur during treatment are different between responders and non-responders and that these differences may be explained by the proliferation-epithelial-to-mesenchymal transition axis through the activity of grainyhead like transcription factor 2. Using tumor purity computed from copy number data, we developed a method for estimating cancer cell-specific methylation to confirm that the association to response reflects DNA methylation in cancer cells. Taken together, these results support the potential for clinical benefit of the addition of bevacizumab to chemotherapy when administered to the correct patients.

Multi-omics analysis reveals epigenetically regulated processes and patient classification in lung adenocarcinoma.

Aberrant DNA methylation is a hallmark of many cancer types. Despite our knowledge of epigenetic and transcriptomic alterations in lung adenocarcinoma (LUAD), we lack robust multi-modal molecular classifications for patient stratification. This is partly because the impact of epigenetic alterations on lung cancer development and progression is still not fully understood. To that end, we identified disease-associated processes under epigenetic regulation in LUAD. We performed a genome-wide expression-methylation Quantitative Trait Loci (emQTL) analysis by integrating DNA methylation and gene expression data from 453 patients in the TCGA cohort. Using a community detection algorithm, we identified distinct communities of CpG-gene associations with diverse biological processes. Interestingly, we identified a community linked to hormone response and lipid metabolism; the identified CpGs in this community were enriched in enhancer regions and binding regions of transcription factors such as FOXA1/2, GRHL2, HNF1B, AR, and ESR1. Furthermore, the CpGs were connected to their associated genes through chromatin interaction loops. These findings suggest that the expression of genes involved in hormone response and lipid metabolism in LUAD is epigenetically regulated through DNA methylation and enhancer-promoter interactions. By applying consensus clustering on the integrated expression-methylation pattern of the emQTL-genes and CpGs linked to hormone response and lipid metabolism, we further identified subclasses of patients with distinct prognoses. This novel patient stratification was validated in an independent patient cohort of 135 patients and showed increased prognostic significance compared to previously defined molecular subtypes.

MRI Assessment of Changes in Tumor Vascularization during Neoadjuvant Anti-Angiogenic Treatment in Locally Advanced Breast Cancer Patients.

Anti-VEGF (vascular endothelial growth factor) treatment improves response rates, but not progression-free or overall survival in advanced breast cancer. It has been suggested that subgroups of patients may benefit from this treatment; however, the effects of adding anti-VEGF treatment to a standard chemotherapy regimen in breast cancer patients are not well studied. Understanding the effects of the anti-vascular treatment on tumor vasculature may provide a selection of patients that can benefit. The aim of this study was to study the vascular effect of bevacizumab using clinical dynamic contrast-enhanced MRI (DCE-MRI). A total of 70 women were randomized to receive either chemotherapy alone or chemotherapy with bevacizumab for 25 weeks. DCE-MRI was performed at baseline and at 12 and 25 weeks, and in addition 25 of 70 patients agreed to participate in an early MRI after one week. Voxel-wise pharmacokinetic analysis was performed using semi-quantitative methods and the extended Tofts model. Vascular architecture was assessed by calculating the fractal dimension of the contrast-enhanced images. Changes during treatment were compared with baseline and between the treatment groups. There was no significant difference in tumor volume at any point; however, DCE-MRI parameters revealed differences in vascular function and vessel architecture. Adding bevacizumab to chemotherapy led to a pronounced reduction in vascular DCE-MRI parameters, indicating decreased vascularity. At 12 and 25 weeks, the difference between the treatment groups is severely reduced.

ABCB1 Amplicon Contains Cyclic AMP Response Element-Driven TRIP6 Gene in Taxane-Resistant MCF-7 Breast Cancer Sublines.

A limited number of studies are devoted to regulating TRIP6 expression in cancer. Hence, we aimed to unveil the regulation of TRIP6 expression in MCF-7 breast cancer cells (with high TRIP6 expression) and taxane-resistant MCF-7 sublines (manifesting even higher TRIP6 expression). We found that TRIP6 transcription is regulated primarily by the cyclic AMP response element (CRE) in hypomethylated proximal promoters in both taxane-sensitive and taxane-resistant MCF-7 cells. Furthermore, in taxane-resistant MCF-7 sublines, TRIP6 co-amplification with the neighboring ABCB1 gene, as witnessed by fluorescence in situ hybridization (FISH), led to TRIP6 overexpression. Ultimately, we found high TRIP6 mRNA levels in progesterone receptor-positive breast cancer and samples resected from premenopausal women.

Transcriptomic pan-cancer analysis using rank-based Bayesian inference.

The analysis of whole genomes of pan-cancer data sets provides a challenge for researchers, and we contribute to the literature concerning the identification of robust subgroups with clear biological interpretation. Specifically, we tackle this unsupervised problem via a novel rank-based Bayesian clustering method. The advantages of our method are the integration and quantification of all uncertainties related to both the input data and the model, the probabilistic interpretation of final results to allow straightforward assessment of the stability of clusters leading to reliable conclusions, and the transparent biological interpretation of the identified clusters since each cluster is characterized by its top-ranked genomic features. We applied our method to RNA-seq data from cancer samples from 12 tumor types from the Cancer Genome Atlas. We identified a robust clustering that mostly reflects tissue of origin but also includes pan-cancer clusters. Importantly, we identified three pan-squamous clusters composed of a mix of lung squamous cell carcinoma, head and neck squamous carcinoma, and bladder cancer, with different biological functions over-represented in the top genes that characterize the three clusters. We also found two novel subtypes of kidney cancer that show different prognosis, and we reproduced known subtypes of breast cancer. Taken together, our method allows the identification of robust and biologically meaningful clusters of pan-cancer samples.

Complex molecular profile of DNA repair genes in epithelial ovarian carcinoma patients with different sensitivity to platinum-based therapy.

Epithelial ovarian carcinoma (EOC) is known for high mortality due to diagnosis at advanced stages and frequent therapy resistance. Previous findings suggested that the DNA repair system is involved in the therapeutic response of cancer patients and DNA repair genes are promising targets for novel therapies. This study aimed to address complex inter-relations among gene expression levels, methylation profiles, and somatic mutations in DNA repair genes and EOC prognosis and therapy resistance status. We found significant associations of DUT expression with the presence of peritoneal metastases in EOC patients. The high-grade serous EOC subtype was enriched with TP53 mutations compared to other subtypes. Furthermore, somatic mutations in XPC and PRKDC were significantly associated with worse overall survival of EOC patients, and higher FAAP20 expression in platinum-resistant than platinum-sensitive patients was observed. We found higher methylation of RAD50 in platinum-resistant than in platinum-sensitive patients. Somatic mutations in BRCA1 and RAD9A were significantly associated with higher RBBP8 methylation in platinum-sensitive compared to platinum-resistant EOC patients. In conclusion, we discovered associations of several candidate genes from the DNA repair pathway with the prognosis and platinum resistance status of EOC patients, which deserve further validation as potential predictive biomarkers.

PHLPP1 regulates CFTR activity and lumen expansion through AMPK.

Complex organ development depends on single lumen formation and its expansion during tubulogenesis. This can be achieved by correct mitotic spindle orientation during cell division, combined with luminal fluid filling that generates hydrostatic pressure. Using a human 3D cell culture model, we have identified two regulators of these processes. We find that pleckstrin homology leucine-rich repeat protein phosphatase (PHLPP) 2 regulates mitotic spindle orientation, and thereby midbody positioning and maintenance of a single lumen. Silencing the sole PHLPP family phosphatase in Drosophila melanogaster, phlpp, resulted in defective spindle orientation in Drosophila neuroblasts. Importantly, cystic fibrosis transmembrane conductance regulator (CFTR) is the main channel regulating fluid transport in this system, stimulated by phosphorylation by protein kinase A and inhibited by the AMP-activated protein kinase AMPK. During lumen expansion, CFTR remains open through the action of PHLPP1, which stops activated AMPK from inhibiting ion transport through CFTR. In the absence of PHLPP1, the restraint on AMPK activity is lost and this tips the balance in the favour of channel closing, resulting in the lack of lumen expansion and accumulation of mucus.

Subtype and cell type specific expression of lncRNAs provide insight into breast cancer.

Long non-coding RNAs (lncRNAs) are involved in breast cancer pathogenesis through chromatin remodeling, transcriptional and post-transcriptional gene regulation. We report robust associations between lncRNA expression and breast cancer clinicopathological features in two population-based cohorts: SCAN-B and TCGA. Using co-expression analysis of lncRNAs with protein coding genes, we discovered three distinct clusters of lncRNAs. In silico cell type deconvolution coupled with single-cell RNA-seq analyses revealed that these three clusters were driven by cell type specific expression of lncRNAs. In one cluster lncRNAs were expressed by cancer cells and were mostly associated with the estrogen signaling pathways. In the two other clusters, lncRNAs were expressed either by immune cells or fibroblasts of the tumor microenvironment. To further investigate the cis-regulatory regions driving lncRNA expression in breast cancer, we identified subtype-specific transcription factor (TF) occupancy at lncRNA promoters. We also integrated lncRNA expression with DNA methylation data to identify long-range regulatory regions for lncRNA which were validated using ChiA-Pet-Pol2 loops. lncRNAs play an important role in shaping the gene regulatory landscape in breast cancer. We provide a detailed subtype and cell type-specific expression of lncRNA, which improves the understanding of underlying transcriptional regulation in breast cancer.

Rank-based Bayesian variable selection for genome-wide transcriptomic analyses.

Variable selection is crucial in high-dimensional omics-based analyses, since it is biologically reasonable to assume only a subset of non-noisy features contributes to the data structures. However, the task is particularly hard in an unsupervised setting, and a priori ad hoc variable selection is still a very frequent approach, despite the evident drawbacks and lack of reproducibility. We propose a Bayesian variable selection approach for rank-based unsupervised transcriptomic analysis. Making use of data rankings instead of the actual continuous measurements increases the robustness of conclusions when compared to classical statistical methods, and embedding variable selection into the inferential tasks allows complete reproducibility. Specifically, we develop a novel extension of the Bayesian Mallows model for variable selection that allows for a full probabilistic analysis, leading to coherent quantification of uncertainties. Simulation studies demonstrate the versatility and robustness of the proposed method in a variety of scenarios, as well as its superiority with respect to several competitors when varying the data dimension or data generating process. We use the novel approach to analyze genome-wide RNAseq gene expression data from ovarian cancer patients: several genes that affect cancer development are correctly detected in a completely unsupervised fashion, showing the usefulness of the method in the context of signature discovery for cancer genomics. Moreover, the possibility to also perform uncertainty quantification plays a key role in the subsequent biological investigation.

DNA methylation changes in response to neoadjuvant chemotherapy are associated with breast cancer survival.

BACKGROUND: Locally advanced breast cancer is a heterogeneous disease with respect to response to neoadjuvant chemotherapy (NACT) and survival. It is currently not possible to accurately predict who will benefit from the specific types of NACT. DNA methylation is an epigenetic mechanism known to play an important role in regulating gene expression and may serve as a biomarker for treatment response and survival. We investigated the potential role of DNA methylation as a prognostic marker for long-term survival (> 5 years) after NACT in breast cancer. METHODS: DNA methylation profiles of pre-treatment (n = 55) and post-treatment (n = 75) biopsies from 83 women with locally advanced breast cancer were investigated using the Illumina HumanMethylation450 BeadChip. The patients received neoadjuvant treatment with epirubicin and/or paclitaxel. Linear mixed models were used to associate DNA methylation to treatment response and survival based on clinical response to NACT (partial response or stable disease) and 5-year survival, respectively. LASSO regression was performed to identify a risk score based on the statistically significant methylation sites and Kaplan-Meier curve analysis was used to estimate survival probabilities using ten years of survival follow-up data. The risk score developed in our discovery cohort was validated in an independent validation cohort consisting of paired pre-treatment and post-treatment biopsies from 85 women with locally advanced breast cancer. Patients included in the validation cohort were treated with either doxorubicin or 5-FU and mitomycin NACT. RESULTS: DNA methylation patterns changed from before to after NACT in 5-year survivors, while no significant changes were observed in non-survivors or related to treatment response. DNA methylation changes included an overall loss of methylation at CpG islands and gain of methylation in non-CpG islands, and these changes affected genes linked to transcription factor activity, cell adhesion and immune functions. A risk score was developed based on four methylation sites which successfully predicted long-term survival in our cohort (p = 0.0034) and in an independent validation cohort (p = 0.049). CONCLUSION: Our results demonstrate that DNA methylation patterns in breast tumors change in response to NACT. These changes in DNA methylation show potential as prognostic biomarkers for breast cancer survival.

Pioneer transcription factors are associated with the modulation of DNA methylation patterns across cancers.

Methylation of cytosines on DNA is a prominent modification associated with gene expression regulation. Aberrant DNA methylation patterns have recurrently been linked to dysregulation of the regulatory program in cancer cells. To shed light on the underlying molecular mechanism driving this process, we hypothesised that aberrant methylation patterns could be controlled by the binding of specific transcription factors (TFs) across cancer types. By combining DNA methylation arrays and gene expression data with TF binding sites (TFBSs), we explored the interplay between TF binding and DNA methylation in 19 cancer types. We performed emQTL (expression-methylation quantitative trait loci) analyses independently in each cancer type and identified 13 TFs whose expression levels are correlated with local DNA methylation patterns around their binding sites in at least 2 cancer types. The 13 TFs are mainly associated with local demethylation and are enriched for pioneer function, suggesting a specific role for these TFs in modulating chromatin structure and transcription in cancer patients. Furthermore, we confirmed that de novo methylation is precluded across cancers at CpGs lying in genomic regions enriched for TF binding signatures associated with SP1, CTCF, NRF1, GABPA, KLF9, and/or YY1. The modulation of DNA methylation associated with TF binding was observed at cis-regulatory regions controlling immune- and cancer-associated pathways, corroborating that the emQTL signals were derived from both cancer and tumor-infiltrating cells. As a case example, we experimentally confirmed that FOXA1 knock-down is associated with higher methylation in regions bound by FOXA1 in breast cancer MCF-7 cells. Finally, we reported physical interactions between FOXA1 with TET1 and TET2 both in an in vitro setup and in vivo at physiological levels in MCF-7 cells, adding further support for FOXA1 attracting TET1 and TET2 to induce local demethylation in cancer cells.

Quantification of Tumor Hypoxia through Unsupervised Modelling of Consumption and Supply Hypoxia MR Imaging in Breast Cancer.

The purpose of the present study is to investigate if consumption and supply hypoxia (CSH) MR-imaging can depict breast cancer hypoxia, using the CSH-method initially developed for prostate cancer. Furthermore, to develop a generalized pan-cancer application of the CSH-method that doesn't require a hypoxia reference standard for training the CSH-parameters. In a cohort of 69 breast cancer patients, we generated, based on the principles of intravoxel incoherent motion modelling, images reflecting cellular density (apparent diffusion coefficient; ADC) and vascular density (perfusion fraction; fp). Combinations of the information in these images were compared to a molecular hypoxia score made from gene expression data, aiming to identify a way to apply the CSH-methodology in breast cancer. Attempts to adapt previously proposed models for prostate cancer included direct transfers and model parameter rescaling. A novel approach, based on rescaling ADC and fp data to give more nuanced response in the relevant physiologic range, was also introduced. The new CSH-method was validated in a prostate cancer cohort with known hypoxia status. The proposed CSH-method gave estimates of hypoxia that was strongly correlated to the molecular hypoxia score in breast cancer, and hypoxia as measured in pathology slices stained with pimonidazole in prostate cancer. The generalized approach to CSH-imaging depicted hypoxia in both breast and prostate cancers and requires no model training. It is easy to implement using readily available technology and encourages further investigation of CSH-imaging in other cancer entities and in other settings, with the goal being to overcome hypoxia-induced resistance to treatment.

Epigenetic alterations at distal enhancers are linked to proliferation in human breast cancer.

Aberrant DNA methylation is an early event in breast carcinogenesis and plays a critical role in regulating gene expression. Here, we perform genome-wide expression-methylation Quantitative Trait Loci (emQTL) analysis through the integration of DNA methylation and gene expression to identify disease-driving pathways under epigenetic control. By grouping the emQTLs using biclustering we identify associations representing important biological processes associated with breast cancer pathogenesis including regulation of proliferation and tumor-infiltrating fibroblasts. We report genome-wide loss of enhancer methylation at binding sites of proliferation-driving transcription factors including CEBP-ß, FOSL1, and FOSL2 with concomitant high expression of proliferation-related genes in aggressive breast tumors as we confirm with scRNA-seq. The identified emQTL-CpGs and genes were found connected through chromatin loops, indicating that proliferation in breast tumors is under epigenetic regulation by DNA methylation. Interestingly, the associations between enhancer methylation and proliferation-related gene expression were also observed within known subtypes of breast cancer, suggesting a common role of epigenetic regulation of proliferation. Taken together, we show that proliferation in breast cancer is linked to loss of methylation at specific enhancers and transcription factor binding and gene activation through chromatin looping.

Multi-Omics Marker Analysis Enables Early Prediction of Breast Tumor Progression.

Ductal carcinoma in situ (DCIS) is a preinvasive form of breast cancer with a highly variable potential of becoming invasive and affecting mortality of the patients. Due to the lack of accurate markers of disease progression, many women with detected DCIS are currently overtreated. To distinguish those DCIS cases who are likely to require therapy from those who should be left untreated, there is a need for robust and predictive biomarkers extracted from molecular or genetic profiles. We developed a supervised machine learning approach that implements multi-omics feature selection and model regularization for the identification of biomarker combinations that could be used to distinguish low-risk DCIS lesions from those with a higher likelihood of progression. To investigate the genetic heterogeneity of disease progression, we applied this approach to 40 pure DCIS and 259 invasive breast cancer (IBC) samples profiled with genome-wide transcriptomics, DNA methylation, and DNA copy number variation. Feature selection using the multi-omics Lasso-regularized algorithm identified both known genes involved in breast cancer development, as well as novel markers for early detection. Even though the gene expression-based model features led to the highest classification accuracy alone, methylation data provided a complementary source of features and improved especially the sensitivity of correctly classifying DCIS cases. We also identified a number of repeatedly misclassified DCIS cases when using either the expression or methylation markers. A small panel of 10 gene markers was able to distinguish DCIS and IBC cases with high accuracy in nested cross-validation (AU-ROC = 0.99). The marker panel was not specific to any of the established breast cancer subtypes, suggesting that the 10-gene signature may provide a subtype-agnostic and cost-effective approach for breast cancer detection and patient stratification. We further confirmed high accuracy of the 10-gene signature in an external validation cohort (AU-ROC = 0.95), profiled using distinct transcriptomic assay, hence demonstrating robustness of the risk signature.

Crosstalk between microRNA expression and DNA methylation drives the hormone-dependent phenotype of breast cancer.

BACKGROUND: Abnormal DNA methylation is observed as an early event in breast carcinogenesis. However, how such alterations arise is still poorly understood. microRNAs (miRNAs) regulate gene expression at the post-transcriptional level and play key roles in various biological processes. Here, we integrate miRNA expression and DNA methylation at CpGs to study how miRNAs may affect the breast cancer methylome and how DNA methylation may regulate miRNA expression. METHODS: miRNA expression and DNA methylation data from two breast cancer cohorts, Oslo2 (n = 297) and The Cancer Genome Atlas (n = 439), were integrated through a correlation approach that we term miRNA-methylation Quantitative Trait Loci (mimQTL) analysis. Hierarchical clustering was used to identify clusters of miRNAs and CpGs that were further characterized through analysis of mRNA/protein expression, clinicopathological features, in silico deconvolution, chromatin state and accessibility, transcription factor binding, and long-range interaction data. RESULTS: Clustering of the significant mimQTLs identified distinct groups of miRNAs and CpGs that reflect important biological processes associated with breast cancer pathogenesis. Notably, two major miRNA clusters were related to immune or fibroblast infiltration, hence identifying miRNAs associated with cells of the tumor microenvironment, while another large cluster was related to estrogen receptor (ER) signaling. Studying the chromatin landscape surrounding CpGs associated with the estrogen signaling cluster, we found that miRNAs from this cluster are likely to be regulated through DNA methylation of enhancers bound by FOXA1, GATA2, and ER-alpha. Further, at the hub of the estrogen cluster, we identified hsa-miR-29c-5p as negatively correlated with the mRNA and protein expression of DNA methyltransferase DNMT3A, a key enzyme regulating DNA methylation. We found deregulation of hsa-miR-29c-5p already present in pre-invasive breast lesions and postulate that hsa-miR-29c-5p may trigger early event abnormal DNA methylation in ER-positive breast cancer. CONCLUSIONS: We describe how miRNA expression and DNA methylation interact and associate with distinct breast cancer phenotypes.

Multi-modal meta-analysis of cancer cell line omics profiles identifies ECHDC1 as a novel breast tumor suppressor.

Molecular and functional profiling of cancer cell lines is subject to laboratory-specific experimental practices and data analysis protocols. The current challenge therefore is how to make an integrated use of the omics profiles of cancer cell lines for reliable biological discoveries. Here, we carried out a systematic analysis of nine types of data modalities using meta-analysis of 53 omics studies across 12 research laboratories for 2,018 cell lines. To account for a relatively low consistency observed for certain data modalities, we developed a robust data integration approach that identifies reproducible signals shared among multiple data modalities and studies. We demonstrated the power of the integrative analyses by identifying a novel driver gene, ECHDC1, with tumor suppressive role validated both in breast cancer cells and patient tumors. The multi-modal meta-analysis approach also identified synthetic lethal partners of cancer drivers, including a co-dependency of PTEN deficient endometrial cancer cells on RNA helicases.

Epigenome-wide gene-age interaction analysis reveals reversed effects of PRODH DNA methylation on survival between young and elderly early-stage NSCLC patients.

DNA methylation changes during aging, but it remains unclear whether the effect of DNA methylation on lung cancer survival varies with age. Such an effect could decrease prediction accuracy and treatment efficacy. We performed a methylation-age interaction analysis using 1,230 early-stage lung adenocarcinoma patients from five cohorts. A Cox proportional hazards model was used to investigate lung adenocarcinoma and squamous cell carcinoma patients for methylation-age interactions, which were further confirmed in a validation phase. We identified one adenocarcinoma-specific CpG probe, cg14326354PRODH, with effects significantly modified by age (HRinteraction = 0.989; 95% CI: 0.986-0.994; P = 9.18×10-7). The effect of low methylation was reversed for young and elderly patients categorized by the boundary of 95% CI standard (HRyoung = 2.44; 95% CI: 1.26-4.72; P = 8.34×10-3; HRelderly = 0.58; 95% CI: 0.42-0.82; P = 1.67×10-3). Moreover, there was an antagonistic interaction between low cg14326354PRODH methylation and elderly age (HRinteraction = 0.21; 95% CI: 0.11-0.40; P = 2.20×10-6). In summary, low methylation of cg14326354PRODH might benefit survival of elderly lung adenocarcinoma patients, providing new insight to age-specific prediction and potential drug targeting.

Independent Validation of Early-Stage Non-Small Cell Lung Cancer Prognostic Scores Incorporating Epigenetic and Transcriptional Biomarkers With Gene-Gene Interactions and Main Effects.

BACKGROUND: DNA methylation and gene expression are promising biomarkers of various cancers, including non-small cell lung cancer (NSCLC). Besides the main effects of biomarkers, the progression of complex diseases is also influenced by gene-gene (G×G) interactions. RESEARCH QUESTION: Would screening the functional capacity of biomarkers on the basis of main effects or interactions, using multiomics data, improve the accuracy of cancer prognosis? STUDY DESIGN AND METHODS: Biomarker screening and model validation were used to construct and validate a prognostic prediction model. NSCLC prognosis-associated biomarkers were identified on the basis of either their main effects or interactions with two types of omics data. A prognostic score incorporating epigenetic and transcriptional biomarkers, as well as clinical information, was independently validated. RESULTS: Twenty-six pairs of biomarkers with G×G interactions and two biomarkers with main effects were significantly associated with NSCLC survival. Compared with a model using clinical information only, the accuracy of the epigenetic and transcriptional biomarker-based prognostic model, measured by area under the receiver operating characteristic curve (AUC), increased by 35.38% (95% CI, 27.09%-42.17%; P = 5.10 × 10-17) and 34.85% (95% CI, 26.33%-41.87%; P = 2.52 × 10-18) for 3- and 5-year survival, respectively, which exhibited a superior predictive ability for NSCLC survival (AUC3 year, 0.88 [95% CI, 0.83-0.93]; and AUC5 year, 0.89 [95% CI, 0.83-0.93]) in an independent Cancer Genome Atlas population. G×G interactions contributed a 65.2% and 91.3% increase in prediction accuracy for 3- and 5-year survival, respectively. INTERPRETATION: The integration of epigenetic and transcriptional biomarkers with main effects and G×G interactions significantly improves the accuracy of prognostic prediction of early-stage NSCLC survival.

Trans-omics biomarker model improves prognostic prediction accuracy for early-stage lung adenocarcinoma.

Limited studies have focused on developing prognostic models with trans-omics biomarkers for early-stage lung adenocarcinoma (LUAD). We performed integrative analysis of clinical information, DNA methylation, and gene expression data using 825 early-stage LUAD patients from 5 cohorts. Ranger algorithm was used to screen prognosis-associated biomarkers, which were confirmed with a validation phase. Clinical and biomarker information was fused using an iCluster plus algorithm, which significantly distinguished patients into high- and low-mortality risk groups (Pdiscovery = 0.01 and Pvalidation = 2.71×10-3). Further, potential functional DNA methylation-gene expression-overall survival pathways were evaluated by causal mediation analysis. The effect of DNA methylation level on LUAD survival was significantly mediated through gene expression level. By adding DNA methylation and gene expression biomarkers to a model of only clinical data, the AUCs of the trans-omics model improved by 18.3% (to 87.2%) and 16.4% (to 85.3%) in discovery and validation phases, respectively. Further, concordance index of the nomogram was 0.81 and 0.77 in discovery and validation phases, respectively. Based on systematic review of published literatures, our model was superior to all existing models for early-stage LUAD. In summary, our trans-omics model may help physicians accurately identify patients with high mortality risk.

Toward Personalized Computer Simulation of Breast Cancer Treatment: A Multiscale Pharmacokinetic and Pharmacodynamic Model Informed by Multitype Patient Data.

The usefulness of mechanistic models to disentangle complex multiscale cancer processes, such as treatment response, has been widely acknowledged. However, a major barrier for multiscale models to predict treatment outcomes in individual patients lies in their initialization and parametrization, which needs to reflect individual cancer characteristics accurately. In this study, we use multitype measurements acquired routinely on a single breast tumor, including histopathology, MRI, and molecular profiling, to personalize parts of a complex multiscale model of breast cancer treated with chemotherapeutic and antiangiogenic agents. The model accounts for drug pharmacokinetics and pharmacodynamics. We developed an open-source computer program that simulates cross-sections of tumors under 12-week therapy regimens and used it to individually reproduce and elucidate treatment outcomes of 4 patients. Two of the tumors did not respond to therapy, and model simulations were used to suggest alternative regimens with improved outcomes dependent on the tumor's individual characteristics. It was determined that more frequent and lower doses of chemotherapy reduce tumor burden in a low proliferative tumor while lower doses of antiangiogenic agents improve drug penetration in a poorly perfused tumor. Furthermore, using this model, we were able to correctly predict the outcome in another patient after 12 weeks of treatment. In summary, our model bridges multitype clinical data to shed light on individual treatment outcomes. SIGNIFICANCE: Mathematical modeling is used to validate possible mechanisms of tumor growth, resistance, and treatment outcome.