scRNA-seq reveals the diversity of the developing cardiac cell lineage and molecular building blocks of the primary pacemaker

The heart is comprised of a variety of specialized cell types that work in unison to maintain blood flow. Here we utilized scRNA-seq analysis to delineate the diversity of cardiac cell types in the zebrafish. With the growing use of the zebrafish to model human heart biology, a deeper insight into its complex cellular composition is critical for a better understanding of heart function, development, and associated malformations. We present a high resolution atlas of zebrafish heart single cells transcriptomics, consisting of over 50 000 cells representing the building blocks of the zebrafish heart at 48 and 72 hpf. We defined 18 discrete cell populations comprising major cell lineages and sublineages of the developing heart. We pinpointed a population of cells likely to be the primary pacemaker and identified the transcriptome profile defining this critical cell type. Our analyses identified two genes, atp1b3b and colec10, which were enriched in the sinoatrial pacemaker cells. CRISPR/Cas9-mediated knockout of these two genes significantly reduced heart rate which is accompanied by arrhythmia or morphological defects, suggesting their novel function in cardiac development and conduction. Additionally, we describe other subpopulations of cardiac cell lineages, including the endothelial and neural cells, whose expression profiles we provide as a resource for further investigations into the cellular and molecular mechanisms of this organ.

[1]  A. Zielezinski,et al.  A ligand-receptor interactome atlas of the zebrafish , 2023, bioRxiv.

[2]  D. Stainier,et al.  Origin and function of activated fibroblast states during zebrafish heart regeneration , 2022, Nature Genetics.

[3]  R. Lovering,et al.  GOing Forward With the Cardiac Conduction System Using Gene Ontology , 2022, Frontiers in Genetics.

[4]  Åsa K. Björklund,et al.  High cardiomyocyte diversity in human early prenatal heart development , 2022, bioRxiv.

[5]  Guoji Guo,et al.  Characterization of the Zebrafish Cell Landscape at Single-Cell Resolution , 2021, Frontiers in Cell and Developmental Biology.

[6]  P. Kohl,et al.  Genomic and physiological analyses of the zebrafish atrioventricular canal reveal molecular building blocks of the secondary pacemaker region , 2021, Cellular and Molecular Life Sciences.

[7]  K. Yutzey,et al.  Periostin-expressing Schwann cells and endoneurial cardiac fibroblasts contribute to sympathetic nerve fasciculation after birth. , 2021, Journal of molecular and cellular cardiology.

[8]  T. Brand,et al.  The Intrinsic Cardiac Nervous System and Its Role in Cardiac Pacemaking and Conduction , 2020, Journal of cardiovascular development and disease.

[9]  Joseph M. Fernandez,et al.  A simple and effective F0 knockout method for rapid screening of behaviour and other complex phenotypes , 2020, bioRxiv.

[10]  K. Ardlie,et al.  Transcriptional and Cellular Diversity of the Human Heart , 2020, bioRxiv.

[11]  James A. Gagnon,et al.  Zebrafish dscaml1 Deficiency Impairs Retinal Patterning and Oculomotor Function , 2020, The Journal of Neuroscience.

[12]  V. Prince,et al.  Neural crest development: insights from the zebrafish , 2020, Developmental dynamics : an official publication of the American Association of Anatomists.

[13]  Mauro J. Muraro,et al.  Single-cell analysis uncovers that metabolic reprogramming by ErbB2 signaling is essential for cardiomyocyte proliferation in the regenerating heart , 2019, eLife.

[14]  Joakim Lundeberg,et al.  A Spatiotemporal Organ-Wide Gene Expression and Cell Atlas of the Developing Human Heart , 2019, Cell.

[15]  A. Charloux,et al.  The Endocrine Function of the Heart: Physiology and Involvements of Natriuretic Peptides and Cyclic Nucleotide Phosphodiesterases in Heart Failure , 2019, Journal of clinical medicine.

[16]  Kornel Labun,et al.  CHOPCHOP v3: expanding the CRISPR web toolbox beyond genome editing , 2019, Nucleic Acids Res..

[17]  S. Sumanas,et al.  Clec14a genetically interacts with Etv2 and Vegf signaling during vasculogenesis and angiogenesis in zebrafish , 2019, BMC Developmental Biology.

[18]  R. Satija,et al.  Normalization and variance stabilization of single-cell RNA-seq data using regularized negative binomial regression , 2019, Genome Biology.

[19]  S. Lakhal-Littleton Mechanisms of cardiac iron homeostasis and their importance to heart function , 2019, Free radical biology & medicine.

[20]  Anushya Muruganujan,et al.  PANTHER version 14: more genomes, a new PANTHER GO-slim and improvements in enrichment analysis tools , 2018, Nucleic Acids Res..

