H3K9 methylation drives resistance to androgen receptor-antagonist therapy in prostate cancer.

Antiandrogen strategies remain the prostate cancer treatment backbone, but drug resistance develops. We show that androgen blockade in prostate cancer leads to derepression of retroelements (REs) followed by a double-stranded RNA (dsRNA)-stimulated interferon response that blocks tumor growth. A forward genetic approach identified H3K9 trimethylation (H3K9me3) as an essential epigenetic adaptation to antiandrogens, which enabled transcriptional silencing of REs that otherwise stimulate interferon signaling and glucocorticoid receptor expression. Elevated expression of terminal H3K9me3 writers was associated with poor patient hormonal therapy outcomes. Forced expression of H3K9me3 writers conferred resistance, whereas inhibiting H3K9-trimethylation writers and readers restored RE expression, blocking antiandrogen resistance. Our work reveals a drug resistance axis that integrates multiple cellular signaling elements and identifies potential pharmacologic vulnerabilities.

SYSMut: decoding the functional significance of rare somatic mutations in cancer.

Current tailored-therapy efforts in cancer are largely focused on a small number of highly recurrently mutated driver genes but therapeutic targeting of these oncogenes remains challenging. However, the vast number of genes mutated infrequently across cancers has received less attention, in part, due to a lack of understanding of their biological significance. We present SYSMut, an extendable systems biology platform that can robustly infer the biologic consequences of somatic mutations by integrating routine multiomics profiles in primary tumors. We establish SYSMut's improved performance vis-à-vis state-of-the-art driver gene identification methodologies by recapitulating the functional impact of known driver genes, while additionally identifying novel functionally impactful mutated genes across 29 cancers. Subsequent application of SYSMut on low-frequency gene mutations in head and neck squamous cell (HNSC) cancers, followed by molecular and pharmacogenetic validation, revealed the lipidogenic network as a novel therapeutic vulnerability in aggressive HNSC cancers. SYSMut is thus a robust scalable framework that enables the discovery of new targetable avenues in cancer.

The Ephrin B2 Receptor Tyrosine Kinase Is a Regulator of Proto-oncogene MYC and Molecular Programs Central to Barrett's Neoplasia.

BACKGROUND & AIMS: Mechanisms contributing to the onset and progression of Barrett's (BE)-associated esophageal adenocarcinoma (EAC) remain elusive. Here, we interrogated the major signaling pathways deregulated early in the development of Barrett's neoplasia. METHODS: Whole-transcriptome RNA sequencing analysis was performed in primary BE, EAC, normal esophageal squamous, and gastric biopsy tissues (n = 89). Select pathway components were confirmed by quantitative polymerase chain reaction in an independent cohort of premalignant and malignant biopsy tissues (n = 885). Functional impact of selected pathway was interrogated using transcriptomic, proteomic, and pharmacogenetic analyses in mammalian esophageal organotypic and patient-derived BE/EAC cell line models, in vitro and/or in vivo. RESULTS: The vast majority of primary BE/EAC tissues and cell line models showed hyperactivation of EphB2 signaling. Transcriptomic/proteomic analyses identified EphB2 as an endogenous binding partner of MYC binding protein 2, and an upstream regulator of c-MYC. Knockdown of EphB2 significantly impeded the viability/proliferation of EAC and BE cells in vitro/in vivo. Activation of EphB2 in normal esophageal squamous 3-dimensional organotypes disrupted epithelial maturation and promoted columnar differentiation programs, notably including MYC. EphB2 and MYC showed selective induction in esophageal submucosal glands with acinar ductal metaplasia, and in a porcine model of BE-like esophageal submucosal gland spheroids. Clinically approved inhibitors of MEK, a protein kinase that regulates MYC, effectively suppressed EAC tumor growth in vivo. CONCLUSIONS: The EphB2 signaling is frequently hyperactivated across the BE-EAC continuum. EphB2 is an upstream regulator of MYC, and activation of EphB2-MYC axis likely precedes BE development. Targeting EphB2/MYC could be a promising therapeutic strategy for this often refractory and aggressive cancer.

The transcriptional repressor BCL11A promotes breast cancer metastasis.

