Glioblastoma evolution and heterogeneity from a 3D whole-tumor perspective.

Treatment failure for the lethal brain tumor glioblastoma (GBM) is attributed to intratumoral heterogeneity and tumor evolution. We utilized 3D neuronavigation during surgical resection to acquire samples representing the whole tumor mapped by 3D spatial coordinates. Integrative tissue and single-cell analysis revealed sources of genomic, epigenomic, and microenvironmental intratumoral heterogeneity and their spatial patterning. By distinguishing tumor-wide molecular features from those with regional specificity, we inferred GBM evolutionary trajectories from neurodevelopmental lineage origins and initiating events such as chromothripsis to emergence of genetic subclones and spatially restricted activation of differential tumor and microenvironmental programs in the core, periphery, and contrast-enhancing regions. Our work depicts GBM evolution and heterogeneity from a 3D whole-tumor perspective, highlights potential therapeutic targets that might circumvent heterogeneity-related failures, and establishes an interactive platform enabling 360° visualization and analysis of 3D spatial patterns for user-selected genes, programs, and other features across whole GBM tumors.

Sequencing Diversity One Cell at a Time.

Single-cell RNA sequencing provides a new approach to an old problem: how to study cellular diversity in complex biological systems. Three studies-Saunders et al., Zeisel et al., and Davie et al.-deploy this technique on an unprecedented scale to reveal transcriptional patterns that distinguish cells in the nervous systems of mice and flies.

Progranulin Deficiency Promotes Circuit-Specific Synaptic Pruning by Microglia via Complement Activation.

Microglia maintain homeostasis in the brain, but whether aberrant microglial activation can cause neurodegeneration remains controversial. Here, we use transcriptome profiling to demonstrate that deficiency in frontotemporal dementia (FTD) gene progranulin (Grn) leads to an age-dependent, progressive upregulation of lysosomal and innate immunity genes, increased complement production, and enhanced synaptic pruning in microglia. During aging, Grn(-/-) mice show profound microglia infiltration and preferential elimination of inhibitory synapses in the ventral thalamus, which lead to hyperexcitability in the thalamocortical circuits and obsessive-compulsive disorder (OCD)-like grooming behaviors. Remarkably, deleting C1qa gene significantly reduces synaptic pruning by Grn(-/-) microglia and mitigates neurodegeneration, behavioral phenotypes, and premature mortality in Grn(-/-) mice. Together, our results uncover a previously unrecognized role of progranulin in suppressing aberrant microglia activation during aging. These results represent an important conceptual advance that complement activation and microglia-mediated synaptic pruning are major drivers, rather than consequences, of neurodegeneration caused by progranulin deficiency.

Molecular identity of human outer radial glia during cortical development.

Radial glia, the neural stem cells of the neocortex, are located in two niches: the ventricular zone and outer subventricular zone. Although outer subventricular zone radial glia may generate the majority of human cortical neurons, their molecular features remain elusive. By analyzing gene expression across single cells, we find that outer radial glia preferentially express genes related to extracellular matrix formation, migration, and stemness, including TNC, PTPRZ1, FAM107A, HOPX, and LIFR. Using dynamic imaging, immunostaining, and clonal analysis, we relate these molecular features to distinctive behaviors of outer radial glia, demonstrate the necessity of STAT3 signaling for their cell cycle progression, and establish their extensive proliferative potential. These results suggest that outer radial glia directly support the subventricular niche through local production of growth factors, potentiation of growth factor signals by extracellular matrix proteins, and activation of self-renewal pathways, thereby enabling the developmental and evolutionary expansion of the human neocortex.

Generation of functional human oligodendrocytes from dermal fibroblasts by direct lineage conversion.

