A human embryonic limb cell atlas resolved in 1 space and time 2

Human limbs emerge during the fourth post-conception week as mesenchymal buds which develop into fully-formed limbs over the subsequent months. Limb development is orchestrated by numerous temporally and spatially restricted gene expression programmes, making congenital alterations in phenotype common. Decades of work with model organisms has outlined the fundamental processes underlying vertebrate limb development, but an in-depth characterisation of this process in humans has yet to be performed. Here we detail the development of the human embryonic limb across space and time, using both single-cell and spatial transcriptomics. We demonstrate extensive diversification of cells, progressing from a restricted number of multipotent progenitors to myriad mature cell states, and identify several novel cell populations, including perineural fibroblasts and multiple distinct mesenchymal states. We uncover two waves of human muscle development, each characterised by different cell states regulated by separate gene expression programmes. We identify musculin (MSC) as a key transcriptional repressor maintaining muscle stem cell identity and validate this by performing MSC knock down in human embryonic myoblasts, which results in significant upregulation of late myogenic genes. Spatially mapping the cell types of the limb across a range of gestational ages demonstrates a clear anatomical segregation between genes linked to brachydactyly and polysyndactyly, and uncovers two transcriptionally and spatially distinct populations of the progress zone, which we term “outer” and “transitional” layers. The latter exhibits a transcriptomic profile similar to that of the chondrocyte lineage, but lacking the key chondrogenic transcription factors SOX5,6 & 9. Finally, we perform scRNA-seq on murine embryonic limbs to facilitate cross-species developmental comparison at single-cell resolution, finding substantial homology between the two species.

[1]  M. Gerstung,et al.  Cell2location maps fine-grained cell types in spatial transcriptomics , 2022, Nature Biotechnology.

[2]  Mingchao Li,et al.  Single-Cell Regulatory Network Inference and Clustering Identifies Cell-Type Specific Expression Pattern of Transcription Factors in Mouse Sciatic Nerve , 2021, Frontiers in Cellular Neuroscience.

[3]  Monika S. Kowalczyk,et al.  Blood and immune development in human fetal bone marrow and Down syndrome , 2021, Nature.

[4]  J. Drouin,et al.  Control of mouse limb initiation and antero-posterior patterning by Meis transcription factors , 2021, Nature Communications.

[5]  Deanne M. Taylor,et al.  A roadmap for the Human Developmental Cell Atlas , 2021, Nature.

[6]  J. Marioni,et al.  Cells of the human intestinal tract mapped across space and time , 2021, Nature.

[7]  S. Teichmann,et al.  MultiMAP: dimensionality reduction and integration of multimodal data , 2021, Genome Biology.

[8]  P. Tschopp,et al.  Assessing evolutionary and developmental transcriptome dynamics in homologous cell types , 2021, bioRxiv.

[9]  S. Mundlos,et al.  Non-coding deletions identify Maenli lncRNA as a limb-specific En1 regulator , 2021, Nature.

[10]  Lia S. Campos,et al.  Mapping the temporal and spatial dynamics of the human endometrium in vivo and in vitro , 2021, Nature Genetics.

[11]  M. Hawkins,et al.  Latent developmental potential to form limb-like skeletal structures in zebrafish , 2018, Cell.

[12]  Lan Zhang,et al.  Apical ectodermal ridge regulates three principal axes of the developing limb , 2020, Journal of Zhejiang University-SCIENCE B.

[13]  R. Haltiwanger,et al.  Canonical Notch ligands and Fringes have distinct effects on NOTCH1 and NOTCH2 , 2020, The Journal of Biological Chemistry.

[14]  A. Zannettino,et al.  HOPX regulates bone marrow-derived mesenchymal stromal cell fate determination via suppression of adipogenic gene pathways , 2020, Scientific Reports.

[15]  Matthew D. Young,et al.  Tumor to normal single cell mRNA comparisons reveal a pan-neuroblastoma cancer cell , 2020, bioRxiv.

