A single-cell transcriptomic atlas of the developing chicken limb

The chicken embryo brings new insights into the evolutionary role of WFDC1 during amniote development

WFDC1 encodes an extracellular matrix protein involved in cell proliferation, migration, and epithelial-mesenchymal transition in disease conditions. Despite this, Wfdc1-null mice display no discernible malformations while cattle bearing a WFDC1 mutation present multiple ocular defects, leaving the role of WFDC1 during embryonic development unclear. To address this, we used the chicken embryo as a model to investigate WFDC1 developmental roles in amniotes. We performed a comparative expression analysis during chicken and mouse development, which revealed expression in ectodermal and mesodermal derivatives, with both conserved and species-specific domains. Conserved expression was observed in the eye, otic vesicle, central and peripheral nervous systems, and neural crest cells. Chicken-specific expression was identified in mesodermal structures, including the notochord, limbs and heart. However, even in the conserved sites like the eyes, WFDC1 localizes to different retinal layers, indicating potential divergence roles in retinal development and function across species. In contrast, WFDC1 expression in the limb buds is specific to chicken, encompassing the distal mesenchyme, interdigital membranes, and blastemas. Functional enrichment analysis links WFDC1 to limb patterning, morphogenesis, and Wnt signaling. The species-specific differences likely stem from evolutionary changes in gene regulation, supported by differences in proximal cis-regulatory elements of the WFDC1 loci between chicken and mouse. The complexity of WFDC1 expression in the chicken embryo, along with WFDC1 regulatory conservation within birds, indicates that this gene may play specific roles in avian development, possibly contributing to features specific to this lineage. Future studies using the chicken model will be valuable in further uncovering the specific roles of WFDC1.

Unravelling early hematoendothelial development through the chick model: Insights and future perspectives

The chicken embryo has been an important model in advancing our understanding of early hematoendothelial development, particularly in the formation of hematopoietic stem cells (HSCs) and the endothelial-to-hematopoietic transition (EHT). The accessibility and ease of manipulation of chicken embryos have made them an invaluable tool for researching development of blood and endothelial cells. Early research using this model provided pivotal insights, demonstrating that intra-embryonic regions, such as the dorsal aorta (DA), are primary sources of HSCs, rather than the yolk sac (YS), as previously believed. The identification of intra-aortic hematopoietic clusters (IAHCs) and the process of EHT in the chicken embryo laid the foundation for similar discoveries in other vertebrate species, including mice and zebrafish. Recent advances in genetic tools, such as transgenic chickens expressing fluorescent proteins, have further enhanced the precision of cell lineage tracing and real-time imaging of dynamic cellular processes. This review highlights both historical contributions and contemporary advancements facilitated by the chicken model, underscoring its continued relevance in developmental biology. By examining key findings and methodological innovations, we aim to demonstrate the importance of the chicken embryo as a model system for understanding hematoendothelial development and its potential for informing therapeutic applications in regenerative medicine and blood disorders. Finally, we will underscore potential applications of the chicken model for comparative and omics-level studies in conjunction with other model systems and what future directions lie ahead.

Characterisation of the Avascular Mesenchyme during Digit Outgrowth

The avascular mesenchyme at the tip of the developing digit contributes to digit outgrowth and patterning, however, it has been poorly characterised. Using newly developed fate mapping approaches, tissue manipulation and single-cell mRNA sequencing data, we explore the transcriptional nature and developmental potential of this tissue. We find that the avascular mesenchyme is essential to normal segmental patterning of the digit and has a distinct transcriptional identity. In addition, we uncover an unexpected relationship between the unspecified tissue of the avascular mesenchyme and the committed phalanx forming region, which controls patterning, but not outgrowth of the digit. This multifaceted approach provides insights into the mechanics and genetic pathways that regulate digit outgrowth and patterning.

ZIC3 shapes digit morphogenesis in avian: Research Note: Zic3 shapes anterior digit morphogenesis in avian diversity

Through the evolutionary diversification of birds, a variety of digit morphologies have evolved as adaptations to diverse environments, with the regression of the anterior digits being a remarkable phenomenon. Transcriptomic analysis of chicken limb development has revealed molecular signatures in the anterior digits, among which Zic family member 3 (Zic3) stands out as a unique highly expressed transcription factor. However, the function of Zic3 in digit development in birds remains unclear. In this study, we investigated the expression pattern of Zic3 across five phylogenetically diverse avian species. Our analysis revealed a consistent Zic3 expression pattern in species from Neoaves and Palaeognathae, except for the ostrich. Ostrich hindlimbs exhibited increased Zic3 expression and an extended expression range from the first toe to the second toe. By locally expanding the expression domain of Zic3 in chicken anterior hindlimb buds we obtained a phenotype similar to that of early ostrich feet, characterized by significant shortening of digits I and II. We further performed mRNA Sequencing (mRNA-Seq) of the Zic3-overexpressed autopods and found that Zic3 inhibited skeletal development through multiple pathways, including the Wnt signaling pathway, ECM-receptor interaction and Focal adhesion. Our results reveal a pivotal role for Zic3 in anterior digits and identify its downstream regulatory mediators in birds.

Distinct gene regulatory dynamics drive skeletogenic cell fate convergence during vertebrate embryogenesis

Cell type repertoires have expanded extensively in metazoan animals, with some clade-specific cells being crucial to evolutionary success. A prime example are the skeletogenic cells of vertebrates. Depending on anatomical location, these cells originate from three different precursor lineages, yet they converge developmentally towards similar cellular phenotypes. Furthermore, their 'skeletogenic competency' arose at distinct evolutionary timepoints, thus questioning to what extent different skeletal body parts rely on truly homologous cell types. Here, we investigate how lineage-specific molecular properties are integrated at the gene regulatory level, to allow for skeletogenic cell fate convergence. Using single-cell functional genomics, we find that distinct transcription factor profiles are inherited from the three precursor states and incorporated at lineage-specific enhancer elements. This lineage-specific regulatory logic suggests that these regionalized skeletogenic cells are distinct cell types, rendering them amenable to individualized selection, to define adaptive morphologies and biomaterial properties in different parts of the vertebrate skeleton.

