Aberrant DNA repair reveals a vulnerability in histone H3.3-mutant brain tumors.

Pediatric high-grade gliomas (pHGG) are devastating and incurable brain tumors with recurrent mutations in histone H3.3. These mutations promote oncogenesis by dysregulating gene expression through alterations of histone modifications. We identify aberrant DNA repair as an independent mechanism, which fosters genome instability in H3.3 mutant pHGG, and opens new therapeutic options. The two most frequent H3.3 mutations in pHGG, K27M and G34R, drive aberrant repair of replication-associated damage by non-homologous end joining (NHEJ). Aberrant NHEJ is mediated by the DNA repair enzyme polynucleotide kinase 3'-phosphatase (PNKP), which shows increased association with mutant H3.3 at damaged replication forks. PNKP sustains the proliferation of cells bearing H3.3 mutations, thus conferring a molecular vulnerability, specific to mutant cells, with potential for therapeutic targeting.

Clinical and Genomic Predictors of Adverse Events in Newly Diagnosed Glioblastoma.

PURPOSE: Adverse clinical events cause significant morbidity in patients with GBM (GBM). We examined whether genomic alterations were associated with AE (AE) in patients with GBM. EXPERIMENTAL DESIGN: We identified adults with histologically confirmed IDH-wild-type GBM with targeted next-generation sequencing (OncoPanel) at Dana Farber Cancer Institute from 2013 to 2019. Seizure at presentation, lymphopenia, thromboembolic events, pseudoprogression, and early progression (within 6 months of diagnosis) were identified as AE. The biologic function of genetic variants was categorized as loss-of-function (LoF), no change in function, or gain-of-function (GoF) using a somatic tumor mutation knowledge base (OncoKB) and consensus protein function predictions. Associations between functional genomic alterations and AE were examined using univariate logistic regressions and multivariable regressions adjusted for additional clinical predictors. RESULTS: Our study included 470 patients diagnosed with GBM who met the study criteria. We focused on 105 genes that had sequencing data available for ≥ 90% of the patients and were altered in ≥10% of the cohort. Following false-discovery rate (FDR) correction and multivariable adjustment, the TP53, RB1, IGF1R, and DIS3 LoF alterations were associated with lower odds of seizures, while EGFR, SMARCA4, GNA11, BRD4, and TCF3 GoF and SETD2 LoF alterations were associated with higher odds of seizures. For all other AE of interest, no significant associations were found with genomic alterations following FDR correction. CONCLUSIONS: Genomic biomarkers based on functional variant analysis of a routine clinical panel may help identify AE in GBM, particularly seizures. Identifying these risk factors could improve the management of patients through better supportive care and consideration of prophylactic therapies.

Genomic sequence context differs between germline and somatic structural variants allowing for their differentiation in tumor samples without paired normals.

There is currently no method to distinguish between germline and somatic structural variants (SVs) in tumor samples that lack a matched normal sample. In this study, we analyzed several features of germline and somatic SVs from a cohort of 974 patients from The Cancer Genome Atlas (TCGA). We identified a total of 21 features that differed significantly between germline and somatic SVs. Several of the germline SV features were associated with each other, as were several of the somatic SV features. We also found that these associations differed between the germline and somatic classes, for example, we found that somatic inversions were more likely to be longer events than their germline counterparts. Using these features we trained a support vector machine (SVM) classifier on 555,849 TCGA SVs to computationally distinguish germline from somatic SVs in the absence of a matched normal. This classifier had an ROC curve AUC of 0.984 when tested on an independent test set of 277,925 TCGA SVs. In this dataset, we achieved a positive predictive value (PPV) of 0.81 for an SV called somatic by the classifier being truly somatic. We further tested the classifier on a separate set of 7,623 SVs from pediatric high-grade gliomas (pHGG). In this non-TCGA cohort, our classifier achieved a PPV of 0.828, showing robust performance across datasets.

Double-strand break repair-associated intragenic deletions and tandem duplications suggest the architecture of the repair replication fork.

