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.

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.

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.

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.

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.

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.

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.

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.

MicroRNA Ratios Distinguish Melanomas from Nevi.

The use of microRNAs as biomarkers has been proposed for many diseases, including the diagnosis of melanoma. Although hundreds of microRNAs have been identified as differentially expressed in melanomas as compared to benign melanocytic lesions, a limited consensus has been achieved across studies, constraining the effective use of these potentially useful markers. In this study, we applied a machine learning-based pipeline to a dataset consisting of genetic features, clinical features, and next-generation microRNA sequencing from micro-dissected formalin-fixed paraffin embedded melanomas and their adjacent benign precursor nevi. We identified patient age and tumor cellularity as variables that frequently confound the measured expression of potentially diagnostic microRNAs. By employing the ratios of microRNAs that were either enriched or depleted in melanoma compared to the nevi as a normalization strategy, we developed a model that classified all the available published cohorts with an area under the receiver operating characteristic curve of 0.98. External validation on an independent cohort classified lesions with 81% sensitivity and 88% specificity and was uninfluenced by the tumor content of the sample or patient age.

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.

Variation among intact tissue samples reveals the core transcriptional features of human CNS cell classes.

It is widely assumed that cells must be physically isolated to study their molecular profiles. However, intact tissue samples naturally exhibit variation in cellular composition, which drives covariation of cell-class-specific molecular features. By analyzing transcriptional covariation in 7,221 intact CNS samples from 840 neurotypical individuals, representing billions of cells, we reveal the core transcriptional identities of major CNS cell classes in humans. By modeling intact CNS transcriptomes as a function of variation in cellular composition, we identify cell-class-specific transcriptional differences in Alzheimer's disease, among brain regions, and between species. Among these, we show that PMP2 is expressed by human but not mouse astrocytes and significantly increases mouse astrocyte size upon ectopic expression in vivo, causing them to more closely resemble their human counterparts. Our work is available as an online resource ( http://oldhamlab.ctec.ucsf.edu/ ) and provides a generalizable strategy for determining the core molecular features of cellular identity in intact biological systems.

A Glial Signature and Wnt7 Signaling Regulate Glioma-Vascular Interactions and Tumor Microenvironment.

Gliomas comprise heterogeneous malignant glial and stromal cells. While blood vessel co-option is a potential mechanism to escape anti-angiogenic therapy, the relevance of glial phenotype in this process is unclear. We show that Olig2+ oligodendrocyte precursor-like glioma cells invade by single-cell vessel co-option and preserve the blood-brain barrier (BBB). Conversely, Olig2-negative glioma cells form dense perivascular collections and promote angiogenesis and BBB breakdown, leading to innate immune cell activation. Experimentally, Olig2 promotes Wnt7b expression, a finding that correlates in human glioma profiling. Targeted Wnt7a/7b deletion or pharmacologic Wnt inhibition blocks Olig2+ glioma single-cell vessel co-option and enhances responses to temozolomide. Finally, Olig2 and Wnt7 become upregulated after anti-VEGF treatment in preclinical models and patients. Thus, glial-encoded pathways regulate distinct glioma-vascular microenvironmental interactions.

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.

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.

Correction: Pleiotropic Mechanisms Indicated for Sex Differences in Autism.

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

Resolving stem and progenitor cells in the adult mouse incisor through gene co-expression analysis.

Investigations into stem cell-fueled renewal of an organ benefit from an inventory of cell type-specific markers and a deep understanding of the cellular diversity within stem cell niches. Using the adult mouse incisor as a model for a continuously renewing organ, we performed an unbiased analysis of gene co-expression relationships to identify modules of co-expressed genes that represent differentiated cells, transit-amplifying cells, and residents of stem cell niches. Through in vivo lineage tracing, we demonstrated the power of this approach by showing that co-expression module members Lrig1 and Igfbp5 define populations of incisor epithelial and mesenchymal stem cells. We further discovered that two adjacent mesenchymal tissues, the periodontium and dental pulp, are maintained by distinct pools of stem cells. These findings reveal novel mechanisms of incisor renewal and illustrate how gene co-expression analysis of intact biological systems can provide insights into the transcriptional basis of cellular identity.

An FAK-YAP-mTOR Signaling Axis Regulates Stem Cell-Based Tissue Renewal in Mice.

Tissue homeostasis requires the production of newly differentiated cells from resident adult stem cells. Central to this process is the expansion of undifferentiated intermediates known as transit-amplifying (TA) cells, but how stem cells are triggered to enter this proliferative TA state remains an important open question. Using the continuously growing mouse incisor as a model of stem cell-based tissue renewal, we found that the transcriptional cofactors YAP and TAZ are required both to maintain TA cell proliferation and to inhibit differentiation. Specifically, we identified a pathway involving activation of integrin α3 in TA cells that signals through an LATS-independent FAK/CDC42/PP1A cascade to control YAP-S397 phosphorylation and nuclear localization. This leads to Rheb expression and potentiates mTOR signaling to drive the proliferation of TA cells. These findings thus reveal a YAP/TAZ signaling mechanism that coordinates stem cell expansion and differentiation during organ renewal.

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.