Deriving a continuum score for group 3 and 4 medulloblastoma tumor samples analyzed via RNA-sequencing or DNA methylation microarray.

Here, we present a protocol for deriving a continuum score for group 3 and 4 medulloblastoma tumor samples analyzed via RNA-sequencing or DNA methylation microarray. We describe steps for utilizing NMF-defined group 3/group 4 metagenes to calculate a continuum score between 0 and 1 that can be projected onto new sample data analyzed via RNA-sequencing. We then detail procedures for reverse engineering a continuum score for samples analyzed via DNA methylation microarray using a random forest classifier.

Medulloblastoma group 3 and 4 tumors comprise a clinically and biologically significant expression continuum reflecting human cerebellar development.

Medulloblastoma is currently subclassified into distinct DNA methylation subgroups/subtypes with particular clinico-molecular features. Using RNA sequencing (RNA-seq) in large, well-annotated cohorts of medulloblastoma, we show that transcriptionally group 3 and group 4 medulloblastomas exist as intermediates on a bipolar continuum between archetypal group 3 and group 4 entities. Continuum position is prognostic, reflecting a propensity for specific DNA copy-number changes, and specific switches in isoform/enhancer usage and RNA editing. Examining single-cell RNA-seq (scRNA-seq) profiles, we show that intratumoral transcriptional heterogeneity along the continuum is limited in a subtype-dependent manner. By integrating with a human scRNA-seq reference atlas, we show that this continuum is mirrored by an equivalent continuum of transcriptional cell types in early fetal cerebellar development. We identify distinct developmental niches for all four major subgroups and link each to a common developmental antecedent. Our findings show a transcriptional continuum arising from oncogenic disruption of highly specific fetal cerebellar cell types, linked to almost every aspect of group 3/group 4 molecular biology and clinico-pathology.