Cell type transcriptome atlas for the planarian Schmidtea mediterranea

Mapping the planarian transcriptome A cell type's transcriptome defines the active genes that control its biology. Two groups used single-cell RNA sequencing to define the transcriptomes for essentially all cell types of a complete animal, the regenerative planarian Schmidtea mediterranea. Because pluripotent stem cells constantly differentiate to rejuvenate any part of the body of this species, all developmental lineages are active in adult animals. Fincher et al. determined the transcriptomes for most, if not all, planarian cell types, including some that were previously unknown. They also identified transition states and genes governing positional information. Plass et al. used single-cell transcriptomics and computational algorithms to reconstruct a lineage tree capturing the developmental progressions from stem to differentiated cells. They could then predict gene programs that are specifically turned on and off along the tree, and they used this approach to study how the cell types behaved during regeneration. These whole-animal transcriptome “atlases” are a powerful way to study metazoan biology. Science, this issue p. eaaq1736, p. eaaq1723 Single-cell RNA sequencing identifies the transcriptomes for most to all cell types of a flatworm. INTRODUCTION The complete sequence of animal genomes has had a transformative impact on biological research. Whereas the genome sequence of an organism contains the information for its development and physiology, the transcriptomes (the sets of actively transcribed genes) of the cell types in an organism define how the genome is used for the unique functions of its cells. Cell number and complexity have historically made the identification of all cell types, much less their transcriptomes, an extreme challenge for most multicellular organisms. Recent advances in single-cell RNA sequencing (SCS) have greatly enhanced the ability to determine cell type transcriptomes, with SCS of thousands of cells readily achievable. RATIONALE We reasoned that it might be possible, given these advances, to determine the transcriptomes of essentially every cell type of a complete organism possessing an unknown number of cell types. The planarian Schmidtea mediterranea, famous for its regeneration ability, is an attractive case study for such an undertaking. Planarians possess a complex anatomy with diverse differentiated cell types, including many found across animals. Furthermore, planarians contain a proliferating cell population called neoblasts that includes pluripotent stem cells. Neoblasts mediate regeneration and constitutive tissue turnover. Consequently, lineage precursors for essentially all differentiated cells types are also present in adults. Finally, planarians constitutively express positional information guiding tissue turnover. Therefore, comprehensive SCS at a single time point (the adult) could allow transcriptome determination for all differentiated cell types and for lineage precursors, and could identify patterning information that guides new cell production and organization. Capturing this information in most organisms would require sampling adults and many transient embryonic stages. RESULTS We used the SCS method Drop-seq to determine the transcriptomes for 66,783 individual cells from adult planarians. We locally saturated cell type coverage by iteratively sequencing distinct body regions and assessing the frequency of known rare cell types in the data. Clustering the cells by shared gene expression grouped cells into broad tissue classes. Subclustering of each broad tissue type in isolation enabled separation of cells into the cell populations constituting each tissue. These analyses enabled the identification of a previously unidentified tissue group and the classification of poorly characterized tissues into their constituent cell types, including numerous previously unknown cell types. Transcriptomes were identified for many rare cell types, including those that exist as rarely as ~10 cells in an animal that has 105 to 106 cells, which suggests that near-to-complete cellular saturation was reached. In addition, transcriptomes for known and novel lineage precursors, from pluripotent stem cell to differentiated cell types, were generated. Precursor transcriptomes identified transcription factors required for maintenance of associated differentiated cells during homeostatic cell turnover. Finally, the data were used to identify genes regionally expressed in muscle, which is the site of planarian patterning gene expression. CONCLUSION We successfully used SCS to generate transcriptomes for most to all cells of a complete organism. This resource provides a wealth of data regarding the cellular site of expression of thousands of conserved genes and the transcriptomes for cell types widely used in animals. These data will inform studies of these genes and cell types broadly, and will provide a resource for the fields of planarian biology and comparative evolutionary biology. This work also provides a template for the generation of cell type transcriptome atlases, which can be applied to a large array of organisms. An atlas of planarian cell type transcriptomes. High-throughput single-cell RNA sequencing of adult planarians reveals a cell type transcriptome atlas that includes rare cell types and many novel cell populations, cellular transition states, and patterning information, as demonstrated by two regionally expressed genes in muscle. The transcriptome of a cell dictates its unique cell type biology. We used single-cell RNA sequencing to determine the transcriptomes for essentially every cell type of a complete animal: the regenerative planarian Schmidtea mediterranea. Planarians contain a diverse array of cell types, possess lineage progenitors for differentiated cells (including pluripotent stem cells), and constitutively express positional information, making them ideal for this undertaking. We generated data for 66,783 cells, defining transcriptomes for known and many previously unknown planarian cell types and for putative transition states between stem and differentiated cells. We also uncovered regionally expressed genes in muscle, which harbors positional information. Identifying the transcriptomes for potentially all cell types for many organisms should be readily attainable and represents a powerful approach to metazoan biology.

