Molecular analysis of cancer genomes in children with Lynch syndrome: Exploring causal associations.

Lynch syndrome (LS) predisposes to cancer in adulthood and is caused by heterozygous germline variants in a mismatch repair (MMR) gene. Recent studies show an increased prevalence of LS among children with cancer, suggesting a causal relationship. For LS-spectrum (LSS) cancers, including high-grade gliomas and colorectal cancer, causality has been supported by typical MMR-related tumor characteristics, but for non-LSS cancers, causality is unclear. We characterized 20 malignant tumors of 18 children with LS, including 16 non-LSS tumors. We investigated second hits, tumor mutational load, mutational signatures and MMR protein expression. In all LSS tumors and three non-LSS tumors, we detected MMR deficiency caused by second hit somatic alterations. Furthermore, these MMR-deficient tumors carried driver variants that likely originated as a consequence of MMR deficiency. However, in 13 non-LSS tumors (81%), a second hit and MMR deficiency were absent, thus a causal link between LS and cancer development in these children is lacking. These findings demonstrate that causality of LS in children with cancer, which can be determined by molecular tumor characterization, seems to be restricted to specific tumor types. Large molecular and epidemiological studies are needed to further refine the tumor spectrum in children with LS.

Transcriptional immunogenomic analysis reveals distinct immunological clusters in paediatric nervous system tumours.

BACKGROUND: Cancer immunotherapies including immune checkpoint inhibitors and Chimeric Antigen Receptor (CAR) T-cell therapy have shown variable response rates in paediatric patients highlighting the need to establish robust biomarkers for patient selection. While the tumour microenvironment in adults has been widely studied to delineate determinants of immune response, the immune composition of paediatric solid tumours remains relatively uncharacterized calling for investigations to identify potential immune biomarkers. METHODS: To inform immunotherapy approaches in paediatric cancers with embryonal origin, we performed an immunogenomic analysis of RNA-seq data from 925 treatment-naïve paediatric nervous system tumours (pedNST) spanning 12 cancer types from three publicly available data sets. RESULTS: Within pedNST, we uncovered four broad immune clusters: Paediatric Inflamed (10%), Myeloid Predominant (30%), Immune Neutral (43%) and Immune Desert (17%). We validated these clusters using immunohistochemistry, methylation immune inference and segmentation analysis of tissue images. We report shared biology of these immune clusters within and across cancer types, and characterization of specific immune cell frequencies as well as T- and B-cell repertoires. We found no associations between immune infiltration levels and tumour mutational burden, although molecular cancer entities were enriched within specific immune clusters. CONCLUSIONS: Given the heterogeneity of immune infiltration within pedNST, our findings suggest personalized immunogenomic profiling is needed to guide selection of immunotherapeutic strategies.

Mapping pediatric brain tumors to their origins in the developing cerebellum.

BACKGROUND: Distinguishing the cellular origins of childhood brain tumors is key for understanding tumor initiation and identifying lineage-restricted, tumor-specific therapeutic targets. Previous strategies to map the cell-of-origin typically involved comparing human tumors to murine embryonal tissues, which is potentially limited due to species-specific differences. The aim of this study was to unravel the cellular origins of the 3 most common pediatric brain tumors, ependymoma, pilocytic astrocytoma, and medulloblastoma, using a developing human cerebellar atlas. METHODS: We used a single-nucleus atlas of the normal developing human cerebellum consisting of 176 645 cells as a reference for an in-depth comparison to 4416 bulk and single-cell transcriptome tumor datasets, using gene set variation analysis, correlation, and single-cell matching techniques. RESULTS: We find that the astroglial cerebellar lineage is potentially the origin for posterior fossa ependymomas. We propose that infratentorial pilocytic astrocytomas originate from the oligodendrocyte lineage and MHC II genes are specifically enriched in these tumors. We confirm that SHH and Group 3/4 medulloblastomas originate from the granule cell and unipolar brush cell lineages. Radiation-induced gliomas stem from cerebellar glial lineages and demonstrate distinct origins from the primary medulloblastoma. We identify tumor genes that are expressed in the cerebellar lineage of origin, and genes that are tumor specific; both gene sets represent promising therapeutic targets for future study. CONCLUSION: Based on our results, individual cells within a tumor may resemble different cell types along a restricted developmental lineage. Therefore, we suggest that tumors can arise from multiple cellular states along the cerebellar "lineage of origin."