[16]  Yvan Saeys,et al.  A scalable SCENIC workflow for single-cell gene regulatory network analysis , 2020, Nature Protocols.

[17]  Yu Zhang,et al.  The changing mouse embryo transcriptome at whole tissue and single-cell resolution , 2020, Nature.

[18]  H. Akiyama,et al.  Induced pluripotent stem cell-derived tenocyte-like cells promote the regeneration of injured tendons in mice , 2020, Scientific Reports.

[19]  Mirjana Efremova,et al.  CellPhoneDB: inferring cell–cell communication from combined expression of multi-subunit ligand–receptor complexes , 2020, Nature Protocols.

[20]  T. Korff,et al.  Loss of the serine protease HTRA1 impairs smooth muscle cells maturation , 2019, Scientific Reports.

[21]  R. Lovering,et al.  Keratin 18 is an integral part of the intermediate filament network in murine skeletal muscle. , 2019, American journal of physiology. Cell physiology.

[22]  A. Uezumi,et al.  Sustained expression of HeyL is critical for the proliferation of muscle stem cells in overloaded muscle , 2019, eLife.

[23]  David McDonald,et al.  Decoding human fetal liver haematopoiesis , 2019, Nature.

[24]  F. Guilak,et al.  Single cell RNA-sequencing reveals cellular heterogeneity and trajectories of lineage specification during murine embryonic limb development , 2019, bioRxiv.

[25]  Oren Parnas,et al.  A single-cell transcriptomic atlas of the developing chicken limb , 2019, bioRxiv.

[26]  Michael J. T. Stubbington,et al.  Single-cell reconstruction of the early maternal–fetal interface in humans , 2018, Nature.

[27]  Hongbo Zhang,et al.  Muscle Stem Cell Immunostaining , 2018, Current protocols in mouse biology.

[28]  Joerg M. Buescher,et al.  A metabolic interplay coordinated by HLX regulates myeloid differentiation and AML through partly overlapping pathways , 2018, Nature Communications.

[29]  S. Ramsey,et al.  FACS-Seq analysis of Pax3-derived cells identifies non-myogenic lineages in the embryonic forelimb , 2018, Scientific reports.

[30]  Jianyun Yan,et al.  Smad4 deficiency impairs chondrocyte hypertrophy via the Runx2 transcription factor in mouse skeletal development , 2018, The Journal of Biological Chemistry.

[31]  Dennis Wolf,et al.  Single-Cell RNA-Seq Reveals the Transcriptional Landscape and Heterogeneity of Aortic Macrophages in Murine Atherosclerosis , 2018, Circulation research.

[32]  Xu Wang,et al.  TFAP2C promotes stemness and chemotherapeutic resistance in colorectal cancer via inactivating hippo signaling pathway , 2018, Journal of experimental & clinical cancer research : CR.

[33]  Megan K. Cox,et al.  Skeletal Characterization of the Fgfr3 Mouse Model of Achondroplasia Using Micro-CT and MRI Volumetric Imaging , 2018, Scientific Reports.

[34]  C. Deng,et al.  Fibroblast Growth Factor Receptor 2 (FGFR2) Mutation Related Syndromic Craniosynostosis , 2017, International journal of biological sciences.

[35]  J. Aerts,et al.  SCENIC: Single-cell regulatory network inference and clustering , 2017, Nature Methods.

[36]  D. Franco,et al.  Pitx2 in Embryonic and Adult Myogenesis , 2017, Front. Cell Dev. Biol..

[37]  Christopher. Simons,et al.  Machine learning with Python , 2017 .

[38]  M. Goldring,et al.  ELF3 modulates type II collagen gene (COL2A1) transcription in chondrocytes by inhibiting SOX9-CBP/p300-driven histone acetyltransferase activity , 2017, Connective tissue research.

[39]  C. Cabrele,et al.  The Id-protein family in developmental and cancer-associated pathways , 2017, Cell Communication and Signaling.