Oligodendrocytes, the myelinating cells of the central nervous system, possess great potential for disease modeling and cell transplantation-based therapies for leukodystrophies. However, caveats to oligodendrocyte differentiation protocols ( Ehrlich et al., 2017; Wang et al., 2013; Douvaras and Fossati, 2015) from human embryonic stem and induced pluripotent stem cells (iPSCs), which include slow and inefficient differentiation, and tumorigenic potential of contaminating undifferentiated pluripotent cells, are major bottlenecks towards their translational utility. Here, we report the rapid generation of human oligodendrocytes by direct lineage conversion of human dermal fibroblasts (HDFs). We show that the combination of the four transcription factors OLIG2, SOX10, ASCL1 and NKX2.2 is sufficient to convert HDFs to induced oligodendrocyte precursor cells (iOPCs). iOPCs resemble human primary and iPSC-derived OPCs based on morphology and transcriptomic analysis. Importantly, iOPCs can differentiate into mature myelinating oligodendrocytes in vitro and in vivo. Finally, iOPCs derived from patients with Pelizaeus Merzbacher disease, a hypomyelinating leukodystrophy caused by mutations in the proteolipid protein 1 (PLP1) gene, showed increased cell death compared with iOPCs from healthy donors. Thus, human iOPCs generated by direct lineage conversion represent an attractive new source for human cell-based disease models and potentially myelinating cell grafts.

Human hippocampal neurogenesis drops sharply in children to undetectable levels in adults.

New neurons continue to be generated in the subgranular zone of the dentate gyrus of the adult mammalian hippocampus. This process has been linked to learning and memory, stress and exercise, and is thought to be altered in neurological disease. In humans, some studies have suggested that hundreds of new neurons are added to the adult dentate gyrus every day, whereas other studies find many fewer putative new neurons. Despite these discrepancies, it is generally believed that the adult human hippocampus continues to generate new neurons. Here we show that a defined population of progenitor cells does not coalesce in the subgranular zone during human fetal or postnatal development. We also find that the number of proliferating progenitors and young neurons in the dentate gyrus declines sharply during the first year of life and only a few isolated young neurons are observed by 7 and 13 years of age. In adult patients with epilepsy and healthy adults (18-77 years; n = 17 post-mortem samples from controls; n = 12 surgical resection samples from patients with epilepsy), young neurons were not detected in the dentate gyrus. In the monkey (Macaca mulatta) hippocampus, proliferation of neurons in the subgranular zone was found in early postnatal life, but this diminished during juvenile development as neurogenesis decreased. We conclude that recruitment of young neurons to the primate hippocampus decreases rapidly during the first years of life, and that neurogenesis in the dentate gyrus does not continue, or is extremely rare, in adult humans. The early decline in hippocampal neurogenesis raises questions about how the function of the dentate gyrus differs between humans and other species in which adult hippocampal neurogenesis is preserved.

Positive Controls in Adults and Children Support That Very Few, If Any, New Neurons Are Born in the Adult Human Hippocampus.

Adult hippocampal neurogenesis was originally discovered in rodents. Subsequent studies identified the adult neural stem cells and found important links between adult neurogenesis and plasticity, behavior, and disease. However, whether new neurons are produced in the human dentate gyrus (DG) during healthy aging is still debated. We and others readily observe proliferating neural progenitors in the infant hippocampus near immature cells expressing doublecortin (DCX), but the number of such cells decreases in children and few, if any, are present in adults. Recent investigations using dual antigen retrieval find many cells stained by DCX antibodies in adult human DG. This has been interpreted as evidence for high rates of adult neurogenesis, even at older ages. However, most of these DCX-labeled cells have mature morphology. Furthermore, studies in the adult human DG have not found a germinal region containing dividing progenitor cells. In this Dual Perspectives article, we show that dual antigen retrieval is not required for the detection of DCX in multiple human brain regions of infants or adults. We review prior studies and present new data showing that DCX is not uniquely expressed by newly born neurons: DCX is present in adult amygdala, entorhinal and parahippocampal cortex neurons despite being absent in the neighboring DG. Analysis of available RNA-sequencing datasets supports the view that DG neurogenesis is rare or absent in the adult human brain. To resolve the conflicting interpretations in humans, it is necessary to identify and visualize dividing neuronal precursors or develop new methods to evaluate the age of a neuron at the single-cell level.

Human-specific regulation of MeCP2 levels in fetal brains by microRNA miR-483-5p.

Proper neurological function in humans requires precise control of levels of the epigenetic regulator methyl CpG-binding protein 2 (MeCP2). MeCP2 protein levels are low in fetal brains, where the predominant MECP2 transcripts have an unusually long 3' untranslated region (UTR). Here, we show that miR-483-5p, an intragenic microRNA of the imprinted IGF2, regulates MeCP2 levels through a human-specific binding site in the MECP2 long 3' UTR. We demonstrate the inverse correlation of miR-483-5p and MeCP2 levels in developing human brains and fibroblasts from Beckwith-Wiedemann syndrome patients. Importantly, expression of miR-483-5p rescues abnormal dendritic spine phenotype of neurons overexpressing human MeCP2. In addition, miR-483-5p modulates the levels of proteins of the MeCP2-interacting corepressor complexes, including HDAC4 and TBL1X. These data provide insight into the role of miR-483-5p in regulating the levels of MeCP2 and interacting proteins during human fetal development.