Single-cell RNA transcriptome uncovers distinct developmental trajectories in the embryonic skeletal muscle of Daheng broiler and Tibetan chicken

Different chicken breeds exhibit distinct muscle phenotypes resulting from selective breeding, but little is known about the molecular mechanisms responsible for this phenotypic difference. Skeletal muscle is composed of a large number of heterogeneous cell populations. Differences in differentiation and interaction of cell populations play a key role in the difference of skeletal muscle phenotype. In the current study, we performed a single-cell RNA sequencing (scRNA-seq) on the leg muscle of Daheng broiler (DH, cultivated breed) and Tibetan chicken (TC, native breed) at embryonic (E) 10, E14 and E18. A comprehensive cell atlas of embryonic chicken skeletal muscle, consisting of 29,579 high-quality cells representing 6 distinct cell types was built. The differentiation trajectory of Myoblasts and fibro-adipogenic progenitors (FAPs) was constructed through pseudotemporal trajectory analysis. Our results revealed the different developmental trajectories and dynamic gene expression profiles in 3 subtypes of myoblasts and 5 FAPs subtypes of the two chicken breeds. Tibetan chicken showed earlier embryonic myogenesis and less myoblasts compared with Daheng broiler. By comparing the switch status and switch time of genes in the two breeds, SNRPG,SNRPE,EIF4EBP1 and HSP90AB1 were considered as potentially critical genes for embryonic myogenesis, and the genes MYOG,MYBPH,APOA1, and MGP played dominant roles in the embryonic adipogenesis. Intercellular interaction networks showed that strong and complex intercellular communication was contained during embryonic skeletal muscle growth and development. These findings revealed the differences of molecular mechanisms in the skeletal muscle development between TC and DH chickens. Our data provide a better understanding of skeletal muscle developmental differences between cultivated and native breeds and valuable information for genetic breeding of chicken.

Single-cell sequencing provides clues about the developmental genetic basis of evolutionary adaptations in syngnathid fishes

Seahorses, pipefishes, and seadragons are fishes from the family Syngnathidae that have evolved extraordinary traits including male pregnancy, elongated snouts, loss of teeth, and dermal bony armor. The developmental genetic and cellular changes that led to the evolution of these traits are largely unknown. Recent syngnathid genome assemblies revealed suggestive gene content differences and provided the opportunity for detailed genetic analyses. We created a single-cell RNA sequencing atlas of Gulf pipefish embryos to understand the developmental basis of four traits: derived head shape, toothlessness, dermal armor, and male pregnancy. We completed marker gene analyses, built genetic networks, and examined the spatial expression of select genes. We identified osteochondrogenic mesenchymal cells in the elongating face that express regulatory genes bmp4, sfrp1a, and prdm16. We found no evidence for tooth primordia cells, and we observed re-deployment of osteoblast genetic networks in developing dermal armor. Finally, we found that epidermal cells expressed nutrient processing and environmental sensing genes, potentially relevant for the brooding environment. The examined pipefish evolutionary innovations are composed of recognizable cell types, suggesting that derived features originate from changes within existing gene networks. Future work addressing syngnathid gene networks across multiple stages and species is essential for understanding how the novelties of these fish evolved.

Innovative Insights into Single-Cell Technologies and Multi-Omics Integration in Livestock and Poultry

In recent years, single-cell RNA sequencing (scRNA-seq) has marked significant strides in livestock and poultry research, especially when integrated with multi-omics approaches. These advancements provide a nuanced view into complex regulatory networks and cellular dynamics. This review outlines the application of scRNA-seq in key species, including poultry, swine, and ruminants, with a focus on outcomes related to cellular heterogeneity, developmental biology, and reproductive mechanisms. We emphasize the synergistic power of combining scRNA-seq with epigenomic, proteomic, and spatial transcriptomic data, enhancing molecular breeding precision, optimizing health management strategies, and refining production traits in livestock and poultry. The integration of these technologies offers a multidimensional approach that not only broadens the scope of data analysis but also provides actionable insights for improving animal health and productivity.

Single-cell transcriptional profiling reveals cell type-specific responses to duck reovirus infection in the Bursa of Fabricius of Cairna moschata

Duck reovirus (DRV) is a universal waterfowl virus that causes significant economic losses in the duck industry. However, the role of the host innate immune response of the Bursa of Fabricius to DRV infection is largely unknown. In the present study, we constructed a single-cell resolution transcriptomic atlas of the Bursa of Fabricius of Cairna moschata after infection with HN10 (a novel DRV). Ten cell-type marker genes were used to annotate the cell type, indicating a high degree of cell heterogeneity in the Bursa of Fabricius. Most of the innate and adaptive immune system-related genes were highly expressed in T cells, B cells, neutrophils, macrophages, and DCs. In the Bursa of Fabricius, the proportions of DCs and macrophages were largely increased by HN10 infection at 14 d, suggesting that DCs and macrophages play important roles in the long-term viral response. Notably, a number of innate and adaptive immune system-related genes were highly expressed at 24 h after HN10 infection, indicating that the Bursa of Fabricius has a very strong immune function even in the early developmental stage. In the immune system, the NOD-like receptor signaling pathway and RIG-I-like receptor signaling pathway were significantly activated at the early stage of HN10 infection, while the Toll-like receptor signaling pathway was significantly activated at the late stage. Enrichment analysis suggested that different immune signaling pathways play roles in specific developmental stages. Our data provide an opportunity to reveal the immune response to DRV infection at the single-cell level.

Limb connective tissue is organized in a continuum of promiscuous fibroblast identities during development

