Wu, Gang

Abstract Although sequencing costs have steadily decreased with advances in technology, they remain high for large scale studies. The design of traditional individual-disease sequencing studies is either case only or cases with relatively few controls, resulting in potential loss of statistical power for discovery of disease associated genes. Here we show that for a given number of sequenced cases, a large control sample size is critical to maximize power for rare variant burden analysis. Furthermore, we have developed an end-to-end workflow based tool (CoCoRV-nf) to facilitate the use of external biobank sequence resources as controls. The modules include consistent variant QC, variant annotation, ancestry population prediction, and gene based burden analysis using summary genotype information, and combined analysis from multiple independent results. The tool supports exomes and genomes from gnomAD and All of Us as controls with preprocessed datasets. We apply the tool in two rare neurological diseases: amyotrophic lateral sclerosis and neuroblastoma. For each disease, two case cohorts are paired with gnomAD and All of Us data, respectively, followed by a combined analysis. Not only did we recapture known genes, but also, we identified new candidate genes for both diseases. By leveraging multiple large external biobank sequence data, we demonstrate the feasibility of using our tool to maximize statistical power to identify new disease predisposition genes.

PURPOSE Recent studies reveal that 5%-18% of children with cancer harbor pathogenic variants in known cancer-predisposing genes. However, DNA damage repair (DDR) genes, which are frequently somatically altered in pediatric tumors, have not been systematically examined as a source of novel cancer-predisposing signals. METHODS To address this gap, we interrogated 189 DDR genes for presence of germline predisposing variants (PV) among 5,993 childhood cancer cases and 14,477 adult noncancer controls (discovery cohort). PV were determined using a tiered approach incorporating ClinVar annotations, InterVar classification, and in silico tools (REVEL, CADD, and MetaSVM). Using logistic and firth regression, we identified genes with PV statistically enriched in the germline of children with tumors and replicated findings among 1,497 additional childhood cancer cases across three independent cohorts. RESULTS Analysis across all cases with cancer revealed enrichment of TP53 PV. Cancer-specific analyses confirmed known associations including germline TP53 PV in adrenocortical carcinoma, high-grade glioma (HGG), and medulloblastoma (MB), PMS2 in HGG and non-Hodgkin lymphoma (NHL), MLH1 in HGG, BRCA2 in NHL, and BARD1 in neuroblastoma. In addition, four novel associations were uncovered, including BRCA1 in ependymoma, SPIDR in HGG, SMC5 in MB, and SMARCAL1 in osteosarcoma (OS). Importantly, the SMARCAL1 :OS association was significant in the discovery (6/230, 2.6%, false discovery rate [FDR] logistic = 0.0189) as well as all three replication cohorts (Childhood Cancer Survivor Study: 8/275, 2.9%; P Fisher < .0001; Cancer Predisposition Syndrome-German Childhood Cancer Registry: 4/135, 3%, P Fisher = .002; Individualized Therapy for Relapsed Malignancies in Childhood: 4/217, 1.8%, P Fisher = .012). The remaining wild-type SMARCAL1 allele was deleted in three of four OS tumors with available data. CONCLUSION Our study confirms the relevance DDR genetic variation in pediatric cancer risk and establishes SMARCAL1 as a novel OS predisposing gene, providing insights into tumor biology and creating opportunities to optimize care for patients with this challenging tumor.

PURPOSE Although cure rates for childhood acute lymphoblastic leukemia (ALL) exceed 90%, ALL remains a leading cause of cancer death in children. Half of relapses arise in children initially classified with standard-risk (SR) disease. MATERIALS AND METHODS To identify genomic determinants of relapse in children with SR ALL, we performed genome and transcriptome sequencing of diagnostic and remission samples of children with SR (n = 1,381) or high-risk B-ALL with favorable cytogenetic features (n = 115) enrolled on Children's Oncology Group trials. We used a case-control study design analyzing 439 patients who relapsed and 1,057 who remained in complete remission for at least 5 years. RESULTS Genomic subtype was associated with relapse, which occurred in approximately 50% of cases of PAX5 -altered ALL (odds ratio [OR], 3.31 [95% CI, 2.17 to 5.03]; P = 3.18 × 10 –8 ). Within high-hyperdiploid ALL, gain of chromosome 10 with disomy of chromosome 7 was associated with favorable outcome (OR, 0.27 [95% CI, 0.17 to 0.42]; P = 8.02 × 10 –10 ; St Jude Children's Research Hospital validation cohort: OR, 0.22 [95% CI, 0.05 to 0.80]; P = .009), and disomy of chromosomes 10 and 17 with gain of chromosome 6 was associated with relapse (OR, 7.16 [95% CI, 2.63 to 21.51]; P = 2.19 × 10 –5 ; validation cohort: OR, 21.32 [95% CI, 3.62 to 119.30]; P = .0004). Genomic alterations were associated with relapse in a subtype-dependent manner, including alterations of INO80 in ETV6::RUNX1 ALL, IKZF1 , and CREBBP in high-hyperdiploid ALL and FHIT in BCR::ABL1 -like ALL. Genomic alterations were also associated with the presence of minimal residual disease, including NRAS and CREBBP in high-hyperdiploid ALL. CONCLUSION Genetic subtype, patterns of aneuploidy, and secondary genomic alterations determine risk of relapse in childhood ALL. Comprehensive genomic analysis is required for optimal risk stratification.