[21]  T. Sauka-Spengler,et al.  Functional Heterogeneity within the Developing Zebrafish Epicardium , 2018, bioRxiv.

[22]  Erica Bresciani,et al.  An efficient dissociation protocol for generation of single cell suspension from zebrafish embryos and larvae , 2018, MethodsX.

[23]  A. van Oudenaarden,et al.  Single-Cell Sequencing of the Healthy and Diseased Heart Reveals Cytoskeleton-Associated Protein 4 as a New Modulator of Fibroblasts Activation , 2018, Circulation.

[24]  M. Giorgi,et al.  Egr-1 Maintains NSC Proliferation and Its Overexpression Counteracts Cell Cycle Exit Triggered by the Withdrawal of Epidermal Growth Factor , 2018, Developmental Neuroscience.

[25]  Rahul C. Deo,et al.  A Rapid Method for Directed Gene Knockout for Screening in G0 Zebrafish. , 2018, Developmental cell.

[26]  Zev J. Gartner,et al.  DoubletFinder: Doublet detection in single-cell RNA sequencing data using artificial nearest neighbors , 2018, bioRxiv.

[27]  J. Bakkers,et al.  Spatially resolved RNA-sequencing of the embryonic heart identifies a role for Wnt/β-catenin signaling in autonomic control of heart rate , 2018, eLife.

[28]  N. Mercader,et al.  Transient fibrosis resolves via fibroblast inactivation in the regenerating zebrafish heart , 2018, Proceedings of the National Academy of Sciences.

[29]  J. Molkentin,et al.  Redefining the identity of cardiac fibroblasts , 2017, Nature Reviews Cardiology.

[30]  P. Beales,et al.  COLEC10 is mutated in 3MC patients and regulates early craniofacial development , 2017, PLoS genetics.

[31]  D. Stainier,et al.  Isl2b regulates anterior second heart field development in zebrafish , 2017, Scientific Reports.

[32]  S. Husson,et al.  Peptidomics of the zebrafish Danio rerio: In search for neuropeptides. , 2017, Journal of proteomics.

[33]  R. Croll,et al.  Zebrafish heart as a model to study the integrative autonomic control of pacemaker function. , 2016, American journal of physiology. Heart and circulatory physiology.

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

[35]  M. Liebling,et al.  Development of the cardiac conduction system in zebrafish. , 2016, Gene expression patterns : GEP.

[36]  M. Robinson,et al.  Maximizing mutagenesis with solubilized CRISPR-Cas9 ribonucleoprotein complexes , 2016, Development.

[37]  S. Abdelilah-Seyfried,et al.  The force within: endocardial development, mechanotransduction and signalling during cardiac morphogenesis , 2016, Development.

[38]  M. Vornanen,et al.  Zebrafish heart as a model for human cardiac electrophysiology , 2016, Channels.

[39]  S. Kuraku,et al.  Evolution of the fish heart by sub/neofunctionalization of an elastin gene , 2016, Nature Communications.

[40]  F. Grosveld,et al.  The Isl1/Ldb1 Complex Orchestrates Genome-wide Chromatin Organization to Instruct Differentiation of Multipotent Cardiac Progenitors. , 2015, Cell stem cell.

[41]  Ann M. Cavanaugh,et al.  Two developmentally distinct populations of neural crest cells contribute to the zebrafish heart. , 2015, Developmental biology.

[42]  Mitchell P. Levesque,et al.  A Floor-Plate Extracellular Protein-Protein Interaction Screen Identifies Draxin as a Secreted Netrin-1 Antagonist. , 2015, Cell Reports.

[43]  Arndt F. Siekmann,et al.  Chemokine-guided angiogenesis directs coronary vasculature formation in zebrafish. , 2015, Developmental cell.

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

[45]  W. Pu,et al.  Cellular origin and developmental program of coronary angiogenesis. , 2015, Circulation research.

[46]  S. Abdelilah-Seyfried,et al.  Large-scale zebrafish embryonic heart dissection for transcriptional analysis. , 2015, Journal of visualized experiments : JoVE.

[47]  T. Huh,et al.  Zebrafish Crip2 Plays a Critical Role in Atrioventricular Valve Development by Downregulating the Expression of ECM Genes in the Endocardial Cushion , 2014, Molecules and cells.

[48]  M. Bronner,et al.  Insights into neural crest development and evolution from genomic analysis , 2013, Genome research.

[49]  H. Takeda,et al.  Neuronal Subtypes Are Specified by the Level of neurod Expression in the Zebrafish Lateral Line , 2013, The Journal of Neuroscience.

[50]  David Staudt,et al.  Uncovering the molecular and cellular mechanisms of heart development using the zebrafish. , 2012, Annual review of genetics.

