Multi-scale systems genomics analysis predicts pathways, cell types, and drug targets involved in normative variation in peri-adolescent human cognition.

An open challenge in human genetics is to better understand the systems-level impact of genotype variation on developmental cognition. To characterize the genetic underpinnings of peri-adolescent cognition, we performed genotype-phenotype and systems analysis for binarized accuracy in nine cognitive tasks from the Philadelphia Neurodevelopmental Cohort (~2,200 individuals of European continental ancestry aged 8-21 years). We report a region of genome-wide significance within the 3' end of the Fibulin-1 gene (P = 4.6 × 10-8), associated with accuracy in nonverbal reasoning, a heritable form of complex reasoning ability. Diffusion tensor imaging data from a subset of these participants identified a significant association of white matter fractional anisotropy with FBLN1 genotypes (P < 0.025); poor performers show an increase in the C and A allele for rs77601382 and rs5765534, respectively, which is associated with increased fractional anisotropy. Integration of published human brain-specific 'omic maps, including single-cell transcriptomes of the developing human brain, shows that FBLN1 demonstrates greatest expression in the fetal brain, as a marker of intermediate progenitor cells, demonstrates negligible expression in the adolescent and adult human brain, and demonstrates increased expression in the brain in schizophrenia. Collectively these findings warrant further study of this gene and genetic locus in cognition, neurodevelopment, and disease. Separately, genotype-pathway analysis identified an enrichment of variants associated with working memory accuracy in pathways related to development and to autonomic nervous system dysfunction. Top-ranking pathway genes include those genetically associated with diseases with working memory deficits, such as schizophrenia and Parkinson's disease. This work advances the "molecules-to-behavior" view of cognition and provides a framework for using systems-level organization of data for other biomedical domains.

PRMT5 inhibition disrupts splicing and stemness in glioblastoma.

Glioblastoma (GBM) is a deadly cancer in which cancer stem cells (CSCs) sustain tumor growth and contribute to therapeutic resistance. Protein arginine methyltransferase 5 (PRMT5) has recently emerged as a promising target in GBM. Using two orthogonal-acting inhibitors of PRMT5 (GSK591 or LLY-283), we show that pharmacological inhibition of PRMT5 suppresses the growth of a cohort of 46 patient-derived GBM stem cell cultures, with the proneural subtype showing greater sensitivity. We show that PRMT5 inhibition causes widespread disruption of splicing across the transcriptome, particularly affecting cell cycle gene products. We identify a GBM splicing signature that correlates with the degree of response to PRMT5 inhibition. Importantly, we demonstrate that LLY-283 is brain-penetrant and significantly prolongs the survival of mice with orthotopic patient-derived xenografts. Collectively, our findings provide a rationale for the clinical development of brain penetrant PRMT5 inhibitors as treatment for GBM.

Gradient of Developmental and Injury Response transcriptional states defines functional vulnerabilities underpinning glioblastoma heterogeneity.

Glioblastomas harbor diverse cell populations, including rare glioblastoma stem cells (GSCs) that drive tumorigenesis. To characterize functional diversity within this population, we performed single-cell RNA sequencing on >69,000 GSCs cultured from the tumors of 26 patients. We observed a high degree of inter- and intra-GSC transcriptional heterogeneity that could not be fully explained by DNA somatic alterations. Instead, we found that GSCs mapped along a transcriptional gradient spanning two cellular states reminiscent of normal neural development and inflammatory wound response. Genome-wide CRISPR-Cas9 dropout screens independently recapitulated this observation, with each state characterized by unique essential genes. Further single-cell RNA sequencing of >56,000 malignant cells from primary tumors found that the majority organize along an orthogonal astrocyte maturation gradient yet retain expression of founder GSC transcriptional programs. We propose that glioblastomas grow out of a fundamental GSC-based neural wound response transcriptional program, which is a promising target for new therapy development.

Wnt activation as a therapeutic strategy in medulloblastoma.

Medulloblastoma (MB) is defined by four molecular subgroups (Wnt, Shh, Group 3, Group 4) with Wnt MB having the most favorable prognosis. Since prior reports have illustrated the antitumorigenic role of Wnt activation in Shh MB, we aimed to assess the effects of activated canonical Wnt signaling in Group 3 and 4 MBs. By using primary patient-derived MB brain tumor-initiating cell (BTIC) lines, we characterize differences in the tumor-initiating capacity of Wnt, Group 3, and Group 4 MB. With single cell RNA-seq technology, we demonstrate the presence of rare Wnt-active cells in non-Wnt MBs, which functionally retain the impaired tumorigenic potential of Wnt MB. In treating MB xenografts with a Wnt agonist, we provide a rational therapeutic option in which the protective effects of Wnt-driven MBs may be augmented in Group 3 and 4 MB and thereby support emerging data for a context-dependent tumor suppressive role for Wnt/ß-catenin signaling.

Medulloblastoma Arises from the Persistence of a Rare and Transient Sox2+ Granule Neuron Precursor.

Medulloblastoma (MB) is a neoplasm linked to dysregulated cerebellar development. Previously, we demonstrated that the Sonic Hedgehog (SHH) subgroup grows hierarchically, with Sox2+ cells at the apex of tumor progression and relapse. To test whether this mechanism is rooted in a normal developmental process, we studied the role of Sox2 in cerebellar development. We find that the external germinal layer (EGL) is derived from embryonic Sox2+ precursors and that the EGL maintains a rare fraction of Sox2+ cells during the first postnatal week. Through lineage tracing and single-cell analysis, we demonstrate that these Sox2+ cells are within the Atoh1+ lineage, contribute extensively to adult granule neurons, and resemble Sox2+ tumor cells. Critically, constitutive activation of the SHH pathway leads to their aberrant persistence in the EGL and rapid tumor onset. We propose that failure to eliminate this rare but potent developmental population is the tumor initiation mechanism in SHH-subgroup MB.

Oscillatory cortical forces promote three dimensional cell intercalations that shape the murine mandibular arch.

Multiple vertebrate embryonic structures such as organ primordia are composed of confluent cells. Although mechanisms that shape tissue sheets are increasingly understood, those which shape a volume of cells remain obscure. Here we show that 3D mesenchymal cell intercalations are essential to shape the mandibular arch of the mouse embryo. Using a genetically encoded vinculin tension sensor that we knock-in to the mouse genome, we show that cortical force oscillations promote these intercalations. Genetic loss- and gain-of-function approaches show that Wnt5a functions as a spatial cue to coordinate cell polarity and cytoskeletal oscillation. These processes diminish tissue rigidity and help cells to overcome the energy barrier to intercalation. YAP/TAZ and PIEZO1 serve as downstream effectors of Wnt5a-mediated actomyosin polarity and cytosolic calcium transients that orient and drive mesenchymal cell intercalations. These findings advance our understanding of how developmental pathways regulate biophysical properties and forces to shape a solid organ primordium.