{'fraction_non_alphanumeric': 0.030936236797496415, 'fraction_numerical': 0.015419220237319077, 'mean_word_length': 5.319942315615987, 'pattern_counts': {'":': 0, '<': 1, '<?xml version=': 0, '>': 2, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 1, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Childhood Cerebellar Tumors Mirror Conserved Fetal Transcriptional Programs

Summary: The study of the origin and development of cerebellar tumors has been hampered by the complexity and heterogeneity of cerebellar cells that change over the course of development. We used single-cell transcriptomics to study >60,000 cells from the developing murine cerebellum, and show that different molecular subgroups of childhood cerebellar tumors mirror the transcription of cells from distinct, temporally restricted cerebellar lineages. Sonic Hedgehog medulloblastoma transcriptionally mirrors the granule cell hierarchy as expected, whereas Group3 MB resemble Nestin+ve stem cells, Group 4 MB resemble unipolar brush cells, and PFA/PFB ependymoma and cerebellar pilocytic astrocytoma resemble the pre-natal gliogenic progenitor cells. Furthermore, single-cell transcriptomics of human childhood cerebellar tumors demonstrates that many bulk tumors contain a mixed population of cells with divergent differentiation. Our data highlight cerebellar tumors as a disorder of early brain development, and provide a proximate explanation for the peak incidence of cerebellar tumors in early childhood.


Introduction
Pediatric brain tumors are most commonly found in the posterior fossa, particularly the cerebellum, where medulloblastoma, ependymoma, and pilocytic astrocytoma are most prevalent. Cerebellar tumors are currently treated using non-specific therapies. They have few somatically mutated driver genes 2,17,[22][23][24][25][26] which has impeded the understanding of their biology and resulted in the development of targeted therapies lagging behind other major cancer types. Nevertheless, childhood cerebellar tumor types are known to have distinct molecular subtypes with different clinical behaviours. 'Medulloblastoma' is now known to comprise four molecularly distinct diseases (subgroups), with further clinical and molecular heterogeneity within each subgroup [1][2][3][4][5] . Sonic Hedgehog (Shh-MB), Group 3 (Grp3-MB), and Group 4 medulloblastoma (Grp4-MB) are thought to originate in the cerebellum [6][7][8][9][10][11] , while Wnt-MB arises from the lower rhombic lip of the developing brain stem 12 . Although Shh-MB is thought to arise from the granule cell precursors (GCPs), careful examination of fully developed Shh-MB revealed populations of cells displaying varying levels of differentiation and growth capacity, mirroring the temporal evolution of the developing GCP hierarchy 13 . It is currently unclear to what extent the other molecular subgroups of medulloblastoma recapitulate a similar developmental hierarchy. Ependymomas are found throughout the central nervous system, but in the cerebellum are thought to be largely limited to two variants: PFA and PFB [14][15][16] , and have been suggested to arise from regional radial glial-like cells [17][18][19] . Molecular subgroups of medulloblastoma and ependymoma have been delineated by transcriptomics as well as patterns of DNA CpG methylation, both of which have been suggested to reflect the cell of origin 20,21 .
The cerebellum is made up of a large variety of cell types, with many undergoing temporally regulated differentiation through defined developmental hierarchies [27][28][29] . GABAergic neurons, including Purkinje cells and a variety of interneurons, arise from the ventricular zone (VZ), while glutamatergic neurons, including those of the cerebellar nuclei (CN), the inner granule cell layer, and the unipolar brush cells (UBCs) arise from the upper rhombic lip (URL). Cerebellar glial cells, including radial glia, astrocytes, and Bergmann glia also arise from stem cells in the VZ that produce a proliferating progenitor still present in the cerebellar cortex after birth [30][31][32] . In the past, transcriptional studies of bulk cerebellar tissue were performed on a complex mixture of GABAergic neurons, glutamatergic neurons, glia, and non-neuronal cells. The use of mixed transcriptomes from normal bulk cerebellum precludes a meaningful comparison to the transcriptome or epigenome of cerebellar tumors. Massive changes during early development, and the relative inaccessibility of the cerebellum inside the skull further complicate the study of the normal developing cerebellum, particularly from human samples. However, the recent development of large-scale single cell RNA sequencing (scRNA-seq) permits the development of a 'cellular scaffold' 53,54 for cerebellar development in which the transcriptomes of distinct hierarchies can be determined at various points in time, and subsequently compared to the transcriptomes of childhood cerebellar tumors. Identification of temporally and lineage restricted cell populations in the developing cerebellum that most closely mirror the transcriptome of cerebellar neoplasms could allow for identification of tumor cells of origin, as transcriptional similarity can be construed as good evidence for the lineage of cellular origin. Although it remains possible that a more differentiated cell could 'de-differentiate,' or that a cell from another lineage hierarchy could 'trans-differentiate' during the process of transformation 32 there is scant evidence for de-differentiation or trans-differentiation in human brain tumors. Furthermore, because of the low mutational burden in cerebellar tumors, it would seem likely that the tumor remains similar to its cell of origin. Benefits of matching cerebellar tumors to their transcriptional best match from cerebellar development include the discovery of developmental checkpoints defective in cerebellar tumors, the development of clinically relevant mouse models, and direct comparison of tumor and normal cell transcriptomes to discover tumor-related changes in transcriptional programs. Insights into the process of transformation from a normal cerebellar cell to a cerebellar cancer cell could lead to the development of novel targeted therapies.

Identification of transcriptional clusters in the developing murine cerebellum
We isolated the mesial cerebellum (E14-P14), or hindbrain (E10-E12) from wild-type mice, and performed scRNA-seq on > 60,000 cells from five embryonal time points and four early postnatal time points (Figure- Conversely non-CNS cells were found across all time points. Several distinct clusters of cells only appear during restricted developmental time points, with many not found in the postnatal period (Extended Figure-4). We conclude that scRNA-seq is able to identify biologically distinct cerebellar cellular populations.
Based on known developmental relationships, and transcriptional similarity, we constructed pseudo-time trajectories for the various cerebellar lineages (Figure-

Cerebellar tumors mirror the transcriptomes of specific embryonic cerebellar cell clusters.
Using a carefully curated list of mouse/human orthologs and an algorithm to deconvolute complex RNA mixtures against a series of single cell type transcriptional profiles, we compared the transcriptomes of bulk human ependymomas (PFA and PFB, 43 tumors) as well as cerebellar pilocytic astrocytomas (C-PA, 10 tumors) to distinct developmental clusters. PFA, PFB, and C-PA are all transcriptionally most similar to the gliogenic progenitors-1 cell cluster (Extended Figure-5), with some similarity to the proliferating VZ progenitors, and to a novel cluster of 'roof plate-like' stem cells which has transcriptional similarity to the developing roof plate (Lmx1a, Msx1, Bmp7) 33,34 . The gliogenic progenitor cluster is observed initially at E12, peaks at E18, and has largely disappeared by P5 (Extended Figure-4). The 'roof plate-like' stem cells are seen much earlier, appearing only between E10 and E12 (Extended Figure-4). This is consistent with PFA, PFB, and C-PA being classified as 'gliomas', and prior publications that suggest that PFA arises from the regional radial glia 19 .
Posterior fossa ependymomas and C-PAs are both histologically and clinically distinct from each other. Re-clustering of the gliogenic progenitors, early VZ radial glia, and roof platestem cells identifies eight distinct transcriptional clusters within this lineage (Extended Figure-7). Both PFA and PFB remain transcriptionally best matched to the same developmental population (Extended Figure-7i). However, C-PA now transcriptionally match to a very distinct sub-cluster of the cerebellar gliogenic progenitor cell cluster, supporting a model in which posterior fossa ependymomas and C-PAs have distinct cells of origin (Extended Figure-7i).
Deconvoluting bulk RNA-seq transcriptomes from human medulloblastomas (Shh, Grp3-MB, and Grp4-MB, total 145 tumors) against transcriptionally defined cerebellar clusters (Figure-3) demonstrated that Shh-MB is most similar to GCP clusters, as supported by prior publications 6,7 . Re-clustering of cells in the GCP lineage revealed additional heterogeneity, with the identification of seven distinct clusters (Extended Figure-8). Comparison of bulk human Shh-MB transcriptomes to these seven GCP cell lineage clusters revealed heterogeneity within Shh-MB (Extended Figure-8c). Shh MBs that transcriptionally resemble later time points in cerebellar development (SHH1 -similar to the Post-natal GCPs 2.1) have a worse prognosis than those that resemble the earlier arising GCPs (SHH2similar to the Post-natal GCPs 1.1 (Extended Figure-8g) (p=0.00442). Furthermore, Shh-β subtype medulloblastomas are more similar to the earlier "SHH2" subset, demonstrating the similarity of this subtype to differentiation states in earlier GCP development (Extended Figure-8j).
The cell of origin for Grp3-MB is not yet definitively known, but has been suggested to be a Nestin +ve cerebellar stem cell 10 . Comparison of Grp3-MB to developmental cerebellar cell clusters revealed a broad resemblance across Grp3-MB to Nestin +ve cerebellar early stem cells (Figure-3). Interestingly, subsets of bulk Grp3-MB transcriptomes also resemble developmental cell clusters in the GCP and UBC lineages, as well as similarity to cerebellar GABAergic interneurons (Figure-3). This multi-lineage differentiation is consistent with a model in which Grp3-MBs contain Nestin +ve early neural stem-like cells that give rise to a variety of differentiated progeny; however this possibility cannot be evaluated properly using bulk RNA-seq data.
Mouse models of Grp4-MB are not available, and the cell of origin for Grp4-MB MB is unknown. Unexpectedly, Grp4-MB transcriptionally best matched to cells of the UBC lineage ( Figure-3). UBCs are glutamatergic interneurons that derive from the URL, which have not been studied as deeply as other types of cerebellar neurons 33 . The cluster that best matched Grp4-MB MB (cluster #8) (Extended Figure-1) is first observed at E14, peaks in prevalence at E18, and has largely disappeared by P0 (Extended Figure-4). Re-clustering of cells in the UBC lineage followed by comparison to the transcriptomes of bulk human Grp4-MB affirms that in all cases, the bulk Grp4-MB MB transcriptome is most similar to the UBC progenitor cell cluster (Extended Figure-9g). These data are consistent with a model in which Grp4-MB arises from a cell in the UBC lineage.

Temporal mirroring of specific embryonic cell clusters.
Different childhood tumors transcriptionally mirror specific clusters in defined cerebellar lineages. Many of these lineages are only detected over a defined period of development, while others persist into adulthood (Extended Figure-4). We compared human tumor transcriptomes to their best-matched developmental cluster as a function of time. Comparison of bulk PFA and PFB transcriptomes to murine single cells in the gliogenic progenitor cell lineage from E10 to P0 revealed a specific match to the E16 gliogenic progenitors ( Figure-4a,b). Similarly, C-PAs revealed a very strong match at both E16 and E18 (Figure- Comparison of Grp4-MBs that more closely mirror E16 to those that mirror E18 demonstrates that Grp4-β is largely restricted to E16 similar tumors, while Grp4-γ is completely restricted to E18 similar tumors (Extended Figure-9r, P=0.00004) 1 . These data suggest that differences between Grp4-β and Grp4-γ could be secondary to their arising at different time points in the UBC lineage 1 .
We did not attempt to temporally position Grp3-MBs using bulk transcriptomics, as Grp3-MB often transcriptionally match more than one cluster in the developing cerebellum. While PFA, PFB, C-PA, and Grp4-MB are all transcriptionally best matched to cell clusters present during fetal development, Shh-MBs are best matched to the GCP lineage in the early postnatal period at P5 (Figure-4g,h) 6 . We conclude that in addition to transcriptional mirroring of specific cell populations in the developing cerebellum, human cerebellar tumor transcriptomes are most similar to specific time points during development, predominantly during fetal life.

Single-cell heterogeneity in cerebellar tumors.
Medulloblastomas exhibit well-characterized intertumoral heterogeneity, as well as geographic, spatial (metastases) and temporal (at recurrence) heterogeneity 1,2,24,35-39 . Cerebellar ependymomas also show marked intertumoral heterogeneity [15][16][17][18][19]26 , while the intertumoral heterogeneity amongst C-PAs is not well characterized 25 . We undertook scRNA-seq of human cerebellar tumors including medulloblastoma (8 patients), PFA ependymoma (4 patients) and C-PA (3 patients). Analysis of clusters from scRNA-seq of human tumors revealed both tumor cell clusters, and non-tumor cell clusters (i.e., endothelium, monocytes/microglia, lymphocytes) (Extended Figure Figure-9h). Two distinct types of UBCs have been described in the mammalian cerebellum 40 . Consistent with this, re-clustering of murine cerebellar cells in the UBC lineage revealed two distinct types of UBCs (Extended Figure-9c). Grp4 MB is predominantly similar to only one of these subtypes (CR +ve UBCs -Extended Figure-9h). Individual cells from the UBC/GCP progenitor cluster simultaneously express both GCP and UBC marker genes, which is not observed in cells committed to a GCP or a UBC fate (Extended Figure-10), but is observed in Grp4-MB cells. These data are consistent with a model in which Grp4-MB arises from a bipotential progenitor cell population (likely in cluster #8) that is capable of giving rise to cells in both the GCP, and the UBC lineages. Cumulatively, across the medulloblastoma subgroups our scRNA-seq data demonstrates high levels of single cell heterogeneity, with evidence of multiple lines of differentiation, and cells at different time-points in the differentiation hierarchy.
Single cell RNA-seq of PFA ependymomas revealed single cell heterogeneity, with each tumor containing clusters that best matched sub-clusters of the gliogenic progenitor cell population or roof plate-like stem cells (Figure-6a, Extended Figure-7). Of note, we did not observe clusters of more differentiated cell types such as Bergmann glia or astrocytes within the PFA ependymomas, but only observed less differentiated cell types. This lack of differentiated cell types within the tumor is unique to PFA among the childhood cerebellar tumors examined in this study. Some single cell transcriptional clusters of human PFA tumors are more similar to single cell clusters from another patient's tumor than they are to other clusters from within the same individual (Extended Figure-7d). C-PA scRNA-seq also demonstrated clusters with transcriptional similarity to the relatively undifferentiated roof plate-like stem cell cluster and the gliogenic progenitors ( Figure-6b). Unlike PFA ependymoma, some C-PA tumor cell clusters demonstrate similarity to more differentiated cell types, such as astrocytes, Bergmann glia, and oligodendrocytes ( Figure-6b).

Discussion
Cumulatively, across the medulloblastoma subgroups (Shh, Grp3-MB, and Group 4), as well as PFA ependymoma and C-PA, our scRNA-seq data demonstrated high levels of single cell heterogeneity, with evidence of multiple lineages of differentiation, and tumor cells matching different time-points in the differentiation hierarchy. Childhood cerebellar tumor transcriptomes demonstrate high levels of similarity to discrete cell populations within the developing cerebellum, supporting a model in which discrete cells of origin have a profound influence on the transcriptome and biology of the observed mature tumors. While the use of cell cycle genes is essential in defining developmental stages of cell populations in the cerebellum (i.e., progenitor cell versus differentiated cell), we excluded cell cycle genes when comparing human tumor cells to murine developmental clusters to avoid spurious comparisons based largely on cell cycle phenotypes (Supplementary Table 3); however cell cycle content and genomic alterations were quantified independently in the human single cell data sets (Extended Figure-10, Supplementary Table 7). The relative biological significance of well separated clusters in the t-SNE plots, versus those that are spread out and partially overlapping is uncertain, and a well-known issue in single cell sequence analysis that will require future exploration.
Many of the normal murine cerebellar cell populations, which are transcriptionally most similar to human cerebellar tumors are only present for a restricted time period during fetal development, or are found only in the very early postnatal period. It is certainly possible that each of the cerebellar tumor types discussed above arises in a particular cell type indicated by the cluster and developmental time point identified by transcriptional matching. However extensive cell lineage tracing using in vivo models of these diseases will be necessary to exclude the possibility that tumors arise in other cell lineages and undergo trans-differentiation during transformation, or arise in cerebellar cells at later time points and then undergo dedifferentiation.
Our data does not reveal substantial populations of differentiated cells within PFA ependymomas, differentiated glia within the medulloblastomas, nor differentiated neuronal cell types in the C-PAs. Distinguishing infiltrating normal cells within pediatric cerebellar tumors is difficult as these tumor types contain very few somatic genetic events that would allow us to discern transformed from non-transformed cells. The normal, non-transformed cell populations that best transcriptionally match cerebellar tumors are only present in utero, or immediately postnatally, and are therefore not present in the brains of children at the time of presentation, and likely could not contaminate their tumors. However, we cannot exclude the possibility that some of our single cell tumor clusters also contain small populations of infiltrating non-transformed differentiated cells.
The ability of these different progenitor cell populations restricted to specific time points to give rise to different types of cerebellar tumors should ideally be tested functionally in vivo, as previously demonstrated 8,12,[18][19] . While there would be great value in comparing human cerebellar tumors to single cell transcriptomes from normal human cerebellar cells from various time points in development, these types of samples are not readily available. The presence of multiple lineages and stages of differentiation within bulk medulloblastoma samples illustrates the difficulty of using the bulk tumor population as a tool to decipher tumor biology, or to develop tumor diagnostics. The absence of differentiated cell types in PFA ependymoma, but not other tumor types, suggests the presence of a differentiation block in PFA ependymomas. Due to space limitations, the current manuscript has focused on the relationship of posterior fossa tumors to cell types in the developing cerebellum. The current single cell RNA-seq datasets from posterior fossa tumors should allow further analyses and insights into both tumor cell autonomous, and infiltrating non-tumor cell biology from these entities. A more complete understanding of the biology and transcriptomes of the specific cerebellar hierarchies identified above, and their developmental timing may allow a better comprehension of cerebellar tumor biology, and promote the subsequent development of novel mouse models, improved tumor diagnostics, and eventually the development of novel rational therapeutics based on the differences between tumor cells and their normal cells of origin.

Patient recruitment.
Participants recruitment was in compliance with Hospital for Sick Children and McGill University Health Centre ethical regulations. Sample size of participants was determined based on consent availability and diagnosis. Available patient characteristic information is described in Supplementary

Animal experiments.
All mouse breeding and procedures were approved by The Centre for Phenogenomics (Toronto). Mated C57BL/6J female mice were dissected in order to collect embryos from the following gestational time points: 10, 12, 14, 16 and 18. C57BL/6J pups were dissected to collect tissue from the following postnatal time points: 0, 5, 7 and 14.
Tissue handling and dissociation.
Fresh tumor tissue was collected at the time of resection. The tumor tissue was mechanically and enzymatically dissociated using a collagenase-based dissociation method as previously reported 41 . Early embryonic hindbrain structures were dissected from the gestational time points E10 and E12. An incision was made between the midbrain and hindbrain boundary, as well as between the prepontine hindbrain and pontine hindbrain, in order to isolate the isthmus, and rhombomeres −1 and −2 at these early time points in development. Late embryonic cerebellar primordia were collected at embryonic time points 14, 16 and 18. All embryonic mouse dissections were performed under a Leica stereoscope with a pair of Moria ultra fine forceps (Fine Science Tools). The tissue was transferred into ice cold Leibovitz's medium, followed by single cell dissociation with the Papain Dissociation System (Worthington Biochemical Corporation). Postnatal cerebella were dissected from the following time points: day 0, 5, 7 and 14. The central nervous system was fully dissected, then embedded in 2% Low melting point agarose. One mid sagittal slice of 300 um was generated using the Leica vibratome 42 . Under the stereoscope the cerebellum was isolated from the slice, followed by immediate single cell dissociation as described above.

RT, amplification and sequencing.
The concentration of the single cell suspension was assessed with a Trypan blue count. Approximately 10,000-14,000 cells per time point were loaded on the Chromium Controller and generated single cell GEMs. GEM-RT, DynaBeads cleanup, PCR amplification and SPRIselect beads cleanup were performed using Chromium Single Cell 3' Gel Bead kit. Indexed single cell libraries were generated using the Chromium Single Cell 3' Library kit and the Chromium i7 Multiplex kit. Size, quality, concentration and purity of the cDNAs and the corresponding 10× library was evaluated by the Agilent 2100 Bioanalyzer system. The 10× libraries were sequenced in the Illumina 2500 sequencing platform.

Alignment of raw reads.
Through 10× CellRanger's pipeline 43 , the raw base call (BCL) files were demultiplexed into FASTQ files. The FASTQ files were aligned to the reference mouse genome GRCm38 (mm10) to generate raw gene-barcode count matrices. When clustering of multiple samples, we aggregated the multiple runs together to normalize on sequencing depth, and recomputed the gene-barcode matrices. Alignment quality control metrics are found in Supplementary Table 1-tab 1.

QC and normalization.
Low-quality cells were identified and removed from the datasets. We considered low-quality cells as cells with <200-300 genes expressed, and cells with high mitochondrial gene content (4 S.D.s above median). We predicted doublets to be cells with relatively high library sizes (4-5 S.D.s above median), and removed them from the analysis. Lowabundance genes were also removed from the datasets (genes expressed in less than 3 cells). The human single cell datasets were normalized using methods adapted from Scran's pipeline 44 . Size factors were computed, and were applied to normalize gene expression across the cells, to produce normalized log-expression values to each dataset individually. Mouse single cell datasets were processed using CellRanger aggregation in order to account for sequencing depth variation, followed by QC and normalization as described above.

Clustering analysis and visualization.
Clustering analysis was performed with both the presence (Figure-1) and absence (Extended Figure-1) of cell cycle related genes obtained from Ensembl's biomart 49 . When implemented, cell cycle genes removal was done prior to highly variable genes (HGVs) detection. HGVs were detected using Seurat's pipeline 45,46 , calculating average expression and dispersion for each gene, diving genes into bins, and computing a z-score for dispersion within each bin. We used a z-score of 0.5 as the cutoff of dispersion, and a bottom cutoff of 0.0125 and a high cutoff of 3.0 for average expression. Linear dimensionality reduction was performed using principal component analysis (PCA), and statistically significant principal components were selected using the elbow and jackstraw methods from Seurat. The clusters of cells were identified by a shared nearest neighbor (SNN) modularity optimization based clustering algorithm from Seurat. We then visualized these clusters using t-SNE, tdistributed stochastic neighbor embedding.

Pseudo-time trajectory analysis.
The barcodes of selected clusters were normalized using Monocle's dPFeature to remove lowly expressed genes and perform PCA analysis on the remaining genes, for significant PC selection 47 . Cells are then grouped using 'density peak' clustering algorithm. Differential gene expression analysis was performed using a generalized linear model (GLM), and the top 1000 genes per cluster were selected. Reverse graph embedding (RGE) was then used to reduce the high-dimensionality data into lower dimensional space and build the trajectory 48 . The structure of the trajectory was plotted into two dimensional space using the DRTree dimensionality reduction algorithm and order the cells in pseudo-time.

Creation of cell-type-specific signatures.
For each cluster identified, the average expression of each gene was calculated. Differential gene expression was performed using Seurat's likelihood ratio test (LRT) method, and we filtered out genes expressed in less than 25% of the cells. The top differentially expressed genes were used as markers to build a signature gene expression matrix. Genes involved in cell proliferation and ribosome biogenesis were obtained from Ensembl's biomart 49 , and omitted from the matrix. Human orthologues of mouse genes were identified, and used to create the final matrix. Complete gene signatures and inputs can be found in Supplementary Table 5 and 6.

Deconvolution analysis.
CIBERSORT was used to perform the deconvolution analysis of the bulk and single-cell RNA-seq tumor data against the mouse clusters 50 . The full transcriptomes of the tumor data were used as the input mixture, and the signature gene input was the mouse cluster expression matrix after removal of all genes that can introduce bias in the deconvolution process including 1400 cell cycle genes, 300 genes associated with ribosome biogenesis, and ~100 mitochondrial and apoptosis-related genes. Quantile normalization was disabled and 100-500 permutations were run. To test CIBERSORT on our datasets, we created synthetic bulk mixtures from the mouse clusters, and we selected known amount of reads from various clusters. CIBERSORT roughly yielded the expected relative abundances. In order to generate reliable input expression profiles, tumor clusters with very low number of cells were discarded from the analysis. To validate our mouse cluster signatures, we obtained human and mouse data of brain cell types from published datasets (Saunders et Table 1-tab 2), and deconvoluted them against our mouse signatures in order to ensure that our expected abundances had similar values to our cell of origin matches.

Bulk RNA-seq human tumor samples.
60 Shh, 40 Group 3, and 45 Group 4 human medulloblastoma bulk RNA-seq samples were obtained from MAGIC, Medulloblastoma Advanced Genomics International Consortium. The raw data was aligned to the reference human genome GRChg38 using STAR to generate raw counts 51 . FPKMs were then obtained from DESeq2 52 . We performed bulk RNAsequencing on 22 PFA, 4 PFB, and 10 C-PA patient samples, and obtained FPKMs using the same strategy. Microarray expression data was obtained from 17 PFB samples.

Relative gene expression panels.
Relative expression of genes within selected clusters was measured by calculating the average of the gene count from the log normalized matrix in each cluster. The averages of specific genes were then scaled between the values of −1.0 -1.0 among the selected clusters in order to reflect higher versus lower expression levels.

Cell cycle analysis of human scRNA-seq tumor samples.
Cell cycle phase specific annotations were acquired from Whitfield et al. 2002 and used to quantify the percentage of phase specific cell cycle genes in each individual cell. The percentages was then normalized from a scale 0-1, after which the ratio between G1/S and G2/M was calculated. Cells with a low ratio values for the G2/M and G1/S ratios were labelled to be in G0 phase, whereas cells with increasing values for G1/S but low ratio values were labelled as cells progressing into G1 phase. Cell with high ratio values for G2/M were labelled to be progressing into S, G2 and M phase respectively.

Data availability statement.
The datasets generated and analysed during the current study are available in the following repositories: BAMs and filtered gene matrices of mouse developmental time points scRNAseq (GSE118068), FASTQs of PFB bulk RNA-seq and microarray expression  Survival curve, corrected for metastatic dissemination and molecular subtype, of SHH-1 (n=15) and SHH-2 (n=45) identified through matching to a re-clustered granule cell lineage (g). p value (p=0.00442) was determined by log rank test and + indicates censored cases. Comparison of additional clinical characteristics including histology (h), sex (i), molecular subtype affiliation (j), and metastatic status (k) of SHH-1 (n=15) and SHH-2 (n=45) patient samples. p values were determined using Fisher exact test.