The Open Pediatric Cancer Project.

Background: In 2019, the Open Pediatric Brain Tumor Atlas (OpenPBTA) was created as a global, collaborative open-science initiative to genomically characterize 1,074 pediatric brain tumors and 22 patient-derived cell lines. Here, we extend the OpenPBTA to create the Open Pediatric Cancer (OpenPedCan) Project, a harmonized open-source multi-omic dataset from 6,112 pediatric cancer patients with 7,096 tumor events across more than 100 histologies. Combined with RNA-Seq from the Genotype-Tissue Expression (GTEx) and The Cancer Genome Atlas (TCGA), OpenPedCan contains nearly 48,000 total biospecimens (24,002 tumor and 23,893 normal specimens). Findings: We utilized Gabriella Miller Kids First (GMKF) workflows to harmonize WGS, WXS, RNA-seq, and Targeted Sequencing datasets to include somatic SNVs, InDels, CNVs, SVs, RNA expression, fusions, and splice variants. We integrated summarized CPTAC whole cell proteomics and phospho-proteomics data, miRNA-Seq data, and have developed a methylation array harmonization workflow to include m-values, beta-vales, and copy number calls. OpenPedCan contains reproducible, dockerized workflows in GitHub, CAVATICA, and Amazon Web Services (AWS) to deliver harmonized and processed data from over 60 scalable modules which can be leveraged both locally and on AWS. The processed data are released in a versioned manner and accessible through CAVATICA or AWS S3 download (from GitHub), and queryable through PedcBioPortal and the NCI's pediatric Molecular Targets Platform. Notably, we have expanded PBTA molecular subtyping to include methylation information to align with the WHO 2021 Central Nervous System Tumor classifications, allowing us to create research- grade integrated diagnoses for these tumors. Conclusions: OpenPedCan data and its reproducible analysis module framework are openly available and can be utilized and/or adapted by researchers to accelerate discovery, validation, and clinical translation.

Everolimus for Children With Recurrent or Progressive Low-Grade Glioma: Results From the Phase II PNOC001 Trial.

PURPOSE: The PNOC001 phase II single-arm trial sought to estimate progression-free survival (PFS) associated with everolimus therapy for progressive/recurrent pediatric low-grade glioma (pLGG) on the basis of phosphatidylinositol 3-kinase (PI3K)/AKT/mammalian target of rapamycin (mTOR) pathway activation as measured by phosphorylated-ribosomal protein S6 and to identify prognostic and predictive biomarkers. PATIENTS AND METHODS: Patients, age 3-21 years, with progressive/recurrent pLGG received everolimus orally, 5 mg/m2 once daily. Frequency of driver gene alterations was compared among independent pLGG cohorts of newly diagnosed and progressive/recurrent patients. PFS at 6 months (primary end point) and median PFS (secondary end point) were estimated for association with everolimus therapy. RESULTS: Between 2012 and 2019, 65 subjects with progressive/recurrent pLGG (median age, 9.6 years; range, 3.0-19.9; 46% female) were enrolled, with a median follow-up of 57.5 months. The 6-month PFS was 67.4% (95% CI, 60.0 to 80.0) and median PFS was 11.1 months (95% CI, 7.6 to 19.8). Hypertriglyceridemia was the most common grade ≥3 adverse event. PI3K/AKT/mTOR pathway activation did not correlate with clinical outcomes (6-month PFS, active 68.4% v nonactive 63.3%; median PFS, active 11.2 months v nonactive 11.1 months; P = .80). Rare/novel KIAA1549::BRAF fusion breakpoints were most frequent in supratentorial midline pilocytic astrocytomas, in patients with progressive/recurrent disease, and correlated with poor clinical outcomes (median PFS, rare/novel KIAA1549::BRAF fusion breakpoints 6.1 months v common KIAA1549::BRAF fusion breakpoints 16.7 months; P < .05). Multivariate analysis confirmed their independent risk factor status for disease progression in PNOC001 and other, independent cohorts. Additionally, rare pathogenic germline variants in homologous recombination genes were identified in 6.8% of PNOC001 patients. CONCLUSION: Everolimus is a well-tolerated therapy for progressive/recurrent pLGGs. Rare/novel KIAA1549::BRAF fusion breakpoints may define biomarkers for progressive disease and should be assessed in future clinical trials.