Connective tissue (CT), which includes tendon and muscle CT, plays critical roles in development, in particular as positional cue provider. Nonetheless, our understanding of fibroblast developmental programs is hampered because fibroblasts are highly heterogeneous and poorly characterized. Combining single-cell RNA-sequencing-based strategies including trajectory inference and in situ hybridization analyses, we address the diversity of fibroblasts and their developmental trajectories during chicken limb fetal development. We show that fibroblasts switch from a positional information to a lineage diversification program at the fetal period onset. Muscle CT and tendon are composed of several fibroblast populations that emerge asynchronously. Once the final muscle pattern is set, transcriptionally close populations are found in neighboring locations in limbs, prefiguring the adult fibroblast layers. We propose that the limb CT is organized in a continuum of promiscuous fibroblast identities, allowing for the robust and efficient connection of muscle to bone and skin.

Mechanisms of regeneration: to what extent do they recapitulate development?

One of the enduring debates in regeneration biology is the degree to which regeneration mirrors development. Recent technical advances, such as single-cell transcriptomics and the broad applicability of CRISPR systems, coupled with new model organisms in research, have led to the exploration of this longstanding concept from a broader perspective. In this Review, I outline the historical parallels between development and regeneration before focusing on recent research that highlights how dissecting the divergence between these processes can uncover previously unreported biological mechanisms. Finally, I discuss how these advances position regeneration as a more dynamic and variable process with expanded possibilities for morphogenesis compared with development. Collectively, these insights into mechanisms that orchestrate morphogenesis may reshape our understanding of the evolution of regeneration, reveal hidden biology activated by injury, and offer non-developmental strategies for restoring lost or damaged organs and tissues.

Single-cell transcriptomic reveals a cell atlas and diversity of chicken amygdala responded to social hierarchy

Amygdala serves as a highly cellular, heterogeneous brain region containing excitatory and inhibitory neurons and is involved in the dopamine and serotoninergic neuron systems. An increasing number of studies have revealed the underpinned mechanism mediating social hierarchy in mammal and vertebrate, however, there are rare studies conducted on how amygdala on social hierarchy in poultry. In this study, we conducted food competition tests and determined the social hierarchy of the rooster. We performed cross-species analysis with mammalian amygdala, and found that cell types of human and rhesus monkeys were more closely related and that of chickens were more distant. We identified 26 clusters and divided them into 10 main clusters, of which GABAergic and glutamatergic neurons were associated with social behaviors. In conclusion, our results provide to serve the developmental studies of the amygdala neuron system and new insights into the underpinned mechanism of social hierarchy in roosters.

Self-organized BMP signaling dynamics underlie the development and evolution of digit segmentation patterns in birds and mammals

Repeating patterns of synovial joints are a highly conserved feature of articulated digits, with variations in joint number and location resulting in diverse digit morphologies and limb functions across the tetrapod clade. During the development of the amniote limb, joints form iteratively within the growing digit ray, as a population of distal progenitors alternately specifies joint and phalanx cell fates to segment the digit into distinct elements. While numerous molecular pathways have been implicated in this fate choice, it remains unclear how they give rise to a repeating pattern. Here, using single-cell RNA sequencing and spatial gene expression profiling, we investigate the transcriptional dynamics of interphalangeal joint specification in vivo. Combined with mathematical modeling, we predict that interactions within the BMP signaling pathway-between the ligand GDF5, the inhibitor NOGGIN, and the intracellular effector pSMAD-result in a self-organizing Turing system that forms periodic joint patterns. Our model is able to recapitulate the spatiotemporal gene expression dynamics observed in vivo, as well as phenocopy digit malformations caused by BMP pathway perturbations. By contrasting in silico simulations with in vivo morphometrics of two morphologically distinct digits, we show how changes in signaling parameters and growth dynamics can result in variations in the size and number of phalanges. Together, our results reveal a self-organizing mechanism that underpins amniote digit segmentation and its evolvability and, more broadly, illustrate how Turing systems based on a single molecular pathway may generate complex repetitive patterns in a wide variety of organisms.

Multi-species atlas resolves an axolotl limb development and regeneration paradox

Humans and other tetrapods are considered to require apical-ectodermal-ridge (AER) cells for limb development, and AER-like cells are suggested to be re-formed to initiate limb regeneration. Paradoxically, the presence of AER in the axolotl, a primary model organism for regeneration, remains controversial. Here, by leveraging a single-cell transcriptomics-based multi-species atlas, composed of axolotl, human, mouse, chicken, and frog cells, we first establish that axolotls contain cells with AER characteristics. Further analyses and spatial transcriptomics reveal that axolotl limbs do not fully re-form AER cells during regeneration. Moreover, the axolotl mesoderm displays part of the AER machinery, revealing a program for limb (re)growth. These results clarify the debate about the axolotl AER and the extent to which the limb developmental program is recapitulated during regeneration.

Distinct patterning responses of wing and leg neuromuscular systems to different preaxial polydactylies

The tetrapod limb has long served as a paradigm to study vertebrate pattern formation and evolutionary diversification. The distal part of the limb, the so-called autopod, is of particular interest in this regard, given the numerous modifications in both its morphology and behavioral motor output. While the underlying alterations in skeletal form have received considerable attention, much less is known about the accompanying changes in the neuromuscular system. However, modifications in the skeleton need to be properly integrated with both muscle and nerve patterns, to result in a fully functional limb. This task is further complicated by the distinct embryonic origins of the three main tissue types involved-skeleton, muscles and nerves-and, accordingly, how they are patterned and connected with one another during development. To evaluate the degree of regulative crosstalk in this complex limb patterning process, here we analyze the developing limb neuromuscular system of Silkie breed chicken. These animals display a preaxial polydactyly, due to a polymorphism in the limb regulatory region of the Sonic Hedgehog gene. Using lightsheet microscopy and 3D-reconstructions, we investigate the neuromuscular patterns of extra digits in Silkie wings and legs, and compare our results to Retinoic Acid-induced polydactylies. Contrary to previous findings, Silkie autopod muscle patterns do not adjust to alterations in the underlying skeletal topology, while nerves show partial responsiveness. We discuss the implications of tissue-specific sensitivities to global limb patterning cues for our understanding of the evolution of novel forms and functions in the distal tetrapod limb.