Radial glia require PDGFD-PDGFRß signalling in human but not mouse neocortex.

Evolutionary expansion of the human neocortex underlies many of our unique mental abilities. This expansion has been attributed to the increased proliferative potential of radial glia (RG; neural stem cells) and their subventricular dispersion from the periventricular niche during neocortical development. Such adaptations may have evolved through gene expression changes in RG. However, whether or how RG gene expression varies between humans and other species is unknown. Here we show that the transcriptional profiles of human and mouse neocortical RG are broadly conserved during neurogenesis, yet diverge for specific signalling pathways. By analysing differential gene co-expression relationships between the species, we demonstrate that the growth factor PDGFD is specifically expressed by RG in human, but not mouse, corticogenesis. We also show that the expression domain of PDGFRß, the cognate receptor for PDGFD, is evolutionarily divergent, with high expression in the germinal region of dorsal human neocortex but not in the mouse. Pharmacological inhibition of PDGFD-PDGFRß signalling in slice culture prevents normal cell cycle progression of neocortical RG in human, but not mouse. Conversely, injection of recombinant PDGFD or ectopic expression of constitutively active PDGFRß in developing mouse neocortex increases the proportion of RG and their subventricular dispersion. These findings highlight the requirement of PDGFD-PDGFRß signalling for human neocortical development and suggest that local production of growth factors by RG supports the expanded germinal region and progenitor heterogeneity of species with large brains.

Correction: Pleiotropic Mechanisms Indicated for Sex Differences in Autism.

[This corrects the article DOI: 10.1371/journal.pgen.1006425.].

Secretagogin is Expressed by Developing Neocortical GABAergic Neurons in Humans but not Mice and Increases Neurite Arbor Size and Complexity.

The neocortex of primates, including humans, contains more abundant and diverse inhibitory neurons compared with rodents, but the molecular foundations of these observations are unknown. Through integrative gene coexpression analysis, we determined a consensus transcriptional profile of GABAergic neurons in mid-gestation human neocortex. By comparing this profile to genes expressed in GABAergic neurons purified from neonatal mouse neocortex, we identified conserved and distinct aspects of gene expression in these cells between the species. We show here that the calcium-binding protein secretagogin (SCGN) is robustly expressed by neocortical GABAergic neurons derived from caudal ganglionic eminences (CGE) and lateral ganglionic eminences during human but not mouse brain development. Through electrophysiological and morphometric analyses, we examined the effects of SCGN expression on GABAergic neuron function and form. Forced expression of SCGN in CGE-derived mouse GABAergic neurons significantly increased total neurite length and arbor complexity following transplantation into mouse neocortex, revealing a molecular pathway that contributes to morphological differences in these cells between rodents and primates.

Pleiotropic Mechanisms Indicated for Sex Differences in Autism.

Sexual dimorphism in common disease is pervasive, including a dramatic male preponderance in autism spectrum disorders (ASDs). Potential genetic explanations include a liability threshold model requiring increased polymorphism risk in females, sex-limited X-chromosome contribution, gene-environment interaction driven by differences in hormonal milieu, risk influenced by genes sex-differentially expressed in early brain development, or contribution from general mechanisms of sexual dimorphism shared with secondary sex characteristics. Utilizing a large single nucleotide polymorphism (SNP) dataset, we identify distinct sex-specific genome-wide significant loci. We investigate genetic hypotheses and find no evidence for increased genetic risk load in females, but evidence for sex heterogeneity on the X chromosome, and contribution of sex-heterogeneous SNPs for anthropometric traits to ASD risk. Thus, our results support pleiotropy between secondary sex characteristic determination and ASDs, providing a biological basis for sex differences in ASDs and implicating non brain-limited mechanisms.

Neurons show distinctive DNA methylation profile and higher interindividual variations compared with non-neurons.