Double-strand break (DSB) repair is associated with a 1000-fold increase in mutations compared to normal replication of the same sequences. In budding yeast, repair of an HO endonuclease-induced DSB at the MATα locus can be repaired by using a homologous, heterochromatic HMR::Kl-URA3 donor harboring a transcriptionally silenced URA3 gene, resulting in a MAT::URA3 (Ura+) repair product where URA3 is expressed. Repair-associated ura3- mutations can be selected by resistance to 5-fluoroorotic acid (FOA). Using this system, we find that a major class of mutations are -1 deletions, almost always in homonucleotide runs, but there are few +1 insertions. In contrast, +1 and -1 insertions in homonucleotide runs are nearly equal among spontaneous mutations. Approximately 10% of repair-associated mutations are interchromosomal template switches (ICTS), even though the K. lactis URA3 sequence embedded in HMR is only 72% identical with S. cerevisiae ura3-52 sequences on a different chromosome. ICTS events begin and end in regions of short microhomology, averaging 7 bp. Long microhomologies are favored, but some ICTS junctions are as short as 2 bp. Both repair-associated intragenic deletions (IDs) and tandem duplications (TDs) are recovered, with junctions sharing short stretches of, on average, 6 bp of microhomology. Intragenic deletions are more than 5 times more frequent than TDs. IDs have a mean length of 60 bp, but, surprisingly there are almost no deletions shorter than 25 bp. In contrast, TDs average only 12 bp. The usage of microhomologies among intragenic deletions is not strongly influenced by the degree of adjacent homeology. Together, these data provide a picture of the structure of the repair replication fork. We suggest that IDs and TDs occur within the migrating D-loop in which DNA polymerase δ copies the template, where the 3' end of a partly copied new DNA strand can dissociate and anneal with a single-stranded region of microhomology that lies either in front or behind the 3' end, within the open structure of a migrating D-loop. Our data suggest that ~100 bp ahead of the polymerase is "open," but that part of the repair replication apparatus remains bound in the 25 bp ahead of the newly copied DNA, preventing annealing. In contrast, the template region behind the polymerase appears to be rapidly reannealed, limiting template switching to a very short region.

Inaugural Results of the Individualized Screening Trial of Innovative Glioblastoma Therapy: A Phase II Platform Trial for Newly Diagnosed Glioblastoma Using Bayesian Adaptive Randomization.

PURPOSE: The Individualized Screening Trial of Innovative Glioblastoma Therapy (INSIGhT) is a phase II platform trial that uses response adaptive randomization and genomic profiling to efficiently identify novel therapies for phase III testing. Three initial experimental arms (abemaciclib [a cyclin-dependent kinase [CDK]4/6 inhibitor], neratinib [an epidermal growth factor receptor [EGFR]/human epidermal growth factor receptor 2 inhibitor], and CC-115 [a deoxyribonucleic acid-dependent protein kinase/mammalian target of rapamycin inhibitor]) were simultaneously evaluated against a common control arm. We report the results for each arm and examine the feasibility and conduct of the adaptive platform design. PATIENTS AND METHODS: Patients with newly diagnosed O6-methylguanine-DNA methyltransferase-unmethylated glioblastoma were eligible if they had tumor genotyping to identify prespecified biomarker subpopulations of dominant glioblastoma signaling pathways (EGFR, phosphatidylinositol 3-kinase, and CDK). Initial random assignment was 1:1:1:1 between control (radiation therapy and temozolomide) and the experimental arms. Subsequent Bayesian adaptive randomization was incorporated on the basis of biomarker-specific progression-free survival (PFS) data. The primary end point was overall survival (OS), and one-sided P values are reported. The trial is registered with ClinicalTrials.gov (identifier: NCT02977780). RESULTS: Two hundred thirty-seven patients were treated (71 control; 73 abemaciclib; 81 neratinib; 12 CC-115) in years 2017-2021. Abemaciclib and neratinib were well tolerated, but CC-115 was associated with ≥ grade 3 treatment-related toxicity in 58% of patients. PFS was significantly longer with abemaciclib (hazard ratio [HR], 0.72; 95% CI, 0.49 to 1.06; one-sided P = .046) and neratinib (HR, 0.72; 95% CI, 0.50 to 1.02; one-sided P = .033) relative to the control arm but there was no PFS benefit with CC-115 (one-sided P = .523). None of the experimental therapies demonstrated a significant OS benefit (P > .05). CONCLUSION: The INSIGhT design enabled efficient simultaneous testing of three experimental agents using a shared control arm and adaptive randomization. Two investigational arms had superior PFS compared with the control arm, but none demonstrated an OS benefit. The INSIGhT design may promote improved and more efficient therapeutic discovery in glioblastoma. New arms have been added to the trial.