A single-cell census of mouse limb development identifies complex spatiotemporal dynamics of skeleton formation

Limb development has long served as a model system for coordinated spatial patterning of progenitor cells. Here, we identify a population of naive limb progenitors and show that they differentiate progressively to form the skeleton in a complex, non-consecutive, three-dimensional pattern. Single-cell RNA sequencing of the developing mouse forelimb identified three progenitor states: naive, proximal, and autopodial, as well as Msx1 as a marker for the naive progenitors. In vivo lineage tracing confirmed this role and localized the naive progenitors to the outer margin of the limb, along the anterior-posterior axis. Sequential pulse-chase experiments showed that the progressive transition of Msx1⁺ naive progenitors into proximal and autopodial progenitors coincides with their differentiation to Sox9⁺ chondroprogenitors, which occurs along all the forming skeletal segments. Indeed, tracking the spatiotemporal sequence of differentiation showed that the skeleton forms progressively in a complex pattern. These findings suggest an alternative model for limb skeleton development.

Single-cell RNA sequencing revealed the liver heterogeneity between egg-laying duck and ceased-laying duck

BACKGROUND

In the late phase of production, ducks untimely cease laying, leading to a lower feed conversion. Liver plays a vital role in the synthesis and transport of yolk materials during egg formation in birds. However, the molecular mechanism of liver in ceased-laying duck is far from clear, higher resolution and deeper analysis is needed. Sing-cell RNA-sequencing of 10 × Genomics platform can help to map the liver single cell gene expression atlas of Shaoxing duck and provide new insights into the liver between egg-laying and ceased-laying ducks.

RESULTS

About 20,000 single cells were profiled and 22 clusters were identified. All the clusters were identified as 6 cell types. The dominant cell type is hepatocyte, accounted for about 60% of all the cells. Of note, the heterogeneity of cells between egg-laying duck and ceased-laying duck mainly occurred in hepatocytes. Cells of cluster 3 and 12 were the unique hepatocyte states of egg-laying ducks, while cells of cluster 0 and 15 were the unique hepatocyte states of ceased-laying ducks. The expression mode of yolk precursor transporters, lipid metabolizing enzymes and fibrinogens were different in hepatocytes between egg-laying duck and ceased-laying duck. APOV1, VTG2, VTG1, APOB, RBP, VTDB and SCD might be activated in egg-laying ducks, while APOA1, APOA4, APOC3, FGB and FGG might be activated in ceased-laying ducks.

CONCLUSIONS

Our study further proofs that APOV1 and APOB play key roles in egg production, rather than APOA1 and APOA4. It is also the first to detect a correlation between the higher expression of APOC3, FGB, FGG and ceased-laying in duck.

Advances in Single-Cell Sequencing Technology and Its Application in Poultry Science

Single-cell sequencing (SCS) uses a single cell as the research material and involves three dimensions: genes, phenotypes and cell biological mechanisms. This type of research can locate target cells, analyze the dynamic changes in the target cells and the relationships between the cells, and pinpoint the molecular mechanism of cell formation. Currently, a common problem faced by animal husbandry scientists is how to apply existing science and technology to promote the production of high-quality livestock and poultry products and to breed livestock for disease resistance; this is also a bottleneck for the sustainable development of animal husbandry. In recent years, although SCS technology has been successfully applied in the fields of medicine and bioscience, its application in poultry science has been rarely reported. With the sustainable development of science and technology and the poultry industry, SCS technology has great potential in the application of poultry science (or animal husbandry). Therefore, it is necessary to review the innovation of SCS technology and its application in poultry science. This article summarizes the current main technical methods of SCS and its application in poultry, which can provide potential references for its future applications in precision breeding, disease prevention and control, immunity, and cell identification.

Cell lineage specification and signalling pathway use during development of the lateral plate mesoderm and forelimb mesenchyme

The lateral plate mesoderm (LPM) is a transient tissue that produces a diverse range of differentiated structures, including the limbs. However, the molecular mechanisms that drive early LPM specification and development are poorly understood. In this study, we use single-cell transcriptomics to define the cell-fate decisions directing LPM specification, subdivision and early initiation of the forelimb mesenchyme in chicken embryos. We establish a transcriptional atlas and global cell-cell signalling interactions in progenitor, transitional and mature cell types throughout the developing forelimb field. During LPM subdivision, somatic and splanchnic LPM fate is achieved through activation of lineage-specific gene modules. During the earliest stages of limb initiation, we identify activation of TWIST1 in the somatic LPM as a putative driver of limb bud epithelial-to-mesenchymal transition. Furthermore, we define a new role for BMP signalling during early limb development, revealing that it is necessary for inducing a somatic LPM fate and initiation of limb outgrowth, potentially through activation of TBX5. Together, these findings provide new insights into the mechanisms underlying LPM development, somatic LPM fate choice and early initiation of the vertebrate limb.

Identification of candidate enhancers controlling the transcriptome during the formation of interphalangeal joints

