Clinical detection of sequence and structural variants in known cancer genes points to viable treatment options for a minority of children with cancer.¹ To increase the number of children who benefit from genomic profiling, gene expression information must be considered alongside mutations.^2,3 Although high expression has been used to nominate drug targets for pediatric cancers,^4,5 its utility has not been evaluated in a systematic way.⁶ We describe a child with a rare sarcoma that was profiled with whole-genome and RNA sequencing (RNA-Seq) techniques. Although the tumor did not harbor DNA mutations targetable by available therapies, incorporation of gene expression information derived from RNA-Seq analysis led to a therapy that produced a significant clinical response. We use this case to describe a framework for inclusion of gene expression into the clinical genomic evaluation of pediatric tumors.

Aim: To provide a comprehensive understanding of gene regulatory networks in the developing human brain and a foundation for interpreting pathogenic deregulation. Materials & methods: We generated reference epigenomes and transcriptomes of dissected brain regions and primary neural progenitor cells (NPCs) derived from cortical and ganglionic eminence tissues of four normal human fetuses. Results: Integration of these data across developmental stages revealed a directional increase in active regulatory states, transcription factor activities and gene transcription with developmental stage. Consistent with differences in their biology, NPCs derived from cortical and ganglionic eminence regions contained common, region specific, and gestational week specific regulatory states. Conclusion: We provide a high-resolution regulatory network for NPCs from different brain regions as a comprehensive reference for future studies.

Clinical detection of sequence and structural variants in known cancer genes points to viable treatment options for a minority of children with cancer.1 To increase the number of children who benefit from genomic profiling, gene expression information must be considered alongside mutations.2,3 Although high expression has been used to nominate drug targets for pediatric cancers,4,5 its utility has not been evaluated in a systematic way.6 We describe a child with a rare sarcoma that was profiled with whole-genome and RNA sequencing (RNA-Seq) techniques. Although the tumor did not harbor DNA mutations targetable by available therapies, incorporation of gene expression information derived from RNA-Seq analysis led to a therapy that produced a significant clinical response. We use this case to describe a framework for inclusion of gene expression into the clinical genomic evaluation of pediatric tumors.

Cervical cancer is the most common cancer affecting sub-Saharan African women and is prevalent among HIV-positive (HIV⁺) individuals. No comprehensive profiling of cancer genomes, transcriptomes or epigenomes has been performed in this population thus far. We characterized 118 tumors from Ugandan patients, of whom 72 were HIV⁺, and performed extended mutation analysis on an additional 89 tumors. We detected human papillomavirus (HPV)-clade-specific differences in tumor DNA methylation, promoter- and enhancer-associated histone marks, gene expression and pathway dysregulation. Changes in histone modification at HPV integration events were correlated with upregulation of nearby genes and endogenous retroviruses.

In this study, we evaluate the impact of whole genome and transcriptome analysis (WGTA) on predictive molecular profiling and histologic diagnosis in a cohort of advanced malignancies. WGTA was used to generate reports including molecular alterations and site/tissue of origin prediction. Two reviewers analyzed genomic reports, clinical history, and tumor pathology. We used National Comprehensive Cancer Network (NCCN) consensus guidelines, Food and Drug Administration (FDA) approvals, and provincially reimbursed treatments to define genomic biomarkers associated with approved targeted therapeutic options (TTOs). Tumor tissue/site of origin was reassessed for most cases using genomic analysis, including a machine learning algorithm (Supervised Cancer Origin Prediction Using Expression [SCOPE]) trained on The Cancer Genome Atlas data. WGTA was performed on 652 cases, including a range of primary tumor types/tumor sites and 15 malignant tumors of uncertain histogenesis (MTUH). At the time WGTA was performed, alterations associated with an approved TTO were identified in 39 (6%) cases; 3 of these were not identified through routine pathology workup. In seven (1%) cases, the pathology workup either failed, was not performed, or gave a different result from the WGTA. Approved TTOs identified by WGTA increased to 103 (16%) when applying 2021 guidelines. The histopathologic diagnosis was reviewed in 389 cases and agreed with the diagnostic consensus after WGTA in 94% of non-MTUH cases (n = 374). The remainder included situations where the morphologic diagnosis was changed based on WGTA and clinical data (0.5%), or where the WGTA was non-contributory (5%). The 15 MTUH were all diagnosed as specific tumor types by WGTA. Tumor board reviews including WGTA agreed with almost all initial predictive molecular profile and histopathologic diagnoses. WGTA was a powerful tool to assign site/tissue of origin in MTUH. Current efforts focus on improving therapeutic predictive power and decreasing cost to enhance use of WGTA data as a routine clinical test.