Mitochondrial DNA mutations in Medulloblastoma.

To date, several studies on genomic events underlying medulloblastoma (MB) biology have expanded our understanding of this tumour entity and led to its division into four groups-WNT, SHH, group 3 (G3) and group 4 (G4). However, there is little information about the relevance of pathogenic mitochondrial DNA (mtDNA) mutations and their consequences across these. In this report, we describe the case of a female patient with MB and a mitochondriopathy, followed by a study of mtDNA variants in MB groups. After being diagnosed with G4 MB, the index patient was treated in line with the HIT 2000 protocol with no indications of relapse after five years. Long-term side effects of treatment were complemented by additional neurological symptoms and elevated lactate levels ten years later, resulting in suspected mitochondrial disease. This was confirmed by identifying a mutation in the MT-TS1 gene which appeared homoplasmic in patient tissue and heteroplasmic in the patient's mother. Motivated by this case, we explored mtDNA mutations across 444 patients from ICGC and HIT cohorts. While there was no statistically significant enrichment of mutations in one MB group, both cohorts encompassed a small group of patients harbouring potentially deleterious mtDNA variants. The case presented here highlights the possible similarities between sequelae caused by MB treatment and neurological symptoms of mitochondrial dysfunction, which may apply to patients across all MB groups. In the context of the current advances in characterising and interpreting mtDNA aberrations, recognising affected patients could enhance our future knowledge regarding the mutations' impact on carcinogenesis and cancer treatment.

Group-specific cellular metabolism in Medulloblastoma.

BACKGROUND: Cancer metabolism influences multiple aspects of tumorigenesis and causes diversity across malignancies. Although comprehensive research has extended our knowledge of molecular subgroups in medulloblastoma (MB), discrete analysis of metabolic heterogeneity is currently lacking. This study seeks to improve our understanding of metabolic phenotypes in MB and their impact on patients' outcomes. METHODS: Data from four independent MB cohorts encompassing 1,288 patients were analysed. We explored metabolic characteristics of 902 patients (ICGC and MAGIC cohorts) on bulk RNA level. Moreover, data from 491 patients (ICGC cohort) were searched for DNA alterations in genes regulating cell metabolism. To determine the role of intratumoral metabolic differences, we examined single-cell RNA-sequencing (scRNA-seq) data from 34 additional patients. Findings on metabolic heterogeneity were correlated to clinical data. RESULTS: Established MB groups exhibit substantial differences in metabolic gene expression. By employing unsupervised analyses, we identified three clusters of group 3 and 4 samples with distinct metabolic features in ICGC and MAGIC cohorts. Analysis of scRNA-seq data confirmed our results of intertumoral heterogeneity underlying the according differences in metabolic gene expression. On DNA level, we discovered clear associations between altered regulatory genes involved in MB development and lipid metabolism. Additionally, we determined the prognostic value of metabolic gene expression in MB and showed that expression of genes involved in metabolism of inositol phosphates and nucleotides correlates with patient survival. CONCLUSION: Our research underlines the biological and clinical relevance of metabolic alterations in MB. Thus, distinct metabolic signatures presented here might be the first step towards future metabolism-targeted therapeutic options.

Outcome of Children and Adolescents With Relapsed/Refractory/Progressive Malignancies Treated With Molecularly Informed Targeted Drugs in the Pediatric Precision Oncology Registry INFORM.