The formation of the synovial joint begins with the visible emergence of a stripe of densely packed mesenchymal cells located between distal ends of the developing skeletal anlagen called the interzone. Recently the transcriptome of the early synovial joint was reported. Knowledge about enhancers would complement these data and lead to a better understanding of the control of gene transcription at the onset of joint development. Using ChIP-sequencing we have mapped the H3-signatures H3K27ac and H3K4me1 to locate regulatory elements specific for the interzone and adjacent phalange, respectively. This one-stage atlas of candidate enhancers (CEs) was used to map the association between these respective joint tissue specific CEs and biological processes. Subsequently, integrative analysis of transcriptomic data and CEs identified new putative regulatory elements of genes expressed in interzone (e.g., GDF5, BMP2 and DACT2) and phalange (e.g., MATN1, HAPLN1 and SNAI1). We also linked such CEs to genes known as crucial in synovial joint hypermobility and osteoarthritis, as well as phalange malformations. These analyses show that the CE atlas can serve as resource for identifying, and as starting point for experimentally validating, putative disease-causing genomic regulatory regions in patients with synovial joint dysfunctions and/or phalange disorders, and enhancer-controlled synovial joint and phalange formation.

Chicken Peripheral Blood Mononuclear Cells Response to Avian Leukosis Virus Subgroup J Infection Assessed by Single-Cell RNA Sequencing

Chicken peripheral blood mononuclear cells (PBMCs) exhibit wide-ranging cell types, but current understanding of their subclasses, immune cell classification, and function is limited and incomplete. Here we performed single-cell RNA sequencing (scRNA-seq) of PBMCs in Avian leukosis virus subgroup J (ALV-J) infected and control chickens at 21 days post infection (DPI) to determine chicken PBMCs subsets and their specific molecular and cellular characteristics. Eight cell populations and their potential marker genes were identified in PBMCs. T cell populations had the strongest response to (ALV-J) infection, based on the detection of the largest number of differentially expressed genes (DEGs), and could be further grouped into four subsets: activated CD4⁺ T cells, Th1-like cells, Th2-like cells, and cytotoxic CD8⁺ T cells. Furthermore, pseudotime analysis results suggested that chicken CD4⁺ T cells could potentially differentiate into Th1-like and Th2-like cells. Moreover, ALV-J infection activated CD4⁺ T cell was probably inclined to differentiate into Th1-like cells. Compared to the control PBMCs, ALV-J infection also had an obvious impact on PBMCs composition. B cells showed inconspicuous response and their numbers decreased in PBMCs from ALV-J infected chicken. Proportions of cytotoxic Th1-like cells and CD8⁺ T cells increased in the T cell population of PBMCs from ALV-J infected chicken, which were potentially key mitigating effectors against ALV-J infection. More importantly, our results provide a rich resource of gene expression profiles of chicken PBMCs subsets for a systems-level understanding of their function in homeostatic condition as well as in response to viral infection.

Towards Tabula Gallus

The Tabula Gallus is a proposed project that aims to create a map of every cell type in the chicken body and chick embryos. Chickens (Gallus gallus) are one of the most recognized model animals that recapitulate the development and physiology of mammals. The Tabula Gallus will generate a compendium of single-cell transcriptome data from Gallus gallus, characterize each cell type, and provide tools for the study of the biology of this species, similar to other ongoing cell atlas projects (Tabula Muris and Tabula Sapiens/Human Cell Atlas for mice and humans, respectively). The Tabula Gallus will potentially become an international collaboration between many researchers. This project will be useful for the basic scientific study of Gallus gallus and other birds (e.g., cell biology, molecular biology, developmental biology, neuroscience, physiology, oncology, virology, behavior, ecology, and evolution). It will eventually be beneficial for a better understanding of human health and diseases.

Advances in single-cell sequencing and its application to musculoskeletal system research

In recent years, single-cell sequencing (SCS) technologies have continued to advance with improved operating procedures and reduced cost, leading to increasing practical adoption among researchers. These emerging technologies have superior abilities to analyse cell heterogeneity at a single-cell level, which have elevated multi-omics research to a higher level. In some fields of research, application of SCS has enabled many valuable discoveries, and musculoskeletal system offers typical examples. This article reviews some major scientific issues and recent advances in musculoskeletal system. In addition, combined with SCS technologies, the research of cell or tissue heterogeneity in limb development and various musculoskeletal system clinical diseases also provides new possibilities for treatment strategies. Finally, this article discusses the challenges and future development potential of SCS and recommends the direction of future applications of SCS to musculoskeletal medicine.

Genome-scale actions of master regulators directing skeletal development

The mammalian skeleton develops through two distinct modes of ossification: intramembranous ossification and endochondral ossification. During the process of skeletal development, SRY-box containing gene 9 (Sox9), runt-related transcription factor 2 (Runx2), and Sp7 work as master transcription factors (TFs) or transcriptional regulators, underlying cell fate specification of the two distinct populations: bone-forming osteoblasts and cartilage-forming chondrocytes. In the past two decades, core transcriptional circuits underlying skeletal development have been identified mainly through mouse genetics and biochemical approaches. Recently emerging next-generation sequencer (NGS)-based studies have provided genome-scale views on the gene regulatory landscape programmed by the master TFs/transcriptional regulators. With particular focus on Sox9, Runx2, and Sp7, this review aims to discuss the gene regulatory landscape in skeletal development, which has been identified by genome-scale data, and provide future perspectives in this field.

Single-cell co-expression analysis reveals that transcriptional modules are shared across cell types in the brain

Gene-gene relationships are commonly measured via the co-variation of gene expression across samples, also known as gene co-expression. Because shared expression patterns are thought to reflect shared function, co-expression networks describe functional relationships between genes, including co-regulation. However, the heterogeneity of cell types in bulk RNA-seq samples creates connections in co-expression networks that potentially obscure co-regulatory modules. The brain initiative cell census network (BICCN) single-cell RNA sequencing (scRNA-seq) datasets provide an unparalleled opportunity to understand how gene-gene relationships shape cell identity. Comparison of the BICCN data (500,000 cells/nuclei across 7 BICCN datasets) with that of bulk RNA-seq networks (2,000 mouse brain samples across 52 studies) reveals a consistent topology reflecting a shared co-regulatory signal. Differential signals between broad cell classes persist in driving variation at finer levels, indicating that convergent regulatory processes affect cell phenotype at multiple scales.

