High-throughput genotyping data are useful for making inferences about human evolutionary history. However, the populations sampled to date are unevenly distributed, and some areas (e.g., South and Central Asia) have rarely been sampled in large-scale studies. To assess human genetic variation more evenly, we sampled 296 individuals from 13 worldwide populations that are not covered by previous studies. By combining these samples with a data set from our laboratory and the HapMap II samples, we assembled a final dataset of ~250,000 SNPs in 850 individuals from 40 populations. With more uniform sampling, the estimate of global genetic differentiation (F(ST)) substantially decreases from ~16% with the HapMap II samples to ~11%. A panel of copy number variations typed in the same populations shows patterns of diversity similar to the SNP data, with highest diversity in African populations. This unique sample collection also permits new inferences about human evolutionary history. The comparison of haplotype variation among populations supports a single out-of-Africa migration event and suggests that the founding population of Eurasia may have been relatively large but isolated from Africans for a period of time. We also found a substantial affinity between populations from central Asia (Kyrgyzstani and Mongolian Buryat) and America, suggesting a central Asian contribution to New World founder populations.

Telomerase reverse transcriptase (TERT) promoter mutations were recently shown to drive telomerase activity in various cancer types, including medulloblastoma. However, the clinical and biological implications of TERT mutations in medulloblastoma have not been described. Hence, we sought to describe these mutations and their impact in a subgroup-specific manner. We analyzed the TERT promoter by direct sequencing and genotyping in 466 medulloblastomas. The mutational distributions were determined according to subgroup affiliation, demographics, and clinical, prognostic, and molecular features. Integrated genomics approaches were used to identify specific somatic copy number alterations in TERT promoter-mutated and wild-type tumors. Overall, TERT promoter mutations were identified in 21 % of medulloblastomas. Strikingly, the highest frequencies of TERT mutations were observed in SHH (83 %; 55/66) and WNT (31 %; 4/13) medulloblastomas derived from adult patients. Group 3 and Group 4 harbored this alteration in <5 % of cases and showed no association with increased patient age. The prognostic implications of these mutations were highly subgroup-specific. TERT mutations identified a subset with good and poor prognosis in SHH and Group 4 tumors, respectively. Monosomy 6 was mostly restricted to WNT tumors without TERT mutations. Hallmark SHH focal copy number aberrations and chromosome 10q deletion were mutually exclusive with TERT mutations within SHH tumors. TERT promoter mutations are the most common recurrent somatic point mutation in medulloblastoma, and are very highly enriched in adult SHH and WNT tumors. TERT mutations define a subset of SHH medulloblastoma with distinct demographics, cytogenetics, and outcomes.

DDX3X encodes a DEAD-box family RNA helicase (DDX3) commonly mutated in medulloblastoma, a highly aggressive cerebellar tumor affecting both children and adults. Despite being implicated in several facets of RNA metabolism, the nature and scope of DDX3's interactions with RNA remain unclear. Here, we show DDX3 collaborates extensively with the translation initiation machinery through direct binding to 5'UTRs of nearly all coding RNAs, specific sites on the 18S rRNA, and multiple components of the translation initiation complex. Impairment of translation initiation is also evident in primary medulloblastomas harboring mutations in DDX3X, further highlighting DDX3's role in this process. Arsenite-induced stress shifts DDX3 binding from the 5'UTR into the coding region of mRNAs concomitant with a general reduction of translation, and both the shift of DDX3 on mRNA and decreased translation are blunted by expression of a catalytically-impaired, medulloblastoma-associated DDX3R534H variant. Furthermore, despite the global repression of translation induced by arsenite, translation is preserved on select genes involved in chromatin organization in DDX3R534H-expressing cells. Thus, DDX3 interacts extensively with RNA and ribosomal machinery to help remodel the translation landscape in response to stress, while cancer-related DDX3 variants adapt this response to selectively preserve translation.

