Kollmeyer, Thomas

ABSTRACT IDH-mutant astrocytomas maintain telomeres through the alternative lengthening of telomeres (ALT) pathway, producing extreme inter-arm telomere length heterogeneity, yet how this heterogeneity shapes structural genome evolution remains unknown. Using Oxford Nanopore long-read sequencing of 20 IDH-mutant astrocytomas, we profiled structural variants (SVs), copy number variants, extrachromosomal DNA (ecDNA) and measured allele-specific telomere lengths from individual long reads. We identified pervasive complex rearrangements, including chromothripsis and foldback events consistent with breakage-fusion-bridge cycles, and widespread ecDNAs. SV breakpoints were enriched at telomeric and centromeric regions regardless of local telomere length, revealing constitutive structural fragility. Arm-level telomere length analysis uncovered a dual-mode model: arms with short telomeres preferentially harbored breakage-associated events, while arms with long ALT-maintained telomeres were enriched for ecDNA and amplification-associated events. These findings identify chromosome-arm-specific telomere length as a determinant of structural genome evolution in ALT-driven tumors.

Abstract Glioblastoma (GBM) diagnosis integrates morphology and tumor molecular alterations. Patients with classic morphology are designated as histologic GBM (hGBM). Patients without classic morphology harboring TERT promoter mutation (TERTp), or epidermal growth factor receptor amplification (EGFRamp), or combined gain of chromosome 7 and loss of 10 (+7/-10) designate molecular GBM (mGBM). We compared radiographic features, genomic alterations, and clinical outcome between hGBM and mGBM. 109 mGBM Mayo Clinic cases diagnosed between 1993-2021 were matched to 101 cGBM by age, sex, year of diagnosis, and TERT promoter (TERTp) status. Clinical and radiographic data abstraction was blinded. Progression-free and overall survival were compared across groups in the subset of patients who received standard of care. Contrast enhancement was present on index MRI (MRI at diagnosis) for 98% (96/98) hGBM and 50% (54/107) mGBM (p<0.0001) patients. mGBM tumors harbored fewer total copy number alterations (median of 16.5 vs 20, p=0.029) and were less likely to have CDKN2A/B homozygous deletion (56% vs 77%, p=0.005). Among the 87 mGBM and 82 hGBM receiving standard of care, median overall survival was longer for mGBM versus hGBM (2.24 vs 1.62 years, p=0.016). mGBM patients with +7/-10 (HR=2.60, p=0.035) or CDKN2A/B (HR=3.61, p=0.002) had poorer overall survival compared to mGBM patients without the corresponding alteration. On index MRI mGBM patients without contrast enhancement had better overall survival compared to hGBM (HR=0.50; p=0.003) whereas there was no survival difference between mGBM with enhancement and hGBM. mGBM with enhancement on index MRI likely represent undersampled, or early / evolving hGBM. mGBM patients with non-enhancing MRI, younger than 60, with long symptom profile including seizures, and some combination of non-anaplastic morphology, TERTp only, no homozygous deletion of CDKN2A/B, and low copy number complexity may represent other IDHwt glioma entities. Caution is advised in counseling, management, discussing prognosis, and inclusion in clinical trials.

Abstract The 2021 WHO classification redefined glioblastoma IDH-wildtype (GBM), integrating histopathological and tumor molecular features. Diffuse astrocytic glioma lacking histologic features of necrosis and microvascular proliferation are classified as (molecular) GBM if the tumor is IDH- and H3-wildtype and has one or more additional tumor molecular alterations: TERT promoter mutation, EGFR amplification, or +7/-10. Due to the absence of histopathological features of classic GBM and presentation without enhancement on MRI in approximately 50% of cases, molecular GBM frequently resemble low-grade IDH-mutant gliomas. We hypothesized that MRI-based deep learning models could distinguish IDH-mutant glioma from molecular GBM. We retrospectively identified 97 individuals with molecular GBM, who were matched to 97 individuals with IDH-mutant glioma by age, sex, year of diagnosis, and enhancement on MRI. Of the 97 IDH-mutant gliomas, 51% were 1p/19q codeleted, 42% grade 2, 47% grade 3, and 10% grade 4. Preoperative postcontrast T1-weighted (T1C) and T2-weighted (T2) MRI sequences were used to train a 3D-DenseNet121 model. The model was validated on an independent prospective set of 150 IDH-mutant gliomas and 71 molecular GBM. The MRI-based deep learning model achieved an area under the receiver operator curve (AUC) of 0.75 (95% CI: 0.68-0.82) in the prospective validation set. Sex-stratified AUCs were 0.72 (95% CI: 0.59-0.85) for females and 0.76 (95% CI: 0.68-0.85) for males. Age stratified AUCs were 0.71 (95% CI: 0.60-0.83) for patients under 55 years of age and 0.64 (95% CI: 0.50-0.79) for patients older than 55. These results suggest that discriminative MRI features of IDH-mutant glioma and molecular GBM exist and can be identified from deep learning models. However, radiographic ambiguity exists and underscores the need for further optimization of the deep learning architecture or by integrating additional non-invasive modalities. For example, future work entails integrating germline variants into the model, including rs55705857, which is a causal variant for IDH-mutant glioma. With further validation and integration of multi-modal data, a noninvasive tool may support more personalized decision-making.