The phenotypes of each breast cancer subtype are defined by their transcriptomes. However, the transcription factors that regulate differential patterns of gene expression that contribute to specific disease outcomes are not well understood. Here, using gene silencing and overexpression approaches, RNA-Seq, and splicing analysis, we report that the transcription factor B-cell leukemia/lymphoma 11A (BCL11A) is highly expressed in triple-negative breast cancer (TNBC) and drives metastatic disease. Moreover, BCL11A promotes cancer cell invasion by suppressing the expression of muscleblind-like splicing regulator 1 (MBNL1), a splicing regulator that suppresses metastasis. This ultimately increases the levels of an alternatively spliced isoform of integrin-α6 (ITGA6), which is associated with worse patient outcomes. These results suggest that BCL11A sustains TNBC cell invasion and metastatic growth by repressing MBNL1-directed splicing of ITGA6 Our findings also indicate that BCL11A lies at the interface of transcription and splicing and promotes aggressive TNBC phenotypes.

Systems Biology Analyses Show Hyperactivation of Transforming Growth Factor-ß and JNK Signaling Pathways in Esophageal Cancer.

BACKGROUND & AIMS: Esophageal adenocarcinoma (EAC) is resistant to standard chemoradiation treatments, and few targeted therapies are available. We used large-scale tissue profiling and pharmacogenetic analyses to identify deregulated signaling pathways in EAC tissues that might be targeted to slow tumor growth or progression. METHODS: We collected 397 biopsy specimens from patients with EAC and nonmalignant Barrett's esophagus (BE), with or without dysplasia. We performed RNA-sequencing analyses and used systems biology approaches to identify pathways that are differentially activated in EAC vs nonmalignant dysplastic tissues; pathway activities were confirmed with immunohistochemistry and quantitative real-time polymerase chain reaction analyses of signaling components in patient tissue samples. Human EAC (FLO-1 and EsoAd1), dysplastic BE (CP-B, CP-C, CP-D), and nondysplastic BE (CP-A) cells were incubated with pharmacologic inhibitors or transfected with small interfering RNAs. We measured effects on proliferation, colony formation, migration, and/or growth of xenograft tumors in nude mice. RESULTS: Comparisons of EAC vs nondysplastic BE tissues showed hyperactivation of transforming growth factor-ß (TGFB) and/or Jun N-terminal kinase (JNK) signaling pathways in more than 80% of EAC samples. Immunohistochemical analyses showed increased nuclear localization of phosphorylated JUN and SMAD proteins in EAC tumor tissues compared with nonmalignant tissues. Genes regulated by the TGFB and JNK pathway were overexpressed specifically in EAC and dysplastic BE. Pharmacologic inhibition or knockdown of TGFB or JNK signaling components in EAC cells (FLO-1 or EsoAd1) significantly reduced cell proliferation, colony formation, cell migration, and/or growth of xenograft tumors in mice in a SMAD4-independent manner. Inhibition of the TGFB pathway in BE cell lines reduced the proliferation of dysplastic, but not nondysplastic, cells. CONCLUSIONS: In a transcriptome analysis of EAC and nondysplastic BE tissues, we found the TGFB and JNK signaling pathways to be hyperactivated in EACs and the genes regulated by these pathways to be overexpressed in EAC and dysplastic BE. Inhibiting these pathways in EAC cells reduces their proliferation, migration, and formation of xenograft tumors. Strategies to block the TGFB and JNK signaling pathways might be developed for treatment of EAC.

Interferon-beta represses cancer stem cell properties in triple-negative breast cancer.

Triple-negative breast cancer (TNBC), the deadliest form of this disease, lacks a targeted therapy. TNBC tumors that fail to respond to chemotherapy are characterized by a repressed IFN/signal transducer and activator of transcription (IFN/STAT) gene signature and are often enriched for cancer stem cells (CSCs). We have found that human mammary epithelial cells that undergo an epithelial-to-mesenchymal transition (EMT) following transformation acquire CSC properties. These mesenchymal/CSCs have a significantly repressed IFN/STAT gene expression signature and an enhanced ability to migrate and form tumor spheres. Treatment with IFN-beta (IFN-ß) led to a less aggressive epithelial/non-CSC-like state, with repressed expression of mesenchymal proteins (VIMENTIN, SLUG), reduced migration and tumor sphere formation, and reexpression of CD24 (a surface marker for non-CSCs), concomitant with an epithelium-like morphology. The CSC-like properties were correlated with high levels of unphosphorylated IFN-stimulated gene factor 3 (U-ISGF3), which was previously linked to resistance to DNA damage. Inhibiting the expression of IRF9 (the DNA-binding component of U-ISGF3) reduced the migration of mesenchymal/CSCs. Here we report a positive translational role for IFN-ß, as gene expression profiling of patient-derived TNBC tumors demonstrates that an IFN-ß metagene signature correlates with improved patient survival, an immune response linked with tumor-infiltrating lymphocytes (TILs), and a repressed CSC metagene signature. Taken together, our findings indicate that repressed IFN signaling in TNBCs with CSC-like properties is due to high levels of U-ISGF3 and that treatment with IFN-ß reduces CSC properties, suggesting a therapeutic strategy to treat drug-resistant, highly aggressive TNBC tumors.