PURPOSE: INFORM is an international pediatric precision oncology registry, prospectively collecting molecular and clinical data of children with recurrent, progressive, or very high-risk malignancies. We have previously identified a subgroup of patients with improved outcomes on the basis of molecular profiling. The present analysis systematically investigates progression-free survival (PFS) and overall survival (OS) of patients receiving matching targeted treatment (MTT) with the most frequently applied drug classes and its correlation with underlying molecular alterations. METHODS: A cohort of 519 patients with relapsed or refractory high-risk malignancies who had completed a follow-up of at least 2 years or shorter in the case of death or loss to follow-up was analyzed. Survival times were compared using the log-rank test. RESULTS: MTT with anaplastic lymphoma kinase (ALK), neurotrophic tyrosine receptor kinase (NTRK), and B-RAF kinase (BRAF) inhibitors showed significantly improved PFS (P = .012) and OS (P = .036) in comparison with conventional treatment or no treatment. However, analysis of the four most commonly applied MTT groups, mitogen-activated protein kinase (MEK- n = 19), cyclin-dependent kinase (CDK- n = 23), other kinase (n = 62), and mammalian-target of rapamycin (mTOR- n = 20) inhibitors, did not reveal differences in PFS or OS compared with conventional treatment or no treatment in patients with similar molecular pathway alterations. We did not observe differences in the type of pathway alterations (eg, copy number alterations, single-nucleotide variants, InDels, gene fusions) addressed by MTT. CONCLUSION: Patients with respective molecular alterations benefit from treatment with ALK, NTRK, and BRAF inhibitors as previously described. No survival benefit was observed with MTT for mutations in the MEK, CDK, other kinase, or mTOR signaling pathways. The noninterventional character of a registry has to be taken into account when interpreting these data and underlines the need for innovative interventional biomarker-driven clinical trials in pediatric oncology.

The potential of PARP as a therapeutic target across pediatric solid malignancies.

BACKGROUND: Pediatric cancer is the leading cause of disease-related death in children and the need for better therapeutic options remains urgent. Due to the limited number of patients, target and drug development for pediatrics is often supplemented by data from studies focused on adult cancers. Recent evidence shows that pediatric cancers possess different vulnerabilities that should be explored independently from adult cancers. METHODS: Using the publicly available Genomics of Drug Sensitivity in Cancer database, we explore therapeutic targets and biomarkers specific to the pediatric solid malignancies Ewing sarcoma, medulloblastoma, neuroblastoma, osteosarcoma, and rhabdomyosarcoma. Results are validated using cell viability assays and high-throughput drug screens are used to identify synergistic combinations. RESULTS: Using published drug screening data, PARP is identified as a drug target of interest across multiple different pediatric malignancies. We validate these findings, and we show that efficacy can be improved when combined with conventional chemotherapeutics, namely topoisomerase inhibitors. Additionally, using gene set enrichment analysis, we identify ribosome biogenesis as a potential biomarker for PARP inhibition in pediatric cancer cell lines. CONCLUSION: Collectively, our results provide evidence to support the further development of PARP inhibition and the combination with TOP1 inhibition as a therapeutic approach in solid pediatric malignancies. Additionally, we propose ribosome biogenesis as a component to PARP inhibitor sensitivity that should be further investigated to help maximize the potential utility of PARP inhibition and combinations across pediatric solid malignancies.

3D genome mapping identifies subgroup-specific chromosome conformations and tumor-dependency genes in ependymoma.

Ependymoma is a tumor of the brain or spinal cord. The two most common and aggressive molecular groups of ependymoma are the supratentorial ZFTA-fusion associated and the posterior fossa ependymoma group A. In both groups, tumors occur mainly in young children and frequently recur after treatment. Although molecular mechanisms underlying these diseases have recently been uncovered, they remain difficult to target and innovative therapeutic approaches are urgently needed. Here, we use genome-wide chromosome conformation capture (Hi-C), complemented with CTCF and H3K27ac ChIP-seq, as well as gene expression and DNA methylation analysis in primary and relapsed ependymoma tumors, to identify chromosomal conformations and regulatory mechanisms associated with aberrant gene expression. In particular, we observe the formation of new topologically associating domains ('neo-TADs') caused by structural variants, group-specific 3D chromatin loops, and the replacement of CTCF insulators by DNA hyper-methylation. Through inhibition experiments, we validate that genes implicated by these 3D genome conformations are essential for the survival of patient-derived ependymoma models in a group-specific manner. Thus, this study extends our ability to reveal tumor-dependency genes by 3D genome conformations even in tumors that lack targetable genetic alterations.