[40]  M. Kmita,et al.  Evolution of Hoxa11 regulation in vertebrates is linked to the pentadactyl state , 2016, Nature.

[41]  R. Atit,et al.  Defining the identity of mouse embryonic dermal fibroblasts , 2016, Genesis.

[42]  Bahubali D. Gane,et al.  Split-hand/feet malformation: A rare syndrome , 2016, Journal of family medicine and primary care.

[43]  J. Lopez-Rios The many lives of SHH in limb development and evolution. , 2016, Seminars in cell & developmental biology.

[44]  A. Zuniga Next generation limb development and evolution: old questions, new perspectives , 2015, Development.

[45]  S. Rho,et al.  WIF1 can effectively co-regulate pro-apoptotic activity through the combination with DKK1. , 2014, Cellular signalling.

[46]  R. Salonen,et al.  Periderm prevents pathological epithelial adhesions during embryogenesis. , 2014, The Journal of clinical investigation.

[47]  Kuo-I Lin,et al.  Transcription Factor ABF-1 Suppresses Plasma Cell Differentiation but Facilitates Memory B Cell Formation , 2014, The Journal of Immunology.

[48]  J. Lecardonnel,et al.  Chondrocytes play a major role in the stimulation of bone growth by thyroid hormone. , 2014, Endocrinology.

[49]  Nan Su,et al.  Role of FGF/FGFR signaling in skeletal development and homeostasis: learning from mouse models , 2014, Bone Research.

[50]  J. Chimal-Monroy,et al.  Molecular Control of Interdigital Cell Death and Cell Differentiation by Retinoic Acid during Digit Development , 2014 .

[51]  Sevan Hopyan,et al.  Formation of proximal and anterior limb skeleton requires early function of Irx3 and Irx5 and is negatively regulated by Shh signaling. , 2014, Developmental cell.

[52]  Y. Kawakami,et al.  Distinct populations within Isl1 lineages contribute to appendicular and facial skeletogenesis through the β-catenin pathway. , 2014, Developmental biology.

[53]  P. Rigby,et al.  Gene regulatory networks and transcriptional mechanisms that control myogenesis. , 2014, Developmental cell.

[54]  Julie Moss,et al.  EMAGE mouse embryo spatial gene expression database: 2014 update , 2013, Nucleic Acids Res..

[55]  Abraham P. Fong,et al.  Genome-wide binding of the basic helix-loop-helix myogenic inhibitor musculin has substantial overlap with MyoD: implications for buffering activity , 2013, Skeletal Muscle.

[56]  A. Jamsheer,et al.  Split-hand/foot malformation - molecular cause and implications in genetic counseling , 2013, Journal of Applied Genetics.

[57]  Fátima Sánchez-Cabo,et al.  Analysis of the DNA-binding profile and function of TALE homeoproteins reveals their specialization and specific interactions with Hox genes/proteins. , 2013, Cell reports.

[58]  Sung-Jan Lin,et al.  Inducible deletion of the Blimp-1 gene in adult epidermis causes granulocyte-dominated chronic skin inflammation in mice , 2013, Proceedings of the National Academy of Sciences.

[59]  J. Chimal-Monroy,et al.  Irx1 and Irx2 Are Coordinately Expressed and Regulated by Retinoic Acid, TGFβ and FGF Signaling during Chick Hindlimb Development , 2013, PloS one.

[60]  H. Ouyang,et al.  Force and scleraxis synergistically promote the commitment of human ES cells derived MSCs to tenocytes , 2012, Scientific Reports.

[61]  Johannes E. Schindelin,et al.  Fiji: an open-source platform for biological-image analysis , 2012, Nature Methods.

[62]  Allison M. Haaning,et al.  Preaxial polydactyly caused by Gli3 haploinsufficiency is rescued by Zic3 loss of function in mice. , 2012, Human molecular genetics.

[63]  H. Yoshikawa,et al.  Nkx3.2 Promotes Primary Chondrogenic Differentiation by Upregulating Col2a1 Transcription , 2012, PloS one.

