As whole-genome sequencing for cancer genome analysis becomes a clinical tool, a full understanding of the variables affecting sequencing analysis output is required. Here using tumour-normal sample pairs from two different types of cancer, chronic lymphocytic leukaemia and medulloblastoma, we conduct a benchmarking exercise within the context of the International Cancer Genome Consortium. We compare sequencing methods, analysis pipelines and validation methods. We show that using PCR-free methods and increasing sequencing depth to ∼ 100 × shows benefits, as long as the tumour:control coverage ratio remains balanced. We observe widely varying mutation call rates and low concordance among analysis pipelines, reflecting the artefact-prone nature of the raw data and lack of standards for dealing with the artefacts. However, we show that, using the benchmark mutation set we have created, many issues are in fact easy to remedy and have an immediate positive impact on mutation detection accuracy.

TP53 mutations confer subgroup specific poor survival for children with medulloblastoma. We hypothesized that WNT activation which is associated with improved survival for such children abrogates TP53 related radioresistance and can be used to sensitize TP53 mutant tumors for radiation. We examined the subgroup-specific role of TP53 mutations in a cohort of 314 patients treated with radiation. TP53 wild-type or mutant human medulloblastoma cell-lines and normal neural stem cells were used to test radioresistance of TP53 mutations and the radiosensitizing effect of WNT activation on tumors and the developing brain. Children with WNT/TP53 mutant medulloblastoma had higher 5-year survival than those with SHH/TP53 mutant tumours (100% and 36.6%±8.7%, respectively (p<0.001)). Introduction of TP53 mutation into medulloblastoma cells induced radioresistance (survival fractions at 2Gy (SF2) of 89%±2% vs. 57.4%±1.8% (p<0.01)). In contrast, β-catenin mutation sensitized TP53 mutant cells to radiation (p<0.05). Lithium, an activator of the WNT pathway, sensitized TP53 mutant medulloblastoma to radiation (SF2 of 43.5%±1.5% in lithium treated cells vs. 56.6±3% (p<0.01)) accompanied by increased number of γH2AX foci. Normal neural stem cells were protected from lithium induced radiation damage (SF2 of 33%±8% for lithium treated cells vs. 27%±3% for untreated controls (p=0.05). Poor survival of patients with TP53 mutant medulloblastoma may be related to radiation resistance. Since constitutive activation of the WNT pathway by lithium sensitizes TP53 mutant medulloblastoma cells and protect normal neural stem cells from radiation, this oral drug may represent an attractive novel therapy for high-risk medulloblastomas.

DNA methylation plays a key role in regulating eukaryotic gene expression. Although mitotically heritable and stable over time, patterns of DNA methylation frequently change in response to cell differentiation, disease and environmental influences. Several methods have been developed to map DNA methylation on a genomic scale. Here, we benchmark four of these approaches by analyzing two human embryonic stem cell lines derived from genetically unrelated embryos and a matched pair of colon tumor and adjacent normal colon tissue obtained from the same donor. Our analysis reveals that methylated DNA immunoprecipitation sequencing (MeDIP-seq), methylated DNA capture by affinity purification (MethylCap-seq), reduced representation bisulfite sequencing (RRBS) and the Infinium HumanMethylation27 assay all produce accurate DNA methylation data. However, these methods differ in their ability to detect differentially methylated regions between pairs of samples. We highlight strengths and weaknesses of the four methods and give practical recommendations for the design of epigenomic case-control studies.

We collated data from 157 unpublished cases of pediatric high-grade glioma and diffuse intrinsic pontine glioma and 20 publicly available datasets in an integrated analysis of >1,000 cases. We identified co-segregating mutations in histone-mutant subgroups including loss of FBXW7 in H3.3G34R/V, TOP3A rearrangements in H3.3K27M, and BCOR mutations in H3.1K27M. Histone wild-type subgroups are refined by the presence of key oncogenic events or methylation profiles more closely resembling lower-grade tumors. Genomic aberrations increase with age, highlighting the infant population as biologically and clinically distinct. Uncommon pathway dysregulation is seen in small subsets of tumors, further defining the molecular diversity of the disease, opening up avenues for biological study and providing a basis for functionally defined future treatment stratification.

Establishment of telomere maintenance mechanisms is a universal step in tumor development to achieve replicative immortality. These processes leave molecular footprints in cancer genomes in the form of altered telomere content and aberrations in telomere composition. To retrieve these telomere characteristics from high-throughput sequencing data the available computational approaches need to be extended and optimized to fully exploit the information provided by large scale cancer genome data sets.

Tuberous sclerosis complex (TSC) is a monogenic disorder caused by mutations in either the TSC1 or TSC2 gene, two key regulators of the mechanistic target of the rapamycin complex pathway. Phenotypically, this leads to growth and formation of hamartomas in several organs, including the brain. Subependymal giant cell astrocytomas (SEGAs) are low-grade brain tumors commonly associated with TSC. Recently, gene expression studies provided evidence that the immune system, the MAPK pathway and extracellular matrix organization play an important role in SEGA development. However, the precise mechanisms behind the gene expression changes in SEGA are still largely unknown, providing a potential role for DNA methylation. We investigated the methylation profile of SEGAs using the Illumina Infinium HumanMethylation450 BeadChip (SEGAs n = 42, periventricular control n = 8). The SEGA methylation profile was enriched for the adaptive immune system, T cell activation, leukocyte mediated immunity, extracellular structure organization and the ERK1 & ERK2 cascade. More interestingly, we identified two subgroups in the SEGA methylation data and show that the differentially expressed genes between the two subgroups are related to the MAPK cascade and adaptive immune response. Overall, this study shows that the immune system, the MAPK pathway and extracellular matrix organization are also affected on DNA methylation level, suggesting that therapeutic intervention on DNA level could be useful for these specific pathways in SEGA. Moreover, we identified two subgroups in SEGA that seem to be driven by changes in the adaptive immune response and MAPK pathway and could potentially hold predictive information on target treatment response.

High throughput screening methods, measuring the sensitivity and resistance of tumor cells to drug treatments have been rapidly evolving. Not only do these screens allow correlating response profiles to tumor genomic features for developing novel predictors of treatment response, but they can also add evidence for therapy decision making in precision oncology. Recent analysis methods developed for either assessing single agents or combination drug efficacies enable quantification of dose-response curves with restricted symmetric fit settings. Here, we introduce iTReX, a user-friendly and interactive Shiny/R application, for both the analysis of mono- and combination therapy responses. The application features an extended version of the drug sensitivity score (DSS) based on the integral of an advanced five-parameter dose-response curve model and a differential DSS for combination therapy profiling. Additionally, iTReX includes modules that visualize drug target interaction networks and support the detection of matches between top therapy hits and the sample omics features to enable the identification of druggable targets and biomarkers. iTReX enables the analysis of various quantitative drug or therapy response readouts (e.g. luminescence, fluorescence microscopy) and multiple treatment strategies (drug treatments, radiation). Using iTReX we validate a cost-effective drug combination screening approach and reveal the application's ability to identify potential sample-specific biomarkers based on drug target interaction networks. The iTReX web application is accessible at https://itrex.kitz-heidelberg.de.

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.

BRAF is a cytoplasmic protein kinase, which activates the MEK-ERK signalling pathway. Deregulation of the pathway is associated with the presence of BRAF mutations in human cancer, the most common being V600E BRAF, although structural rearrangements, which remove N-terminal regulatory sequences, have also been reported. RAF-MEK-ERK signalling is normally thought to occur in the cytoplasm of the cell. However, in an investigation of BRAF localisation using fluorescence microscopy combined with subcellular fractionation of Green Fluorescent Protein (GFP)-tagged proteins expressed in NIH3T3 cells, surprisingly, we detected N-terminally truncated BRAF (ΔBRAF) in both nuclear and cytoplasmic compartments. In contrast, ΔCRAF and full-length, wild-type BRAF (WTBRAF) were detected at lower levels in the nucleus while full-length V600EBRAF was virtually excluded from this compartment. Similar results were obtained using ΔBRAF tagged with the hormone-binding domain of the oestrogen receptor (hbER) and with the KIAA1549-ΔBRAF translocation mutant found in human pilocytic astrocytomas. Here we show that GFP-ΔBRAF nuclear translocation does not involve a canonical Nuclear Localisation Signal (NLS), but is suppressed by N-terminal sequences. Nuclear GFP-ΔBRAF retains MEK/ERK activating potential and is associated with the accumulation of phosphorylated MEK and ERK in the nucleus. In contrast, full-length GFP-WTBRAF and GFP-V600EBRAF are associated with the accumulation of phosphorylated ERK but not phosphorylated MEK in the nucleus. These data have implications for cancers bearing single nucleotide variants or N-terminal deleted structural variants of BRAF.

Chromothripsis and chromoanasynthesis are catastrophic events leading to clustered genomic rearrangements. Whole-genome sequencing revealed frequent complex genomic rearrangements (n = 16/26) in brain tumors developing in mice deficient for factors involved in homologous-recombination-repair or non-homologous-end-joining. Catastrophic events were tightly linked to Myc/Mycn amplification, with increased DNA damage and inefficient apoptotic response already observable at early postnatal stages. Inhibition of repair processes and comparison of the mouse tumors with human medulloblastomas (n = 68) and glioblastomas (n = 32) identified chromothripsis as associated with MYC/MYCN gains and with DNA repair deficiencies, pointing towards therapeutic opportunities to target DNA repair defects in tumors with complex genomic rearrangements.

Diffuse intrinsic pontine glioma (DIPG) is an aggressive childhood tumor of the brainstem with currently no curative treatment available. The vast majority of DIPGs carry a histone H3 mutation leading to a lysine 27-to-methionine exchange (H3K27M). We engineered human induced pluripotent stem cells (iPSCs) to carry an inducible H3.3-K27M allele in the endogenous locus and studied the effects of the mutation in different disease-relevant neural cell types. H3.3-K27M upregulated bivalent promoter-associated developmental genes, producing diverse outcomes in different cell types. While being fatal for iPSCs, H3.3-K27M increased proliferation in neural stem cells (NSCs) and to a lesser extent in oligodendrocyte progenitor cells (OPCs). Only NSCs gave rise to tumors upon induction of H3.3-K27M and TP53 inactivation in an orthotopic xenograft model recapitulating human DIPGs. In NSCs, H3.3-K27M leads to maintained expression of stemness and proliferative genes and a premature activation of OPC programs that together may cause tumor initiation.

Embryonal rhabdomyosarcoma (ERMS) of the uterus has recently been shown to frequently harbor DICER1 mutations. Interestingly, only rare cases of extrauterine DICER1-associated ERMS, mostly located in the genitourinary tract, have been reported to date. Our goal was to study clinicopathologic and molecular profiles of DICER1-mutant (DICER1-mut) and DICER1-wild type (DICER1-wt) ERMS in a cohort of genitourinary tumors. We collected a cohort of 17 ERMS including nine uterine (four uterine corpus and five cervix), one vaginal, and seven urinary tract tumors. DNA sequencing revealed mutations of DICER1 in 9/9 uterine ERMS. All other ERMS of our cohort were DICER1-wt. The median age at diagnosis of patients with DICER1-mut and DICER1-wt ERMS was 36 years and 5 years, respectively. Limited follow-up data (available for 15/17 patients) suggested that DICER1-mut ERMS might show a less aggressive clinical course than DICER1-wt ERMS. Histological features only observed in DICER1-mut ERMS were cartilaginous nodules (6/9 DICER1-mut ERMS), in one case accompanied by foci of ossification. Recurrent mutations identified in both DICER1-mut and DICER1-wt ERMS affected KRAS, NRAS, and TP53. Copy number analysis revealed similar structural variations with frequent gains on chromosomes 2, 3, and 8, independent of DICER1 mutation status. Unsupervised hierarchical clustering of array-based whole-genome DNA methylation data of our study cohort together with an extended methylation data set including different RMS subtypes from genitourinary and extra-genitourinary locations (n = 102), revealed a distinct cluster for DICER1-mut ERMS. Such tumors clearly segregated from the clusters of DICER1-wt ERMS, alveolar RMS, and MYOD1-mutant spindle cell and sclerosing RMS. Only one tumor, previously diagnosed as ERMS arising in the maxilla of a 6-year-old boy clustered with DICER1-mut ERMS of the uterus. Subsequent sequencing analysis identified two DICER1 mutations in the latter case. Our results suggest that DICER1-mut ERMS might qualify as a distinct subtype in future classifications of RMS.

Histone 3.3 (H3.3) hotspot mutations in bone tumors occur in the vast majority of giant cell tumors of bone (GCTBs; 96%), chondroblastomas (95%) and in a few cases of osteosarcomas. However, clinical presentation, histopathological features, and additional molecular characteristics of H3.3 mutant osteosarcomas are largely unknown.

The identification of a heritable tumor predisposition often leads to changes in management and increased surveillance of individuals who are at risk; however, for many rare entities, our knowledge of heritable predisposition is incomplete.

Pediatric central nervous system (CNS) tumors represent the most common cause of cancer-related death in children aged 0-14 years. They differ from their adult counterparts, showing extensive clinical and molecular heterogeneity as well as a challenging histopathological spectrum that often impairs accurate diagnosis. Here, we use DNA methylation-based CNS tumor classification in combination with copy number, RNA-seq, and ChIP-seq analysis to characterize a newly identified CNS tumor type. In addition, we report histology, patient characteristics, and survival data in this tumor type. We describe a biologically distinct pediatric CNS tumor type (n = 31 cases) that is characterized by focal high-level amplification and resultant overexpression of either PLAGL1 or PLAGL2, and an absence of recurrent genetic alterations characteristic of other pediatric CNS tumor types. Both genes act as transcription factors for a regulatory subset of imprinted genes (IGs), components of the Wnt/β-Catenin pathway, and the potential drug targets RET and CYP2W1, which are also specifically overexpressed in this tumor type. A derived PLAGL-specific gene expression signature indicates dysregulation of imprinting control and differentiation/development. These tumors occurred throughout the neuroaxis including the cerebral hemispheres, cerebellum, and brainstem, and were predominantly composed of primitive embryonal-like cells lacking robust expression of markers of glial or neuronal differentiation (e.g., GFAP, OLIG2, and synaptophysin). Tumors with PLAGL1 amplification were typically diagnosed during adolescence (median age 10.5 years), whereas those with PLAGL2 amplification were diagnosed during early childhood (median age 2 years). The 10-year overall survival was 66% for PLAGL1-amplified tumors, 25% for PLAGL2-amplified tumors, 18% for male patients, and 82% for female patients. In summary, we describe a new type of biologically distinct CNS tumor characterized by PLAGL1/2 amplification that occurs predominantly in infants and toddlers (PLAGL2) or adolescents (PLAGL1) which we consider best classified as a CNS embryonal tumor and which is associated with intermediate survival. The cell of origin and optimal treatment strategies remain to be defined.

The diagnosis of sinonasal tumors is challenging due to a heterogeneous spectrum of various differential diagnoses as well as poorly defined, disputed entities such as sinonasal undifferentiated carcinomas (SNUCs). In this study, we apply a machine learning algorithm based on DNA methylation patterns to classify sinonasal tumors with clinical-grade reliability. We further show that sinonasal tumors with SNUC morphology are not as undifferentiated as their current terminology suggests but rather reassigned to four distinct molecular classes defined by epigenetic, mutational and proteomic profiles. This includes two classes with neuroendocrine differentiation, characterized by IDH2 or SMARCA4/ARID1A mutations with an overall favorable clinical course, one class composed of highly aggressive SMARCB1-deficient carcinomas and another class with tumors that represent potentially previously misclassified adenoid cystic carcinomas. Our findings can aid in improving the diagnostic classification of sinonasal tumors and could help to change the current perception of SNUCs.

Although pediatric low-grade gliomas (pLGG) are the most common pediatric brain tumors, patient-derived cell lines reflecting pLGG biology in culture are scarce. This also applies to the most common pLGG subtype pilocytic astrocytoma (PA). Conventional cell culture approaches adapted from higher-grade tumors fail in PA due to oncogene-induced senescence (OIS) driving tumor cells into arrest. Here, we describe a PA modeling workflow using the Simian Virus large T antigen (SV40-TAg) to circumvent OIS.

Pediatric high-grade gliomas of the subclass MYCN (HGG-MYCN) are highly aggressive tumors frequently carrying MYCN amplifications, TP53 mutations, or both alterations. Due to their rarity, such tumors have only recently been identified as a distinct entity, and biological as well as clinical characteristics have not been addressed specifically. To gain insights into tumorigenesis and molecular profiles of these tumors, and to ultimately suggest alternative treatment options, we generated a genetically engineered mouse model by breeding hGFAP-cre::Trp53Fl/Fl::lsl-MYCN mice. All mice developed aggressive forebrain tumors early in their lifetime that mimic human HGG-MYCN regarding histology, DNA methylation, and gene expression. Single-cell RNA sequencing revealed a high intratumoral heterogeneity with neuronal and oligodendroglial lineage signatures. High-throughput drug screening using both mouse and human tumor cells finally indicated high efficacy of Doxorubicin, Irinotecan, and Etoposide as possible therapy options that children with HGG-MYCN might benefit from.

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.

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.