Multiomic neuropathology improves diagnostic accuracy in pediatric neuro-oncology.

The large diversity of central nervous system (CNS) tumor types in children and adolescents results in disparate patient outcomes and renders accurate diagnosis challenging. In this study, we prospectively integrated DNA methylation profiling and targeted gene panel sequencing with blinded neuropathological reference diagnostics for a population-based cohort of more than 1,200 newly diagnosed pediatric patients with CNS tumors, to assess their utility in routine neuropathology. We show that the multi-omic integration increased diagnostic accuracy in a substantial proportion of patients through annotation to a refining DNA methylation class (50%), detection of diagnostic or therapeutically relevant genetic alterations (47%) or identification of cancer predisposition syndromes (10%). Discrepant results by neuropathological WHO-based and DNA methylation-based classification (30%) were enriched in histological high-grade gliomas, implicating relevance for current clinical patient management in 5% of all patients. Follow-up (median 2.5 years) suggests improved survival for patients with histological high-grade gliomas displaying lower-grade molecular profiles. These results provide preliminary evidence of the utility of integrating multi-omics in neuropathology for pediatric neuro-oncology.

Comprehensive analysis of mutational signatures reveals distinct patterns and molecular processes across 27 pediatric cancers.

Analysis of mutational signatures can reveal underlying molecular mechanisms of the processes that have imprinted the somatic mutations found in cancer genomes. Here, we analyze single base substitutions and small insertions and deletions in pediatric cancers encompassing 785 whole-genome sequenced tumors from 27 molecularly defined cancer subtypes. We identified only a small number of mutational signatures active in pediatric cancers, compared with previously analyzed adult cancers. Further, we report a significant difference in the proportion of pediatric tumors showing homologous recombination repair defect signatures compared with previous analyses. In pediatric leukemias, we identified an indel signature, not previously reported, characterized by long insertions in nonrepeat regions, affecting mainly intronic and intergenic regions, but also exons of known cancer genes. We provide a systematic overview of COSMIC v.3 mutational signatures active across pediatric cancers, which is highly relevant for understanding tumor biology and enabling future research in defining biomarkers of treatment response.

Revana: a comprehensive tool for regulatory variant analysis and visualization of cancer genomes.

MOTIVATION: As non-coding driver mutations move more into the focus of cancer research, a comprehensive and easy-to-use software solution for regulatory variant analysis and data visualization is highly relevant. The interpretation of regulatory variants in large tumor genome cohorts requires specialized analysis and visualization of multiple layers of data, including for example breakpoints of structural variants, enhancer elements and additional available gene locus annotation, in the context of changes in gene expression. RESULTS: We introduce a user-friendly tool, Revana (REgulatory Variant ANAlysis), that can aggregate and visually represent regulatory variants from cancer genomes in a gene-centric manner. It requires whole-genome and RNA sequencing data of a cohort of tumor samples and creates interactive HTML reports summarizing the most important regulatory events. AVAILABILITY AND IMPLEMENTATION: Revana is implemented in R and JavaScript. It is available for download as an R package under . Sample results can be viewed under and a short walkthrough is available under . SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Drug sensitivity profiling of 3D tumor tissue cultures in the pediatric precision oncology program INFORM.

