Drucker, Kristen

Abstract BACKGROUND IDH-mutant gliomas are diffusely infiltrating, slow-growing brain tumors that inevitably undergo malignant progression to grade IV astrocytoma - a highly aggressive and uniformly fatal disease with dismal prognosis. Despite extensive research, the molecular mechanisms underlying this transformation remain poorly understood. While alterations in canonical driver genes such as CDKN2A, PTEN, and PDGFRA have been implicated, the broader landscape of genetic drivers involved in IDH-mutant glioma progression is largely uncharacterized. In contrast to many other malignancies, IDH-mutant gliomas acquire relatively few somatic point mutations during progression - typically in the range of 30 to 40 - suggesting that alternative mechanisms may be driving tumor evolution. Notably, these tumors exhibit widespread chromosomal instability, resulting in a heterogeneous and highly prevalent pattern of copy number alterations (CNAs). This observation raises the possibility that gene dosage changes due to CNAs play a critical role in tumor progression. MATERIAL AND METHODS To investigate this hypothesis, we performed integrative analysis of IDH-mutant glioma cohorts from the Mayo Clinic and TCGA datasets (~1,500 patients), identifying twelve highly recurrent CNAs significantly associated with patient outcomes. We refer to these as malignant IDH-associated CNAs (mIDH CNAs). We hypothesized that these recurrent CNAs harbor driver genes that confer selective growth advantages, thereby contributing to their recurrence and prognostic relevance. To systematically assess the functional impact of genes within these regions, we developed a novel in vivo CRISPR-based screening platform, termed CRISPR-Knock Out and Activation Linked Assay (CRISPR-KOALA). This approach utilizes stereotaxic delivery of lentiviral gRNA libraries - targeting mouse orthologs of ~1,800 genes within recurrently deleted regions and ~1,200 genes within amplified regions - into the brains of a genetically engineered mouse model of IDH-mutant glioma (rs557mut;IDHmut;P53KO;AtrxKO;Cas9GFP: Yanchus et. al Science 2022). RESULTS Mice injected with mIDH CNA-targeting libraries developed tumors significantly more rapidly than control cohorts. Deconvolution of the CRISPR-KOALA screen identified several putative driver genes, including regulators of the NOTCH signaling pathway, lipid metabolism modulators such as THEM6, and the NMDA receptor modulator GRINA - implicating dysregulation of metabolism, immune response, and glutamatergic signaling in the malignant transformation of IDH-mutant gliomas. CONCLUSION Together, these findings provide critical insights into the genetic mechanisms driving the progression of IDH-mutant gliomas and identify new candidate genes and pathways for potential therapeutic targeting.

Background: Tumefactive demyelination (TD) is a rare variant of multiple sclerosis (MS) characterized by tumor-like lesions that often require aggressive management. Genome-wide association studies (GWAS) identified variants associated with MS; similar analyses in TD are lacking. Objective: A GWAS was performed to identify variants associated with TD. Methods: The case–control study included 142 TD cases and 293 controls. TD patients were required to have a demyelinating event and magnetic resonance imaging (MRI) showing one or more lesions. Controls were patients without a neurologic or systemic inflammatory disease or cancer. Logistic regression was used to compare cases versus controls for each variant; age, sex, and principal components were included as covariates. A p -value threshold of 5 × 10 −8 was GWAS significant and 5 × 10 −6 nominally significant. A polygenic risk score (PRS) was compared across TD and controls. Results: Variants on chromosome 14 (rs117797734, p = 2.06 × 10 −11 , odds ratio (OR) = 13.14) and chromosome 6 (most significant rs6936540, p = 5.5 × 10 −7 , OR = 2.61) near DCBLD1 were significant. Seven non-MHC and two MHC variants associated with MS were associated with TD. The PRS was significantly higher in TD versus controls. Conclusion: We identified novel regions associated with TD, demonstrating the importance of performing GWAS in homogeneous subtypes of MS. Further validation and functional experiments are necessary.

