CELLO: a longitudinal data analysis toolbox untangling cancer evolution

Genetic and microenvironmental evolution of colorectal liver metastases under chemotherapy

Drug resistance limits the efficacy of chemotherapy for colorectal cancer liver metastasis (CRLM). However, the evolution of CRLM during drug treatment remains poorly elucidated. Multi-omics and treatment response data from 115 samples of 49 patients with CRLM undergoing bevacizumab (BVZ)-based chemotherapy show little difference in genomic alterations in 92% of cases, while remarkable differences are observed at the transcriptomic level. By decoupling intrinsic and acquired resistance, we find that hepatocyte and myeloid cell infiltration contribute to 38.5% and 23.1% of acquired resistance, respectively. Importantly, SMAD4 mutations and chr20q copy-number gain are associated with intrinsic chemoresistance. Gene interference experiments suggest that SMAD4R361H/C mutations confer BVZ and 5-fluorouracil (5-FU) resistance through STAT3 signaling. Notably, supplementing BVZ and 5-FU with the STAT3 inhibitor GB201 restores therapeutic efficacy in SMAD4R361H/C cancer cells. Our study uncovers the evolutionary dynamics of CRLM and its microenvironment during treatment and offers strategies to overcome drug resistance.

Identifying predictors of glioma evolution from longitudinal sequencing

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.

MGMT genomic rearrangements contribute to chemotherapy resistance in gliomas

Temozolomide (TMZ) is an oral alkylating agent used for the treatment of glioblastoma and is now becoming a chemotherapeutic option in patients diagnosed with high-risk low-grade gliomas. The O-6-methylguanine-DNA methyltransferase (MGMT) is responsible for the direct repair of the main TMZ-induced toxic DNA adduct, the O6-Methylguanine lesion. MGMT promoter hypermethylation is currently the only known biomarker for TMZ response in glioblastoma patients. Here we show that a subset of recurrent gliomas carries MGMT genomic rearrangements that lead to MGMT overexpression, independently from changes in its promoter methylation. By leveraging the CRISPR/Cas9 technology we generated some of these MGMT rearrangements in glioma cells and demonstrated that the MGMT genomic rearrangements contribute to TMZ resistance both in vitro and in vivo. Lastly, we showed that such fusions can be detected in tumor-derived exosomes and could potentially represent an early detection marker of tumor recurrence in a subset of patients treated with TMZ.