The international precision oncology program INFORM enrolls relapsed/refractory pediatric cancer patients for comprehensive molecular analysis. We report a two-year pilot study implementing ex vivo drug sensitivity profiling (DSP) using a library of 75-78 clinically relevant drugs. We included 132 viable tumor samples from 35 pediatric oncology centers in seven countries. DSP was conducted on multicellular fresh tumor tissue spheroid cultures in 384-well plates with an overall mean processing time of three weeks. In 89 cases (67%), sufficient viable tissue was received; 69 (78%) passed internal quality controls. The DSP results matched the identified molecular targets, including BRAF, ALK, MET, and TP53 status. Drug vulnerabilities were identified in 80% of cases lacking actionable (very) high-evidence molecular events, adding value to the molecular data. Striking parallels between clinical courses and the DSP results were observed in selected patients. Overall, DSP in clinical real-time is feasible in international multicenter precision oncology programs.

Heterozygous BRCA1 and BRCA2 and Mismatch Repair Gene Pathogenic Variants in Children and Adolescents With Cancer.

BACKGROUND: Genetic predisposition is has been identified as a cause of cancer, yet little is known about the role of adult cancer predisposition syndromes in childhood cancer. We examined the extent to which heterozygous pathogenic germline variants in BRCA1, BRCA2, PALB2, ATM, CHEK2, MSH2, MSH6, MLH1, and PMS2 contribute to cancer risk in children and adolescents. METHODS: We conducted a meta-analysis of 11 studies that incorporated comprehensive germline testing for children and adolescents with cancer. ClinVar pathogenic or likely pathogenic variants (PVs) in genes of interest were compared with 2 control groups. Results were validated in a cohort of mainly European patients and controls. We employed the Proxy External Controls Association Test to account for different pipelines. RESULTS: Among 3975 children and adolescents with cancer, statistically significant associations with cancer risk were observed for PVs in BRCA1 and 2 (26 PVs vs 63 PVs among 27 501 controls, odds ratio = 2.78, 95% confidence interval = 1.69 to 4.45; P < .001) and mismatch repair genes (19 PVs vs 14 PVs among 27 501 controls, odds ratio = 7.33, 95% confidence interval = 3.64 to 14.82; P <.001). Associations were seen in brain and other solid tumors but not in hematologic neoplasms. We confirmed similar findings in 1664 pediatric cancer patients primarily of European descent. CONCLUSION: These data suggest that heterozygous PVs in BRCA1 and 2 and mismatch repair genes contribute with reduced penetrance to cancer risk in children and adolescents. No changes to predictive genetic testing and surveillance recommendations are required.

FOXR2 Is an Epigenetically Regulated Pan-Cancer Oncogene That Activates ETS Transcriptional Circuits.

Forkhead box R2 (FOXR2) is a forkhead transcription factor located on the X chromosome whose expression is normally restricted to the testis. In this study, we performed a pan-cancer analysis of FOXR2 activation across more than 10,000 adult and pediatric cancer samples and found FOXR2 to be aberrantly upregulated in 70% of all cancer types and 8% of all individual tumors. The majority of tumors (78%) aberrantly expressed FOXR2 through a previously undescribed epigenetic mechanism that involves hypomethylation of a novel promoter, which was functionally validated as necessary for FOXR2 expression and proliferation in FOXR2-expressing cancer cells. FOXR2 promoted tumor growth across multiple cancer lineages and co-opted ETS family transcription circuits across cancers. Taken together, this study identifies FOXR2 as a potent and ubiquitous oncogene that is epigenetically activated across the majority of human cancers. The identification of hijacking of ETS transcription circuits by FOXR2 extends the mechanisms known to active ETS transcription factors and highlights how transcription factor families cooperate to enhance tumorigenesis. SIGNIFICANCE: This work identifies a novel promoter that drives aberrant FOXR2 expression and delineates FOXR2 as a pan-cancer oncogene that specifically activates ETS transcriptional circuits across human cancers. See related commentary by Liu and Northcott, p. 2977.

MiR-1248: a new prognostic biomarker able to identify supratentorial hemispheric pediatric low-grade gliomas patients associated with progression.