[51]  Arie O. Verkerk,et al.  Identification and Functional Characterization of Cardiac Pacemaker Cells in Zebrafish , 2012, PloS one.

[52]  Guangchuang Yu,et al.  clusterProfiler: an R package for comparing biological themes among gene clusters. , 2012, Omics : a journal of integrative biology.

[53]  D. Raible,et al.  Bmps and Id2a Act Upstream of Twist1 To Restrict Ectomesenchyme Potential of the Cranial Neural Crest , 2012, PLoS genetics.

[54]  S. Sumanas,et al.  Hedgehog signaling is required for differentiation of endocardial progenitors in zebrafish. , 2012, Developmental biology.

[55]  Jonathan J. Sager,et al.  Hypokinesia and Reduced Dopamine Levels in Zebrafish Lacking β- and γ1-Synucleins* , 2011, The Journal of Biological Chemistry.

[56]  M. Shattock,et al.  The Inhibitory Effect of Phospholemman on the Sodium Pump Requires Its Palmitoylation* , 2011, The Journal of Biological Chemistry.

[57]  Joshua D. Wythe,et al.  A Slit/miR-218/Robo regulatory loop is required during heart tube formation in zebrafish , 2011, Development.

[58]  M. Howell,et al.  SNW1 Is a Critical Regulator of Spatial BMP Activity, Neural Plate Border Formation, and Neural Crest Specification in Vertebrate Embryos , 2011, PLoS biology.

[59]  G. A. Wilkinson,et al.  Nogo-B receptor is essential for angiogenesis in zebrafish via Akt pathway. , 2010, Blood.

[60]  Annemarie H Meijer,et al.  Macrophage-specific gene functions in Spi1-directed innate immunity. , 2010, Blood.

[61]  M. Veldman,et al.  Tuba1a gene expression is regulated by KLF6/7 and is necessary for CNS development and regeneration in zebrafish , 2010, Molecular and Cellular Neuroscience.

[62]  M. Liebling,et al.  Zebrafish cardiac enhancer trap lines: New tools for in vivo studies of cardiovascular development and disease , 2010, Developmental dynamics : an official publication of the American Association of Anatomists.

[63]  L. Clijsters,et al.  Distinct phases of cardiomyocyte differentiation regulate growth of the zebrafish heart , 2009, Development.

[64]  E. Conway,et al.  Role of the 2 zebrafish survivin genes in vasculo-angiogenesis, neurogenesis, cardiogenesis and hematopoiesis , 2009, BMC Developmental Biology.

[65]  Shuo Lin,et al.  Discovery and Characterization of Novel Vascular and Hematopoietic Genes Downstream of Etsrp in Zebrafish , 2009, PloS one.

[66]  Herwig Baier,et al.  Genetic and Physiologic Dissection of the Vertebrate Cardiac Conduction System , 2008, PLoS biology.

[67]  F. Serluca Development of the proepicardial organ in the zebrafish. , 2008, Developmental biology.

[68]  S. Abdelilah-Seyfried,et al.  Asymmetric Involution of the Myocardial Field Drives Heart Tube Formation in Zebrafish , 2008, Circulation research.

[69]  J. Bakkers,et al.  Early Endocardial Morphogenesis Requires Scl/Tal1 , 2007, PLoS genetics.

[70]  M. Samuels,et al.  The brain-heart connection. , 2007, Circulation.

[71]  D. Stainier,et al.  A Mutation in Zebrafish hmgcr1b Reveals a Role for Isoprenoids in Vertebrate Heart-Tube Formation , 2007, Current Biology.

[72]  Calum A MacRae,et al.  Notch1b and neuregulin are required for specification of central cardiac conduction tissue , 2006, Development.

[73]  L. Zon,et al.  A Critical Role for Calponin 2 in Vascular Development* , 2006, Journal of Biological Chemistry.

[74]  B. Paw,et al.  Analysis of thrombocyte development in CD41-GFP transgenic zebrafish. , 2005, Blood.

[75]  Heather Verkade,et al.  Genetic and cellular analyses of zebrafish atrioventricular cushion and valve development , 2005, Development.

[76]  Shuo Lin,et al.  15-zinc finger protein Bloody Fingers is required for zebrafish morphogenetic movements during neurulation. , 2005, Developmental biology.

[77]  D. Stainier,et al.  Microarray analysis of zebrafish cloche mutant using amplified cDNA and identification of potential downstream target genes , 2005, Developmental dynamics : an official publication of the American Association of Anatomists.

[78]  Robert H. Anderson,et al.  Lineage and Morphogenetic Analysis of the Cardiac Valves , 2004, Circulation research.

[79]  P. Byers,et al.  Rare autosomal recessive cardiac valvular form of Ehlers-Danlos syndrome results from mutations in the COL1A2 gene that activate the nonsense-mediated RNA decay pathway. , 2004, American journal of human genetics.