Novel recurrently mutated genes in African American colon cancers.

We used whole-exome and targeted sequencing to characterize somatic mutations in 103 colorectal cancers (CRC) from African Americans, identifying 20 new genes as significantly mutated in CRC. Resequencing 129 Caucasian derived CRCs confirmed a 15-gene set as a preferential target for mutations in African American CRCs. Two predominant genes, ephrin type A receptor 6 (EPHA6) and folliculin (FLCN), with mutations exclusive to African American CRCs, are by genetic and biological criteria highly likely African American CRC driver genes. These previously unsuspected differences in the mutational landscapes of CRCs arising among individuals of different ethnicities have potential to impact on broader disparities in cancer behaviors.

Brief-exposure to preoperative bevacizumab reveals a TGF-ß signature predictive of response in HER2-negative breast cancers.

To best define biomarkers of response, and to shed insight on mechanism of action of certain clinically important agents for early breast cancer, we used a brief-exposure paradigm in the preoperative setting to study transcriptional changes in patient tumors that occur with one dose of therapy prior to combination chemotherapy. Tumor biopsies from breast cancer patients enrolled in two preoperative clinical trials were obtained at baseline and after one dose of bevacizumab (HER2-negative), trastuzumab (HER2-positive) or nab-paclitaxel, followed by treatment with combination chemo-biologic therapy. RNA-Sequencing based PAM50 subtyping at baseline of 46 HER2-negative patients revealed a strong association between the basal-like subtype and pathologic complete response (pCR) to chemotherapy plus bevacizumab (p ≤ 0.0027), but did not provide sufficient specificity to predict response. However, a single dose of bevacizumab resulted in down-regulation of a well-characterized TGF-ß activity signature in every single breast tumor that achieved pCR (p ≤ 0.004). The TGF-ß signature was confirmed to be a tumor-specific read-out of the canonical TGF-ß pathway using pSMAD2 (p ≤ 0.04), with predictive power unique to brief-exposure to bevacizumab (p ≤ 0.016), but not trastuzumab or nab-paclitaxel. Down-regulation of TGF-ß activity was associated with reduction in tumor hypoxia by transcription and protein levels, suggesting therapy-induced disruption of an autocrine-loop between tumor stroma and malignant cells. Modulation of the TGF-ß pathway upon brief-exposure to bevacizumab may provide an early functional readout of pCR to preoperative anti-angiogenic therapy in HER2-negative breast cancer, thus providing additional avenues for exploration in both preclinical and clinical settings with these agents.

Biomarkers for Predicting Response to Anti-HER2 Agents.

The HER2 receptor is amplified or overexpressed in approximately 20% of all breast cancers, but despite significant efforts of the clinical research community and a growing number of anti-HER2 agents, a significant number of patients with HER2-positive breast cancer either progress or suffer disease relapse within 5-10 years. The development of robust biomarkers that predict response to anti-HER2 agents is therefore an important clinical need to prevent overtreatment and to enable earlier assignment of patients to more optimal therapies. Here we review some of the recent advances in the field by focusing on pathways mediating resistance to anti-HER2 therapies, and the role of the immune system and cancer stem cells in therapy response. We also review preoperative treatment strategies and research paradigms that show promise in identifying novel biomarkers of response while also enabling the delineation of the mechanisms underlying clinical benefit from anti-HER2 therapies.

A microRNA panel to discriminate carcinomas from high-grade intraepithelial neoplasms in colonoscopy biopsy tissue.