Histological grading of choroid plexus tumors (CPTs) remains the best prognostic tool to distinguish between aggressive choroid plexus carcinoma (CPC) and the more benign choroid plexus papilloma (CPP) or atypical choroid plexus papilloma (aCPP); however, these distinctions can be challenging. Standard treatment of CPC is very aggressive and often leads to severe damage to the young child's brain. Therefore, it is crucial to distinguish between CPC and less aggressive entities (CPP or aCPP) to avoid unnecessary exposure of the young patient to neurotoxic therapy. To better stratify CPTs, we utilized DNA methylation (DNAm) to identify prognostic epigenetic biomarkers for CPCs.

Study of the origin and development of cerebellar tumours has been hampered by the complexity and heterogeneity of cerebellar cells that change over the course of development. Here we use single-cell transcriptomics to study more than 60,000 cells from the developing mouse cerebellum and show that different molecular subgroups of childhood cerebellar tumours mirror the transcription of cells from distinct, temporally restricted cerebellar lineages. The Sonic Hedgehog medulloblastoma subgroup transcriptionally mirrors the granule cell hierarchy as expected, while group 3 medulloblastoma resembles Nestin+ stem cells, group 4 medulloblastoma resembles unipolar brush cells, and PFA/PFB ependymoma and cerebellar pilocytic astrocytoma resemble the prenatal gliogenic progenitor cells. Furthermore, single-cell transcriptomics of human childhood cerebellar tumours demonstrates that many bulk tumours contain a mixed population of cells with divergent differentiation. Our data highlight cerebellar tumours as a disorder of early brain development and provide a proximate explanation for the peak incidence of cerebellar tumours in early childhood.

TP53 missense mutations leading to the expression of mutant p53 oncoproteins are frequent driver events during tumorigenesis. p53 mutants promote tumor growth, metastasis and chemoresistance by affecting fundamental cellular pathways and functions. Here, we demonstrate that p53 mutants modify structure and function of the Golgi apparatus, culminating in the increased release of a pro-malignant secretome by tumor cells and primary fibroblasts from patients with Li-Fraumeni cancer predisposition syndrome. Mechanistically, interacting with the hypoxia responsive factor HIF1α, mutant p53 induces the expression of miR-30d, which in turn causes tubulo-vesiculation of the Golgi apparatus, leading to enhanced vesicular trafficking and secretion. The mut-p53/HIF1α/miR-30d axis potentiates the release of soluble factors and the deposition and remodeling of the ECM, affecting mechano-signaling and stromal cells activation within the tumor microenvironment, thereby enhancing tumor growth and metastatic colonization.

We present an extensive assessment of mutation burden through sequencing analysis of >81,000 tumors from pediatric and adult patients, including tumors with hypermutation caused by chemotherapy, carcinogens, or germline alterations. Hypermutation was detected in tumor types not previously associated with high mutation burden. Replication repair deficiency was a major contributing factor. We uncovered new driver mutations in the replication-repair-associated DNA polymerases and a distinct impact of microsatellite instability and replication repair deficiency on the scale of mutation load. Unbiased clustering, based on mutational context, revealed clinically relevant subgroups regardless of the tumors' tissue of origin, highlighting similarities in evolutionary dynamics leading to hypermutation. Mutagens, such as UV light, were implicated in unexpected cancers, including sarcomas and lung tumors. The order of mutational signatures identified previous treatment and germline replication repair deficiency, which improved management of patients and families. These data will inform tumor classification, genetic testing, and clinical trial design.

We conducted integrative somatic-germline analyses by deeply sequencing 864 cancer-associated genes, complete genomes and transcriptomes for 300 mostly previously treated children and adolescents/young adults with cancer of poor prognosis or with rare tumors enrolled in the SickKids Cancer Sequencing (KiCS) program. Clinically actionable variants were identified in 56% of patients. Improved diagnostic accuracy led to modified management in a subset. Therapeutically targetable variants (54% of patients) were of unanticipated timing and type, with over 20% derived from the germline. Corroborating mutational signatures (SBS3/BRCAness) in patients with germline homologous recombination defects demonstrates the potential utility of PARP inhibitors. Mutational burden was significantly elevated in 9% of patients. Sequential sampling identified changes in therapeutically targetable drivers in over one-third of patients, suggesting benefit from rebiopsy for genomic analysis at the time of relapse. Comprehensive cancer genomic profiling is useful at multiple points in the care trajectory for children and adolescents/young adults with cancer, supporting its integration into early clinical management.