Abstract IDH-wildtype diffuse astrocytic gliomas with lower-grade histological features alongside specific genomic alterations (e.g., TERT promoter mutation, EGFR amplification, gain Chr7/loss Chr10) are presently diagnosed as glioblastoma, IDH-wildtype (GBM) under WHO 2021 criteria. Previous analyses categorized these into “early/evolving” or “surgically undersampled” GBMs. However, it remains unclear whether these tumors share identical clinical and prognostic features with histologic GBM. Data from IDH-wt molecular GBMs across six institutions were aggregated. Astrocytic IDH-wt tumors of any grade that lacked necrosis or microvascular proliferation on initial histopathology yet exhibited molecular alterations consistent with WHO 2021 GBM diagnosis, were included. Recursive partitioning analysis (RPA) was used to categorize patients into risk groups based on overall survival (OS). Among the 264 patients included, the median age was 60 years (IQR: 51-67), with 162 (61%) being male. Additionally, 210 (80%) were diagnosed since 2016, and 216 (82%) received chemoradiation plus adjuvant temozolomide (GBM-SOC). The median follow-up period was 3.9 years (95% CI: 3.4, 4.6). Identified univariate risk factors for worse survival included older age, biopsy-only status, EGFR amplification, absence of initial adjuvant treatment, low post-operative KPS, and tumor histologic grade of 3. RPA divided patients into four distinct survival groups (p<0.0001). The “best” risk group comprised lower grade (<=2) tumors (n=58, median OS (mOS): 40 months (95% CI: 33, NA)). The “second best” group included Grade 3 patients receiving GBM-SOC without steroids during initial RT (n=139, mOS: 25 months (95% CI: 22, 28)). The “second worst” group encompassed Grade 3 patients who were on steroids during initial RT while receiving GBM-SOC (n=38, mOS: 15 months (95% CI: 13, 20)). The “worst” group contained Grade 3 tumors that had not undergone SOC (n=29, mOS: 9 months (95% CI: 4, 26)). These data show that histologic features and treatment are strong indicators of patient outcomes.

Abstract Telomere length variation is a hallmark of genomic instability, however, its role in shaping the structural variant (SV) landscape of diffuse astrocytoma remains poorly understood. In this study, we employed Oxford Nanopore long-read sequencing to investigate the association between SVs and chromosome arm-specific telomere length in a clinical patient cohort consisting of n = 20 individuals diagnosed with an IDH mutant astrocytoma. Oxford Nanopore long-read sequencing was performed to a depth of ~20-40× coverage. SVs were identified using Sniffles, copy number alterations were assessed using DNAcopy and chromosome arm-specific telomere length was estimated with Telogator. Statistical analyses were conducted to assess associations between telomere dynamics, structural variation, and copy number changes. Astrocytoma genomes exhibited a broad spectrum of structural variants, including deletions and insertions, with significant enrichment in the short arms of acrocentric chromosomes. A subset of samples demonstrated markedly shorter telomeres relative to the overall cohort, accompanied by an increased burden of structural variants, particularly duplications and inversions, suggesting enhanced genomic instability. Notably, a prominent accumulation of deletion and insertion SVs was observed near chromosome ends and centromeres, highlighting the (sub)telomeric and peri(centromeric) regions as SV hotspots. Statistical modeling revealed a significant association between telomere length and arm-specific variation, including telomeric deletions, inversions, and copy number losses. In summary, our findings reveal a previously underappreciated link between SVs and telomere length in astrocytomas. Genomic rearrangements caused by telomere shortening may represent critical drivers of genome instability and could serve as novel biomarkers for disease progression.

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.

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.