[64]  G. Narkis,et al.  Isl1 and Ldb Co‐regulators of transcription are essential early determinants of mouse limb development , 2012, Developmental dynamics : an official publication of the American Association of Anatomists.

[65]  R. Behringer,et al.  Regionalized Twist1 activity in the forelimb bud drives the morphogenesis of the proximal and preaxial skeleton. , 2012, Developmental biology.

[66]  Jianbo Wang,et al.  Hierarchical Interactions of Homeodomain and Forkhead Transcription Factors in Regulating Odontogenic Gene Expression* , 2011, The Journal of Biological Chemistry.

[67]  M. Khokha,et al.  Cooperative activity of noggin and gremlin 1 in axial skeleton development , 2011, Development.

[68]  J. Cobb,et al.  Shox2 function couples neural, muscular and skeletal development in the proximal forelimb. , 2011, Developmental biology.

[69]  Gaël Varoquaux,et al.  Scikit-learn: Machine Learning in Python , 2011, J. Mach. Learn. Res..

[70]  P. Scambler,et al.  SHOX interacts with the chondrogenic transcription factors SOX5 and SOX6 to activate the aggrecan enhancer. , 2011, Human molecular genetics.

[71]  P. Kondaiah,et al.  Gene expression profile of epithelial cells and mesenchymal cells derived from limbal explant culture , 2010, Molecular vision.

[72]  J. David,et al.  Glucocorticoids suppress bone formation by attenuating osteoblast differentiation via the monomeric glucocorticoid receptor. , 2010, Cell metabolism.

[73]  R. Paro,et al.  Distinct Roles of Hand2 in Initiating Polarity and Posterior Shh Expression during the Onset of Mouse Limb Bud Development , 2010, PLoS genetics.

[74]  S. Gay,et al.  The transcription factor Fra-2 regulates the production of extracellular matrix in systemic sclerosis. , 2010, Arthritis and rheumatism.

[75]  J. Montero,et al.  Expression of Id2 in the developing limb is associated with zones of active BMP signaling and marks the regions of growth and differentiation of the developing digits. , 2009, The International journal of developmental biology.

[76]  S. Mundlos,et al.  A mutation in Ihh that causes digit abnormalities alters its signalling capacity and range , 2009, Nature.

[77]  C. Tabin,et al.  Fgf-dependent Etv4/5 activity is required for posterior restriction of Sonic Hedgehog and promoting outgrowth of the vertebrate limb. , 2009, Developmental cell.

[78]  H. Westphal,et al.  LIM homeobox transcription factors integrate signaling events that control three-dimensional limb patterning and growth , 2009, Development.

[79]  Naoko Koyano-Nakagawa,et al.  Sall genes regulate region-specific morphogenesis in the mouse limb by modulating Hox activities , 2009, Development.

[80]  Matthew Towers,et al.  Growing models of vertebrate limb development , 2009, Development.

[81]  T. Kadesch,et al.  Inhibition of myogenesis by Notch: Evidence for multiple pathways , 2009, Journal of cellular physiology.

[82]  G. Weinmaster,et al.  The many facets of Notch ligands , 2008, Oncogene.

[83]  C. Tickle,et al.  Integration of growth and specification in chick wing digit-patterning , 2008, Nature.

[84]  B. Wainwright,et al.  Expression of the NET family member Zfp503 is regulated by hedgehog and BMP signaling in the limb , 2008, Developmental dynamics : an official publication of the American Association of Anatomists.

[85]  C. Ponting,et al.  Mutations in BMP4 cause eye, brain, and digit developmental anomalies: overlap between the BMP4 and hedgehog signaling pathways. , 2008, American journal of human genetics.

[86]  B. Mankoo,et al.  A comparative analysis of Meox1 and Meox2 in the developing somites and limbs of the chick embryo. , 2007, The International journal of developmental biology.