Ovarian cancer presents as an aggressive, advanced stage cancer with widespread metastases that depend primarily on multicellular spheroids in the peritoneal fluid. To identify new druggable pathways related to metastatic progression and spheroid formation, we integrated microRNA and mRNA sequencing data from 293 tumors from The Cancer Genome Atlas (TCGA) ovarian cancer cohort. We identified miR-509-3p as a clinically significant microRNA that is more abundant in patients with favorable survival in both the TCGA cohort (P = 2.3E-3), and, by in situ hybridization (ISH), in an independent cohort of 157 tumors (P < 1.0E-3). We found that miR-509-3p attenuated migration and disrupted multi-cellular spheroids in HEYA8, OVCAR8, SKOV3, OVCAR3, OVCAR4 and OVCAR5 cell lines. Consistent with disrupted spheroid formation, in TCGA data miR-509-3p's most strongly anti-correlated predicted targets were enriched in components of the extracellular matrix (ECM). We validated the Hippo pathway effector YAP1 as a direct miR-509-3p target. We showed that siRNA to YAP1 replicated 90% of miR-509-3p-mediated migration attenuation in OVCAR8, which contained high levels of YAP1 protein, but not in the other cell lines, in which levels of this protein were moderate to low. Our data suggest that the miR-509-3p/YAP1 axis may be a new druggable target in cancers with high YAP1, and we propose that therapeutically targeting the miR-509-3p/YAP1/ECM axis may disrupt early steps in multi-cellular spheroid formation, and so inhibit metastasis in epithelial ovarian cancer and potentially in other cancers.

Tumor recurrence is a leading cause of cancer mortality. Therapies for recurrent disease may fail, at least in part, because the genomic alterations driving the growth of recurrences are distinct from those in the initial tumor. To explore this hypothesis, we sequenced the exomes of 23 initial low-grade gliomas and recurrent tumors resected from the same patients. In 43% of cases, at least half of the mutations in the initial tumor were undetected at recurrence, including driver mutations in TP53, ATRX, SMARCA4, and BRAF; this suggests that recurrent tumors are often seeded by cells derived from the initial tumor at a very early stage of their evolution. Notably, tumors from 6 of 10 patients treated with the chemotherapeutic drug temozolomide (TMZ) followed an alternative evolutionary path to high-grade glioma. At recurrence, these tumors were hypermutated and harbored driver mutations in the RB (retinoblastoma) and Akt-mTOR (mammalian target of rapamycin) pathways that bore the signature of TMZ-induced mutagenesis.

Chromosomal rearrangements are initiating events in acute lymphoblastic leukaemia (ALL). Here using RNA sequencing of 560 ALL cases, we identify rearrangements between MEF2D (myocyte enhancer factor 2D) and five genes (BCL9, CSF1R, DAZAP1, HNRNPUL1 and SS18) in 22 B progenitor ALL (B-ALL) cases with a distinct gene expression profile, the most common of which is MEF2D-BCL9. Examination of an extended cohort of 1,164 B-ALL cases identified 30 cases with MEF2D rearrangements, which include an additional fusion partner, FOXJ2; thus, MEF2D-rearranged cases comprise 5.3% of cases lacking recurring alterations. MEF2D-rearranged ALL is characterized by a distinct immunophenotype, DNA copy number alterations at the rearrangement sites, older diagnosis age and poor outcome. The rearrangements result in enhanced MEF2D transcriptional activity, lymphoid transformation, activation of HDAC9 expression and sensitive to histone deacetylase inhibitor treatment. Thus, MEF2D-rearranged ALL represents a distinct form of high-risk leukaemia, for which new therapeutic approaches should be considered.

Current therapies for medulloblastoma, a highly malignant childhood brain tumour, impose debilitating effects on the developing child, and highlight the need for molecularly targeted treatments with reduced toxicity. Previous studies have been unable to identify the full spectrum of driver genes and molecular processes that operate in medulloblastoma subgroups. Here we analyse the somatic landscape across 491 sequenced medulloblastoma samples and the molecular heterogeneity among 1,256 epigenetically analysed cases, and identify subgroup-specific driver alterations that include previously undiscovered actionable targets. Driver mutations were confidently assigned to most patients belonging to Group 3 and Group 4 medulloblastoma subgroups, greatly enhancing previous knowledge. New molecular subtypes were differentially enriched for specific driver events, including hotspot in-frame insertions that target KBTBD4 and 'enhancer hijacking' events that activate PRDM6. Thus, the application of integrative genomics to an extensive cohort of clinical samples derived from a single childhood cancer entity revealed a series of cancer genes and biologically relevant subtype diversity that represent attractive therapeutic targets for the treatment of patients with medulloblastoma.

Sonic hedgehog medulloblastoma encompasses a clinically and molecularly diverse group of cancers of the developing central nervous system. Here, we use unbiased sequencing of the transcriptome across a large cohort of 250 tumors to reveal differences among molecular subtypes of the disease, and demonstrate the previously unappreciated importance of non-coding RNA transcripts. We identify alterations within the cAMP dependent pathway (GNAS, PRKAR1A) which converge on GLI2 activity and show that 18% of tumors have a genetic event that directly targets the abundance and/or stability of MYCN. Furthermore, we discover an extensive network of fusions in focally amplified regions encompassing GLI2, and several loss-of-function fusions in tumor suppressor genes PTCH1, SUFU and NCOR1. Molecular convergence on a subset of genes by nucleotide variants, copy number aberrations, and gene fusions highlight the key roles of specific pathways in the pathogenesis of Sonic hedgehog medulloblastoma and open up opportunities for therapeutic intervention.