Cancer aneuploidies are shaped primarily by effects on tumour fitness.

Aneuploidies-whole-chromosome or whole-arm imbalances-are the most prevalent alteration in cancer genomes1,2. However, it is still debated whether their prevalence is due to selection or ease of generation as passenger events1,2. Here we developed a method, BISCUT, that identifies loci subject to fitness advantages or disadvantages by interrogating length distributions of telomere- or centromere-bounded copy-number events. These loci were significantly enriched for known cancer driver genes, including genes not detected through analysis of focal copy-number events, and were often lineage specific. BISCUT identified the helicase-encoding gene WRN as a haploinsufficient tumour-suppressor gene on chromosome 8p, which is supported by several lines of evidence. We also formally quantified the role of selection and mechanical biases in driving aneuploidy, finding that rates of arm-level copy-number alterations are most highly correlated with their effects on cellular fitness1,2. These results provide insight into the driving forces behind aneuploidy and its contribution to tumorigenesis.

ALK Amplification and Rearrangements Are Recurrent Targetable Events in Congenital and Adult Glioblastoma.

PURPOSE: Anaplastic lymphoma kinase (ALK) aberrations have been identified in pediatric-type infant gliomas, but their occurrence across age groups, functional effects, and treatment response has not been broadly established. EXPERIMENTAL DESIGN: We performed a comprehensive analysis of ALK expression and genomic aberrations in both newly generated and retrospective data from 371 glioblastomas (156 adult, 205 infant/pediatric, and 10 congenital) with in vitro and in vivo validation of aberrations. RESULTS: ALK aberrations at the protein or genomic level were detected in 12% of gliomas (45/371) in a wide age range (0-80 years). Recurrent as well as novel ALK fusions (LRRFIP1-ALK, DCTN1-ALK, PRKD3-ALK) were present in 50% (5/10) of congenital/infant, 1.4% (3/205) of pediatric, and 1.9% (3/156) of adult GBMs. ALK fusions were present as the only candidate driver in congenital/infant GBMs and were sometimes focally amplified. In contrast, adult ALK fusions co-occurred with other oncogenic drivers. No activating ALK mutations were identified in any age group. Novel and recurrent ALK rearrangements promoted STAT3 and ERK1/2 pathways and transformation in vitro and in vivo. ALK-fused GBM cellular and mouse models were responsive to ALK inhibitors, including in patient cells derived from a congenital GBM. Relevant to the treatment of infant gliomas, we showed that ALK protein appears minimally expressed in the forebrain at perinatal stages, and no gross effects on perinatal brain development were seen in pregnant mice treated with the ALK inhibitor ceritinib. CONCLUSIONS: These findings support use of brain-penetrant ALK inhibitors in clinical trials across infant, pediatric, and adult GBMs. See related commentary by Mack and Bertrand, p. 2567.

Investigative needle core biopsies for multi-omics in Glioblastoma.