Approximately 20 TP53 retrogenes exist in the African and Asian elephant genomes (Loxodonta Africana, Elephas Maximus) in addition to a conserved TP53 gene that encodes a full-length protein. Elephant TP53-RETROGENE 9 (TP53-R9) encodes a p53 protein (p53-R9) that is truncated in the middle of the canonical DNA binding domain. This C-terminally truncated p53 retrogene protein lacks the nuclear localization signals and oligomerization domain of its full-length counterpart. When expressed in human osteosarcoma cells (U2OS), p53-R9 binds to Tid1, the chaperone protein responsible for mitochondrial translocation of human p53 in response to cellular stress. Tid1 expression is required for p53-R9-induced apoptosis. At the mitochondria, p53-R9 binds to the pro-apoptotic BCL-2 family member Bax, which leads to caspase activation, cytochrome c release, and cell death. Our data show, for the first time, that expression of this truncated elephant p53 retrogene protein induces apoptosis in human cancer cells. Understanding the molecular mechanism by which the additional elephant TP53 retrogenes function may provide evolutionary insight that can be utilized for the development of therapeutics to treat human cancers.

Cancers arising from germline DNA mismatch repair deficiency or polymerase proofreading deficiency (MMRD and PPD) in children harbour the highest mutational and microsatellite insertion-deletion (MS-indel) burden in humans. MMRD and PPD cancers are commonly lethal due to the inherent resistance to chemo-irradiation. Although immune checkpoint inhibitors (ICIs) have failed to benefit children in previous studies, we hypothesized that hypermutation caused by MMRD and PPD will improve outcomes following ICI treatment in these patients. Using an international consortium registry study, we report on the ICI treatment of 45 progressive or recurrent tumors from 38 patients. Durable objective responses were observed in most patients, culminating in a 3 year survival of 41.4%. High mutation burden predicted response for ultra-hypermutant cancers (>100 mutations per Mb) enriched for combined MMRD + PPD, while MS-indels predicted response in MMRD tumors with lower mutation burden (10-100 mutations per Mb). Furthermore, both mechanisms were associated with increased immune infiltration even in 'immunologically cold' tumors such as gliomas, contributing to the favorable response. Pseudo-progression (flare) was common and was associated with immune activation in the tumor microenvironment and systemically. Furthermore, patients with flare who continued ICI treatment achieved durable responses. This study demonstrates improved survival for patients with tumors not previously known to respond to ICI treatment, including central nervous system and synchronous cancers, and identifies the dual roles of mutation burden and MS-indels in predicting sustained response to immunotherapy.

Functional heterogeneity within tumors presents a significant therapeutic challenge. Here we show that quiescent, therapy-resistant Sox2(+) cells propagate sonic hedgehog subgroup medulloblastoma by a mechanism that mirrors a neurogenic program. Rare Sox2(+) cells produce rapidly cycling doublecortin(+) progenitors that, together with their postmitotic progeny expressing NeuN, comprise tumor bulk. Sox2(+) cells are enriched following anti-mitotic chemotherapy and Smoothened inhibition, creating a reservoir for tumor regrowth. Lineage traces from Sox2(+) cells increase following treatment, suggesting that this population is responsible for relapse. Targeting Sox2(+) cells with the antineoplastic mithramycin abrogated tumor growth. Addressing functional heterogeneity and eliminating Sox2(+) cells presents a promising therapeutic paradigm for treatment of sonic hedgehog subgroup medulloblastoma.

Choroid plexus carcinomas (CPCs) are poorly understood and frequently lethal brain tumors with few treatment options. Using a mouse model of the disease and a large cohort of human CPCs, we performed a cross-species, genome-wide search for oncogenes within syntenic regions of chromosome gain. TAF12, NFYC, and RAD54L co-located on human chromosome 1p32-35.3 and mouse chromosome 4qD1-D3 were identified as oncogenes that are gained in tumors in both species and required for disease initiation and progression. TAF12 and NFYC are transcription factors that regulate the epigenome, whereas RAD54L plays a central role in DNA repair. Our data identify a group of concurrently gained oncogenes that cooperate in the formation of CPC and reveal potential avenues for therapy.