OpenPBTA: The Open Pediatric Brain Tumor Atlas.

Pediatric brain and spinal cancers are collectively the leading disease-related cause of death in children; thus, we urgently need curative therapeutic strategies for these tumors. To accelerate such discoveries, the Children's Brain Tumor Network (CBTN) and Pacific Pediatric Neuro-Oncology Consortium (PNOC) created a systematic process for tumor biobanking, model generation, and sequencing with immediate access to harmonized data. We leverage these data to establish OpenPBTA, an open collaborative project with over 40 scalable analysis modules that genomically characterize 1,074 pediatric brain tumors. Transcriptomic classification reveals universal TP53 dysregulation in mismatch repair-deficient hypermutant high-grade gliomas and TP53 loss as a significant marker for poor overall survival in ependymomas and H3 K28-mutant diffuse midline gliomas. Already being actively applied to other pediatric cancers and PNOC molecular tumor board decision-making, OpenPBTA is an invaluable resource to the pediatric oncology community.

Radiomics for characterization of the glioma immune microenvironment.

Increasing evidence suggests that besides mutational and molecular alterations, the immune component of the tumor microenvironment also substantially impacts tumor behavior and complicates treatment response, particularly to immunotherapies. Although the standard method for characterizing tumor immune profile is through performing integrated genomic analysis on tissue biopsies, the dynamic change in the immune composition of the tumor microenvironment makes this approach not feasible, especially for brain tumors. Radiomics is a rapidly growing field that uses advanced imaging techniques and computational algorithms to extract numerous quantitative features from medical images. Recent advances in machine learning methods are facilitating biological validation of radiomic signatures and allowing them to "mine" for a variety of significant correlates, including genetic, immunologic, and histologic data. Radiomics has the potential to be used as a non-invasive approach to predict the presence and density of immune cells within the microenvironment, as well as to assess the expression of immune-related genes and pathways. This information can be essential for patient stratification, informing treatment decisions and predicting patients' response to immunotherapies. This is particularly important for tumors with difficult surgical access such as gliomas. In this review, we provide an overview of the glioma microenvironment, describe novel approaches for clustering patients based on their tumor immune profile, and discuss the latest progress on utilization of radiomics for immune profiling of glioma based on current literature.

Analysis of matched primary and recurrent BRCA1/2 mutation-associated tumors identifies recurrence-specific drivers.

Recurrence is a major cause of death among BRCA1/2 mutation carriers with breast (BrCa) and ovarian cancers (OvCa). Herein we perform multi-omic sequencing on 67 paired primary and recurrent BrCa and OvCa from 27 BRCA1/2 mutation carriers to identify potential recurrence-specific drivers. PARP1 amplifications are identified in recurrences (False Discovery Rate q = 0.05), and PARP1 is significantly overexpressed across primary BrCa and recurrent BrCa and OvCa, independent of amplification status. RNA sequencing analysis finds two BRCA2 isoforms, BRCA2-201/Long and BRCA2-001/Short, respectively predicted to be sensitive and insensitive to nonsense-mediated decay. BRCA2-001/Short is expressed more frequently in recurrences and associated with reduced overall survival in breast cancer (87 vs. 121 months; Hazard Ratio = 2.5 [1.18-5.5]). Loss of heterozygosity (LOH) status is discordant in 25% of patient's primary and recurrent tumors, with switching between both LOH and lack of LOH found. Our study reveals multiple potential drivers of recurrent disease in BRCA1/2 mutation-associated cancer, improving our understanding of tumor evolution and suggesting potential biomarkers.

Upfront Biology-Guided Therapy in Diffuse Intrinsic Pontine Glioma: Therapeutic, Molecular, and Biomarker Outcomes from PNOC003.

