Integrative inference of brain cell similarities and differences from single-cell genomics

Defining cell types requires integrating diverse measurements from multiple experiments and biological contexts. Recent technological developments in single-cell analysis have enabled high-throughput profiling of gene expression, epigenetic regulation, and spatial relationships amongst cells in complex tissues, but computational approaches that deliver a sensitive and specific joint analysis of these datasets are lacking. We developed LIGER, an algorithm that delineates shared and dataset-specific features of cell identity, allowing flexible modeling of highly heterogeneous single-cell datasets. We demonstrated its broad utility by applying it to four diverse and challenging analyses of human and mouse brain cells. First, we defined both cell-type-specific and sexually dimorphic gene expression in the mouse bed nucleus of the stria terminalis, an anatomically complex brain region that plays important roles in sex-specific behaviors. Second, we analyzed gene expression in the substantia nigra of seven postmortem human subjects, comparing cell states in specific donors, and relating cell types to those in the mouse. Third, we jointly leveraged in situ gene expression and scRNA-seq data to spatially locate fine subtypes of cells present in the mouse frontal cortex. Finally, we integrated mouse cortical scRNA-seq profiles with single-cell DNA methylation signatures, revealing mechanisms of cell-type-specific gene regulation. Integrative analyses using the LIGER algorithm promise to accelerate single-cell investigations of cell-type definition, gene regulation, and disease states.

[1]  R. Gorski,et al.  Sex difference in the bed nucleus of the stria terminalis of the human brain , 1990, The Journal of comparative neurology.

[2]  R. Gorski,et al.  Sex differences in subregions of the medial nucleus of the amygdala and the bed nucleus of the stria terminalis of the rat , 1992, Brain Research.

[3]  G. Fishell,et al.  The caudal ganglionic eminence is a source of distinct cortical and subcortical cell populations , 2002, Nature Neuroscience.

[4]  Xin Liu,et al.  Document clustering based on non-negative matrix factorization , 2003, SIGIR.

[5]  S. Amir,et al.  A Circadian Rhythm in the Expression of PERIOD2 Protein Reveals a Novel SCN-Controlled Oscillator in the Oval Nucleus of the Bed Nucleus of the Stria Terminalis , 2004, The Journal of Neuroscience.

[6]  D. Haussler,et al.  Ultraconserved Elements in the Human Genome , 2004, Science.

[7]  Larry W Swanson,et al.  Projections from bed nuclei of the stria terminalis, posterior division: Implications for cerebral hemisphere regulation of defensive and reproductive behaviors , 2004, The Journal of comparative neurology.

[8]  P. Greengard,et al.  Localization of DARPP-32 immunoreactive neurons in the bed nucleus of the stria terminalis and central nucleus of the amygdala: co-distribution with axons containing tyrosine hydroxylase, vasoactive intestinal polypeptide, and calcitonin gene-related peptide , 2004, Experimental Brain Research.

[9]  E. Maywood,et al.  Developmental and reproductive performance in circadian mutant mice. , 2006, Human reproduction.

[10]  Cheng Li,et al.  Adjusting batch effects in microarray expression data using empirical Bayes methods. , 2007, Biostatistics.

[11]  Allan R. Jones,et al.  Genome-wide atlas of gene expression in the adult mouse brain , 2007, Nature.

[12]  J. Rubenstein,et al.  Arx Is a Direct Target of Dlx2 and Thereby Contributes to the Tangential Migration of GABAergic Interneurons , 2008, The Journal of Neuroscience.

[13]  M. Götz,et al.  Progeny of Olig2-Expressing Progenitors in the Gray and White Matter of the Adult Mouse Cerebral Cortex , 2008, The Journal of Neuroscience.

[14]  Israel Steinfeld,et al.  BMC Bioinformatics BioMed Central , 2008 .

[15]  S. Greenberg,et al.  Cerebral Amyloid Angiopathy: A Systematic Review , 2011, Journal of clinical neurology.

[16]  Matko Bosnjak,et al.  REVIGO Summarizes and Visualizes Long Lists of Gene Ontology Terms , 2011, PloS one.

[17]  Ilona B. Bruinsma,et al.  Small Heat Shock Proteins Induce a Cerebral Inflammatory Reaction , 2011, The Journal of Neuroscience.