Spatial heterogeneity of transcriptional and genetic markers between physically isolated biopsies of a single tumor poses major barriers to the identification of biomarkers and the development of targeted therapies that will be effective against the entire tumor. We analyzed the spatial heterogeneity of multiregional biopsies from 35 patients, using a combination of transcriptomic and genomic profiles. Medulloblastomas (MBs), but not high-grade gliomas (HGGs), demonstrated spatially homogeneous transcriptomes, which allowed for accurate subgrouping of tumors from a single biopsy. Conversely, somatic mutations that affect genes suitable for targeted therapeutics demonstrated high levels of spatial heterogeneity in MB, malignant glioma, and renal cell carcinoma (RCC). Actionable targets found in a single MB biopsy were seldom clonal across the entire tumor, which brings the efficacy of monotherapies against a single target into question. Clinical trials of targeted therapies for MB should first ensure the spatially ubiquitous nature of the target mutation.

We have optimized point mutation knock-ins into zebrafish genomic sites using clustered regularly interspaced palindromic repeats (CRISPR)/Cas9 reagents and single-stranded oligodeoxynucleotides. The efficiency of knock-ins was assessed by a novel application of allele-specific polymerase chain reaction and confirmed by high-throughput sequencing. Anti-sense asymmetric oligo design was found to be the most successful optimization strategy. However, cut site proximity to the mutation and phosphorothioate oligo modifications also greatly improved knock-in efficiency. A previously unrecognized risk of off-target trans knock-ins was identified that we obviated through the development of a workflow for correct knock-in detection. Together these strategies greatly facilitate the study of human genetic diseases in zebrafish, with additional applicability to enhance CRISPR-based approaches in other animal model systems.

In cancer, recurrent somatic single-nucleotide variants-which are rare in most paediatric cancers-are confined largely to protein-coding genes1-3. Here we report highly recurrent hotspot mutations (r.3A>G) of U1 spliceosomal small nuclear RNAs (snRNAs) in about 50% of Sonic hedgehog (SHH) medulloblastomas. These mutations were not present across other subgroups of medulloblastoma, and we identified these hotspot mutations in U1 snRNA in only <0.1% of 2,442 cancers, across 36 other tumour types. The mutations occur in 97% of adults (subtype SHHδ) and 25% of adolescents (subtype SHHα) with SHH medulloblastoma, but are largely absent from SHH medulloblastoma in infants. The U1 snRNA mutations occur in the 5' splice-site binding region, and snRNA-mutant tumours have significantly disrupted RNA splicing and an excess of 5' cryptic splicing events. Alternative splicing mediated by mutant U1 snRNA inactivates tumour-suppressor genes (PTCH1) and activates oncogenes (GLI2 and CCND2), and represents a target for therapy. These U1 snRNA mutations provide an example of highly recurrent and tissue-specific mutations of a non-protein-coding gene in cancer.

Loss of p53 promotes vascular endothelial growth factor (VEGF)-A up-regulation and the angiogenic potential of cancer cells. We investigated TP53 somatic mutations in 110 primary gastric adenocarcinomas of two retrospective metastatic series including 48 patients treated with second-line Ramucirumab/Paclitaxel and 62 patients who received first-line chemotherapy with Cisplatin or Oxaliplatin plus 5-Fluorouracil. Missense mutations were classified by tumor protein p53 (TP53) mutant-specific residual transcriptional activity scores (TP53RTAS) and used to stratify patients into two groups: transcriptionally TP53Active and TP53Inactive. The primary endpoint was overall survival (OS). An additional analysis was addressed to measure VEGF/VEGF receptor 2 (VEGFR2) expression levels in relation to the TP53RTAS. In the Ramucirumab/Paclitaxel group, 29/48 (60.4%) patients had TP53 mutations. Ten patients with TP53Inactive mutations showed better OS than carriers of other TP53 mutations. This effect was retained in the multivariate model analysis (Hazard Ratio = 0.29, 95% confidence interval = 0.17-0.85, p = 0.02). In the chemotherapy group, 41/62 (66%) patients had TP53 mutations, and the 11 carriers of TP53Inactive mutations showed the worst OS (Hazard Ratio = 2.64, 95% confidence interval = 1.17-5.95, p = 0.02). VEGF-A mRNA expression levels were significantly increased in TP53Inactive cases. Further studies are warranted to explore the effect of TP53Inactive mutations in different anti-cancer regimens. This information would lead to new tailored therapy strategies for this lethal disease.