BACKGROUND: Pediatric low-grade gliomas (pLGGs), particularly incompletely resected supratentorial tumours, can undergo progression after surgery. However to date, there are no predictive biomarkers for progression. Here, we aimed to identify pLGG-specific microRNA signatures and evaluate their value as a prognostic tool. METHODS: We identified and validated supratentorial incompletey resected pLGG-specific microRNAs in independent cohorts from four European Pediatric Neuro-Oncology Centres. RESULTS: These microRNAs demonstrated high accuracy in differentiating patients with or without progression. Specifically, incompletely resected supratentorial pLGGs with disease progression showed significantly higher miR-1248 combined with lower miR-376a-3p and miR-888-5p levels than tumours without progression. A significant (p < 0.001) prognostic performance for miR-1248 was reported with an area under the curve (AUC) of 1.00. We also highlighted a critical oncogenic role for miR-1248 in gliomas tumours. Indeed, high miR-1248 levels maintain low its validated target genes (CDKN1A (p21)/FRK/SPOP/VHL/MTAP) and consequently sustain the activation of oncogenic pathways. CONCLUSIONS: Altogether, we provide a novel molecular biomarker able to successfully identify pLGG patients associated with disease progression that could support the clinicians in the decision-making strategy, advancing personalized medicine.

The genomic landscape of pediatric renal cell carcinomas.

Pediatric renal cell carcinomas (RCC) differ from their adult counterparts not only in histologic subtypes but also in clinical characteristics and outcome. However, the underlying biology is still largely unclear. For this reason, we performed whole-exome and transcriptome sequencing analyses on a cohort of 25 pediatric RCC patients with various histologic subtypes, including 10 MiT family translocation (MiT) and 10 papillary RCCs. In this cohort of pediatric RCC, we find only limited genomic overlap with adult RCC, even within the same histologic subtype. Recurrent somatic mutations in genes not previously reported in RCC were detected, such as in CCDC168, PLEKHA1, VWF, and MAP3K9. Our papillary pediatric RCCs, which represent the largest cohort to date with comprehensive molecular profiling in this age group, appeared as a distinct genomic subtype differing in terms of gene mutations and gene expression patterns not only from MiT-RCC but also from their adult counterparts.

Predictive modeling of resistance to SMO inhibition in a patient-derived orthotopic xenograft model of SHH medulloblastoma.

Background: Inhibition of the sonic hedgehog (SHH) pathway with Smoothened (SMO) inhibitors is a promising treatment strategy in SHH-activated medulloblastoma, especially in adult patients. However, the problem is that tumors frequently acquire resistance to the treatment. To understand the underlying resistance mechanisms and to find ways to overcome the resistance, preclinical models that became resistant to SMO inhibition are needed. Methods: To induce SMO inhibitor resistant tumors, we have treated a patient-derived xenograft (PDX) model of SHH medulloblastoma, sensitive to SMO inhibition, with 20 mg/kg Sonidegib using an intermitted treatment schedule. Vehicle-treated and resistant models were subjected to whole-genome and RNA sequencing for molecular characterization and target engagement. In vitro drug screens (76 drugs) were performed using Sonidegib-sensitive and -resistant lines to find other drugs to target the resistant lines. One of the top hits was then validated in vivo. Results: Nine independent Sonidegib-resistant PDX lines were generated. Molecular characterization of the resistant models showed that eight models developed missense mutations in SMO and one gained an inactivating point mutation in MEGF8, which acts downstream of SMO as a repressor in the SHH pathway. The in vitro drug screen with Sonidegib-sensitive and -resistant lines identified good efficacy for Selinexor in the resistant line. Indeed, in vivo treatment with Selinexor revealed that it is more effective in resistant than in sensitive models. Conclusions: We report the first human SMO inhibitor resistant medulloblastoma PDX models, which can be used for further preclinical experiments to develop the best strategies to overcome the resistance to SMO inhibitors in patients.