PURPOSE: PNOC003 is a multicenter precision medicine trial for children and young adults with newly diagnosed diffuse intrinsic pontine glioma (DIPG). PATIENTS AND METHODS: Patients (3-25 years) were enrolled on the basis of imaging consistent with DIPG. Biopsy tissue was collected for whole-exome and mRNA sequencing. After radiotherapy (RT), patients were assigned up to four FDA-approved drugs based on molecular tumor board recommendations. H3K27M-mutant circulating tumor DNA (ctDNA) was longitudinally measured. Tumor tissue and matched primary cell lines were characterized using whole-genome sequencing and DNA methylation profiling. When applicable, results were verified in an independent cohort from the Children's Brain Tumor Network (CBTN). RESULTS: Of 38 patients enrolled, 28 patients (median 6 years, 10 females) were reviewed by the molecular tumor board. Of those, 19 followed treatment recommendations. Median overall survival (OS) was 13.1 months [95% confidence interval (CI), 11.2-18.4] with no difference between patients who followed recommendations and those who did not. H3K27M-mutant ctDNA was detected at baseline in 60% of cases tested and associated with response to RT and survival. Eleven cell lines were established, showing 100% fidelity of key somatic driver gene alterations in the primary tumor. In H3K27-altered DIPGs, TP53 mutations were associated with worse OS (TP53mut 11.1 mo; 95% CI, 8.7-14; TP53wt 13.3 mo; 95% CI, 11.8-NA; P = 3.4e-2), genome instability (P = 3.1e-3), and RT resistance (P = 6.4e-4). The CBTN cohort confirmed an association between TP53 mutation status, genome instability, and clinical outcome. CONCLUSIONS: Upfront treatment-naïve biopsy provides insight into clinically relevant molecular alterations and prognostic biomarkers for H3K27-altered DIPGs.

PTEN Loss and BRCA1 Promoter Hypermethylation Negatively Predict for Immunogenicity in BRCA-Deficient Ovarian Cancer.

PURPOSE: Ovarian cancers can exhibit a prominent immune infiltrate, but clinical trials have not demonstrated substantive response rates to immune checkpoint blockade monotherapy. We aimed to understand genomic features associated with immunogenicity in BRCA1/2 mutation-associated cancers. MATERIALS AND METHODS: Using the Cancer Genome Atlas whole-exome sequencing, methylation, and expression data, we analyzed 66 ovarian cancers with either germline or somatic loss of BRCA1/2 and whole-exome sequencing, immunohistochemistry, and CyTOF in 20 ovarian cancers with germline BRCA1/2 pathogenic variants from Penn. RESULTS: We found two groups of BRCA1/2 ovarian cancers differing in their immunogenicity: (1) 37 tumors significantly enriched for PTEN loss (11, 30%) and BRCA1 promoter-hypermethylated (10, 27%; P = .0016) and (2) PTEN wild-type (28 of 29 tumors) cancers, with the latter group having longer overall survival (OS; P = .0186, median OS not reached v median OS = 66.1 months). BRCA1/2-mutant PTEN loss and BRCA1 promoter-hypermethylated cancers were characterized by the decreased composition of lymphocytes estimated by gene expression (P = .0030), cytolytic index (P = .034), and cytokine expression but higher homologous recombination deficiency scores (P = .00013). Large-scale state transitions were the primary discriminating feature (P = .001); neither mutational burden nor neoantigen burden could explain differences in immunogenicity. In Penn tumors, PTEN loss and high homologous recombination deficiency cancers exhibited fewer CD3+ (P = .05), CD8+ (P = .012), and FOXP3+ (P = .0087) T cells; decreased PRF1 expression (P = .041); and lower immune costimulatory and inhibitory molecule expression. CONCLUSION: Our study suggests that within ovarian cancers with genetic loss of BRCA1/2 are two subsets exhibiting differential immunogenicity, with lower levels associated with PTEN loss and BRCA hypermethylation. These genomic features of BRCA1/2-associated ovarian cancers may inform considerations around how to optimally deploy immune checkpoint inhibitors in the clinic.