[80]  S. Ekker,et al.  Syndecan-2 is essential for angiogenic sprouting during zebrafish development. , 2004, Blood.

[81]  S. Korsching,et al.  Zebrafish beta tubulin 1 expression is limited to the nervous system throughout development, and in the adult brain is restricted to a subset of proliferative regions. , 2004, Gene expression patterns : GEP.

[82]  Didier Y. R. Stainier,et al.  Mutation of weak atrium/atrial myosin heavy chain disrupts atrial function and influences ventricular morphogenesis in zebrafish , 2003, Development.

[83]  R. Geisler,et al.  lockjaw encodes a zebrafish tfap2a required for early neural crest development , 2003, Development.

[84]  J. Volff,et al.  Functional divergence of two zebrafish midkine growth factors following fish-specific gene duplication. , 2003, Genome research.

[85]  A. Yamanaka,et al.  Characterization of a family of endogenous neuropeptide ligands for the G protein-coupled receptors GPR7 and GPR8 , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[86]  Didier Y. R. Stainier,et al.  Cardiac troponin T is essential in sarcomere assembly and cardiac contractility , 2002, Nature Genetics.

[87]  A. Meng,et al.  Zebrafish sox9b is an early neural crest marker , 2002, Development Genes and Evolution.

[88]  P. Olsson,et al.  Generating Transparent Zebrafish: A Refined Method to Improve Detection of Gene Expression During Embryonic Development , 2001, Marine Biotechnology.

[89]  A. Hudspeth,et al.  Expression and phylogeny of claudins in vertebrate primordia , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[90]  E. Olson,et al.  Calsarcins, a novel family of sarcomeric calcineurin-binding proteins. , 2000, Proceedings of the National Academy of Sciences of the United States of America.

[91]  L. Brodin,et al.  Zebrafish Genes for Neuropeptide Y and Peptide YY Reveal Origin by Chromosome Duplication from an Ancestral Gene Linked to the Homeobox Cluster , 2000, Journal of neurochemistry.

[92]  C. Kimmel,et al.  lazarus is a novel pbx gene that globally mediates hox gene function in zebrafish. , 2000, Molecular cell.

[93]  S. Horne,et al.  Restricted expression of cardiac myosin genes reveals regulated aspects of heart tube assembly in zebrafish. , 1999, Developmental biology.

[94]  H. Yamazaki,et al.  Molecular Cloning of a Novel Human Collectin from Liver (CL-L1)* , 1999, The Journal of Biological Chemistry.

[95]  D. Stainier,et al.  Patterning during organogenesis: genetic analysis of cardiac chamber formation. , 1999, Seminars in cell & developmental biology.

[96]  A. Brownlie,et al.  Positional cloning of the zebrafish sauternes gene: a model for congenital sideroblastic anaemia , 1998, Nature Genetics.

[97]  B. Göttgens,et al.  The SCL gene specifies haemangioblast development from early mesoderm , 1998, The EMBO journal.

[98]  A. Amores,et al.  The cloche and spadetail genes differentially affect hematopoiesis and vasculogenesis. , 1998, Developmental biology.

[99]  D. Stainier,et al.  Isolation of the zebrafish homologues for the tie‐1 and tie‐2 endothelium‐specific receptor tyrosine kinases , 1998, Developmental dynamics : an official publication of the American Association of Anatomists.

[100]  W H Lamers,et al.  Development of the cardiac conduction system. , 1998, Circulation research.

[101]  J. Campos-Ortega,et al.  A zebrafish Id homologue and its pattern of expression during embryogenesis , 1997, Mechanisms of Development.

[102]  Cheol‐Hee Kim,et al.  Zebrafish elav/HuC homologue as a very early neuronal marker , 1996, Neuroscience Letters.

[103]  P. Gros,et al.  Identification of the Mammalian Na,K-ATPase β3 Subunit* , 1996, The Journal of Biological Chemistry.

[104]  R. Markwald,et al.  Formation and early morphogenesis of endocardial endothelial precursor cells and the role of endoderm. , 1996, Developmental biology.

[105]  A. McDonough,et al.  Subcellular distribution of sodium pump isoform subunits in mammalian cardiac myocytes. , 1996, The American journal of physiology.

[106]  M. Schachner,et al.  Tenascin-C mRNA is expressed in cranial neural crest cells, in some placodal derivatives, and in discrete domains of the embryonic zebrafish brain. , 1995, Journal of neurobiology.

[107]  M. Fishman,et al.  Cardiovascular development in the zebrafish. I. Myocardial fate map and heart tube formation. , 1993, Development.

[108]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .

[109]  A. M. Spanier,et al.  Structural and functional characteristics of cardiac myocytes. , 1985, Basic research in cardiology.