Time-sequenced transcriptomes of developing distal mouse limb buds: A comparative tissue layer analysis

BACKGROUND

The development of the amniote limb has been an important model system to study patterning mechanisms and morphogenesis. For proper growth and patterning, it requires the interaction between the distal sub-apical mesenchyme and the apical ectodermal ridge (AER) that involve the separate implementation of coordinated and tissue-specific genetic programs.

RESULTS

Here, we produce and analyze the transcriptomes of both distal limb mesenchymal progenitors and the overlying ectodermal cells, following time-coursed dissections that cover from limb bud initiation to fully patterned limbs. The comparison of transcriptomes within each layer as well as between layers over time, allowed the identification of specific transcriptional signatures for each of the developmental stages. Special attention was given to the identification of genes whose transcription dynamics suggest a previously unnoticed role in the context of limb development and also to signaling pathways enriched between layers.

CONCLUSION

We interpret the transcriptomic data in light of the known development pattern and we conclude that a major transcriptional transition occurs in distal limb buds between E9.5 and E10.5, coincident with the switch from an early phase continuation of the signature of trunk progenitors, related to the initial proximo distal specification, to a late intrinsic phase of development.

ETS-dependent enhancers for endothelial-specific expression of serum/glucocorticoid-regulated kinase 1 during mouse embryo development

Serum/glucocorticoid-regulated kinase 1 (SGK1) is predominantly expressed in endothelial cells of mouse embryos, and Sgk1 null mice show embryonic lethality due to impaired vascular formation. However, how the SGK1 expression is controlled in developing vasculature remains unknown. In this study, we first identified a proximal endothelial enhancer through lacZ reporter mouse analyses. The mouse Sgk1 proximal enhancer was narrowed down to the 5' region of the major transcription initiation site, while a human corresponding region possessed relatively weak activity. We then searched for distal enhancer candidates using in silico analyses of publicly available databases for DNase accessibility, RNA polymerase association and chromatin modification. A region approximately 500 kb distant from the human SGK1 gene was conserved in the mouse, and the mouse and human genomic fragments drove transcription restricted to embryonic endothelial cells. Minimal fragments of both proximal and distal enhancers had consensus binding elements for the ETS transcription factors, which were essential for the responsiveness to ERG, FLI1 and ETS1 proteins in luciferase assays and the endothelial lacZ reporter expression in mouse embryos. These results suggest that endothelial SGK1 expression in embryonic vasculature is maintained through at least two ETS-regulated enhancers located in the proximal and distal regions.

Assessing evolutionary and developmental transcriptome dynamics in homologous cell types

Background

During development, complex organ patterns emerge through the precise temporal and spatial specification of different cell types. On an evolutionary timescale, these patterns can change, resulting in morphological diversification. It is generally believed that homologous anatomical structures are built—largely—by homologous cell types. However, whether a common evolutionary origin of such cell types is always reflected in the conservation of their intrinsic transcriptional specification programs is less clear.

Results

Here, we developed a user‐friendly bioinformatics workflow to detect gene co‐expression modules and test for their conservation across developmental stages and species boundaries. Using a paradigm of morphological diversification, the tetrapod limb, and single‐cell RNA‐sequencing data from two distantly related species, chicken and mouse, we assessed the transcriptional dynamics of homologous cell types during embryonic patterning. With mouse limb data as reference, we identified 19 gene co‐expression modules with varying tissue or cell type‐restricted activities. Testing for co‐expression conservation revealed modules with high evolutionary turnover, while others seemed maintained—to different degrees, in module make‐up, density or connectivity—over developmental and evolutionary timescales.

Conclusions

We present an approach to identify evolutionary and developmental dynamics in gene co‐expression modules during patterning‐relevant stages of homologous cell type specification using single‐cell RNA‐sequencing data.

We present an approach to identify evolutionary and developmental dynamics in gene co‐expression modules during patterning‐relevant stages of homologous cell type specification using single‐cell RNA‐sequencing data.

Extracellular matrix gene expression signatures as cell type and cell state identifiers

Transcriptomic signatures based on cellular mRNA expression profiles can be used to categorize cell types and states. Yet whether different functional groups of genes perform better or worse in this process remains largely unexplored. Here we test the core matrisome - that is, all genes coding for structural proteins of the extracellular matrix - for its ability to delineate distinct cell types in embryonic single-cell RNA-sequencing (scRNA-seq) data. We show that even though expressed core matrisome genes correspond to less than 2% of an entire cellular transcriptome, their RNA expression levels suffice to recapitulate essential aspects of cell type-specific clustering. Notably, using scRNA-seq data from the embryonic limb, we demonstrate that core matrisome gene expression outperforms random gene subsets of similar sizes and can match and exceed the predictive power of transcription factors. While transcription factor signatures generally perform better in predicting cell types at early stages of chicken and mouse limb development, i.e., when cells are less differentiated, the information content of the core matrisome signature increases in more differentiated cells. Moreover, using cross-species analyses, we show that these cell type-specific signatures are evolutionarily conserved. Our findings suggest that each cell type produces its own unique extracellular matrix, or matreotype, which becomes progressively more refined and cell type-specific as embryonic tissues mature.

Single-cell transcriptomics of the Drosophila wing disc reveals instructive epithelium-to-myoblast interactions

In both vertebrates and invertebrates, generating a functional appendage requires interactions between ectoderm-derived epithelia and mesoderm-derived cells. To investigate such interactions, we used single-cell transcriptomics to generate a temporal cell atlas of the Drosophila wing disc from two developmental time points. Using these data, we visualized gene expression using a multilayered model of the wing disc and cataloged ligand–receptor pairs that could mediate signaling between epithelial cells and adult muscle precursors (AMPs). We found that localized expression of the fibroblast growth factor ligands, Thisbe and Pyramus, in the disc epithelium regulates the number and location of the AMPs. In addition, Hedgehog ligand from the epithelium activates a specific transcriptional program within adjacent AMP cells, defined by AMP-specific targets Neurotactin and midline, that is critical for proper formation of direct flight muscles. More generally, our annotated temporal cell atlas provides an organ-wide view of potential cell–cell interactions between epithelial and myogenic cells.