[87]  Denis Duboule,et al.  The role of Hox genes during vertebrate limb development. , 2007, Current opinion in genetics & development.

[88]  Lewis Wolpert,et al.  Rethinking the proximodistal axis of the vertebrate limb in the molecular era. , 2007, Genes & development.

[89]  P. Ernfors,et al.  Emergence of the sensory nervous system as defined by Foxs1 expression. , 2007, Differentiation; research in biological diversity.

[90]  L. Piccio,et al.  Cutting Edge: TREM-2 Attenuates Macrophage Activation1 , 2006, The Journal of Immunology.

[91]  J. Cidlowski,et al.  Glucocorticoid receptor isoforms generate transcription specificity. , 2006, Trends in cell biology.

[92]  R. Klein,et al.  TEF-1 and C/EBPβ are major p38α MAPK-regulated transcription factors in proliferating cardiomyocytes , 2006 .

[93]  K. Anderson,et al.  Signaling from Smo to Ci/Gli: conservation and divergence of Hedgehog pathways from Drosophila to vertebrates , 2006, Development.

[94]  D. Duboule,et al.  Control of Hoxd genes' collinearity during early limb development. , 2006, Developmental cell.

[95]  A. Berdal,et al.  Expression pattern of Dlx3 during cell differentiation in mineralized tissues. , 2005, Bone.

[96]  M. Soares,et al.  Pax9 and Jagged1 act downstream of Gli3 in vertebrate limb development , 2005, Mechanisms of Development.

[97]  Jamie M. Verheyden,et al.  Conditional inactivation of Fgfr1 in mouse defines its role in limb bud establishment, outgrowth and digit patterning , 2005, Development.

[98]  J. Korving,et al.  Function and regulation of Alx4 in limb development: complex genetic interactions with Gli3 and Shh. , 2005, Developmental biology.

[99]  M. Massett,et al.  GIT1 Mediates HDAC5 Activation by Angiotensin II in Vascular Smooth Muscle Cells , 2008, Arteriosclerosis, thrombosis, and vascular biology.

[100]  Wei Liu,et al.  Bmp4 in limb bud mesoderm regulates digit pattern by controlling AER development. , 2004, Developmental biology.

[101]  D. Simon,et al.  Transcription Factor CHF1/Hey2 Regulates Neointimal Formation In Vivo and Vascular Smooth Muscle Proliferation and Migration In Vitro , 2004, Arteriosclerosis, thrombosis, and vascular biology.

[102]  M. Justice,et al.  A twist code determines the onset of osteoblast differentiation. , 2004, Developmental cell.

[103]  B. Christ,et al.  Experimental analysis of the origin of the wing musculature in avian embryos , 1977, Anatomy and Embryology.

[104]  J. Charron,et al.  Role of Plk2 (Snk) in Mouse Development and Cell Proliferation , 2003, Molecular and Cellular Biology.

[105]  A. Wilkie Why study human limb malformations? , 2003, Journal of anatomy.

[106]  M. Buckingham,et al.  The formation of skeletal muscle: from somite to limb , 2003, Journal of anatomy.

[107]  J. Drouin,et al.  Pitx1 and Pitx2 are required for development of hindlimb buds , 2003, Development.

[108]  C. Marcelle,et al.  FGFR4 signaling is a necessary step in limb muscle differentiation. , 2002, Development.

[109]  G. Martin,et al.  Functions of FGF signalling from the apical ectodermal ridge in limb development , 2002, Nature.

[110]  G. Merlo,et al.  Mouse model of split hand/foot malformation type I , 2002, Genesis.

[111]  D. Prockop,et al.  Targeted disruption of Col11a2 produces a mild cartilage phenotype in transgenic mice: Comparison with the human disorder otospondylomegaepiphyseal dysplasia (OSMED) , 2001, Developmental dynamics : an official publication of the American Association of Anatomists.

[112]  P. Gruss,et al.  Expression pattern of Irx1 and Irx2 during mouse digit development , 2001, Mechanisms of Development.