Evolution of delayed resistance to immunotherapy in a melanoma responder.

Despite initial responses1-3, most melanoma patients develop resistance4 to immune checkpoint blockade (ICB). To understand the evolution of resistance, we studied 37 tumor samples over 9 years from a patient with metastatic melanoma with complete clinical response to ICB followed by delayed recurrence and death. Phylogenetic analysis revealed co-evolution of seven lineages with multiple convergent, but independent resistance-associated alterations. All recurrent tumors emerged from a lineage characterized by loss of chromosome 15q, with post-treatment clones acquiring additional genomic driver events. Deconvolution of bulk RNA sequencing and highly multiplexed immunofluorescence (t-CyCIF) revealed differences in immune composition among different lineages. Imaging revealed a vasculogenic mimicry phenotype in NGFRhi tumor cells with high PD-L1 expression in close proximity to immune cells. Rapid autopsy demonstrated two distinct NGFR spatial patterns with high polarity and proximity to immune cells in subcutaneous tumors versus a diffuse spatial pattern in lung tumors, suggesting different roles of this neural-crest-like program in different tumor microenvironments. Broadly, this study establishes a high-resolution map of the evolutionary dynamics of resistance to ICB, characterizes a de-differentiated neural-crest tumor population in melanoma immunotherapy resistance and describes site-specific differences in tumor-immune interactions via longitudinal analysis of a patient with melanoma with an unusual clinical course.

Association of Filaggrin Loss-of-Function Variants With Race in Children With Atopic Dermatitis.

IMPORTANCE: Atopic dermatitis (AD) is a common chronic illness that has been associated with variation in the filaggrin gene (FLG). Four variants are most often evaluated. OBJECTIVES: To comprehensively describe and compare results from targeted sequencing of FLG loss-of-function (LoF) variants in children of African and European ancestry and the association of these variants with onset and persistence of AD. DESIGN, SETTING, AND PARTICIPANTS: This prospective US cohort study assessed the genetic subcohort of the Pediatric Eczema Elective Registry (PEER). Children with mild to moderate AD were included in the analysis. Massively parallel sequencing (MPS) was used to focus on FLG LoF variation in white and African American children. Patients were enrolled from June 2005 through July 2017. Data were analyzed from January 25 through May 10, 2019. MAIN OUTCOMES AND MEASURES: Associations of FLG LoF variation with white and African American ancestry and with the risk and persistence of AD. RESULTS: A total of 741 children were included in the analysis (394 [53.2%] female and 347 [46.8%] male; mean [SD] age at onset, 1.97 [2.72] years); of these, 394 (53.2%) were white, 326 (44.0%) were African American, and 21 (2.8%) were of other ancestries. Using MPS technology, 23 FLG LoF variants were found in children with AD. The prevalence of FLG LoF variants was 177 participants (23.9%) in the full cohort, 124 white participants (31.5%), and 50 African American participants (15.3%). The odds ratio for carrying any FLG LoF variant in a white child compared with an African American child with AD was 2.44 (95% CI, 1.76-3.39). Some FLG LoF variants are only found in children of a specific ancestry (eg, p.S3316* and p.R826* were not seen in white patients). Children with an FLG LoF were more likely to have persistent AD (odds ratio, 0.67; 95% CI, 0.56-0.80). CONCLUSIONS AND RELEVANCE: The FLG LoF variants in a US cohort of children with mild to moderate AD differ significantly by race and their association with the persistence of AD. Conventional testing of the 4 frequently evaluated variants is inadequate. Any planned genetic diagnostic test for AD based on FLG LoF variants must be inclusive and not rely on the most frequently studied variants.

Genomic Signatures Predict the Immunogenicity of BRCA-Deficient Breast Cancer.