An atlas of neural crest lineages along the posterior developing zebrafish at single-cell resolution

Neural crest cells (NCCs) are vertebrate stem cells that give rise to various cell types throughout the developing body in early life. Here, we utilized single-cell transcriptomic analyses to delineate NCC-derivatives along the posterior developing vertebrate, zebrafish, during the late embryonic to early larval stage, a period when NCCs are actively differentiating into distinct cellular lineages. We identified several major NCC/NCC-derived cell-types including mesenchyme, neural crest, neural, neuronal, glial, and pigment, from which we resolved over three dozen cellular subtypes. We dissected gene expression signatures of pigment progenitors delineating into chromatophore lineages, mesenchyme cells, and enteric NCCs transforming into enteric neurons. Global analysis of NCC derivatives revealed they were demarcated by combinatorial hox gene codes, with distinct profiles within neuronal cells. From these analyses, we present a comprehensive cell-type atlas that can be utilized as a valuable resource for further mechanistic and evolutionary investigations of NCC differentiation.

Turning Nature's own processes into design strategies for living bone implant biomanufacturing: a decade of Developmental Engineering

A decade after the term developmental engineering (DE) was coined to indicate the use of developmental processes as blueprints for the design and development of engineered living implants, a myriad of proof-of-concept studies demonstrate the potential of this approach in small animal models. This review provides an overview of DE work, focusing on applications in bone regeneration. Enabling technologies allow to quantify the distance between in vitro processes and their developmental counterpart, as well as to design strategies to reduce that distance. By embedding Nature's robust mechanisms of action in engineered constructs, predictive large animal data and subsequent positive clinical outcomes can be gradually achieved. To this end, the development of next generation biofabrication technologies should provide the necessary scale and precision for robust living bone implant biomanufacturing.

Dissecting human embryonic skeletal stem cell ontogeny by single-cell transcriptomic and functional analyses

Human skeletal stem cells (SSCs) have been discovered in fetal and adult long bones. However, the spatiotemporal ontogeny of human embryonic SSCs during early skeletogenesis remains elusive. Here we map the transcriptional landscape of human limb buds and embryonic long bones at single-cell resolution to address this fundamental question. We found remarkable heterogeneity within human limb bud mesenchyme and epithelium, and aligned them along the proximal–distal and anterior–posterior axes using known marker genes. Osteo-chondrogenic progenitors first appeared in the core limb bud mesenchyme, which give rise to multiple populations of stem/progenitor cells in embryonic long bones undergoing endochondral ossification. Importantly, a perichondrial embryonic skeletal stem/progenitor cell (eSSPC) subset was identified, which could self-renew and generate the osteochondral lineage cells, but not adipocytes or hematopoietic stroma. eSSPCs are marked by the adhesion molecule CADM1 and highly enriched with FOXP1/2 transcriptional network. Interestingly, neural crest-derived cells with similar phenotypic markers and transcriptional networks were also found in the sagittal suture of human embryonic calvaria. Taken together, this study revealed the cellular heterogeneity and lineage hierarchy during human embryonic skeletogenesis, and identified distinct skeletal stem/progenitor cells that orchestrate endochondral and intramembranous ossification.

Regulation of vertebrate forelimb development and wing reduction in the flightless emu

The vertebrate limb is a dynamic structure which has evolved into many diverse forms to facilitate complex behavioral adaptations. The principle molecular and cellular processes that underlie development of the vertebrate limb are well characterized. However, how these processes are altered to drive differential limb development between vertebrates is less well understood. Several vertebrate models are being utilized to determine the developmental basis of differential limb morphogenesis, though these typically focus on later patterning of the established limb bud and may not represent the complete developmental trajectory. Particularly, heterochronic limb development can occur prior to limb outgrowth and patterning but receives little attention. This review summarizes the genetic regulation of vertebrate forelimb diversity, with particular focus on wing reduction in the flightless emu as a model for examining limb heterochrony. These studies highlight that wing reduction is complex, with heterochronic cellular and genetic events influencing the major stages of limb development. Together, these studies provide a broader picture of how different limb morphologies may be established during development.

Atlas of Musculoskeletal Stem Cells with the Soft and Hard Tissue Differentiation Architecture

Although being of utmost importance for human health and mobility, stem cell identity and hierarchical organization of musculoskeletal progenitors remain largely unexplored. Here, cells from E10.5, E12.5, and E15.5 murine limbs are analyzed by high throughput single-cell RNA sequencing to illustrate the cellular architecture during limb development. Single-cell transcriptional profiling demonstrates the identity and differentiation architecture of musculoskeletal stem cells (MSSC), soft and hard tissue progenitors through expression pattern of musculoskeletal markers (scleraxis [Scx], Hoxd13, Sox9, and Col1a1). This is confirmed by genetic in vivo lineage tracing. Moreover, single-cell analyses of Scx knockout mice tissues illustrates that Scx regulates MSSC self-renewal and proliferation potential. A high-throughput and low-cost multi-tissues RNA sequencing strategy further provides evidence that musculoskeletal system tissues, including muscle, bone, meniscus, and cartilage, are all abnormally developed in Scx knockout mice. These results establish the presence of an indispensable limb Scx+Hoxd13+ MSSC population and their differentiation into soft tissue progenitors (Scx+Col1a1+) and hard tissue progenitors (Scx+Sox9+). Collectively, this study paves the way for systematically decoding the complex molecular mechanisms and cellular programs of musculoskeletal tissues morphogenesis in limb development and regeneration.