[113]  J. Deng,et al.  Haploinsufficiency of Sox9 results in defective cartilage primordia and premature skeletal mineralization , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[114]  C. Valencia,et al.  Differential tissue growth and patterns of cell death in mouse limb autopod morphogenesis , 2001, Developmental dynamics : an official publication of the American Association of Anatomists.

[115]  M. Vekemans,et al.  JAGGED1 Gene Expression During Human Embryogenesis Elucidates the Wide Phenotypic Spectrum of Alagille Syndrome , 2000, Hepatology.

[116]  I. S. Kathiriya,et al.  Role of dHAND in the anterior-posterior polarization of the limb bud: implications for the Sonic hedgehog pathway. , 2000, Development.

[117]  A. Munnich,et al.  Saethre-Chotzen mutations cause TWIST protein degradation or impaired nuclear location. , 2000, Human molecular genetics.

[118]  E. Zackai,et al.  Mutations in the human TWIST gene , 2000, Human mutation.

[119]  H. Brunner,et al.  Autosomal recessive disorder otospondylomegaepiphyseal dysplasia is associated with loss-of-function mutations in the COL11A2 gene. , 2000, American journal of human genetics.

[120]  K. Lyons,et al.  The type I BMP receptor BMPRIB is required for chondrogenesis in the mouse limb. , 2000, Development.

[121]  Andrew P. McMahon,et al.  Signal relay by BMP antagonism controls the SHH/FGF4 feedback loop in vertebrate limb buds , 1999, Nature.

[122]  P. Beachy,et al.  Sonic hedgehog protein signals not as a hydrolytic enzyme but as an apparent ligand for patched. , 1999, Proceedings of the National Academy of Sciences of the United States of America.

[123]  Steven J. M. Jones,et al.  A Wnt5a pathway underlies outgrowth of multiple structures in the vertebrate embryo. , 1999, Development.

[124]  E. Olson,et al.  MyoR: a muscle-restricted basic helix-loop-helix transcription factor that antagonizes the actions of MyoD. , 1999, Proceedings of the National Academy of Sciences of the United States of America.

[125]  J. V. van Dongen,et al.  Characterization of ABF-1, a Novel Basic Helix-Loop-Helix Transcription Factor Expressed in Activated B Lymphocytes , 1998, Molecular and Cellular Biology.

[126]  A. McMahon,et al.  Noggin, cartilage morphogenesis, and joint formation in the mammalian skeleton. , 1998, Science.

[127]  Raphael Kopan,et al.  Notch-1 signalling requires ligand-induced proteolytic release of intracellular domain , 1998, Nature.

[128]  R. Beddington,et al.  Msg1 and Mrg1, founding members of a gene family, show distinct patterns of gene expression during mouse embryogenesis , 1998, Mechanisms of Development.

[129]  P. Leder,et al.  Fibroblast growth factor receptor 2 (FGFR2)-mediated reciprocal regulation loop between FGF8 and FGF10 is essential for limb induction. , 1998, Development.

[130]  A. Rosenthal,et al.  Pseudoautosomal deletions encompassing a novel homeobox gene cause growth failure in idiopathic short stature and Turner syndrome , 1997, Nature Genetics.

[131]  C. Tabin,et al.  Sonic hedgehog mediates the polarizing activity of the ZPA , 1993, Cell.

[132]  H. Kondoh,et al.  Defects of embryonic organogenesis resulting from targeted disruption of the N-myc gene in the mouse. , 1993, Development.

[133]  Denis Duboule,et al.  Coordinate expression of the murine Hox-5 complex homoeobox-containing genes during limb pattern formation , 1989, Nature.

[134]  M. Ontell,et al.  The organogenesis of murine striated muscle: a cytoarchitectural study. , 1984, The American journal of anatomy.

[135]  J. W. Saunders The proximo-distal sequence of origin of the parts of the chick wing and the role of the ectoderm. , 1948, The Journal of experimental zoology.