[18]  William Stafford Noble,et al.  FIMO: scanning for occurrences of a given motif , 2011, Bioinform..

[19]  D. Ron,et al.  Selective Inhibition of a Regulatory Subunit of Protein Phosphatase 1 Restores Proteostasis , 2011, Science.

[20]  G. Fishell,et al.  Three groups of interneurons account for nearly 100% of neocortical GABAergic neurons , 2011, Developmental neurobiology.

[21]  Masahiko Watanabe,et al.  Three Types of Neurochemical Projection from the Bed Nucleus of the Stria Terminalis to the Ventral Tegmental Area in Adult Mice , 2012, The Journal of Neuroscience.

[22]  M. Capogna,et al.  Neurogliaform cells of amygdala: a source of slow phasic inhibition in the basolateral complex , 2012, The Journal of physiology.

[23]  M. Chopp,et al.  Neuroprotective and neurorestorative effects of thymosin β4 treatment following experimental traumatic brain injury , 2012, Annals of the New York Academy of Sciences.

[24]  Data production leads,et al.  An integrated encyclopedia of DNA elements in the human genome , 2012 .

[25]  ENCODEConsortium,et al.  An Integrated Encyclopedia of DNA Elements in the Human Genome , 2012, Nature.

[26]  Tetsuro Izumi,et al.  Modular Genetic Control of Sexually Dimorphic Behaviors , 2012, Cell.

[27]  James H. Marshel,et al.  Diverging neural pathways assemble a behavioural state from separable features in anxiety , 2013, Nature.

[28]  J. Relton,et al.  Colony-stimulating factor 1 receptor (CSF1R) signaling in injured neurons facilitates protection and survival , 2013, The Journal of experimental medicine.

[29]  S. Dudoit,et al.  Normalization of RNA-seq data using factor analysis of control genes or samples , 2014, Nature Biotechnology.

[30]  Haesun Park,et al.  Algorithms for nonnegative matrix and tensor factorizations: a unified view based on block coordinate descent framework , 2014, J. Glob. Optim..

[31]  Hongkui Zeng,et al.  Correlated gene expression and target specificity demonstrate excitatory projection neuron diversity. , 2015, Cerebral cortex.

[32]  Terrence J. Sejnowski,et al.  Epigenomic Signatures of Neuronal Diversity in the Mammalian Brain , 2015, Neuron.

[33]  Allon M. Klein,et al.  Droplet Barcoding for Single-Cell Transcriptomics Applied to Embryonic Stem Cells , 2015, Cell.

[34]  Evan Z. Macosko,et al.  Highly Parallel Genome-wide Expression Profiling of Individual Cells Using Nanoliter Droplets , 2015, Cell.

[35]  Marcello Pelillo,et al.  Context aware nonnegative matrix factorization clustering , 2016, 2016 23rd International Conference on Pattern Recognition (ICPR).

[36]  Samuel L. Wolock,et al.  A Single-Cell Transcriptomic Map of the Human and Mouse Pancreas Reveals Inter- and Intra-cell Population Structure. , 2016, Cell systems.

[37]  G. Wagner,et al.  The origin and evolution of cell types , 2016, Nature Reviews Genetics.

[38]  N. Shah,et al.  Genetic dissection of neural circuits underlying sexually dimorphic social behaviours , 2016, Philosophical Transactions of the Royal Society B: Biological Sciences.

[39]  X. Zhuang,et al.  RNA Imaging with Multiplexed Error-Robust Fluorescence In Situ Hybridization (MERFISH). , 2016, Methods in enzymology.

[40]  Ahmet F. Coskun,et al.  Dense transcript profiling in single cells by image correlation decoding , 2016, Nature Methods.

[41]  Christof Koch,et al.  Adult Mouse Cortical Cell Taxonomy by Single Cell Transcriptomics , 2016, Nature Neuroscience.

[42]  M. Ronaghi,et al.  Neuronal subtypes and diversity revealed by single-nucleus RNA sequencing of the human brain , 2016, Science.

[43]  Allan R. Jones,et al.  Comprehensive cellular‐resolution atlas of the adult human brain , 2016, The Journal of comparative neurology.

[44]  Harrison W. Gabel,et al.  DNA methylation in the gene body influences MeCP2-mediated gene repression , 2016, Proceedings of the National Academy of Sciences.

