Targeting oncogenic pathways holds promise for brain tumor treatment, but inhibition of Sonic Hedgehog (SHH) signaling has failed in SHH-driven medulloblastoma. Cellular diversity within tumors and reduced lineage commitment can undermine targeted therapy by increasing the probability of treatment-resistant populations. Using single-cell RNA-seq and lineage tracing, we analyzed cellular diversity in medulloblastomas in transgenic, medulloblastoma-prone mice, and responses to the SHH-pathway inhibitor vismodegib. In untreated tumors, we find expected stromal cells and tumor-derived cells showing either a spectrum of neural progenitor-differentiation states or glial and stem cell markers. Vismodegib reduces the proliferative population and increases differentiation. However, specific cell types in vismodegib-treated tumors remain proliferative, showing either persistent SHH-pathway activation or stem cell characteristics. Our data show that even in tumors with a single pathway-activating mutation, diverse mechanisms drive tumor growth. This diversity confers early resistance to targeted inhibitor therapy, demonstrating the need to target multiple pathways simultaneously.

Background Medulloblastoma (MB) is a heterogeneous disease in which neoplastic cells and associated immune cells contribute to disease progression. We aimed to determine the influence of neoplastic and immune cell diversity on MB biology in patient samples and animal models. Methods To better characterize cellular heterogeneity in MB we used single-cell RNA sequencing, immunohistochemistry and deconvolution of transcriptomic data to profile neoplastic and immune populations in patient samples and animal models across childhood MB subgroups. Results Neoplastic cells cluster primarily according to individual sample of origin which is influenced by chromosomal copy number variance. Harmony alignment reveals novel MB subgroup/subtype-associated subpopulations that recapitulate neurodevelopmental processes, including photoreceptor and glutamatergic neuron-like cells in molecular subgroups GP3 and GP4, and a specific nodule-associated neuronally-differentiated subpopulation in subgroup molecular SHH. We definitively chart the spectrum of MB immune cell infiltrates, which include subpopulations that recapitulate developmentally-related neuron-pruning and antigen presenting myeloid cells. MB cellular diversity matching human samples is mirrored in subgroup-specific mouse models of MB. Conclusions These findings provide a clearer understanding of the diverse neoplastic and immune cell subpopulations that constitute the MB microenvironment.

Noradrenergic signaling plays a well-established role in promoting the stress response. Here we identify a subpopulation of noradrenergic neurons, defined by developmental expression of Hoxb1, that has a unique role in modulating stress-related behavior. Using an intersectional chemogenetic strategy, in combination with behavioral and physiological analyses, we show that activation of Hoxb1-noradrenergic (Hoxb1-NE) neurons decreases anxiety-like behavior and promotes an active coping strategy in response to acute stressors. In addition, we use cerebral blood volume-weighted functional magnetic resonance imaging to show that chemoactivation of Hoxb1-NE neurons results in reduced activity in stress-related brain regions, including the bed nucleus of the stria terminalis, amygdala, and locus coeruleus. Thus, the actions of Hoxb1-NE neurons are distinct from the well-documented functions of the locus coeruleus in promoting the stress response, demonstrating that the noradrenergic system contains multiple functionally distinct subpopulations.

Previous genetic and public health research in the Pakistani population has focused on the role of consanguinity in increasing recessive disease risk, but little is known about its recent population history or the effects of endogamy. Here, we investigate fine-scale population structure, history and consanguinity patterns using genotype chip data from 2,200 British Pakistanis. We reveal strong recent population structure driven by the biraderi social stratification system. We find that all subgroups have had low recent effective population sizes (Ne), with some showing a decrease 15‒20 generations ago that has resulted in extensive identity-by-descent sharing and homozygosity, increasing the risk of recessive disorders. Our results from two orthogonal methods (one using machine learning and the other coalescent-based) suggest that the detailed reporting of parental relatedness for mothers in the cohort under-represents the true levels of consanguinity. These results demonstrate the impact of cultural practices on population structure and genomic diversity in Pakistanis, and have important implications for medical genetic studies.

Medulloblastoma (MB) is a heterogeneous disease in which neoplastic cells and associated immune cells contribute to disease progression. To better understand cellular heterogeneity in MB we profile neoplastic and immune populations in childhood MB samples using single-cell RNA sequencing, immunohistochemistry and deconvolution of transcriptomic data. Neoplastic cells cluster primarily according to individual sample of origin which is in part due to the effect of chromosomal copy number gains and losses. Harmony alignment reveals novel MB subgroup/subtype-associated subpopulations that recapitulate neurodevelopmental processes and are associated with clinical outcomes. We identify discrete photoreceptor-like cells in MB subgroups GP3 and GP4 and nodule-associated neuronally-differentiated cells in subgroup SHH. MB immune infiltrates consist of both developmentally-related neuron-pruning and antigen presenting myeloid cells. We show that this MB cellular diversity is recapitulated in genetically engineered mouse subgroup-specific models of MB. These findings advance our understanding of both the neoplastic and immune landscape of MB.