[1]  L. Hyman,et al.  A histochemical study of intracellular digestion in triclad flatworms , 1925 .

[2]  P. Reddien,et al.  Smed-βcatenin-1 Is Required for Anteroposterior Blastema Polarity in Planarian Regeneration , 2008, Science.

[3]  Alejandro Sánchez Alvarado,et al.  Allometric scaling and proportion regulation in the freshwater planarian Schmidtea mediterranea , 2003, Developmental dynamics : an official publication of the American Association of Anatomists.

[4]  Rachel H. Roberts-Galbraith,et al.  A functional genomics screen in planarians reveals regulators of whole-brain regeneration , 2016, eLife.

[5]  K. Yoshimoto,et al.  Reconstruction of dopaminergic neural network and locomotion function in planarian regenerates , 2007, Developmental neurobiology.

[6]  F. Cebrià,et al.  Planarian homologs of netrin and netrin receptor are required for proper regeneration of the central nervous system and the maintenance of nervous system architecture , 2005, Development.

[7]  Jochen C. Rink,et al.  The maintenance and regeneration of the planarian excretory system are regulated by EGFR signaling , 2011, Development.

[8]  S. Lapan,et al.  Eye Absence Does Not Regulate Planarian Stem Cells during Eye Regeneration. , 2017, Developmental cell.

[9]  Andrew Smith Genome sequence of the nematode C-elegans: A platform for investigating biology , 1998 .

[10]  J. Rink,et al.  Reactivating head regrowth in a regeneration-deficient planarian species , 2013, Nature.

[11]  A. Aboobaker,et al.  The TALE Class Homeobox Gene Smed-prep Defines the Anterior Compartment for Head Regeneration , 2010, PLoS genetics.

[12]  A. Sánchez Alvarado,et al.  Double-stranded RNA specifically disrupts gene expression during planarian regeneration. , 1999, Proceedings of the National Academy of Sciences of the United States of America.

[13]  P. Reddien,et al.  Fundamentals of planarian regeneration. , 2004, Annual review of cell and developmental biology.

[14]  F. Cebrià,et al.  Evolution of the EGFR pathway in Metazoa and its diversification in the planarian Schmidtea mediterranea , 2016, Scientific Reports.

[15]  Adam L. Bermange,et al.  BMP signaling regulates the dorsal planarian midline and is needed for asymmetric regeneration , 2007, Development.

[16]  K. Yoshimoto,et al.  Characterization of tyramine β-hydroxylase in planarian Dugesia japonica: Cloning and expression , 2008, Neurochemistry International.

[17]  P. Reddien Specialized progenitors and regeneration , 2013, Development.

[18]  K. Nishimura,et al.  Identification of glutamic acid decarboxylase gene and distribution of GABAergic nervous system in the planarian Dugesia japonica , 2008, Neuroscience.

[19]  P. Reddien,et al.  Planarian Epidermal Stem Cells Respond to Positional Cues to Promote Cell-Type Diversity. , 2017, Developmental cell.