OBJECTIVE: It is a challenge to differentiate invasive carcinomas from high-grade intraepithelial neoplasms in colonoscopy biopsy tissues. In this study, microRNA profiles were evaluated in the transformation of colorectal carcinogenesis to discover new molecular markers for identifying a carcinoma in colonoscopy biopsy tissues where the presence of stromal invasion cells is not detectable by microscopic analysis. METHODS: The expression of 723 human microRNAs was measured in laser capture microdissected epithelial tumours from 133 snap-frozen surgical colorectal specimens. Three well-known classification algorithms were used to derive candidate biomarkers for discriminating carcinomas from adenomas. Quantitative reverse-transcriptase PCR was then used to validate the candidates in an independent cohort of macrodissected formalin-fixed paraffin-embedded colorectal tissue samples from 91 surgical resections. The biomarkers were applied to differentiate carcinomas from high-grade intraepithelial neoplasms in 58 colonoscopy biopsy tissue samples with stromal invasion cells undetectable by microscopy. RESULTS: One classifier of 14 microRNAs was identified with a prediction accuracy of 94.1% for discriminating carcinomas from adenomas. In formalin-fixed paraffin-embedded surgical tissue samples, a combination of miR-375, miR-424 and miR-92a yielded an accuracy of 94% (AUC=0.968) in discriminating carcinomas from adenomas. This combination has been applied to differentiate carcinomas from high-grade intraepithelial neoplasms in colonoscopy biopsy tissues with an accuracy of 89% (AUC=0.918). CONCLUSIONS: This study has found a microRNA panel that accurately discriminates carcinomas from high-grade intraepithelial neoplasms in colonoscopy biopsy tissues. This microRNA panel has considerable clinical value in the early diagnosis and optimal surgical decision-making of colorectal cancer.

Effective normalization for copy number variation detection from whole genome sequencing.

BACKGROUND: Whole genome sequencing enables a high resolution view of the human genome and provides unique insights into genome structure at an unprecedented scale. There have been a number of tools to infer copy number variation in the genome. These tools, while validated, also include a number of parameters that are configurable to genome data being analyzed. These algorithms allow for normalization to account for individual and population-specific effects on individual genome CNV estimates but the impact of these changes on the estimated CNVs is not well characterized. We evaluate in detail the effect of normalization methodologies in two CNV algorithms FREEC and CNV-seq using whole genome sequencing data from 8 individuals spanning four populations. METHODS: We apply FREEC and CNV-seq to a sequencing data set consisting of 8 genomes. We use multiple configurations corresponding to different read-count normalization methodologies in FREEC, and statistically characterize the concordance of the CNV calls between FREEC configurations and the analogous output from CNV-seq. The normalization methodologies evaluated in FREEC are: GC content, mappability and control genome. We further stratify the concordance analysis within genic, non-genic, and a collection of validated variant regions. RESULTS: The GC content normalization methodology generates the highest number of altered copy number regions. Both mappability and control genome normalization reduce the total number and length of copy number regions. Mappability normalization yields Jaccard indices in the 0.07 - 0.3 range, whereas using a control genome normalization yields Jaccard index values around 0.4 with normalization based on GC content. The most critical impact of using mappability as a normalization factor is substantial reduction of deletion CNV calls. The output of another method based on control genome normalization, CNV-seq, resulted in comparable CNV call profiles, and substantial agreement in variable gene and CNV region calls. CONCLUSIONS: Choice of read-count normalization methodology has a substantial effect on CNV calls and the use of genomic mappability or an appropriately chosen control genome can optimize the output of CNV analysis.

Multi-scale Pan-cancer Integrative Analyses Identify the STAT3-VSIR Axis as a Key Immunosuppressive Mechanism in Head and Neck Cancer.

PURPOSE: VSIR is a novel immune checkpoint protein whose expression on tumor cells across cancers remains largely uncharacterized. Here we purposed to decode the pan-cancer biologic and clinical significance of VSIR overexpression in the tumor compartment. EXPERIMENTAL DESIGN: We performed multi-omics integrative analyses of 9,735 tumor samples to identify cancers with non-leukocytic expression of VSIR (VSIR High), followed by association with overall survival and immune cell infiltration levels. Orthogonal assessments of VSIR protein expression and lymphocytic infiltration were performed using quantitative immunofluorescence (QIF). RESULTS: Integrative modeling identified a subset of cancer types as being enriched for VSIR High tumors. VSIR High tumors were associated with significantly poorer overall survival in immunogenic ovarian serous adenocarcinoma (SA) and oral cavity squamous cell carcinoma (SCC). QIF assessments in an independent validation cohort confirmed overexpression of VSIR as being associated with poorer overall survival within immunogenic oral cavity SCC. VSIR overexpression was associated with lower CD4 helper T-cell infiltration in both ovarian SA and oral cavity SCC, but did not impact CD8 T-cell infiltration. VSIR overexpressing tumors in both cancer types exhibited significantly higher STAT3 signaling activity. Pharmacologic inhibition of STAT3 signaling resulted in dose-dependent reduction of VSIR expression in ovarian SA and oral cavity SCC cells. CONCLUSIONS: The STAT3-VSIR axis is a potentially significant immunomodulatory mechanism in oral cavity and ovarian cancers, whose activation is associated with poorer survival and an immune microenvironment marked by decreased CD4 helper T-cell activity. The role of VSIR as a tumor-intrinsic modulator of resistance to immunotherapy warrants further exploration.