Transcriptomic similarities and differences between the limb bud AER and unique carapacial ridge of turtle embryos

Evolutionary innovation may arise via major departures from an ancestral condition. Turtle shell morphogenesis depends on a unique structure known as the carapacial ridge (CR). This lateral tissue protrusion in turtle embryos exhibits similar properties as the apical ectodermal ridge (AER)-a well-known molecular signaling center involved in limb development. Still, how the CR influences shell morphogenesis is not entirely clear. The present study aimed to describe the CR transcriptome shortly before ribs were halted within its mesenchyme, as required for shell development. Analyses exposed that the mesenchymal marker VIM was one of the most highly co-expressed genes and numerous appendage formation genes were situated within the core of CR and AER co-expression networks. However, there were tissue-specific differences in the activity of these genes. For instance, WNT5A was most frequently assigned to appendage-related annotations of the CR network core, but not in the AER. Several homeobox transcription factors known to regulate limb bud patterning exhibited their highest expression levels in the AER, but were underexpressed in the CR. The results of this study corroborate that novel body plans often originate via alterations of pre-existing genetic networks. Altogether, this exploratory study enhances the groundwork for future experiments on the molecular underpinnings of turtle shell development and evolution.

Single cell RNA-sequencing reveals cellular heterogeneity and trajectories of lineage specification during murine embryonic limb development

The coordinated spatial and temporal regulation of gene expression in the murine hindlimb determines the identity of mesenchymal progenitors and the development of diversity of musculoskeletal tissues they form. Hindlimb development has historically been studied with lineage tracing of individual genes selected a priori, or at the bulk tissue level, which does not allow for the determination of single cell transcriptional programs yielding mature cell types and tissues. To identify the cellular trajectories of lineage specification during limb bud development, we used single cell mRNA sequencing (scRNA-seq) to profile the developing murine hindlimb between embryonic days (E)11.5-E18.5. We found cell type heterogeneity at all time points, and the expected cell types that form the mouse hindlimb. In addition, we used RNA fluorescence in situ hybridization (FISH) to examine the spatial locations of cell types and cell trajectories to understand the ancestral continuum of cell maturation. This data provides a resource for the transcriptional program of hindlimb development that will support future studies of musculoskeletal development and generate hypotheses for tissue regeneration.

Biochemical characteristics of the chondrocyte-enriched SNORC protein and its transcriptional regulation by SOX9

Snorc (Small NOvel Rich in Cartilage) has been identified as a chondrocyte-specific gene in the mouse. Yet little is known about the SNORC protein biochemical properties, and mechanistically how the gene is regulated transcriptionally in a tissue-specific manner. The goals of the present study were to shed light on those important aspects. The chondrocyte nature of Snorc expression was confirmed in mouse and rat tissues, in differentiated (day 7) ATDC5, and in RCS cells where it was constitutive. Topological mapping and biochemical analysis brought experimental evidences that SNORC is a type I protein carrying a chondroitin sulfate (CS) attached to serine 44. The anomalous migration of SNORC on SDS-PAGE was due to its primary polypeptide features, suggesting no additional post-translational modifications apart from the CS glycosaminoglycan. A highly conserved SOX9-binding enhancer located in intron 1 was necessary to drive transcription of Snorc in the mouse, rat, and human. The enhancer was active independently of orientation and whether located in a heterologous promoter or intron. Crispr-mediated inactivation of the enhancer in RCS cells caused reduction of Snorc. Transgenic mice carrying the intronic multimerized enhancer drove high expression of a βGeo reporter in chondrocytes, but not in the hypertrophic zone. Altogether these data confirmed the chondrocyte-specific nature of Snorc and revealed dependency on the intronic enhancer binding of SOX9 for transcription.

A dot-stripe Turing model of joint patterning in the tetrapod limb

Iterative joints are a hallmark of the tetrapod limb, and their positioning is a key step during limb development. Although the molecular regulation of joint formation is well studied, it remains unclear what controls the location, number and orientation (i.e. the pattern) of joints within each digit. Here, we propose the dot-stripe mechanism for joint patterning, comprising two coupled Turing systems inspired by published gene expression patterns. Our model can explain normal joint morphology in wild-type limbs, hyperphalangy in cetacean flippers, mutant phenotypes with misoriented joints and suggests a reinterpretation of the polydactylous Ichthyosaur fins as a polygonal joint lattice. By formulating a generic dot-stripe model, describing joint patterns rather than molecular joint markers, we demonstrate that the insights from the model should apply regardless of the biological specifics of the underlying mechanism, thus providing a unifying framework to interrogate joint patterning in the tetrapod limb.

Insights into Gene Regulatory Networks in Chondrocytes

Chondrogenesis is a key developmental process that molds the framework of our body and generates the skeletal tissues by coupling with osteogenesis. The developmental processes are well-coordinated by spatiotemporal gene expressions, which are hardwired with gene regulatory elements. Those elements exist as thousands of modules of DNA sequences on the genome. Transcription factors function as key regulatory proteins by binding to regulatory elements and recruiting cofactors. Over the past 30 years, extensive attempts have been made to identify gene regulatory mechanisms in chondrogenesis, mainly through biochemical approaches and genetics. More recently, newly developed next-generation sequencers (NGS) have identified thousands of gene regulatory elements on a genome scale, and provided novel insights into the multiple layers of gene regulatory mechanisms, including the modes of actions of transcription factors, post-translational histone modifications, chromatin accessibility, the concept of pioneer factors, and three-dimensional chromatin architecture. In this review, we summarize the studies that have improved our understanding of the gene regulatory mechanisms in chondrogenesis, from the historical studies to the more recent works using NGS. Finally, we consider the future perspectives, including efforts to improve our understanding of the gene regulatory landscape in chondrogenesis and potential applications to the treatment of chondrocyte-related diseases.