Abstract Background The CCDC26 germline variant rs55705857 is causal for development of IDH mutant (IDHmut) adult glioma. However, ~60% of IDHmut patients do not carry the rs55705857 risk allele. We aimed to identify variants associated with developing IDHmut glioma among patients that do not have the rs55705857 risk allele and to further understand development of IDHwt glioma. Methods We used three datasets that included 1216 IDHmut and 1442 IDHwt glioma patients and a case–case design to perform genome-wide association (GWAS) analyses comparing IDHmut versus IDHwt glioma. Analyses were performed overall and stratified by rs55705857 genotype and sex. Multivariable logistic regression and regression trees were used to develop models to predict IDH status using germline variants, age, and contrast enhancement on MRI. Results Three regions were identified comparing IDHmut versus IDHwt: rs55705857 (meta P-value [P] = 1.35 × 10−43), PHLDB1 (rs7125115, P = 3.46 × 10−17), and D2HGDH (rs71430382, P = 2.43 × 10−12). When analyzing only patients that do not have the rs55705857 risk allele, PHLDB1 (rs7125115, P = 1.73 × 10−13) and D2HGDH (rs71430382, P = 8.86 × 10−10) were identified. Among patients who have the rs55705857 risk allele, four variants between ROBO1 and ROBO2 (rs4680975, P = 4.65 × 10−7) increased the likelihood of having an IDHwt tumor. Tumor expression of ROBO1 was associated with rs4680975 genotype in IDHwt patients that have the rs55705857 risk allele (P = 0.034). Multivariable logistic analysis demonstrated that rs55705857, rs71430382 (D2HGDH), and age predicted IDH mutation status. Conclusions To understand the development of adult glioma, we demonstrate that D2HGDH and PHLDB1 should be prioritized for functional studies in IDHmut tumors. The ROBO1 region warrants further investigation in IDHwt tumors.

Isocitrate dehydrogenase (IDH)–mutant low-grade gliomas are slow-growing brain tumors that frequently progress to aggressive high-grade gliomas that have dismal outcomes. In a recent study, Wu and colleagues provide critical insights into the mechanisms underlying malignant progression by analyzing single-cell gene expression and chromatin accessibility across different tumor grades. Their findings support a two-phase model: In early stages, tumors are primarily driven by oligodendrocyte precursor–like cells and epigenetic alterations that silence tumor suppressors like CDKN2A and activate oncogenes such as PDGFRA. As the disease advances, the tumors become sustained by more proliferative neural precursor–like cells, in which genetic alterations, including PDGFRA, MYCN, and CDK4 amplifications and CDKN2A/B deletion, drive tumor progression. The study further highlights a dynamic regulation of IFN signaling during progression. In low-grade IDH-mutant gliomas, IFN responses are suppressed through epigenetic hypermethylation, which can be reversed with DNA methyltransferase 1 inhibitors or IDH inhibitors, leading to reactivation of the IFN pathway. In contrast, higher grade gliomas evade IFN signaling through genetic deletions of IFN genes. These findings emphasize a broader epigenetic-to-genetic shift in oncogenic regulation that drives glioma progression, provide a valuable framework for understanding the transition from indolent tumors to lethal malignancies, and have implications for therapy and clinical management.

Abstract Low-grade gliomas (LGG) are slowly growing brain cancers that ultimately undergo a malignant progression to lethal secondary glioblastomas with a dismal prognosis. Despite intensive research, very little progress has been made in our understanding of this transformation. Hence, outside of known alternations in well-characterized driver genes CDKN2A, PTEN and PDGFRA, genetic drivers that underlie LGG transformation are largely unknown. Contrary to other types of cancer, LGGs exhibit very few driver mutations (30-40 mutations) as they progress. Interestingly, LGGs also exhibit a highly prevalent and heterogenous pattern of copy-number-alterations (CNAs) suggesting that they might be driven by changes in gene dosages that result from widespread chromosomal instability observed within these tumors. Through analysis of the human Mayo Clinic and TCGA LGG datasets (~1500 patients), we have identified eleven highly recurrent CNAs affecting LGG patients, many of which correlate with patient survival. We hypothesized that these CNAs contain driver genes that confer a growth advantage resulting in the highly recurrent CNA pattern and prognostic associations. To investigate the tumor suppressive potential of the ~1800 genes located in recurrently deleted regions and the oncogenic potential of the ~1200 genes located in recurrently amplified regions, our lab developed a novel in vivo ‘CRISPR-Knock Out and Activation Linked Assay’ (CRISPR-KOALA). Stereotaxic delivery of arm-level CNA lentiviral gRNA libraries targeting the mouse orthologs of genes in CNAs into the brain of our animal model (rs557mut;IDHmut;P53KO;AtrxKO;Cas9GFP: Yanchus et. al Science 2022) resulted in rapid formation of tumors compared to animals injected with control libraries. Deconvoluting the screen revealed enrichment of immune modulators such as Them6 as well as a modulator of NMDA receptors (Grina) implicating immune system and glutamate dysregulation in LGG malignant transformation. Taken together, our study provides insights into mechanisms underlying LGG progression which may serve to inform intervention strategies for this lethal disease.