Abstract BACKGROUND Medulloblastoma (MB) is a highly malignant childhood cerebellar tumor comprising a collection of molecularly and clinically distinct entities. Group 3 and Group 4-MB subgroups originate from overlapping progenitor pools of the developing rhombic lip when neuronal differentiation hierarchies are subverted through somatic alterations. Mutations suspected to deregulate the activity of chromatin-modifying genes are pervasive in Group 3 and Group 4-MB. However, the molecular consequences of these alterations remain largely undefined. METHODS We characterized the chromatin landscape of fifty-two primary MBs using chromatin immunoprecipitation followed by sequencing (ChIP-seq) for six histone modifications and performed multi-modal integration with sample-matched somatic lesions, genome-wide DNA methylation, and transcriptomic profiles. RESULTS Analysis of differential chromatin states across MB subgroups using a Hidden Markov Model revealed a significant enrichment of the bivalent enhancer state (EnhBiv; marked by coincident H3K4me1 and H3K27me3) in Group 4-MB which localized to promoters of neurodevelopmental genes and transcription factors. To functionally dissect the pathogenetic significance of the enriched EnhBiv state in Group 4-MB, we sought to identify key regulators that confer this state. Integrative bioinformatics coupled with CUT&RUN analysis revealed significant enrichment of KDM2B, a histone lysine demethylase overexpressed in Group 4-MB, in gene promoters marked by the EnhBiv state. CRISPR gene targeting or targeted protein degradation of KDM2B selectively suppressed growth of Group 3 and Group 4-MB models in vitro and in vivo. Mechanistically, inhibition of KDM2B disrupted chromatin bivalency by blocking PRC2 complex recruitment and reducing H3K27me3 deposition, profoundly derepressing PRC2 targets involved in neuronal differentiation. CONCLUSIONS Our comprehensive characterization of the MB chromatin landscape in a large cohort of primary tumors provides unprecedented insights into the epigenetic basis of MB subgroups, implicating novel mechanisms of tumorigenesis, and nominating KDM2B as a selective dependency in Group 3/4-MB that should be prioritized as an attractive therapeutic target.