[45]  D. Fitzpatrick,et al.  Opportunities and challenges in modeling human brain disorders in transgenic primates , 2016, Nature Neuroscience.

[46]  George Michailidis,et al.  A non-negative matrix factorization method for detecting modules in heterogeneous omics multi-modal data , 2015, Bioinform..

[47]  Conor Fitzpatrick,et al.  Nuclear RNA-seq of single neurons reveals molecular signatures of activation , 2016, Nature communications.

[48]  V. Gradinaru,et al.  Engineered AAVs for efficient noninvasive gene delivery to the central and peripheral nervous systems , 2017, Nature Neuroscience.

[49]  Grace X. Y. Zheng,et al.  Massively parallel digital transcriptional profiling of single cells , 2016, Nature Communications.

[50]  Hiromasa Funato,et al.  Moxd1 Is a Marker for Sexual Dimorphism in the Medial Preoptic Area, Bed Nucleus of the Stria Terminalis and Medial Amygdala , 2017, Front. Neuroanat..

[51]  I. Amit,et al.  A Unique Microglia Type Associated with Restricting Development of Alzheimer’s Disease , 2017, Cell.

[52]  Comprehensive cellular‐resolution atlas of the adult human brain , 2017, The Journal of comparative neurology.

[53]  Zhaolan Zhou,et al.  The Crucial Role of DNA Methylation and MeCP2 in Neuronal Function , 2017, Genes.

[54]  Allon M. Klein,et al.  Single-Cell Analysis of Experience-Dependent Transcriptomic States in Mouse Visual Cortex , 2017, Nature Neuroscience.

[55]  J. C. Love,et al.  Seq-Well: portable, low-cost RNA sequencing of single cells at high throughput , 2017, Nature Methods.

[56]  Aviv Regev,et al.  Massively-parallel single nucleus RNA-seq with DroNc-seq , 2017, Nature Methods.

[57]  Daniel R. Berger,et al.  Cell diversity and network dynamics in photosensitive human brain organoids , 2017, Nature.

[58]  Justin P Sandoval,et al.  Single-cell methylomes identify neuronal subtypes and regulatory elements in mammalian cortex , 2017, Science.

[59]  Jonathan A. Bernstein,et al.  Assembly of functionally integrated human forebrain spheroids , 2017, Nature.

[60]  Andrew C. Adey,et al.  Highly scalable generation of DNA methylation profiles in single cells , 2018, Nature Biotechnology.

[61]  Christof Koch,et al.  Conserved cell types with divergent features between human and mouse cortex , 2018, bioRxiv.

[62]  Tyler H. Garvin,et al.  Ultraconserved Enhancers Are Required for Normal Development , 2018, Cell.

[63]  Paul Hoffman,et al.  Integrating single-cell transcriptomic data across different conditions, technologies, and species , 2018, Nature Biotechnology.

[64]  Evan Z. Macosko,et al.  Molecular Diversity and Specializations among the Cells of the Adult Mouse Brain , 2018, Cell.

[65]  Laleh Haghverdi,et al.  Batch effects in single-cell RNA-sequencing data are corrected by matching mutual nearest neighbors , 2018, Nature Biotechnology.

[66]  Lars E. Borm,et al.  Molecular Architecture of the Mouse Nervous System , 2018, Cell.

[67]  M. Kreutz,et al.  The Role of Activity-Dependent DNA Demethylation in the Adult Brain and in Neurological Disorders , 2018, Front. Mol. Neurosci..

[68]  David J. Arenillas,et al.  JASPAR 2018: update of the open-access database of transcription factor binding profiles and its web framework , 2017, Nucleic acids research.

[69]  William E. Allen,et al.  Three-dimensional intact-tissue sequencing of single-cell transcriptional states , 2018, Science.

[70]  William S. DeWitt,et al.  A Single-Cell Atlas of In Vivo Mammalian Chromatin Accessibility , 2018, Cell.

[71]  Richard A. Muscat,et al.  Single-cell profiling of the developing mouse brain and spinal cord with split-pool barcoding , 2018, Science.

[72]  Tracy M. Yamawaki,et al.  Evolution of pallium, hippocampus, and cortical cell types revealed by single-cell transcriptomics in reptiles , 2018, Science.

[73]  Brain and Neurological Disorders , .