Sexually dimorphic radiogenomic models identify distinct imaging and biological pathways that are prognostic of overall survival in glioblastoma.

BACKGROUND: Recent epidemiological studies have suggested that sexual dimorphism influences treatment response and prognostic outcome in glioblastoma (GBM). To this end, we sought to (i) identify distinct sex-specific radiomic phenotypes-from tumor subcompartments (peritumoral edema, enhancing tumor, and necrotic core) using pretreatment MRI scans-that are prognostic of overall survival (OS) in GBMs, and (ii) investigate radiogenomic associations of the MRI-based phenotypes with corresponding transcriptomic data, to identify the signaling pathways that drive sex-specific tumor biology and treatment response in GBM. METHODS: In a retrospective setting, 313 GBM patients (male = 196, female = 117) were curated from multiple institutions for radiomic analysis, where 130 were used for training and independently validated on a cohort of 183 patients. For the radiogenomic analysis, 147 GBM patients (male = 94, female = 53) were used, with 125 patients in training and 22 cases for independent validation. RESULTS: Cox regression models of radiomic features from gadolinium T1-weighted MRI allowed for developing more precise prognostic models, when trained separately on male and female cohorts. Our radiogenomic analysis revealed higher expression of Laws energy features that capture spots and ripple-like patterns (representative of increased heterogeneity) from the enhancing tumor region, as well as aggressive biological processes of cell adhesion and angiogenesis to be more enriched in the "high-risk" group of poor OS in the male population. In contrast, higher expressions of Laws energy features (which detect levels and edges) from the necrotic core with significant involvement of immune related signaling pathways was observed in the "low-risk" group of the female population. CONCLUSIONS: Sexually dimorphic radiogenomic models could help risk-stratify GBM patients for personalized treatment decisions.

ICAM1 initiates CTC cluster formation and trans-endothelial migration in lung metastasis of breast cancer.

Circulating tumor cell (CTC) clusters mediate metastasis at a higher efficiency and are associated with lower overall survival in breast cancer compared to single cells. Combining single-cell RNA sequencing and protein analyses, here we report the profiles of primary tumor cells and lung metastases of triple-negative breast cancer (TNBC). ICAM1 expression increases by 200-fold in the lung metastases of three TNBC patient-derived xenografts (PDXs). Depletion of ICAM1 abrogates lung colonization of TNBC cells by inhibiting homotypic tumor cell-tumor cell cluster formation. Machine learning-based algorithms and mutagenesis analyses identify ICAM1 regions responsible for homophilic ICAM1-ICAM1 interactions, thereby directing homotypic tumor cell clustering, as well as heterotypic tumor-endothelial adhesion for trans-endothelial migration. Moreover, ICAM1 promotes metastasis by activating cellular pathways related to cell cycle and stemness. Finally, blocking ICAM1 interactions significantly inhibits CTC cluster formation, tumor cell transendothelial migration, and lung metastasis. Therefore, ICAM1 can serve as a novel therapeutic target for metastasis initiation of TNBC.

Understanding breast cancer disparities-a multi-scale challenge.

Despite convergence of overall breast cancer incidence rates between European American (EA) and African American (AA) women, disparities in mortality persist. The factors contributing to differences in mortality rates across population groups remain controversial and range from population genetics to sociodemographic influences. This review explores the complex multi-factorial nature of tumor-intrinsic and -extrinsic factors that impact the biology and clinical outcomes of breast cancer patients. In addition to summarizing the current state of breast cancer disparities research, we also motivate the development of integrative multi-scale approaches involving interdisciplinary teams to tackle this complex clinical challenge.

Untangling a complex web: Computational analyses of tumor molecular profiles to decode driver mechanisms.