Introduction: ECOG ACRIN E1910 is a randomized phase III trial that showed that adults with newly diagnosed BCR::ABL1 negative acute lymphoblastic leukemia (ALL) who become MRD negative (&lt;0.01%) after induction chemo who receive blinatumomab with conventional chemotherapy (chemo) have improved survival compared with those who received chemo only (Litzow et al, Blood (2022) 140: Supplement 2, LBA-1). However, not all pts were able to receive all four planned cycles of blinatumomab in consolidation. In this report we assessed outcomes of pts in the blinatumomab arm of the trial who received all 4 cycles of blinatumomab compared to those who received 1-2 cycles. Methods: Patient between the ages of 30 and 70 with newly diagnosed BCR::ABL1 negative B-lineage ALL were enrolled and initially received 2.5 months of combination induction chemo utilizing a BFM-like regimen adapted from the E2993/UKALLXII clinical trial, with pegaspargase (after induction only for patients &gt;=55 years of age) and addition of rituximab for CD20 positive patients. After remission induction (step 1) pts in morphologic complete remission (CR/CRi) received an intensification course of high dose methotrexate with pegaspargase for CNS prophylaxis (step 2). At the conclusion of step 2, remission and MRD status were determined centrally by 6-color flow cytometry with MRD negativity defined as &lt;0.01%. In the primary analysis subset, MRD negative patients were randomized to receive an additional 4cycles of consolidation chemo or 2 cycles of blinatumomab at 28 mcg/day (flat dose) for 28 days each cycle followed by 3 cycles of consolidation chemo, another 4-week cycle of blinatumomab (3rd cycle of blinatumomab) followed by an additional cycle of chemo and then a 4th cycle of blinatumomab (step 3). Following completion of consolidation chemo +/- blinatumomab, pts were given 2.5 years of POMP maintenance therapy timed from the start of the intensification cycle (step 4). Estimates of OS were calculated using the Kaplan-Meier method. Comparison of OS between treatment arms were conducted using the stratified log-rank test and Cox proportional-hazards model with age, CD20 status, rituximab use, and whether pts intend to receive HSCT or not as stratification factors. As the number of cycles needed to optimize outcome has not been studied, an unplanned landmark analysis was used to compare OS of pts in the blinatumomab arm who received all 4 cycles o those who received 1-2 cycles. Time 0 was chosen as 9 months post step 3 randomization (the time that patients were supposed to complete 4 cycles of blina). Four pts who received blinatumomab but died within 9 months post step 3 randomization were therefore excluded from this analysis. Results: The study was activated on December 17, 2013 and closed to enrollment on October 15, 2019. 772 pts were screened for the trial and 488 were enrolled on step 1 induction therapy. The median age of the pts was 51 years (range 30-70). 224 MRD negative pts were randomized, 112 pts to each arm. In the blinatumomab arm 12 pts received 1 cycle (11%), 32 pts received 2 cycles (29%), 4 pts received 3 cycles (4%) and 63 pts received 4 cycles (57%). The OS of pts who received 1-2 cycles of blinatumomab was compared to control (Figure 1) The difference was not significant (hazard ratio 0.62, 95% CI 0.28 to 1.34, p=0.22). Fig 2 compares the survival of those who received 1-2 cycles to those who received 4 cycles ( (HR: 0.39, 95% CI 0.12 to 1.16, p=0.076). Conclusion: The addition of blinatumomab to consolidation chemotherapy resulted in a significantly better overall survival compared to chemotherapy alone in pts with newly diagnosed B-lineage ALL who were MRD negative after intensification. This represents a new standard of care for pts aged 30-70 years with BCR::ABL1 negative ALL. The optimal dose and number of cycles is however unknown. In this unplanned subgroup analysis we demonstrate that a survival benefit can only be confirmed in patients who receive the intended 4 cycles of blinatumomab during consolidation. Further studies will be needed to determine if fewer cycles can provide a similar benefit.

Abstract Medulloblastoma (MB) is a highly malignant childhood brain tumor comprising a collection of molecularly and clinically distinct entities. Mutations in chromatin-modifying genes distribute across subgroups and account for nearly half of all MBs. However, the molecular consequences of chromatin modifier alterations remain poorly understood. Herein, we characterize the chromatin landscape of fifty-two primary MBs representing each MB subgroup using chromatin immunoprecipitation followed by sequencing (ChIP-seq) for six histone marks and perform multi-layered integration with somatic genetics, genome-wide DNA methylation, and transcriptomic datasets extracted from the same tumors. Analysis of differential chromatin states across MB subgroups identified significant enrichment of a bivalent enhancer state (EnhBiv) in Group 4 MB. Group 4-specific EnhBivs mapped to promoters of key neurodevelopmental genes and transcription factors and was associated with increased DNA methylation and reduced expression. In accordance with higher EnhBivs coverage in Group 4 MB, we found that KDM2B, a H3K4me3 and H3K36me2 demethylase, is highly expressed in Group 4 MBs and patient-derived xenografts (PDXs). KDM2B knockout mediated by CRISPR targeting suppressed growth of both Group 3 MB cell lines in vitro and Group 4 MB PDXs in vivo but had no effect on control SHH-MB cells, suggesting that KDM2B is a novel selective dependency in MB. Knockout of KDM2B in MB cells induced profound upregulation of PRC2 targets, genes regulated by bivalent histone marks, as well as genes involved in neuronal differentiation. Further biochemical and molecular studies demonstrated that KDM2B interacts with PRC2 complex subunits and occupies gene promoters associated with Group 4-specific EnhBiv chromatin. In summary, our comprehensive characterization of the MB chromatin landscape in a large cohort of primary tumors provides novel insights into the epigenetic basis of MB and implicates KDM2B as a novel Group 4 MB dependency that should be prioritized as a therapeutic target in this subgroup.

Abstract Lineage-ambiguous leukemias are high-risk malignancies of poorly understood genetic basis. Here, we describe a distinct subgroup of acute leukemia with expression of myeloid, T lymphoid, and stem cell markers driven by aberrant allele-specific deregulation of BCL11B, a master transcription factor responsible for thymic T-lineage commitment and specification. Mechanistically, this deregulation was driven by chromosomal rearrangements that juxtapose BCL11B to superenhancers active in hematopoietic progenitors, or focal amplifications that generate a superenhancer from a noncoding element distal to BCL11B. Chromatin conformation analyses demonstrated long-range interactions of rearranged enhancers with the expressed BCL11B allele and association of BCL11B with activated hematopoietic progenitor cell cis-regulatory elements, suggesting BCL11B is aberrantly co-opted into a gene regulatory network that drives transformation by maintaining a progenitor state. These data support a role for ectopic BCL11B expression in primitive hematopoietic cells mediated by enhancer hijacking as an oncogenic driver of human lineage-ambiguous leukemia. Significance: Lineage-ambiguous leukemias pose significant diagnostic and therapeutic challenges due to a poorly understood molecular and cellular basis. We identify oncogenic deregulation of BCL11B driven by diverse structural alterations, including de novo superenhancer generation, as the driving feature of a subset of lineage-ambiguous leukemias that transcend current diagnostic boundaries. This article is highlighted in the In This Issue feature, p. 2659