St. Jude Survivorship Portal: sharing and analyzing large clinical and genomic datasets from pediatric cancer survivors.

Childhood cancer survivorship studies generate comprehensive datasets comprising demographic, diagnosis, treatment, outcome, and genomic data from survivors. To broadly share this data, we created the St. Jude Survivorship Portal (https://survivorship.stjude.cloud), the first data portal for sharing, analyzing, and visualizing pediatric cancer survivorship data. Over 1,600 phenotypic variables and 400 million genetic variants from over 7,700 childhood cancer survivors can be explored on this free, open-access portal. Summary statistics of variables are computed on-the-fly and visualized through interactive and customizable charts. Survivor cohorts can be customized and/or divided into groups for comparative analysis. Users can also seamlessly perform cumulative incidence and regression analyses on the stored survivorship data. Using the portal, we explored the ototoxic effects of platinum-based chemotherapy, uncovered a novel association between mental health, age, and limb amputation, and discovered a novel haplotype in MAGI3 strongly associated with cardiomyopathy specifically in survivors of African ancestry.

High-grade glioma in infants and young children is histologically, molecularly, and clinically diverse: Results from the SJYC07 trial and institutional experience.

BACKGROUND: High-grade gliomas (HGG) in young children pose a challenge due to favorable but unpredictable outcomes. While retrospective studies broadened our understanding of tumor biology, prospective data is lacking. METHODS: A cohort of children with histologically diagnosed HGG from the SJYC07 trial was augmented with nonprotocol patients with HGG treated at St. Jude Children's Research Hospital from November 2007 to December 2020. DNA methylome profiling and whole genome, whole exome, and RNA sequencing were performed. These data were integrated with histopathology to yield an integrated diagnosis. Clinical characteristics and preoperative imaging were analyzed. RESULTS: Fifty-six children (0.0-4.4 years) were identified. Integrated analysis split the cohort into four categories: infant-type hemispheric glioma (IHG), HGG, low-grade glioma (LGG), and other-central nervous system (CNS) tumors. IHG was the most prevalent (n = 22), occurred in the youngest patients (median age = 0.4 years), and commonly harbored receptor tyrosine kinase gene fusions (7 ALK, 2 ROS1, 3 NTRK1/2/3, 4 MET). The 5-year event-free (EFS) and overall survival (OS) for IHG was 53.13% (95%CI: 35.52-79.47) and 90.91% (95%CI: 79.66-100.00) vs. 0.0% and 16.67% (95%CI: 2.78-99.74%) for HGG (p = 0.0043, p = 0.00013). EFS and OS were not different between IHG and LGG (p = 0.95, p = 0.43). Imaging review showed IHGs are associated with circumscribed margins (p = 0.0047), hemispheric location (p = 0.0010), and intratumoral hemorrhage (p = 0.0149). CONCLUSIONS: HGG in young children is heterogeneous and best defined by integrating histopathological and molecular features. Patients with IHG have relatively good outcomes, yet they endure significant deficits, making them good candidates for therapy de-escalation and trials of molecular targeted therapy.

Molecular Classification Improves Therapeutic Options for Infants and Young Children With Medulloblastoma.

Medulloblastoma in infants and young children is a major challenge to treat because craniospinal irradiation (CSI), a cornerstone of therapy for older children, is disproportionately damaging to very young children. As a result, trials have attempted to delay, omit, and replace this therapy. Although success has been limited, the approach has not been a complete failure. In fact, this approach has cured a significant number of children with medulloblastoma. However, many children have endured intensive regimens of chemotherapy only to experience relapse and undergo salvage treatment with CSI, often at higher doses and with worse morbidity than they would have initially experienced. Recent advancements in molecular diagnostics have proven that response to therapy is biologically driven. Medulloblastoma in infants and young children is divided into 2 molecular groups: Sonic Hedgehog (SHH) and group 3 (G3). Both are chemotherapy-sensitive, but only the SHH medulloblastomas are reliably cured with chemotherapy alone. Moreover, SHH can be molecularly parsed into 2 groups: SHH-1 and SHH-2, with SHH-2 showing higher cure rates with less intensive chemotherapy and SHH-1 requiring more intensive regimens. G3 medulloblastoma, on the other hand, has a near universal relapse rate after chemotherapy-only regimens. This predictability represents a significant breakthrough and affords oncologists the ability to properly risk-stratify therapy in such a way that the most curative and least toxic therapy is selected. This review examines the treatment of medulloblastoma in infants and young children, discusses the molecular advancements, and proposes how to use this information to structure the future management of this disease.

ppBAM: ProteinPaint BAM track for read alignment visualization and variant genotyping.

SUMMARY: ProteinPaint BAM track (ppBAM) is designed to assist variant review for cancer research and clinical genomics. With performant server-side computing and rendering, ppBAM supports on-the-fly variant genotyping of thousands of reads using Smith-Waterman alignment. To better visualize support for complex variants, reads are realigned against the mutated reference sequence using ClustalO. ppBAM also supports the BAM slicing API of the NCI Genomic Data Commons (GDC) portal, letting researchers conveniently examine genomic details of vast amounts of cancer sequencing data and reinterpret variant calls. AVAILABILITY AND IMPLEMENTATION: BAM track examples, tutorial, and GDC file access links are available at https://proteinpaint.stjude.org/bam/. Source code is available at https://github.com/stjude/proteinpaint.

The genomic landscape of pediatric acute lymphoblastic leukemia.

Acute lymphoblastic leukemia (ALL) is the most common childhood cancer. Here, using whole-genome, exome and transcriptome sequencing of 2,754 childhood patients with ALL, we find that, despite a generally low mutation burden, ALL cases harbor a median of four putative somatic driver alterations per sample, with 376 putative driver genes identified varying in prevalence across ALL subtypes. Most samples harbor at least one rare gene alteration, including 70 putative cancer driver genes associated with ubiquitination, SUMOylation, noncoding transcripts and other functions. In hyperdiploid B-ALL, chromosomal gains are acquired early and synchronously before ultraviolet-induced mutation. By contrast, ultraviolet-induced mutations precede chromosomal gains in B-ALL cases with intrachromosomal amplification of chromosome 21. We also demonstrate the prognostic significance of genetic alterations within subtypes. Intriguingly, DUX4- and KMT2A-rearranged subtypes separate into CEBPA/FLT3- or NFATC4-expressing subgroups with potential clinical implications. Together, these results deepen understanding of the ALL genomic landscape and associated outcomes.

Molecular classification and outcome of children with rare CNS embryonal tumors: results from St. Jude Children's Research Hospital including the multi-center SJYC07 and SJMB03 clinical trials.

Methylation profiling has radically transformed our understanding of tumors previously called central nervous system primitive neuro-ectodermal tumors (CNS-PNET). While this marks a momentous step toward defining key differences, reclassification has thrown treatment into disarray. To shed light on response to therapy and guide clinical decision-making, we report outcomes and molecular features of children with CNS-PNETs from two multi-center risk-adapted studies (SJMB03 for patients ≥ 3 years; SJYC07 for patients < 3 years) complemented by a non-protocol institutional cohort. Seventy patients who had a histological diagnosis of CNS-PNET or CNS embryonal tumor from one of the new categories that has supplanted CNS-PNET were included. This cohort was molecularly characterized by DNA methylation profiling (n = 70), whole-exome sequencing (n = 53), RNA sequencing (n = 20), and germline sequencing (n = 28). Clinical characteristics were detailed, and treatment was divided into craniospinal irradiation (CSI)-containing (SJMB03 and SJMB03-like) and CSI-sparing therapy (SJYC07 and SJYC07-like). When the cohort was analyzed in its entirety, no differences were observed in the 5-year survival rates even when CSI-containing therapy was compared to CSI-sparing therapy. However, when analyzed by DNA methylation molecular grouping, significant survival differences were observed, and treatment particulars provided suggestions of therapeutic response. Patients with CNS neuroblastoma with FOXR2 activation (CNS-NB-FOXR2) had a 5-year event-free survival (EFS)/overall survival (OS) of 66.7% ± 19.2%/83.3% ± 15.2%, and CIC rearranged sarcoma (CNS-SARC-CIC) had a 5-year EFS/OS both of 57.1% ± 18.7% with most receiving regimens that contained radiation (focal or CSI) and multidrug chemotherapy. Patients with high-grade neuroepithelial tumor with BCOR alteration (HGNET-BCOR) had abysmal responses to upfront chemotherapy-only regimens (5-year EFS = 0%), but survival extended with salvage radiation after progression [5-year OS = 53.6% ± 20.1%]. Patients with embryonal tumor with multilayered rosettes (ETMR) or high-grade glioma/glioblastoma multiforme (HGG/GBM) did not respond favorably to any modality (5-year EFS/OS = 10.7 ± 5.8%/17.9 ± 7.2%, and 10% ± 9.0%/10% ± 9.0%, respectively). As an accompaniment, we have assembled this data onto an interactive website to allow users to probe and query the cases. By reporting on a carefully matched clinical and molecular cohort, we provide the needed insight for future clinical management.

A polygenic score for acute vaso-occlusive pain in pediatric sickle cell disease.

Individuals with monogenic disorders can experience variable phenotypes that are influenced by genetic variation. To investigate this in sickle cell disease (SCD), we performed whole-genome sequencing (WGS) of 722 individuals with hemoglobin HbSS or HbSß0-thalassemia from Baylor College of Medicine and from the St. Jude Children's Research Hospital Sickle Cell Clinical Research and Intervention Program (SCCRIP) longitudinal cohort study. We developed pipelines to identify genetic variants that modulate sickle hemoglobin polymerization in red blood cells and combined these with pain-associated variants to build a polygenic score (PGS) for acute vaso-occlusive pain (VOP). Overall, we interrogated the α-thalassemia deletion -α3.7 and 133 candidate single-nucleotide polymorphisms (SNPs) across 66 genes for associations with VOP in 327 SCCRIP participants followed longitudinally over 6 years. Twenty-one SNPs in 9 loci were associated with VOP, including 3 (BCL11A, MYB, and the ß-like globin gene cluster) that regulate erythrocyte fetal hemoglobin (HbF) levels and 6 (COMT, TBC1D1, KCNJ6, FAAH, NR3C1, and IL1A) that were associated previously with various pain syndromes. An unweighted PGS integrating all 21 SNPs was associated with the VOP event rate (estimate, 0.35; standard error, 0.04; P = 5.9 × 10-14) and VOP event occurrence (estimate, 0.42; standard error, 0.06; P = 4.1 × 10-13). These associations were stronger than those of any single locus. Our findings provide insights into the genetic modulation of VOP in children with SCD. More generally, we demonstrate the utility of WGS for investigating genetic contributions to the variable expression of SCD-associated morbidities.

St. Jude Cloud: A Pediatric Cancer Genomic Data-Sharing Ecosystem.

Effective data sharing is key to accelerating research to improve diagnostic precision, treatment efficacy, and long-term survival in pediatric cancer and other childhood catastrophic diseases. We present St. Jude Cloud (https://www.stjude.cloud), a cloud-based data-sharing ecosystem for accessing, analyzing, and visualizing genomic data from >10,000 pediatric patients with cancer and long-term survivors, and >800 pediatric sickle cell patients. Harmonized genomic data totaling 1.25 petabytes are freely available, including 12,104 whole genomes, 7,697 whole exomes, and 2,202 transcriptomes. The resource is expanding rapidly, with regular data uploads from St. Jude's prospective clinical genomics programs. Three interconnected apps within the ecosystem-Genomics Platform, Pediatric Cancer Knowledgebase, and Visualization Community-enable simultaneously performing advanced data analysis in the cloud and enhancing the Pediatric Cancer knowledgebase. We demonstrate the value of the ecosystem through use cases that classify 135 pediatric cancer subtypes by gene expression profiling and map mutational signatures across 35 pediatric cancer subtypes. SIGNIFICANCE: To advance research and treatment of pediatric cancer, we developed St. Jude Cloud, a data-sharing ecosystem for accessing >1.2 petabytes of raw genomic data from >10,000 pediatric patients and survivors, innovative analysis workflows, integrative multiomics visualizations, and a knowledgebase of published data contributed by the global pediatric cancer community.This article is highlighted in the In This Issue feature, p. 995.

Exploration of Coding and Non-coding Variants in Cancer Using GenomePaint.

GenomePaint (https://genomepaint.stjude.cloud/) is an interactive visualization platform for whole-genome, whole-exome, transcriptome, and epigenomic data of tumor samples. Its design captures the inter-relatedness between DNA variations and RNA expression, supporting in-depth exploration of both individual cancer genomes and full cohorts. Regulatory non-coding variants can be inspected and analyzed along with coding variants, and their functional impact further explored by examining 3D genome data from cancer cell lines. Further, GenomePaint correlates mutation and expression patterns with patient outcomes, and supports custom data upload. We used GenomePaint to unveil aberrant splicing that disrupts the RING domain of CREBBP, discover cis activation of the MYC oncogene by duplication of the NOTCH1-MYC enhancer in B-lineage acute lymphoblastic leukemia, and explore the inter- and intra-tumor heterogeneity at EGFR in adult glioblastomas. These examples demonstrate that deep multi-omics exploration of individual cancer genomes enabled by GenomePaint can lead to biological insights for follow-up validation.

Therapy-induced mutations drive the genomic landscape of relapsed acute lymphoblastic leukemia.

To study the mechanisms of relapse in acute lymphoblastic leukemia (ALL), we performed whole-genome sequencing of 103 diagnosis-relapse-germline trios and ultra-deep sequencing of 208 serial samples in 16 patients. Relapse-specific somatic alterations were enriched in 12 genes (NR3C1, NR3C2, TP53, NT5C2, FPGS, CREBBP, MSH2, MSH6, PMS2, WHSC1, PRPS1, and PRPS2) involved in drug response. Their prevalence was 17% in very early relapse (<9 months from diagnosis), 65% in early relapse (9-36 months), and 32% in late relapse (>36 months) groups. Convergent evolution, in which multiple subclones harbor mutations in the same drug resistance gene, was observed in 6 relapses and confirmed by single-cell sequencing in 1 case. Mathematical modeling and mutational signature analysis indicated that early relapse resistance acquisition was frequently a 2-step process in which a persistent clone survived initial therapy and later acquired bona fide resistance mutations during therapy. In contrast, very early relapses arose from preexisting resistant clone(s). Two novel relapse-specific mutational signatures, one of which was caused by thiopurine treatment based on in vitro drug exposure experiments, were identified in early and late relapses but were absent from 2540 pan-cancer diagnosis samples and 129 non-ALL relapses. The novel signatures were detected in 27% of relapsed ALLs and were responsible for 46% of acquired resistance mutations in NT5C2, PRPS1, NR3C1, and TP53. These results suggest that chemotherapy-induced drug resistance mutations facilitate a subset of pediatric ALL relapses.

Genomic subtyping and therapeutic targeting of acute erythroleukemia.

Acute erythroid leukemia (AEL) is a high-risk leukemia of poorly understood genetic basis, with controversy regarding diagnosis in the spectrum of myelodysplasia and myeloid leukemia. We compared genomic features of 159 childhood and adult AEL cases with non-AEL myeloid disorders and defined five age-related subgroups with distinct transcriptional profiles: adult, TP53 mutated; NPM1 mutated; KMT2A mutated/rearranged; adult, DDX41 mutated; and pediatric, NUP98 rearranged. Genomic features influenced outcome, with NPM1 mutations and HOXB9 overexpression being associated with a favorable prognosis and TP53, FLT3 or RB1 alterations associated with poor survival. Targetable signaling mutations were present in 45% of cases and included recurrent mutations of ALK and NTRK1, the latter of which drives erythroid leukemogenesis sensitive to TRK inhibition. This genomic landscape of AEL provides the framework for accurate diagnosis and risk stratification of this disease, and the rationale for testing targeted therapies in this high-risk leukemia.

Pan-cancer genome and transcriptome analyses of 1,699 paediatric leukaemias and solid tumours.

Analysis of molecular aberrations across multiple cancer types, known as pan-cancer analysis, identifies commonalities and differences in key biological processes that are dysregulated in cancer cells from diverse lineages. Pan-cancer analyses have been performed for adult but not paediatric cancers, which commonly occur in developing mesodermic rather than adult epithelial tissues. Here we present a pan-cancer study of somatic alterations, including single nucleotide variants, small insertions or deletions, structural variations, copy number alterations, gene fusions and internal tandem duplications in 1,699 paediatric leukaemias and solid tumours across six histotypes, with whole-genome, whole-exome and transcriptome sequencing data processed under a uniform analytical framework. We report 142 driver genes in paediatric cancers, of which only 45% match those found in adult pan-cancer studies; copy number alterations and structural variants constituted the majority (62%) of events. Eleven genome-wide mutational signatures were identified, including one attributed to ultraviolet-light exposure in eight aneuploid leukaemias. Transcription of the mutant allele was detectable for 34% of protein-coding mutations, and 20% exhibited allele-specific expression. These data provide a comprehensive genomic architecture for paediatric cancers and emphasize the need for paediatric cancer-specific development of precision therapies.

The genomic landscape of pediatric and young adult T-lineage acute lymphoblastic leukemia.

Genetic alterations that activate NOTCH1 signaling and T cell transcription factors, coupled with inactivation of the INK4/ARF tumor suppressors, are hallmarks of T-lineage acute lymphoblastic leukemia (T-ALL), but detailed genome-wide sequencing of large T-ALL cohorts has not been carried out. Using integrated genomic analysis of 264 T-ALL cases, we identified 106 putative driver genes, half of which had not previously been described in childhood T-ALL (for example, CCND3, CTCF, MYB, SMARCA4, ZFP36L2 and MYCN). We describe new mechanisms of coding and noncoding alteration and identify ten recurrently altered pathways, with associations between mutated genes and pathways, and stage or subtype of T-ALL. For example, NRAS/FLT3 mutations were associated with immature T-ALL, JAK3/STAT5B mutations in HOXA1 deregulated ALL, PTPN2 mutations in TLX1 deregulated T-ALL, and PIK3R1/PTEN mutations in TAL1 deregulated ALL, which suggests that different signaling pathways have distinct roles according to maturational stage. This genomic landscape provides a logical framework for the development of faithful genetic models and new therapeutic approaches.