We aimed to evaluate the potential of radiomics as an imaging biomarker for glioblastoma (GBM) patients and explore the molecular rationale behind radiomics using a radio-genomics approach. A total of 144 primary GBM patients were included in this study (training cohort). Using multi-parametric MR images, radiomics features were extracted from multi-habitats of the tumor. We applied Cox-LASSO algorithm to build a survival prediction model, which we validated using an independent validation cohort. GBM patients were consensus clustered to reveal inherent phenotypic subtypes. GBM patients were successfully stratified by the radiomics risk score, a weighted sum of radiomics features, corroborating the potential of radiomics as a prognostic biomarker. Using consensus clustering, we identified three distinct subtypes which significantly differed in the prognosis (“heterogenous enhancing”, “rim-enhancing necrotic”, and “cystic” subtypes). Transcriptomic traits enriched in individual subtypes were in accordance with imaging phenotypes summarized by radiomics. For example, rim-enhancing necrotic subtype was well described by radiomics profiling (T2 autocorrelation and flat shape) and highlighted by the inflammatory genomic signatures, which well correlated to its phenotypic peculiarity (necrosis). This study showed that imaging subtypes derived from radiomics successfully recapitulated the genomic underpinnings of GBMs and thereby confirmed the feasibility of radiomics as an imaging biomarker for GBM patients with comprehensible biologic annotation.
Clonal evolution drives cancer progression and therapeutic resistance. Recent studies have revealed divergent longitudinal trajectories in gliomas, but early molecular features steering posttreatment cancer evolution remain unclear. Here, we collected sequencing and clinical data of initial-recurrent tumor pairs from 544 adult diffuse gliomas and performed multivariate analysis to identify early molecular predictors of tumor evolution in three diffuse glioma subtypes. We found that CDKN2A deletion at initial diagnosis preceded tumor necrosis and microvascular proliferation that occur at later stages of IDH-mutant glioma. Ki67 expression at diagnosis was positively correlated with acquiring hypermutation at recurrence in the IDH–wild-type glioma. In all glioma subtypes, MYC gain or MYC- target activation at diagnosis was associated with treatment-induced hypermutation at recurrence. To predict glioma evolution, we constructed CELLO2 (Cancer EvoLution for LOngitudinal data version 2), a machine learning model integrating features at diagnosis to forecast hypermutation and progression after treatment. CELLO2 successfully stratified patients into subgroups with distinct prognoses and identified a high-risk patient group featured by MYC gain with worse post-progression survival, from the low-grade IDH-mutant-noncodel subtype. We then performed chronic temozolomide-induction experiments in glioma cell lines and isogenic patient-derived gliomaspheres and demonstrated that MYC drives temozolomide resistance by promoting hypermutation. Mechanistically, we demonstrated that, by binding to open chromatin and transcriptionally active genomic regions, c-MYC increases the vulnerability of key mismatch repair genes to treatment-induced mutagenesis, thus triggering hypermutation. This study reveals early predictors of cancer evolution under therapy and provides a resource for precision oncology targeting cancer dynamics in diffuse gliomas.
Glioblastoma (GBM) is the most lethal brain cancer with a dismal prognosis. Stem-like GBM cells (GSCs) are a major driver of GBM propagation and recurrence; thus, understanding the molecular mechanisms that promote GSCs may lead to effective therapeutic approaches. Through in vitro clonogenic growth-based assays, we determined mitogenic activities of the ligand molecules that are implicated in neural development. We have identified that semaphorin 3A (Sema3A), originally known as an axon guidance molecule in the CNS, promotes clonogenic growth of GBM cells but not normal neural progenitor cells (NPCs). Mechanistically, Sema3A binds to its receptor neuropilin-1 (NRP1) and facilitates an interaction between NRP1 and TGF-β receptor 1 (TGF-βR1), which in turn leads to activation of canonical TGF-β signaling in both GSCs and NPCs. TGF-β signaling enhances self-renewal and survival of GBM tumors through induction of key stem cell factors, but it evokes cytostatic responses in NPCs. Blockage of the Sema3A/NRP1 axis via shRNA-mediated knockdown of Sema3A or NRP1 impeded clonogenic growth and TGF-β pathway activity in GSCs and inhibited tumor growth in vivo. Taken together, these findings suggest that the Sema3A/NRP1/TGF-βR1 signaling axis is a critical regulator of GSC propagation and a potential therapeutic target for GBM.
Purpose Even though pazopanib, a multitargeted tyrosine kinase inhibitor, has been approved for refractory soft tissue sarcoma (STS), little is known about the molecular determinants of the response to pazopanib. We performed integrative molecular characterization to identify potential predictors of pazopanib efficacy.Materials and Methods We obtained fresh pre-treatment tumor tissue from 35 patients with advanced STS receiving pazopanib-based treatment. Among those, 18 (51.4%) received pazopanib monotherapy, and the remaining 17 (48.6%) received pazopanib in combination with durvalumab, programmed death-ligand 1 blockade. Whole-exome and transcriptome sequencing were performed for each tumor and patient germline DNA.Results Of the 35 patients receiving pazopanib-based treatment, nine achieved a partial response (PR), resulting in an objective response rate (ORR) of 27.3%, and the median progression-free survival (PFS) was 6.0 months. Patients with CDK4 amplification (copy ratio tumor to normal > 2) exhibited shorter PFS (3.7 vs. 7.9 months, p=2.09×10–4) and a poorer response (ORR; 0% vs. 33.3%) compared to those without a gene amplification (copy ratio ≤ 2). Moreover, non-responders demonstrated transcriptional activation of CDK4 via DNA amplification, resulting in cell cycle activation. In the durvalumab combination cohort, seven of the 17 patients (41.2%) achieved a PR, and gene expression analysis revealed that durvalumab responders exhibited high immune/stromal cell infiltration, mainly comprising natural killer cells, compared to non-responders as well as increased expression of CD19, a B-cell marker.Conclusion Despite the limitation of heterogeneity in the study population and treatment, we identified possible molecular predictors of pazopanib efficacy that can be employed in future clinical trials aimed at evaluating therapeutic strategies.