Genome-scale studies focusing on molecular profiling of cancers across tissue types have revealed a plethora of aberrations across the genomic, transcriptomic, and epigenomic scales. The significant molecular heterogeneity across individual tumors even within the same tissue context complicates decoding the key etiologic mechanisms of this disease. Furthermore, it is increasingly likely that biologic mechanisms underlying the pathobiology of cancer involve multiple molecular entities interacting across functional scales. This has motivated the development of computational approaches that integrate molecular measurements with prior biological knowledge in increasingly intricate ways to enable the discovery of driver genomic aberrations across cancers. Here, we review diverse methodological approaches that have powered significant advances in our understanding of the genomic underpinnings of cancer at the cohort and at the individual tumor scales. We outline the key advances and challenges in the computational discovery of cancer mechanisms while motivating the development of systems biology approaches to comprehensively decode the biologic drivers of this complex disease.

LIN9 and NEK2 Are Core Regulators of Mitotic Fidelity That Can Be Therapeutically Targeted to Overcome Taxane Resistance.

A significant therapeutic challenge for patients with cancer is resistance to chemotherapies such as taxanes. Overexpression of LIN9, a transcriptional regulator of cell-cycle progression, occurs in 65% of patients with triple-negative breast cancer (TNBC), a disease commonly treated with these drugs. Here, we report that LIN9 is further elevated with acquisition of taxane resistance. Inhibiting LIN9 genetically or by suppressing its expression with a global BET inhibitor restored taxane sensitivity by inducing mitotic progression errors and apoptosis. While sustained LIN9 is necessary to maintain taxane resistance, there are no inhibitors that directly repress its function. Hence, we sought to discover a druggable downstream transcriptional target of LIN9. Using a computational approach, we identified NIMA-related kinase 2 (NEK2), a regulator of centrosome separation that is also elevated in taxane-resistant cells. High expression of NEK2 was predictive of low survival rates in patients who had residual disease following treatment with taxanes plus an anthracycline, suggesting a role for this kinase in modulating taxane sensitivity. Like LIN9, genetic or pharmacologic blockade of NEK2 activity in the presence of paclitaxel synergistically induced mitotic abnormalities in nearly 100% of cells and completely restored sensitivity to paclitaxel, in vitro. In addition, suppressing NEK2 activity with two distinct small molecules potentiated taxane response in multiple in vivo models of TNBC, including a patient-derived xenograft, without inducing toxicity. These data demonstrate that the LIN9/NEK2 pathway is a therapeutically targetable mediator of taxane resistance that can be leveraged to improve response to this core chemotherapy. SIGNIFICANCE: Resistance to chemotherapy is a major hurdle for treating patients with cancer. Combining NEK2 inhibitors with taxanes may be a viable approach for improving patient outcomes by enhancing mitotic defects induced by taxanes alone.

Radiogenomic-Based Survival Risk Stratification of Tumor Habitat on Gd-T1w MRI Is Associated with Biological Processes in Glioblastoma.

PURPOSE: To (i) create a survival risk score using radiomic features from the tumor habitat on routine MRI to predict progression-free survival (PFS) in glioblastoma and (ii) obtain a biological basis for these prognostic radiomic features, by studying their radiogenomic associations with molecular signaling pathways. EXPERIMENTAL DESIGN: Two hundred three patients with pretreatment Gd-T1w, T2w, T2w-FLAIR MRI were obtained from 3 cohorts: The Cancer Imaging Archive (TCIA; n = 130), Ivy GAP (n = 32), and Cleveland Clinic (n = 41). Gene-expression profiles of corresponding patients were obtained for TCIA cohort. For every study, following expert segmentation of tumor subcompartments (necrotic core, enhancing tumor, peritumoral edema), 936 3D radiomic features were extracted from each subcompartment across all MRI protocols. Using Cox regression model, radiomic risk score (RRS) was developed for every protocol to predict PFS on the training cohort (n = 130) and evaluated on the holdout cohort (n = 73). Further, Gene Ontology and single-sample gene set enrichment analysis were used to identify specific molecular signaling pathway networks associated with RRS features. RESULTS: Twenty-five radiomic features from the tumor habitat yielded the RRS. A combination of RRS with clinical (age and gender) and molecular features (MGMT and IDH status) resulted in a concordance index of 0.81 (P < 0.0001) on training and 0.84 (P = 0.03) on the test set. Radiogenomic analysis revealed associations of RRS features with signaling pathways for cell differentiation, cell adhesion, and angiogenesis, which contribute to chemoresistance in GBM. CONCLUSIONS: Our findings suggest that prognostic radiomic features from routine Gd-T1w MRI may also be significantly associated with key biological processes that affect response to chemotherapy in GBM.