PURPOSE: Breast cancers with BRCA1/2 alterations have a relatively high mutational load, suggesting that immune checkpoint blockade may be a potential treatment option. However, the degree of immune cell infiltration varies widely, and molecular features contributing to this variability remain unknown. EXPERIMENTAL DESIGN: We hypothesized that genomic signatures might predict immunogenicity in BRCA1/2 breast cancers. Using The Cancer Genome Atlas (TCGA) genomic data, we compared breast cancers with (89) and without (770) either germline or somatic BRCA1/2 alterations. We also studied 35 breast cancers with germline BRCA1/2 mutations from Penn using WES and IHC. RESULTS: We found that homologous recombination deficiency (HRD) scores were negatively associated with expression-based immune indices [cytolytic index (P = 0.04), immune ESTIMATE (P = 0.002), type II IFN signaling (P = 0.002)] despite being associated with a higher mutational/neoantigen burden, in BRCA1/2 mutant breast cancers. Further, absence of allele-specific loss of heterozygosity (LOH negative; P = 0.01) or subclonality (P = 0.003) of germline and somatic BRCA1/2 mutations, respectively, predicted for heightened cytolytic activity. Gene set analysis found that multiple innate and adaptive immune pathways that converge on NF-κB may contribute to this heightened immunogenicity. IHC of Penn breast cancers demonstrated increased CD45+ (P = 0.039) and CD8+ infiltrates (P = 0.037) and increased PDL1 expression (P = 0.012) in HRD-low or LOH-negative cancers. Triple-negative cancers with low HRD had far greater CD8+ T cells (P = 0.0011) and Perforin 1 expression (P = 0.014) compared with hormone receptor-positive HRD-high cancers. CONCLUSIONS: HRD scores and hormone receptor subtype are predictive of immunogenicity in BRCA1/2 breast cancers and may inform the design of optimal immune therapeutic strategies.

A single dose of neoadjuvant PD-1 blockade predicts clinical outcomes in resectable melanoma.

Immunologic responses to anti-PD-1 therapy in melanoma patients occur rapidly with pharmacodynamic T cell responses detectable in blood by 3 weeks. It is unclear, however, whether these early blood-based observations translate to the tumor microenvironment. We conducted a study of neoadjuvant/adjuvant anti-PD-1 therapy in stage III/IV melanoma. We hypothesized that immune reinvigoration in the tumor would be detectable at 3 weeks and that this response would correlate with disease-free survival. We identified a rapid and potent anti-tumor response, with 8 of 27 patients experiencing a complete or major pathological response after a single dose of anti-PD-1, all of whom remain disease free. These rapid pathologic and clinical responses were associated with accumulation of exhausted CD8 T cells in the tumor at 3 weeks, with reinvigoration in the blood observed as early as 1 week. Transcriptional analysis demonstrated a pretreatment immune signature (neoadjuvant response signature) that was associated with clinical benefit. In contrast, patients with disease recurrence displayed mechanisms of resistance including immune suppression, mutational escape, and/or tumor evolution. Neoadjuvant anti-PD-1 treatment is effective in high-risk resectable stage III/IV melanoma. Pathological response and immunological analyses after a single neoadjuvant dose can be used to predict clinical outcome and to dissect underlying mechanisms in checkpoint blockade.

Uncommon Filaggrin Variants Are Associated with Persistent Atopic Dermatitis in African Americans.

Atopic dermatitis (AD) is a common illness that has been associated with filaggrin gene (FLG) loss of function (LoF) variation. In African Americans, a group that commonly has AD and has not been well studied, FLG LoF variation is rarely found. Our objective was to use massively parallel sequencing to evaluate FLG LoF variation in children of African ancestry to evaluate the association between FLG LoF variation and AD and AD persistence. We studied 262 African American children with AD. Nine unique FLG exon 3 LoF variants were identified for an overall minor variant frequency of 6.30% (95% confidence interval [CI] = 4.37-8.73). The most common variants were p.R501X (1.72%, 95% CI = 0.79-3.24), p.S3316X (1.34%, 95% CI = 0.54-2.73), and p.R826X (0.95%, 95% CI = 0.31-2.2). Over an average follow-up period of 96.4 (95% CI = 92.0-100.8) months, African American children with FLG LoF were less likely to be symptom free (odds ratio = 0.36, 95% CI = 0.14-0.89, P = 0.027) compared with a FLG wild-type child. In contrast to previous reports, uncommon FLG LoF variants in African American children exist and are associated with AD and more persistent AD. In contrast to Europeans, no FLG LoF variants predominate in African American children. Properly determining FLG LoF status requires advanced sequencing techniques.