Epigenome information in mammalian brain cells reflects their developmental history, neuronal activity, and environmental exposures. Studying the epigenetic modifications present in neuronal cells is critical to a more complete understanding of the role of the genome in brain functions. We performed comprehensive DNA methylation analysis in neuronal and non-neuronal nuclei obtained from the human prefrontal cortex. Neuronal nuclei manifest qualitatively and quantitatively distinctive DNA methylation patterns, including relative global hypomethylation, differential enrichment of transcription-factor binding sites, and higher methylation of genes expressed in astrocytes. Non-neuronal nuclei showed indistinguishable DNA methylation patterns from bulk cortex and higher methylation of synaptic transmission-related genes compared with neuronal nuclei. We also found higher variation in DNA methylation in neuronal nuclei, suggesting that neuronal cells have more potential ability to change their epigenetic status in response to developmental and environmental conditions compared with non-neuronal cells in the central nervous system.

Deconstructing intratumoral heterogeneity through multiomic and multiscale analysis of serial sections.

Tumors may contain billions of cells including distinct malignant clones and nonmalignant cell types. Clarifying the evolutionary histories, prevalence, and defining molecular features of these cells is essential for improving clinical outcomes, since intratumoral heterogeneity provides fuel for acquired resistance to targeted therapies. Here we present a statistically motivated strategy for deconstructing intratumoral heterogeneity through multiomic and multiscale analysis of serial tumor sections (MOMA). By combining deep sampling of IDH-mutant astrocytomas with integrative analysis of single-nucleotide variants, copy-number variants, and gene expression, we reconstruct and validate the phylogenies, spatial distributions, and transcriptional profiles of distinct malignant clones. By genotyping nuclei analyzed by single-nucleus RNA-seq for truncal mutations, we further show that commonly used algorithms for identifying cancer cells from single-cell transcriptomes may be inaccurate. We also demonstrate that correlating gene expression with tumor purity in bulk samples can reveal optimal markers of malignant cells and use this approach to identify a core set of genes that is consistently expressed by astrocytoma truncal clones, including AKR1C3, whose expression is associated with poor outcomes in several types of cancer. In summary, MOMA provides a robust and flexible strategy for precisely deconstructing intratumoral heterogeneity and clarifying the core molecular properties of distinct cellular populations in solid tumors.

Deconstructing Intratumoral Heterogeneity through Multiomic and Multiscale Analysis of Serial Sections.

Tumors may contain billions of cells, including distinct malignant clones and nonmalignant cell types. Clarifying the evolutionary histories, prevalence, and defining molecular features of these cells is essential for improving clinical outcomes, since intratumoral heterogeneity provides fuel for acquired resistance to targeted therapies. Here we present a statistically motivated strategy for deconstructing intratumoral heterogeneity through multiomic and multiscale analysis of serial tumor sections (MOMA). By combining deep sampling of IDH-mutant astrocytomas with integrative analysis of single-nucleotide variants, copy-number variants, and gene expression, we reconstruct and validate the phylogenies, spatial distributions, and transcriptional profiles of distinct malignant clones. By genotyping nuclei analyzed by single-nucleus RNA-seq for truncal mutations, we further show that commonly used algorithms for identifying cancer cells from single-cell transcriptomes may be inaccurate. We also demonstrate that correlating gene expression with tumor purity in bulk samples can reveal optimal markers of malignant cells and use this approach to identify a core set of genes that are consistently expressed by astrocytoma truncal clones, including AKR1C3, whose expression is associated with poor outcomes in several types of cancer. In summary, MOMA provides a robust and flexible strategy for precisely deconstructing intratumoral heterogeneity and clarifying the core molecular properties of distinct cellular populations in solid tumors.

Expression profiling of Aldh1l1-precursors in the developing spinal cord reveals glial lineage-specific genes and direct Sox9-Nfe2l1 interactions.

Developmental regulation of gliogenesis in the mammalian CNS is incompletely understood, in part due to a limited repertoire of lineage-specific genes. We used Aldh1l1-GFP as a marker for gliogenic radial glia and later-stage precursors of developing astrocytes and performed gene expression profiling of these cells. We then used this dataset to identify candidate transcription factors that may serve as glial markers or regulators of glial fate. Our analysis generated a database of developmental stage-related markers of Aldh1l1+ cells between murine embryonic day 13.5-18.5. Using these data we identify the bZIP transcription factor Nfe2l1 and demonstrate that it promotes glial fate under direct Sox9 regulatory control. Thus, this dataset represents a resource for identifying novel regulators of glial development.