[20]  Oleg Simakov,et al.  Expression of secreted Wnt pathway components reveals unexpected complexity of the planarian amputation response. , 2010, Developmental biology.

[21]  P. Reddien,et al.  Ingestion of bacterially expressed double-stranded RNA inhibits gene expression in planarians , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[22]  S. Ishii Electron microscopic observations on the Planarian tissues II. The intestine. , 1965, Fukushima journal of medical science.

[23]  Peter W. Reddien,et al.  Genetic regulators of a pluripotent adult stem cell system in planarians identified by RNAi and clonal analysis. , 2012, Cell stem cell.

[24]  P. Reddien,et al.  Two FGFRL-Wnt circuits organize the planarian anteroposterior axis , 2016, eLife.

[25]  T. Gojobori,et al.  Dissecting planarian central nervous system regeneration by the expression of neural‐specific genes , 2002, Development, growth & differentiation.

[26]  A. Regev,et al.  A Generic and Cell-Type-Specific Wound Response Precedes Regeneration in Planarians. , 2015, Developmental cell.

[27]  S. Lapan,et al.  Transcriptome analysis of the planarian eye identifies ovo as a specific regulator of eye regeneration. , 2012, Cell reports.

[28]  A. Sánchez Alvarado,et al.  Bromodeoxyuridine specifically labels the regenerative stem cells of planarians. , 2000, Developmental biology.

[29]  S. Orkin,et al.  Mapping the Mouse Cell Atlas by Microwell-Seq , 2018, Cell.

[30]  Cole Trapnell,et al.  The dynamics and regulators of cell fate decisions are revealed by pseudotemporal ordering of single cells , 2014, Nature Biotechnology.

[31]  A. Sánchez Alvarado,et al.  Molecular analysis of stem cells and their descendants during cell turnover and regeneration in the planarian Schmidtea mediterranea. , 2008, Cell stem cell.

[32]  Kentaro Kato,et al.  Molecular Cloning of Bone Morphogenetic Protein (BMP) Gene from the Planarian Dugesia japonica , 1998 .

[33]  Tetsutaro Hayashi,et al.  Isolation of planarian X‐ray‐sensitive stem cells by fluorescence‐activated cell sorting , 2006, Development, growth & differentiation.

[34]  Antoine H. F. M. Peters,et al.  A Bruno-like gene is required for stem cell maintenance in planarians. , 2006, Developmental cell.

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

[36]  Bret J. Pearson,et al.  A mex3 homolog is required for differentiation during planarian stem cell lineage development , 2015, eLife.

[37]  S. Lapan,et al.  Neoblast Specialization in Regeneration of the Planarian Schmidtea mediterranea , 2014, Stem cell reports.

[38]  F. Cebrià,et al.  Morphogenesis defects are associated with abnormal nervous system regeneration following roboA RNAi in planarians , 2007, Development.

[39]  I. Amit,et al.  Massively Parallel Single-Cell RNA-Seq for Marker-Free Decomposition of Tissues into Cell Types , 2014, Science.

[40]  P. Reddien,et al.  Planarian regeneration involves distinct stem cell responses to wounds and tissue absence. , 2010, Developmental biology.

[41]  Sylvain W. Lapan,et al.  dlx and sp6-9 Control Optic Cup Regeneration in a Prototypic Eye , 2011, PLoS genetics.

[42]  B. Pearson,et al.  Light-induced depigmentation in planarians models the pathophysiology of acute porphyrias , 2016, eLife.

[43]  Melainia McClain,et al.  Stem cells and fluid flow drive cyst formation in an invertebrate excretory organ , 2015, eLife.

[44]  P. Reddien,et al.  Single-cell analysis reveals functionally distinct classes within the planarian stem cell compartment. , 2014, Cell stem cell.

[45]  International Human Genome Sequencing Consortium Initial sequencing and analysis of the human genome , 2001, Nature.