Glioblastoma (GBM) is a primary brain cancer with an abysmal prognosis and few effective therapies. The ability to investigate the tumor microenvironment before and during treatment would greatly enhance both understanding of disease response and progression, as well as the delivery and impact of therapeutics. Stereotactic biopsies are a routine surgical procedure performed primarily for diagnostic histopathologic purposes. The role of investigative biopsies - tissue sampling for the purpose of understanding tumor microenvironmental responses to treatment using integrated multi-modal molecular analyses ('Multi-omics") has yet to be defined. Secondly, it is unknown whether comparatively small tissue samples from brain biopsies can yield sufficient information with such methods. Here we adapt stereotactic needle core biopsy tissue in two separate patients. In the first patient with recurrent GBM we performed highly resolved multi-omics analysis methods including single cell RNA sequencing, spatial-transcriptomics, metabolomics, proteomics, phosphoproteomics, T-cell clonotype analysis, and MHC Class I immunopeptidomics from biopsy tissue that was obtained from a single procedure. In a second patient we analyzed multi-regional core biopsies to decipher spatial and genomic variance. We also investigated the utility of stereotactic biopsies as a method for generating patient derived xenograft models in a separate patient cohort. Dataset integration across modalities showed good correspondence between spatial modalities, highlighted immune cell associated metabolic pathways and revealed poor correlation between RNA expression and the tumor MHC Class I immunopeptidome. In conclusion, stereotactic needle biopsy cores are of sufficient quality to generate multi-omics data, provide data rich insight into a patient's disease process and tumor immune microenvironment and can be of value in evaluating treatment responses. One sentence summary: Integrative multi-omics analysis of stereotactic needle core biopsies in glioblastoma.

Tangent normalization for somatic copy-number inference in cancer genome analysis.

MOTIVATION: Somatic copy-number alterations (SCNAs) play an important role in cancer development. Systematic noise in sequencing and array data present a significant challenge to the inference of SCNAs for cancer genome analyses. As part of The Cancer Genome Atlas, the Broad Institute Genome Characterization Center developed the Tangent normalization method to generate copy-number profiles using data from single-nucleotide polymorphism (SNP) arrays and whole-exome sequencing (WES) technologies for over 10 000 pairs of tumors and matched normal samples. Here, we describe the Tangent method, which uses a unique linear combination of normal samples as a reference for each tumor sample, to subtract systematic errors that vary across samples. We also describe a modification of Tangent, called Pseudo-Tangent, which enables denoising through comparisons between tumor profiles when few normal samples are available. RESULTS: Tangent normalization substantially increases signal-to-noise ratios (SNRs) compared to conventional normalization methods in both SNP array and WES analyses. Tangent and Pseudo-Tangent normalizations improve the SNR by reducing noise with minimal effect on signal and exceed the contribution of other steps in the analysis such as choice of segmentation algorithm. Tangent and Pseudo-Tangent are broadly applicable and enable more accurate inference of SCNAs from DNA sequencing and array data. AVAILABILITY AND IMPLEMENTATION: Tangent is available at https://github.com/broadinstitute/tangent and as a Docker image (https://hub.docker.com/r/broadinstitute/tangent). Tangent is also the normalization method for the copy-number pipeline in Genome Analysis Toolkit 4 (GATK4). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

VRK1 as a synthetic lethal target in VRK2 promoter-methylated cancers of the nervous system.

Collateral lethality occurs when loss of a gene/protein renders cancer cells dependent on its remaining paralog. Combining genome-scale CRISPR/Cas9 loss-of-function screens with RNA sequencing in over 900 cancer cell lines, we found that cancers of nervous system lineage, including adult and pediatric gliomas and neuroblastomas, required the nuclear kinase vaccinia-related kinase 1 (VRK1) for their survival in vivo. VRK1 dependency was inversely correlated with expression of its paralog VRK2. VRK2 knockout sensitized cells to VRK1 loss, and conversely, VRK2 overexpression increased cell fitness in the setting of VRK1 loss. DNA methylation of the VRK2 promoter was associated with low VRK2 expression in human neuroblastomas and adult and pediatric gliomas. Mechanistically, depletion of VRK1 reduced barrier-to-autointegration factor phosphorylation during mitosis, resulting in DNA damage and apoptosis. Together, these studies identify VRK1 as a synthetic lethal target in VRK2 promoter-methylated adult and pediatric gliomas and neuroblastomas.

Analysis of germline-driven ancestry-associated gene expression in cancers.

Differential mRNA expression between ancestry groups can be explained by both genetic and environmental factors. We outline a computational workflow to determine the extent to which germline genetic variation explains cancer-specific molecular differences across ancestry groups. Using multi-omics datasets from The Cancer Genome Atlas (TCGA), we enumerate ancestry-informative markers colocalized with cancer-type-specific expression quantitative trait loci (e-QTLs) at ancestry-associated genes. This approach is generalizable to other settings with paired germline genotyping and mRNA expression data for a multi-ethnic cohort. For complete details on the use and execution of this protocol, please refer to Carrot-Zhang et al. (2020), Robertson et al. (2021), and Sayaman et al. (2021).

Patterns of structural variation define prostate cancer across disease states.

The complex genomic landscape of prostate cancer evolves across disease states under therapeutic pressure directed toward inhibiting androgen receptor (AR) signaling. While significantly altered genes in prostate cancer have been extensively defined, there have been fewer systematic analyses of how structural variation shapes the genomic landscape of this disease across disease states. We uniformly characterized structural alterations across 531 localized and 143 metastatic prostate cancers profiled by whole genome sequencing, 125 metastatic samples of which were also profiled via whole transcriptome sequencing. We observed distinct significantly recurrent breakpoints in localized and metastatic castration-resistant prostate cancers (mCRPC), with pervasive alterations in noncoding regions flanking the AR, MYC, FOXA1, and LSAMP genes enriched in mCRPC and TMPRSS2-ERG rearrangements enriched in localized prostate cancer. We defined 9 subclasses of mCRPC based on signatures of structural variation, each associated with distinct genetic features and clinical outcomes. Our results comprehensively define patterns of structural variation in prostate cancer and identify clinically actionable subgroups based on whole genome profiling.

Structural variants shape driver combinations and outcomes in pediatric high-grade glioma.

We analyzed the contributions of structural variants (SVs) to gliomagenesis across 179 pediatric high-grade gliomas (pHGGs). The most recurrent SVs targeted MYC isoforms and receptor tyrosine kinases (RTKs), including an SV amplifying a MYC enhancer in 12% of diffuse midline gliomas (DMG), indicating an underappreciated role for MYC in pHGG. SV signature analysis revealed that tumors with simple signatures were TP53 wild type (TP53WT) but showed alterations in TP53 pathway members PPM1D and MDM4. Complex signatures were associated with direct aberrations in TP53, CDKN2A and RB1 early in tumor evolution and with later-occurring extrachromosomal amplicons. All pHGGs exhibited at least one simple-SV signature, but complex-SV signatures were primarily restricted to subsets of H3.3K27M DMGs and hemispheric pHGGs. Importantly, DMGs with complex-SV signatures were associated with shorter overall survival independent of histone mutation and TP53 status. These data provide insight into the impact of SVs on gliomagenesis and the mechanisms that shape them.

The genomic landscape of dysembryoplastic neuroepithelial tumours and a comprehensive analysis of recurrent cases.

AIMS: Dysembryoplastic neuroepithelial tumour (DNT) is a glioneuronal tumour that is challenging to diagnose, with a wide spectrum of histological features. Three histopathological patterns have been described: specific DNTs (both the simple form and the complex form) comprising the specific glioneuronal element, and also the non-specific/diffuse form which lacks it, and has unclear phenotype-genotype correlations with numerous differential diagnoses. METHODS: We used targeted methods (immunohistochemistry, fluorescence in situ hybridisation and targeted sequencing) and large-scale genomic methodologies including DNA methylation profiling to perform an integrative analysis to better characterise a large retrospective cohort of 82 DNTs, enriched for tumours that showed progression on imaging. RESULTS: We confirmed that specific DNTs are characterised by a single driver event with a high frequency of FGFR1 variants. However, a subset of DNA methylation-confirmed DNTs harbour alternative genomic alterations to FGFR1 duplication/mutation. We also demonstrated that a subset of DNTs sharing the same FGFR1 alterations can show in situ progression. In contrast to the specific forms, "non-specific/diffuse DNTs" corresponded to a heterogeneous molecular group encompassing diverse, newly-described, molecularly distinct entities. CONCLUSIONS: Specific DNT is a homogeneous group of tumours sharing characteristics of paediatric low-grade gliomas: a quiet genome with a recurrent genomic alteration in the RAS-MAPK signalling pathway, a distinct DNA methylation profile and a good prognosis but showing progression in some cases. The "non-specific/diffuse DNTs" subgroup encompasses various recently described histomolecular entities, such as PLNTY and diffuse astrocytoma, MYB or MYBL1 altered.