Diffuse intrinsic pontine glioma (DIPG) and midline high-grade glioma (mHGG) are lethal childhood brain tumors. Spatial genomic heterogeneity has been well-described in adult HGG but has not been comprehensively characterized in pediatric HGG. We performed whole exome sequencing on 38-matched primary, contiguous, and metastatic tumor sites from eight children with DIPG (n = 7) or mHGG (n = 1) collected using a unique MRI-guided autopsy protocol. Validation was performed using Sanger sequencing, Droplet Digital polymerase-chain reaction, immunohistochemistry, and fluorescent in-situ hybridization.

microRNA-34A is a critical component of the p53 network and expression of miR- 34A is down-regulated by promoter hypermethylation or focal deletions in numerous human cancers. Although miR-34A deregulation may be an important driver in cancer, the endogenous role of this microRNA in cellular homeostasis is not well characterized. To address this knowledge gap, we aimed to determine the transcriptional landscape of the miR-34A-p53 axis in non-transformed cells. Using primary skin-derived fibroblast cell lines from patients who developed childhood cancers, and who harbor either germline TP53 mutations or are TP53 wild type, we sought to characterize the transcriptional response to miR-34A modulation. Through transcriptome-wide RNA-Sequencing, we show for the first time that in human non- transformed cells harboring TP53 mutations, miR-34A functions in a noncanonical manner to influence noncoding RNA networks, including RNA components of the minor (U12) spliceosome, as well as TP53-dependent and independent epigenetic pathways. miR- 34A-regulated transcripts include known cell cycle mediators and abrogation of miR-34A leads to a TP53-dependent increase in the fraction of cells in G2/M. Collectively, these results provide a framework for understanding the endogenous role of the miR-34A signaling axis and identify novel transcripts and pathways regulated by the essential miR-34A-p53 tumor suppressor network.

Li-Fraumeni syndrome (LFS) is a hereditary cancer predisposition syndrome associated with germline TP53 pathogenic variants. Here, we perform whole-genome sequence (WGS) analysis of tumors from 22 patients with TP53 germline pathogenic variants. We observe somatic mutations affecting Wnt, PI3K/AKT signaling, epigenetic modifiers and homologous recombination genes as well as mutational signatures associated with prior chemotherapy. We identify near-ubiquitous early loss of heterozygosity of TP53, with gain of the mutant allele. This occurs earlier in these tumors compared to tumors with somatic TP53 mutations, suggesting the timing of this mark may distinguish germline from somatic TP53 mutations. Phylogenetic trees of tumor evolution, reconstructed from bulk and multi-region WGS, reveal that LFS tumors exhibit comparatively limited heterogeneity. Overall, our study delineates early copy number gains of mutant TP53 as a characteristic mutational process in LFS tumorigenesis, likely arising years prior to tumor diagnosis.

Radiation-induced high-grade gliomas (RIGs) are an incurable late complication of cranial radiation therapy. We performed DNA methylation profiling, RNA-seq, and DNA sequencing on 32 RIG tumors and an in vitro drug screen in two RIG cell lines. We report that based on DNA methylation, RIGs cluster primarily with the pediatric receptor tyrosine kinase I high-grade glioma subtype. Common copy-number alterations include Chromosome (Ch.) 1p loss/1q gain, and Ch. 13q and Ch. 14q loss; focal alterations include PDGFRA and CDK4 gain and CDKN2A and BCOR loss. Transcriptomically, RIGs comprise a stem-like subgroup with lesser mutation burden and Ch. 1p loss and a pro-inflammatory subgroup with greater mutation burden and depleted DNA repair gene expression. Chromothripsis in several RIG samples is associated with extrachromosomal circular DNA-mediated amplification of PDGFRA and CDK4. Drug screening suggests microtubule inhibitors/stabilizers, DNA-damaging agents, MEK inhibition, and, in the inflammatory subgroup, proteasome inhibitors, as potentially effective therapies.