[46]  Bret J. Pearson,et al.  FOX and ETS family transcription factors regulate the pigment cell lineage in planarians , 2017, Development.

[47]  D. Mccormick Sequence the Human Genome , 1986, Bio/Technology.

[48]  Tingxi Guo,et al.  nanos function is essential for development and regeneration of planarian germ cells , 2007, Proceedings of the National Academy of Sciences.

[49]  P. Kharchenko,et al.  Bayesian approach to single-cell differential expression analysis , 2014, Nature Methods.

[50]  K. J. Pedersen Some features of the fine structure and histochemistry of planarian subepidermal gland cells , 1959, Zeitschrift für Zellforschung und Mikroskopische Anatomie.

[51]  Andrew C. Adey,et al.  Single-Cell Transcriptional Profiling of a Multicellular Organism , 2017 .

[52]  Rona S. Gertner,et al.  Single cell RNA Seq reveals dynamic paracrine control of cellular variation , 2014, Nature.

[53]  P. Reddien,et al.  Orthogonal muscle fibers have different instructive roles in planarian regeneration , 2017, Nature.

[54]  David J. Forsthoefel,et al.  RNA interference by feeding in vitro–synthesized double‐stranded RNA to planarians: Methodology and dynamics , 2013, Developmental dynamics : an official publication of the American Association of Anatomists.

[55]  P. Reddien,et al.  Muscle Cells Provide Instructions for Planarian Regeneration , 2013, Cell reports.

[56]  Wallace F. Marshall,et al.  Centrosome Loss in the Evolution of Planarians , 2012, Science.

[57]  L. Gentile,et al.  Stem cell-dependent formation of a functional anterior regeneration pole in planarians requires Zic and Forkhead transcription factors. , 2014, Developmental biology.

[58]  K. J. Pedersen Studies on the nature of planarian connective tissue , 1961, Zeitschrift für Zellforschung und Mikroskopische Anatomie.

[59]  G. Giribet,et al.  Spiralian Phylogeny Informs the Evolution of Microscopic Lineages , 2015, Current Biology.

[60]  Charles H. Yoon,et al.  Dissecting the multicellular ecosystem of metastatic melanoma by single-cell RNA-seq , 2016, Science.

[61]  G. Bell,et al.  A regulatory program for excretory system regeneration in planarians , 2011, Development.

[62]  S. McKinney,et al.  Selective amputation of the pharynx identifies a FoxA-dependent regeneration program in planaria , 2014, eLife.

[63]  K. Yoshimoto,et al.  Identification and distribution of tryptophan hydroxylase (TPH)-positive neurons in the planarian Dugesia japonica , 2007, Neuroscience Research.

[64]  David J. Forsthoefel,et al.  An RNAi screen reveals intestinal regulators of branching morphogenesis, differentiation, and stem cell proliferation in planarians. , 2012, Developmental cell.

[65]  P. Newmark,et al.  In situ hybridization protocol for enhanced detection of gene expression in the planarian Schmidtea mediterranea , 2013, BMC Developmental Biology.

[66]  Tingxi Guo,et al.  The use of lectins as markers for differentiated secretory cells in planarians , 2010, Developmental dynamics : an official publication of the American Association of Anatomists.

[67]  Peter W. Reddien,et al.  Polarized notum Activation at Wounds Inhibits Wnt Function to Promote Planarian Head Regeneration , 2011, Science.

[68]  Teresa Adell,et al.  Smed-Evi/Wntless is required for β-catenin-dependent and -independent processes during planarian regeneration , 2009, Development.

[69]  J. Berg Genome sequence of the nematode C. elegans: a platform for investigating biology. , 1998, Science.

[70]  E. Saló,et al.  The BMP pathway is essential for re-specification and maintenance of the dorsoventral axis in regenerating and intact planarians. , 2007, Developmental biology.