Abstract Acquired tumor alterations are part of the WHO diagnosis of glioma. However, germline predisposition to these alterations is unknown. GWAS in European and Asian populations identified germline variants in 30 regions, many in or near genes that also acquire alterations, suggesting that interactions between inherited and acquired genetics are involved in glioma development. To evaluate such interactions, we used 3686 glioma cases and 1889 controls. A case-control design and multinomial logistic regression were used to identify germline variants associated with development of TERT-promoter mutant gliomas, generating two interesting observations. First, TERT germline variants were associated with predisposition to IDHwt glioma that have a TERT promoter mutation, but not with IDH-mutant codeleted gliomas that have a TERT promoter mutation. Second, two novel regions were observed in this more homogeneous subgroup of TERT mutated tumors: variants in TXNDC11 (OR=2.66, p=5.7x10-8) and RUVBL2 (OR=3.1, p=5.8x10-8). High tumor expression of TXNDC11 is associated with poor prognosis in glioma and RUVBL2 interacts with TERT to affect telomerase activity. Similar analyses are being performed on additional clinically relevant acquired alterations. In parallel, case-case GWAS analyses were performed to evaluate germline variants associated with other tumor alterations. Comparing IDH-mutant versus IDHwt stratified by rs55705857 genotype, the most significant variants were PHLDB1 (p=1.2x10-13) and D2HGDH (p=2.1x10-10), highlighting the importance of these two variants independent of CCDC26 for the development of IDH-mutant glioma. As a parallel approach, we utilized recursive partitioning and regression trees (rpart) to classify IDH mutation status using the 30 known variants. The top four branches of the regression tree included CCDC26, D2HGDH, PHLDB1, and MAML2, as well as two Asian variants (RAB27A, CYP4F12). This highlights the importance of utilizing variants discovered in non-European populations in analyses. Overall, understanding these interactions are necessary to develop therapeutics, stratify clinical trials, and to develop relevant tumor models.