Bulk and mosaic deletions of Egfr reveal regionally defined gliogenesis in the developing mouse forebrain.

The epidermal growth factor receptor (EGFR) plays a role in cell proliferation and differentiation during healthy development and tumor growth; however, its requirement for brain development remains unclear. Here we used a conditional mouse allele for Egfr to examine its contributions to perinatal forebrain development at the tissue level. Subtractive bulk ventral and dorsal forebrain deletions of Egfr uncovered significant and permanent decreases in oligodendrogenesis and myelination in the cortex and corpus callosum. Additionally, an increase in astrogenesis or reactive astrocytes in effected regions was evident in response to cortical scarring. Sparse deletion using mosaic analysis with double markers (MADM) surprisingly revealed a regional requirement for EGFR in rostrodorsal, but not ventrocaudal glial lineages including both astrocytes and oligodendrocytes. The EGFR-independent ventral glial progenitors may compensate for the missing EGFR-dependent dorsal glia in the bulk Egfr-deleted forebrain, potentially exposing a regenerative population of gliogenic progenitors in the mouse forebrain.

Is my network module preserved and reproducible?

In many applications, one is interested in determining which of the properties of a network module change across conditions. For example, to validate the existence of a module, it is desirable to show that it is reproducible (or preserved) in an independent test network. Here we study several types of network preservation statistics that do not require a module assignment in the test network. We distinguish network preservation statistics by the type of the underlying network. Some preservation statistics are defined for a general network (defined by an adjacency matrix) while others are only defined for a correlation network (constructed on the basis of pairwise correlations between numeric variables). Our applications show that the correlation structure facilitates the definition of particularly powerful module preservation statistics. We illustrate that evaluating module preservation is in general different from evaluating cluster preservation. We find that it is advantageous to aggregate multiple preservation statistics into summary preservation statistics. We illustrate the use of these methods in six gene co-expression network applications including 1) preservation of cholesterol biosynthesis pathway in mouse tissues, 2) comparison of human and chimpanzee brain networks, 3) preservation of selected KEGG pathways between human and chimpanzee brain networks, 4) sex differences in human cortical networks, 5) sex differences in mouse liver networks. While we find no evidence for sex specific modules in human cortical networks, we find that several human cortical modules are less preserved in chimpanzees. In particular, apoptosis genes are differentially co-expressed between humans and chimpanzees. Our simulation studies and applications show that module preservation statistics are useful for studying differences between the modular structure of networks. Data, R software and accompanying tutorials can be downloaded from the following webpage: http://www.genetics.ucla.edu/labs/horvath/CoexpressionNetwork/ModulePreservation.

Nests of dividing neuroblasts sustain interneuron production for the developing human brain.

The human cortex contains inhibitory interneurons derived from the medial ganglionic eminence (MGE), a germinal zone in the embryonic ventral forebrain. How this germinal zone generates sufficient interneurons for the human brain remains unclear. We found that the human MGE (hMGE) contains nests of proliferative neuroblasts with ultrastructural and transcriptomic features that distinguish them from other progenitors in the hMGE. When dissociated hMGE cells are transplanted into the neonatal mouse brain, they reform into nests containing proliferating neuroblasts that generate young neurons that migrate extensively into the mouse forebrain and mature into different subtypes of functional interneurons. Together, these results indicate that the nest organization and sustained proliferation of neuroblasts in the hMGE provide a mechanism for the extended production of interneurons for the human forebrain.

Diagnostic blood RNA profiles for human acute spinal cord injury.

Diagnosis of spinal cord injury (SCI) severity at the ultra-acute stage is of great importance for emergency clinical care of patients as well as for potential enrollment into clinical trials. The lack of a diagnostic biomarker for SCI has played a major role in the poor results of clinical trials. We analyzed global gene expression in peripheral white blood cells during the acute injury phase and identified 197 genes whose expression changed after SCI compared with healthy and trauma controls and in direct relation to SCI severity. Unsupervised coexpression network analysis identified several gene modules that predicted injury severity (AIS grades) with an overall accuracy of 72.7% and included signatures of immune cell subtypes. Specifically, for complete SCIs (AIS A), ROC analysis showed impressive specificity and sensitivity (AUC: 0.865). Similar precision was also shown for AIS D SCIs (AUC: 0.938). Our findings indicate that global transcriptomic changes in peripheral blood cells have diagnostic and potentially prognostic value for SCI severity.