BRCA locus-specific loss of heterozygosity in germline BRCA1 and BRCA2 carriers.

Complete loss of BRCA1 or BRCA2 function is associated with sensitivity to DNA damaging agents. However, not all BRCA1 and BRCA2 germline mutation-associated tumors respond. Herein we report analyses of 160 BRCA1 and BRCA2 germline mutation-associated breast and ovarian tumors. Retention of the normal BRCA1 or BRCA2 allele (absence of locus-specific loss of heterozygosity (LOH)) is observed in 7% of BRCA1 ovarian, 16% of BRCA2 ovarian, 10% of BRCA1 breast, and 46% of BRCA2 breast tumors. These tumors have equivalent homologous recombination deficiency scores to sporadic tumors, significantly lower than scores in tumors with locus-specific LOH (ovarian, P = 0.0004; breast P < 0.0001, two-tailed Student's t-test). Absence of locus-specific LOH is associated with decreased overall survival in ovarian cancer patients treated with platinum chemotherapy (P = 0.01, log-rank test). Locus-specific LOH may be a clinically useful biomarker to predict primary resistance to DNA damaging agents in patients with germline BRCA1 and BRCA2 mutations.Most tumours associated with germline BRCA1/BRCA2 loss of function mutations respond to DNA damaging agents, however, some do not. Herein, the authors identify that a subset of breast/ovarian tumors retain a normal allele, which is associated with decreased overall survival after DNA damage-inducing platinum chemotherapy.

Identification of secreted proteins that reflect autophagy dynamics within tumor cells.

Macroautophagy, a catabolic process of cellular self-digestion, is an important tumor cell survival mechanism and a potential target in antineoplastic therapies. Recent discoveries have implicated autophagy in the cellular secretory process, but potential roles of autophagy-mediated secretion in modifying the tumor microenvironment are poorly understood. Furthermore, efforts to inhibit autophagy in clinical trials have been hampered by suboptimal methods to quantitatively measure tumor autophagy levels. Here, we leveraged the autophagy-based involvement in cellular secretion to identify shed proteins associated with autophagy levels in melanoma. The secretome of low-autophagy WM793 melanoma cells was compared to its highly autophagic metastatic derivative, 1205Lu in physiological 3-dimensional cell culture using quantitative proteomics. These comparisons identified candidate autophagy biomarkers IL1B (interleukin 1, ß), CXCL8 (chemokine (C-X-C motif) ligand 8), LIF (leukemia inhibitory factor), FAM3C (family with sequence similarity 3, member C), and DKK3 (dickkopf WNT signaling pathway inhibitor 3) with known roles in inflammation and tumorigenesis, and these proteins were subsequently shown to be elevated in supernatants of an independent panel of high-autophagy melanoma cell lines. Secretion levels of these proteins increased when low-autophagy melanoma cells were treated with the autophagy-inducing tat-BECN1 (Beclin 1) peptide and decreased when ATG7 (autophagy-related 7) was silenced in high-autophagy cells, thereby supporting a mechanistic link between these secreted proteins and autophagy. In addition, serum from metastatic melanoma patients with high tumor autophagy levels exhibited higher levels of these proteins than serum from patients with low-autophagy tumors. These results suggest that autophagy-related secretion affects the tumor microenvironment and measurement of autophagy-associated secreted proteins in plasma and possibly in tumors can serve as surrogates for intracellular autophagy dynamics in tumor cells.