It is unclear why medulloblastoma patients receiving similar treatments experience different outcomes. Transcriptomic profiling identified subgroups with different prognoses, but in each subgroup, individuals remain at risk of incurable recurrence. To investigate why similar-appearing tumors produce variable outcomes, we analyzed medulloblastomas triggered in transgenic mice by a common driver mutation expressed at different points in brain development. We genetically engineered mice to express oncogenic SmoM2, starting in multipotent glio-neuronal stem cells, or committed neural progenitors. Both groups developed medulloblastomas with similar transcriptomic profiles. We compared medulloblastoma progression, radiosensitivity, and cellular heterogeneity, determined by single-cell transcriptomic analysis (scRNA-seq). Stem cell-triggered medulloblastomas progressed faster, contained more OLIG2-expressing stem-like cells, and consistently showed radioresistance. In contrast, progenitor-triggered MBs progressed slower, down-regulated stem-like cells and were curable with radiation. Progenitor-triggered medulloblastomas also contained more diverse stromal populations, with more Ccr2+ macrophages and fewer Igf1+ microglia, indicating that developmental events affected the subsequent tumor microenvironment. Reduced mTORC1 activity in M-Smo tumors suggests that differential Igf1 contributed to differences in phenotype. Developmental events in tumorigenesis that were obscure in transcriptomic profiles thus remained cryptic determinants of tumor composition and outcome. Precise understanding of medulloblastoma pathogenesis and prognosis requires supplementing transcriptomic/methylomic studies with analyses that resolve cellular heterogeneity.

Patients with medulloblastoma are typically treated with a narrow range of therapies, but may experience widely divergent outcomes; 80-90% become long-term survivors while 20% develop incurable recurrence. Transcriptomic profiling has identified four subgroups with different recurrence risks, but outcomes remain variable for individual patients within each subgroup. To gain new insight into why patients with similar-appearing tumors have variable outcomes, we examined how the timing of tumor initiation effects medulloblastomas triggered by a single, common driver mutation. We genetically-engineered mice to express an oncogenic Smo allele starting early in development in the broad lineage of neural stem cells, or later, in the more committed lineage of cerebellar granule neuron progenitors. Both groups developed medulloblastomas and no other tumors. We compared medulloblastoma progression, response to therapy, gene expression profile and cellular heterogeneity, determined by single cell transcriptomic analysis (scRNA-seq). The average transcriptomic profiles of the tumors were similar. However, stem cell-triggered medulloblastomas progressed faster, contained more OLIG2-expressing tumor stem cells, and consistently showed radioresistance. In contrast, progenitor-triggered MBs progressed slower, lost stem cell character over time and were radiosensitive. Progenitor-triggered medulloblastomas also contained more diverse stromal populations, including tumor-associated macrophages, indicating that the timing of oncogenesis affected the subsequent interactions between the tumor and microenvironment. Our findings show that developmental events in tumorigenesis may be impossible to infer from transcriptomic profile, but while remaining cryptic can nevertheless influence tumor composition and the outcome of therapy. Precise understanding of medulloblastoma pathogenesis and prognosis requires supplementing transcriptomic data with biomarkers of cellular heterogeneity.

The therapeutic potential of CDK4/6 inhibitors for brain tumors has been limited by recurrence. To address recurrence, we tested a nanoparticle formulation of CDK4/6 inhibitor palbociclib (POx-Palbo) in mice genetically-engineered to develop SHH-driven medulloblastoma, alone or in combination with specific agents suggested by our analysis. Nanoparticle encapsulation reduced palbociclib toxicity, enabled parenteral administration, improved CNS pharmacokinetics, and extended mouse survival, but recurrence persisted. scRNA-seq identified up-regulation of glutamate transporter Slc1a2 and down-regulation of diverse ribosomal genes in proliferating medulloblastoma cells in POx-Palbo–treated mice, suggesting mTORC1 signaling suppression, subsequently confirmed by decreased 4EBP1 phosphorylation. Combining POx-Palbo with the mTORC1 inhibitor sapanisertib produced mutually enhancing effects and prolonged mouse survival compared to either agent alone, contrasting markedly with other tested drug combinations. Our data show the potential of nanoparticle formulation and scRNA-seq analysis of resistance to improve brain tumor treatment and identify POx-Palbo + Sapanisertib as effective combinatorial therapy for SHH medulloblastoma.