[71]  Phillip A. Newmark,et al.  PRMT5 and the role of symmetrical dimethylarginine in chromatoid bodies of planarian stem cells. , 2012, Development.

[72]  E. Saló,et al.  Expression pattern of the expanded noggin gene family in the planarian Schmidtea mediterranea. , 2009, Gene expression patterns : GEP.

[73]  A. Regev,et al.  Spatial reconstruction of single-cell gene expression data , 2015 .

[74]  Y. Kitamura,et al.  Analysis of motor function modulated by cholinergic neurons in planarian dugesia japonica , 2010, Neuroscience.

[75]  Kerstin Bartscherer,et al.  A Pitx transcription factor controls the establishment and maintenance of the serotonergic lineage in planarians , 2013, Development.

[76]  Christian P. Petersen,et al.  Wnt/Notum spatial feedback inhibition controls neoblast differentiation to regulate reversible growth of the planarian brain , 2015, Development.

[77]  Irving E. Wang,et al.  Clonogenic Neoblasts Are Pluripotent Adult Stem Cells That Underlie Planarian Regeneration , 2011, Science.

[78]  Peter W. Reddien,et al.  SMEDWI-2 Is a PIWI-Like Protein That Regulates Planarian Stem Cells , 2005, Science.

[79]  L. Hyman,et al.  The Invertebrates: Platyhelminthes And Rhynchocoela:The Acoelomate Bilateria (Volume-2) , 1951 .

[80]  Andrew McDavid,et al.  Data exploration, quality control and testing in single-cell qPCR-based gene expression experiments , 2012, Bioinform..

[81]  Jochen C. Rink,et al.  β-Catenin Defines Head Versus Tail Identity During Planarian Regeneration and Homeostasis , 2008, Science.

[82]  Ignacio Maeso,et al.  Noggin and Noggin-Like Genes Control Dorsoventral Axis Regeneration in Planarians , 2011, Current Biology.

[83]  J. Collins,et al.  Genome-Wide Analyses Reveal a Role for Peptide Hormones in Planarian Germline Development , 2010, PLoS biology.

[84]  S. Lapan,et al.  A forkhead Transcription Factor Is Wound-Induced at the Planarian Midline and Required for Anterior Pole Regeneration , 2014, PLoS genetics.

[85]  P. Reddien,et al.  A wound-induced Wnt expression program controls planarian regeneration polarity , 2009, Proceedings of the National Academy of Sciences.

[86]  J. Baguñá,et al.  Quantitative analysis of cell types during growth, degrowth and regeneration in the planarians Dugesia mediterranea and Dugesia tigrina , 1981, Hydrobiologia.

[87]  S. McKinney,et al.  Egr-5 is a post-mitotic regulator of planarian epidermal differentiation , 2015, eLife.

[88]  Irving E. Wang,et al.  Hedgehog signaling regulates gene expression in planarian glia , 2016, eLife.

[89]  F. Cebrià Planarian Body-Wall Muscle: Regeneration and Function beyond a Simple Skeletal Support , 2016, Front. Cell Dev. Biol..

[90]  Jochen C. Rink,et al.  Planarian Hh Signaling Regulates Regeneration Polarity and Links Hh Pathway Evolution to Cilia , 2009, Science.

[91]  P. Reddien,et al.  Constitutive gene expression and the specification of tissue identity in adult planarian biology. , 2011, Trends in genetics : TIG.

[92]  Christian P. Petersen,et al.  zic-1 Expression in Planarian Neoblasts after Injury Controls Anterior Pole Regeneration , 2014, PLoS genetics.

[93]  Ning Ma,et al.  BLAST+: architecture and applications , 2009, BMC Bioinformatics.

[94]  L. Hyman,et al.  The Invertebrates, Vol. II: Platyhelminthes and Rhynchocoela , 1951 .

[95]  Peter W. Reddien,et al.  A Bmp/Admp Regulatory Circuit Controls Maintenance and Regeneration of Dorsal-Ventral Polarity in Planarians , 2011, Current Biology.