Abstract Histologic features, presence of CDKN2A/B homozygous deletion and methylation subtype are strong determinants of outcome and of therapy in IDH-mutant gliomas. We were interested in learning if there are additional markers of outcome in these tumors. To address this question, we assembled a cohort of 963 IDH-mutant gliomas with overall survival data (median follow-up alive 2.7 years, median follow-up deceased 3.9 years). Array-based copy number data, targeted NGS, germline genotyping data and Illumina EPIC array methylation data were available on 872, 226, 542, and 401 patients, respectively. Morphologic type (oligodendroglioma versus astrocytoma, grade (2/3 vs 4) and methylation subtype (A_IDH_LG vs A_IDH_HG) were all significantly associated with survival (p <0.001). There was no significant difference in survival between grade 2 vs 3 IDH-mutant astrocytomas (p=0.53 adjusting for age and sex). However, subjects with the risk allele (G) for rs55705857 had a better overall survival than subjects with the non-risk allele (A) in IDH-mutant oligodendroglioma (AG/GG vs AA, median overall survival 21.6 vs 15.7 years, p=0.048 after adjusting for grade, age, and sex). Conversely, rs55705857 was associated with a markedly worse survival in grade 2 IDH-mutant astrocytomas (AG/GG vs AA, median overall survival 8.6 vs 19.6 years, p=0.007 after adjusting for age and sex). Despite similar sample sizes and number of deaths, there was no significant difference in overall survival by rs55705857 genotype in grade 3 or 4 IDH-mutant astrocytomas. Upon further functional analyses, TLX1 and GALNT17 were found to be differentially methylated in IDH1-mutant astrocytomas and oligodendrogliomas, respectively, when stratified by rs55705857 risk allele. This is the first observation of a germline characteristic (e.g. rs55705857) predicting survival in low-grade IDH-mutant astrocytoma. While this observation requires validation, it suggests that germline genotype may be practice changing and inform results of trials such as CATNON and CODEL.

Abstract INTRODUCTION Prior to 2021, (classic) glioblastoma (GBM) required absence of IDH mutation and histologic features of hypercellularity, nuclear atypia, microvascular proliferation (MVP), and necrosis. Absence of histologic features can be a result of undersampling, with molecular, radiologic, and clinical features consistent with GBM. Accordingly, WHO 2021 defined molecular glioblastoma as diffuse glioma IDHwt, H3wt, with one or more of MVP, necrosis, TERT promoter mutation, EGFR gene amplification (EGFRamp), and/or +7/−10. We hypothesized that there is clinical and radiographic heterogeneity among tumors defined as molecular GBM. METHODS We identified 121 patients that met the criteria for molecular GBM. 121 Classic GBM, defined as IDHwt with hypercellularity, MVP or necrosis, were matched to molecular GBM cases by age, sex, and diagnosis date. All cases were reviewed for relevant histologic features and molecular features. Examined clinical variables included: first symptom, time from first symptom to date of pathological diagnosis, date of first MRI, tumor location, radiographic features of enhancement, diffusion or mass effect, time to progression and overall survival. Clinical, radiologic, and molecular features were compared between classic and molecular GBM. RESULTS Currently,160 cases have been reviewed. Preliminary results show that there is a longer time from first symptom to diagnosis for molecular GBM. Molecular GBM is significantly more likely to be discovered coincidentally, and without enhancement or radiographic necrosis. Classic GBM is significantly more likely to present with focal neurologic deficit as first symptom, regional mass effect on index imaging and have monosomy 10. CONCLUSION While pathologic undersampling accounts for a significant proportion of molecular GBM, our data suggest that there are significant clinical, radiographic, and outcome differences between classical and molecular GBM. Final analysis of the entire cohort will help better define this clinical spectrum.

The CRISPR-Cas9 system has enabled researchers to precisely modify/edit the sequence of a genome. A typical editing experiment consists of two steps: (1) editing cultured cells; (2) cell cloning and selection of clones with and without intended edit, presumed to be isogenic. The application of CRISPR-Cas9 system may result in off-target edits, whereas cloning will reveal culture-acquired mutations. We analyzed the extent of the former and the latter by whole genome sequencing in three experiments involving separate genomic loci and conducted by three independent laboratories. In all experiments we hardly found any off-target edits, whereas detecting hundreds to thousands of single nucleotide mutations unique to each clone after relatively short culture of 10–20 passages. Notably, clones also differed in copy number alterations (CNAs) that were several kb to several mb in size and represented the largest source of genomic divergence among clones. We suggest that screening of clones for mutations and CNAs acquired in culture is a necessary step to allow correct interpretation of DNA editing experiments. Furthermore, since culture associated mutations are inevitable, we propose that experiments involving derivation of clonal lines should compare a mix of multiple unedited lines and a mix of multiple edited lines.

Abstract BACKGROUND The topologically associated domain (TAD) on 8q24 surrounding MYC is important for many types of cancer. The region is also important in the formation of IDH-mutant glioma; the risk allele rs55705857 being associated with the development of these tumors. The rs55705857 risk allele significantly increased the formation of tumors in an IDH-mutant mouse model. Cerebral organoids have been found to recapitulate the early development of the brain. We hypothesized that the rs55705857 risk allele may alter the phenotype of cerebral organoids. METHODS Isogenic induced pluripotent stem cells (iPSCs) were developed with (n=3) and without (n=2) the rs55705857 risk allele from the parental cell line PGP-1 (GM23338). These cells were cultured and differentiated into unpatterned cerebral organoids and observed throughout differentiation. Organoids were fixed, sectioned, and stained by H&E and immunofluorescence. RESULTS Organoids containing the risk allele were smaller throughout development, starting at the embryoid body stage and continuing up to over 75 days in culture, compared to two isogenic non-risk lines. At 4 weeks the risk allele containing organoids also contained a larger number of small rosette-like neuroepithelia, instead of the more continuous neuroepithelia observed in the non-risk allele organoids. At 75 days, the organoids with the risk allele maintained a larger proportion of Sox2 positive cells, which are enriched in the tight clusters of small cells with round nuclei. CONCLUSIONS The risk allele rs55705857 in the MYC TAD at 8q24 is altering the development of unpatterned cerebral organoids. This difference is visible as early as embryoid body development and persists through neuroepithelial development. Further studies of these differences may help our understanding of how rs55705857 accelerates IDH-mutant glioma development.