An interdigit signalling centre instructs coordinate phalanx-joint formation governed by 5'Hoxd-Gli3 antagonism

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.

The BAF chromatin complex component SMARCC1 does not mediate GLI transcriptional repression of Hedgehog target genes in limb buds

Transcriptional responses to the Hedgehog (HH) signaling pathway are primarily modulated by GLI repression in the mouse limb. Previous studies suggested a role for the BAF chromatin remodeling complex in mediating GLI repression. Consistent with this possibility, the core BAF complex protein SMARCC1 is present at most active limb enhancers including the majority of GLI enhancers. However, in contrast to GLI repression which reduces chromatin accessibility, SMARCC1 maintains chromatin accessibility at most enhancers, including those bound by GLI. Moreover, SMARCC1 binding at GLI-regulated enhancers occurs independently of GLI3. Consistent with previous studies, some individual GLI target genes are mis-regulated in Smarcc1 conditional knockouts, though most GLI target genes are unaffected. Moreover, SMARCC1 is not necessary for mediating constitutive GLI repression in HH mutant limb buds. We conclude that SMARCC1 does not mediate GLI3 repression, which we propose utilizes alternative chromatin remodeling complexes.

The BAF chromatin complex component SMARCC1 does not mediate GLI transcriptional repression of Hedgehog target genes in limb buds

Transcriptional responses to the Hedgehog (HH) signaling pathway are primarily modulated by GLI repression in the mouse limb. Previous studies suggested a role for the BAF chromatin remodeling complex in mediating GLI repression. Consistent with this possibility, the core BAF complex protein SMARCC1 is present at most active limb enhancers including the majority of GLI enhancers. However, in contrast to GLI repression which reduces chromatin accessibility, SMARCC1 maintains chromatin accessibility at most enhancers, including those bound by GLI. Moreover, SMARCC1 binding at GLI-regulated enhancers occurs independently of GLI3. Consistent with previous studies, some individual GLI target genes are mis-regulated in Smarcc1 conditional knockouts, though most GLI target genes are unaffected. Moreover, SMARCC1 is not necessary for mediating constitutive GLI repression in HH mutant limb buds. We conclude that SMARCC1 does not mediate GLI3 repression, which we propose utilizes alternative chromatin remodeling complexes.

Computational model of articular cartilage regeneration induced by scaffold implantation in vivo

Computational models allow to explain phenomena that cannot be observed through an animal model, such as the strain and stress states which can highly influence regeneration of the tissue. For this purpose, we have developed a simulation tool to determine the mechanical conditions provided by the polymeric scaffold. The computational model considered the articular cartilage, the subchondral bone, and the scaffold. All materials were modeled as poroelastic, and the cartilage had linear-elastic oriented collagen fibers. This model was able to explain the remodeling process that subchondral bone goes through, and how the scaffold allowed the conditions for cartilage regeneration. These results suggest that the use of scaffolds might lead the cartilaginous tissue growth in vivo by providing a better mechanical environment. Moreover, the developed computational model demonstrated to be useful as a tool prior experimental in vivo studies, by predicting the possible outcome of newly proposed treatments allowing to discard approaches that might not bring good results.

Digits in a dish: An in vitro system to assess the molecular genetics of hand/foot development at single-cell resolution

In vitro models allow for the study of developmental processes outside of the embryo. To gain access to the cells mediating digit and joint development, we identified a unique property of undifferentiated mesenchyme isolated from the distal early autopod to autonomously re-assemble forming multiple autopod structures including: digits, interdigital tissues, joints, muscles and tendons. Single-cell transcriptomic analysis of these developing structures revealed distinct cell clusters that express canonical markers of distal limb development including: Col2a1, Col10a1, and Sp7 (phalanx formation), Thbs2 and Col1a1 (perichondrium), Gdf5, Wnt5a, and Jun (joint interzone), Aldh1a2 and Msx1 (interdigital tissues), Myod1 (muscle progenitors), Prg4 (articular perichondrium/articular cartilage), and Scx and Tnmd (tenocytes/tendons). Analysis of the gene expression patterns for these signature genes indicates that developmental timing and tissue-specific localization were also recapitulated in a manner similar to the initiation and maturation of the developing murine autopod. Finally, the in vitro digit system also recapitulates congenital malformations associated with genetic mutations as in vitro cultures of Hoxa13 mutant mesenchyme produced defects present in Hoxa13 mutant autopods including digit fusions, reduced phalangeal segment numbers, and poor mesenchymal condensation. These findings demonstrate the robustness of the in vitro digit system to recapitulate digit and joint development. As an in vitro model of murine digit and joint development, this innovative system will provide access to the developing limb tissues facilitating studies to discern how digit and articular joint formation is initiated and how undifferentiated mesenchyme is patterned to establish individual digit morphologies. The in vitro digit system also provides a platform to rapidly evaluate treatments aimed at stimulating the repair or regeneration of mammalian digits impacted by congenital malformation, injury, or disease.

Air-breathing behavior underlies the cell death in limbs of Rana pirica tadpoles

Amphibians shape their limbs by differential outgrowth of digits and interdigital regions. In contrast, amniotes employ cell death, an additional developmental system, to determine the final shape of limbs. Previous work has shown that high oxygen availability is correlated with the induction of cell death in developing limbs. Given the diversity of life histories of amphibians, it is conceivable that some amphibians are exposed to a high-oxygen environment during the tadpole phase and exhibit cell death in their limbs. Here, we examined whether air-breathing behavior underlies the cell death in limbs of aquatic tadpoles of the frog species Rana pirica. Our experimental approach revealed that R. pirica tadpoles exhibit cell death in their limbs that is likely to be induced by oxidative stress associated with their frequent air-breathing behavior.

Sonic hedgehog is not a limb morphogen but acts as a trigger to specify all digits in mice

Limb patterning by Sonic hedgehog (Shh), via either graded spatial or temporal signal integration, is a paradigm for "morphogen" function, yet how Shh instructs distinct digit identities remains controversial. Here, we bypass the Shh requirement in cell survival during outgrowth and demonstrate that a transient, early Shh pulse is both necessary and sufficient for normal mouse limb development. Shh response is only short range and is limited to the Shh-expressing zone during this time window. Shh patterns digits 1-3, anterior to this zone, by an indirect mechanism rather than direct spatial or temporal signal integration. Using a genetic relay-signaling assay, we discover that Shh also specifies digit 1/thumb (thought to be exclusively Shh independent) indirectly, and this finding implicates Shh in a unique regulatory hierarchy for digit 1 evolutionary adaptations such as opposable thumbs. This study illuminates Shh as a trigger for an indirect downstream network that becomes rapidly self-sustaining, with mechanistic relevance for limb development, regeneration, and evolution.

Downregulation of Grem1 expression in the distal limb mesoderm is a necessary precondition for phalanx development

BACKGROUND

The phalanges are the final skeletal elements to form in the vertebrate limb and their identity is regulated by signaling at the phalanx forming region (PFR) located at the tip of the developing digit ray. Here, we seek to explore the relationship between PFR activity and phalanx morphogenesis, which define the most distal limb skeletal elements, and signals associated with termination of limb outgrowth.

RESULTS

As Grem1 is extinguished in the distal chick limb mesoderm, the chondrogenesis marker Aggrecan is up-regulated in the metatarsals and phalanges. Fate mapping confirms that subridge mesoderm cells contribute to the metatarsal and phalanges when subridge Grem1 is down-regulated. Grem1 overexpression specifically blocks chick phalanx development by inhibiting PFR activity. PFR activity and digit development are also disrupted following overexpression of a Gli3 repressor, which results in Grem1 expression in the distal limb and downregulation of Bmpr1b.

CONCLUSIONS

Based on expression and fate mapping studies, we propose that downregulation of Grem1 in the distal limb marks the transition from metatarsal to phalanx development. This suggests that downregulation of Grem1 in the distal limb mesoderm is necessary for phalanx development. Grem1 downregulation allows for full PFR activity and phalanx progenitor cell commitment to digit fate.

Mechanical feedback defines organizing centers to drive digit emergence

During embryonic development, digits gradually emerge in a periodic pattern. Although genetic evidence indicates that digit formation results from a self-organizing process, the underlying mechanisms are still unclear. Here, we find that convergent-extension tissue flows driven by active stresses underlie digit formation. These active stresses simultaneously shape cartilage condensations and lead to the emergence of a compressive stress region that promotes high activin/p-SMAD/SOX9 expression, thereby defining digit-organizing centers via a mechanical feedback. In Wnt5a mutants, such mechanical feedback is disrupted due to the loss of active stresses, organizing centers do not emerge, and digit formation is precluded. Thus, digit emergence does not result solely from molecular interactions, as was previously thought, but requires a mechanical feedback that ensures continuous coupling between phalanx specification and elongation. Our work, which links mechanical and molecular signals, provides a mechanistic context for the emergence of organizing centers that may underlie various developmental processes.

BMP signaling and skeletal development in fibrodysplasia ossificans progressiva (FOP)

Fibrodysplasia ossificans progressiva (FOP) is an ultra-rare genetic disease caused by increased BMP pathway signaling due to mutation of ACVR1, a bone morphogenetic protein (BMP) type 1 receptor. The primary clinical manifestation of FOP is extra-skeletal bone formation (heterotopic ossification) within soft connective tissues. However, the underlying ACVR1 mutation additionally alters skeletal bone development and nearly all people born with FOP have bilateral malformation of the great toes as well as other skeletal malformations at diverse anatomic sites. The specific mechanisms through which ACVR1 mutations and altered BMP pathway signaling in FOP influence skeletal bone formation during development remain to be elucidated; however, recent investigations are providing a clearer understanding of the molecular and developmental processes associated with ACVR1-regulated skeletal formation.

Rethinking positional information and digit identity: The role of late interdigit signaling

Seminal work from John Fallon's lab has illuminated how digit identity determination involves ongoing late regulation and occurs progressively during phalanx formation. Complementary genetic analyses in mice and several papers in this special issue have begun to flesh out how interdigit signaling accomplishes this, but major questions remain unaddressed, including how uncommitted progenitors from which phalanges arise are maintained, and what factors set limits on digit extension and phalanx number, particularly in mammals. This review summarizes what has been learned in the two decades since control of digit identity by late interdigit signals was first identified and what remains poorly understood, and will hopefully spark renewed interest in a process that is critical to evolutionary limb adaptations but nevertheless remains enigmatic.

Genetic basis for an evolutionary shift from ancestral preaxial to postaxial limb polarity in non-urodele vertebrates

In most tetrapod vertebrates, limb skeletal progenitors condense with postaxial dominance. Posterior elements (such as ulna and fibula) appear prior to their anterior counterparts (radius and tibia), followed by digit-appearance order with continuing postaxial polarity. The only exceptions are urodele amphibians (salamanders), whose limb elements develop with preaxial polarity and who are also notable for their unique ability to regenerate complete limbs as adults. The mechanistic basis for this preaxial dominance has remained an enigma and has even been proposed to relate to the acquisition of novel genes involved in regeneration. However, recent fossil evidence suggests that preaxial polarity represents an ancestral rather than derived state. Here, we report that 5'Hoxd (Hoxd11-d13) gene deletion in mouse is atavistic and uncovers an underlying preaxial polarity in mammalian limb formation. We demonstrate this shift from postaxial to preaxial dominance in mouse results from excess Gli3 repressor (Gli3R) activity due to the loss of 5'Hoxd-Gli3 antagonism and is associated with cell-cycle changes promoting precocious cell-cycle exit in the anterior limb bud. We further show that Gli3 knockdown in axolotl results in a shift to postaxial dominant limb skeleton formation, as well as expanded paddle-shaped limb-bud morphology and ensuing polydactyly. Evolutionary changes in Gli3R activity level, which also played a key role in the fin-to-limb transition, appear to be fundamental to the shift from preaxial to postaxial polarity in formation of the tetrapod limb skeleton.

Selection on phalanx development in the evolution of the bird wing

The frameshift hypothesis is a widely-accepted model of bird wing evolution. This hypothesis postulates a shift in positional values, or molecular-developmental identity, that caused a change in digit phenotype. The hypothesis synthesised developmental and palaeontological data on wing digit homology. The 'most anterior digit' (MAD) hypothesis presents an alternative view based on changes in transcriptional regulation in the limb. The molecular evidence for both hypotheses is that the most anterior digit expresses Hoxd13 but not Hoxd11 and Hoxd12. This digit I 'signature' is thought to characterise all amniotes. Here, we studied Hoxd expression patterns in a phylogenetic sample of 18 amniotes. Instead of a conserved molecular signature in digit I, we find wide variation of Hoxd11, Hoxd12 and Hoxd13 expression in digit I. Patterns of apoptosis, and Sox9 expression, a marker of the phalanx-forming region, suggest that phalanges were lost from wing digit IV because of early arrest of the phalanx-forming region followed by cell death. Finally, we show that multiple amniote lineages lost phalanges with no frameshift. Our findings suggest that the bird wing evolved by targeted loss of phalanges under selection. Consistent with our view, some recent phylogenies based on dinosaur fossils eliminate the need to postulate a frameshift in the first place. We suggest that the phenotype of the Archaeopteryx lithographica wing is also consistent with phalanx loss. More broadly, our results support a gradualist model of evolution based on tinkering with developmental gene expression.

Over-expression of the bottlenose dolphin Hoxd13 gene in zebrafish provides new insights into the cetacean flipper formation

Cetaceans have evolved elongated soft-tissue flipper with digits made of hyperphalangy. Cetaceans were found to have 2-3 more alanine residues in Hoxd13 than other mammals, which were suggested to be related to their flipper. However, how Hoxd13 regulates other genes and induces hyperphalangy in cetaceans remain poorly understood. Here, we overexpressed the bottlenose dolphin Hoxd13 in zebrafish (Danio rerio). Combined with transcriptome data and evolutionary analyses, our results revealed that the Wingless/Integrated (Wnt) and Hedgehog signaling pathways and multiple genes might regulate hyperphalangy development in cetaceans. Meanwhile, the Notch and mitogen-activated protein kinase (Mapk) signaling pathways and Fibroblast growth factor receptor 1 (Fgfr1) are probably correlated with interdigital tissues retained in the cetacean flipper. In conclusion, this is the first study to use a transgenic zebrafish to explore the molecular evolution of Hoxd13 in cetaceans, and it provides new insights into cetacean flipper formation.

Dbx2 regulation in limbs suggests interTAD sharing of enhancers

BACKGROUND

During tetrapod limb development, the HOXA13 and HOXD13 transcription factors are critical for the emergence and organization of the autopod, the most distal aspect where digits will develop. Since previous work had suggested that the Dbx2 gene is a target of these factors, we set up to analyze in detail this potential regulatory interaction.

RESULTS

We show that HOX13 proteins bind to mammalian-specific sequences at the vicinity of the Dbx2 locus that have enhancer activity in developing digits. However, the functional inactivation of the DBX2 protein did not elicit any particular phenotype related to Hox genes inactivation in digits, suggesting either redundant or compensatory mechanisms. We report that the neighboring Nell2 and Ano6 genes are also expressed in distal limb buds and are in part controlled by the same Dbx2 enhancers despite being localized into two different topologically associating domains (TADs) flanking the Dbx2 locus.

CONCLUSIONS

We conclude that Hoxa13 and Hoxd genes cooperatively activate Dbx2 expression in developing digits through binding to mammalian specific regulatory sequences in the Dbx2 neighborhood. Furthermore, these enhancers can overcome TAD boundaries in either direction to co-regulate a set of genes located in distinct chromatin domains.

Morphospace analysis leads to an evo-devo model of digit patterning

Biological forms occupy a constrained portion of theoretical morphospaces. Developmental models accounting for empirical morphospaces are necessary to achieve a better understanding of this phenomenon. We analyzed the phalangeal formulas (PFs) in lizards and relatives' hands by comparing them with a set of simulated PFs that compose a theoretical morphospace. We detected that: (1) the empirical morphospace is severely limited in size, (2) the PFs comply with two properties of phalangeal count per digit, namely the ordering rule (DI ≤ DII ≤ DIII ≤ DIV ≥ DV), and the contiguity relationship (neighbor digits differ on average in one phalanx), (3) the totality of the PFs can be categorized into four categories of hands aligned along a feasibility gradient. We also reconstructed the evolution of PFs and found a stepwise trajectory from the plesiomorphic PF towards reduced conditions. Finally, we propose a developmental model as the generative mechanism behind the PFs. It is consistent with the bulk of evidence managed and involves an ordered digit primordia initialization timed with periodic signals of joint formation coming from digit tips. Our approach is also useful to address the study of other meristic sequences in nature such as dental, floral, and branchial formulas.

WNT5A-Ca2+-CaN-NFAT signalling plays a permissive role during cartilage differentiation in embryonic chick digit development

During digit development, the correct balance of chondrogenic signals ensures the recruitment of undifferentiated cells into the cartilage lineage or the maintenance of cells at the undifferentiated stage. WNT/β catenin maintains the pool of progenitor cells, whereas TGFβ signalling promotes cartilage differentiation by inducing Sox9 expression. Moreover, WNT5A promotes the degradation of β catenin during mouse limb development. Although these mechanisms are well established, it is still unknown whether the signalling pathway downstream WNT5A is also involved in early chondrogenesis during digit formation. Thus, the aim of this study was to determine the role of WNT5A during the recruitment of progenitor cells during digit development. Our results showed that WNT5A activated calcium (Ca²⁺) release in the undifferentiated region during digit development. Further, the blockade of Ca²⁺ release or calcineurin (CaN) or nuclear factor of activated T-cells (NFAT) functions resulted in an inhibition of cartilage differentiation. Together, our results demonstrate that non canonical WNT5A-Ca²⁺-CaN-NFAT signalling plays a key role during embryonic digit development in vivo promoting the competence for chondrogenic signals and also acts as a permissive factor for chondrogenesis independently of cell death mechanisms.

The Developmental Phenotype of the Great Toe in Fibrodysplasia Ossificans Progressiva

Fibrodysplasia ossificans progressiva (FOP) is a rare genetic disorder in which extensive heterotopic ossification (HO) begins to form during early childhood and progresses throughout life. Although HO does not occur during embryonic development, children who carry the ACVR1^R206H mutation that causes most cases of FOP characteristically exhibit malformation of their great toes at birth, indicating that the mutation acts during embryonic development to alter skeletal formation. Despite the high prevalence of the great toe malformation in the FOP population, it has received relatively little attention due to its clinically benign nature. In this study, we examined radiographs from a cohort of 41 FOP patients ranging from 2 months to 48 years of age to provide a detailed analysis of the developmental features, progression, and variability of the great toe malformation of FOP, which include absent skeletal structures, malformed epiphyses, ectopic ossification centers, malformed first metatarsals and phalangeal fusion.

Dysregulated BMP signaling through ACVR1 impairs digit joint development in fibrodysplasia ossificans progressiva (FOP)

The development of joints in the mammalian skeleton depends on the precise regulation of multiple interacting signaling pathways including the bone morphogenetic protein (BMP) pathway, a key regulator of joint development, digit patterning, skeletal growth, and chondrogenesis. Mutations in the BMP receptor ACVR1 cause the rare genetic disease fibrodysplasia ossificans progressiva (FOP) in which extensive and progressive extra-skeletal bone forms in soft connective tissues after birth. These mutations, which enhance BMP-pSmad1/5 pathway activity to induce ectopic bone, also affect skeletal development. FOP can be diagnosed at birth by symmetric, characteristic malformations of the great toes (first digits) that are associated with decreased joint mobility, shortened digit length, and absent, fused, and/or malformed phalanges. To elucidate the role of ACVR1-mediated BMP signaling in digit skeletal development, we used an Acvr1^R206H/+;Prrx1-Cre knock-in mouse model that mimics the first digit phenotype of human FOP. We have determined that the effects of increased Acvr1-mediated signaling by the Acvr1^R206H mutation are not limited to the first digit but alter BMP signaling, Gdf5+ joint progenitor cell localization, and joint development in a manner that differently affects individual digits during embryogenesis. The Acvr1^R206H mutation leads to delayed and disrupted joint specification and cleavage in the digits and alters the development of cartilage and endochondral ossification at sites of joint morphogenesis. These findings demonstrate an important role for ACVR1-mediated BMP signaling in the regulation of joint and skeletal formation, show a direct link between failure to restrict BMP signaling in the digit joint interzone and failure of joint cleavage at the presumptive interzone, and implicate impaired, digit-specific joint development as the proximal cause of digit malformation in FOP.

Regulation of the limb shape during the development of the Chinese softshell turtles

Interdigital cell death is an important mechanism employed by amniotes to shape their limbs; inhibiting this process leads to the formation of webbed fingers, as seen in bats and ducks. The Chinese softshell turtle Pelodiscus sinensis (Reptilia: Testudines: Trionychidae) has a distinctive limb morphology: the anterior side of the limbs has partially webbed fingers with claw‐like protrusions, while the posterior fingers are completely enclosed in webbings. Here, P. sinensis embryos were investigated to gain insights on the evolution of limb‐shaping mechanisms in amniotes. We found cell death and cell senescence in their interdigital webbings. Spatial or temporal modulation of these processes were correlated with the appearance of indentations in the webbings, but not a complete regression of this tissue. No differences in interdigital cell proliferation were found. In subsequent stages, differential growth of the finger cartilages led to a major difference in limb shape. While no asymmetry in bone morphogenetic protein signaling was evident during interdigital cell death stages, some components of this pathway were expressed exclusively in the clawed digit tips, which also had earlier ossification. In addition, a delay and/or truncation in the chondrogenesis of the posterior digits was found in comparison with the anterior digits of P. sinensis, and also when compared with the previously published pattern of digit skeletogenesis of turtles without posterior webbings. In conclusion, modulation of cell death, as well as a heterochrony in digit chondrogenesis, may contribute to the formation of the unique limbs of the Chinese softshell turtles.

Cell death and senescence shape the interdigital webbings of Pelodiscus sinensis.

Delayed chondrogenesis/ossification and truncated tips are found in posterior digits, as well as differential expression of bone morphogenetic proteins and Msh homeobox 1 transcription factors.

The tissues and regulatory pattern of limb chondrogenesis

Initial limb chondrogenesis offers the first differentiated tissues that resemble the mature skeletal anatomy. It is a developmental progression of three tissues. The limb begins with undifferentiated mesenchyme-1, some of which differentiates into condensations-2, and this tissue then transforms into cartilage-3. Each tissue is identified by physical characteristics of cell density, shape, and extracellular matrix composition. Tissue specific regimes of gene regulation underlie the diagnostic physical and chemical properties of these three tissues. These three tissue based regimes co-exist amid a background of other gene regulatory regimes within the same tissues and time-frame of limb development. The bio-molecular indicators of gene regulation reveal six identifiable patterns. Three of these patterns describe the unique bio-molecular indicators of each of the three tissues. A fourth pattern shares bio-molecular indicators between condensation and cartilage. Finally, a fifth pattern is composed of bio-molecular indicators that are found in undifferentiated mesenchyme prior to any condensation differentiation, then these bio-molecular indicators are upregulated in condensations and downregulated in undifferentiated mesenchyme. The undifferentiated mesenchyme that remains in between the condensations and cartilage, the interdigit, contains a unique set of bio-molecular indicators that exhibit dynamic behaviour during chondrogenesis and therefore argue for its own inclusion as a tissue in its own right and for more study into this process of differentiation.

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.

The formation of the thumb requires direct modulation of Gli3 transcription by Hoxa13

In the tetrapod limb, the digits (fingers or toes) are the elements most subject to morphological diversification in response to functional adaptations. However, despite their functional importance, the mechanisms controlling digit morphology remain poorly understood. Here we have focused on understanding the special morphology of the thumb (digit 1), the acquisition of which was an important adaptation of the human hand. To this end, we have studied the limbs of the Hoxa13 mouse mutant that specifically fail to form digit 1. We show that, consistent with the role of Hoxa13 in Hoxd transcriptional regulation, the expression of Hoxd13 in Hoxa13 mutant limbs does not extend into the presumptive digit 1 territory, which is therefore devoid of distal Hox transcripts, a circumstance that can explain its agenesis. The loss of Hoxd13 expression, exclusively in digit 1 territory, correlates with increased Gli3 repressor activity, a Hoxd negative regulator, resulting from increased Gli3 transcription that, in turn, is due to the release from the negative modulation exerted by Hox13 paralogs on Gli3 regulatory sequences. Our results indicate that Hoxa13 acts hierarchically to initiate the formation of digit 1 by reducing Gli3 transcription and by enabling expansion of the 5'Hoxd second expression phase, thereby establishing anterior-posterior asymmetry in the handplate. Our work uncovers a mutual antagonism between Gli3 and Hox13 paralogs that has important implications for Hox and Gli3 gene regulation in the context of development and evolution.

The benefits differential equations bring to limb development

Systems biology is a large field, offering a number of advantages to a variety of biological disciplines. In limb development, differential-equation based models can provide insightful hypotheses about the gene/protein interactions and tissue differentiation events that form the core of limb development research. Differential equations are like any other communicative tool, with misuse and limitations that can come along with their advantages. Every theory should be critically analyzed to best ascertain whether they reflect the reality in biology as well they claim. Differential equation-based models have consistent features which researchers have drawn upon to aid in more realistic descriptions and hypotheses. Nine features are described that highlight these trade-offs. The advantages range from more detailed descriptions of gene interactions and their consequence and the capacity to model robustness to the incorporation of tissue size and shape. The drawbacks come with the added complication that additional genes and signaling pathways that require additional terms within the mathematical model. They also come in the translation between the mathematical terms of the model, values and matrices, to the real world of genes, proteins, and tissues that constitute limb development. A critical analysis is necessary to ensure that these models effectively expand the understanding of the origins of a diversity of limb anatomy, from evolution to teratology. This article is categorized under: Vertebrate Organogenesis > Musculoskeletal and Vascular Gene Expression and Transcriptional Hierarchies > Regulatory Mechanisms Establishment of Spatial and Temporal Patterns > Repeating Patterns and Lateral Inhibition.

Mechanisms of synovial joint and articular cartilage development

Articular cartilage is formed at the end of epiphyses in the synovial joint cavity and permanently contributes to the smooth movement of synovial joints. Most skeletal elements develop from transient cartilage by a biological process known as endochondral ossification. Accumulating evidence indicates that articular and growth plate cartilage are derived from different cell sources and that different molecules and signaling pathways regulate these two kinds of cartilage. As the first sign of joint development, the interzone emerges at the presumptive joint site within a pre-cartilage tissue. After that, joint cavitation occurs in the center of the interzone, and the cells in the interzone and its surroundings gradually form articular cartilage and the synovial joint. During joint development, the interzone cells continuously migrate out to the epiphyseal cartilage and the surrounding cells influx into the joint region. These complicated phenomena are regulated by various molecules and signaling pathways, including GDF5, Wnt, IHH, PTHrP, BMP, TGF-β, and FGF. Here, we summarize current literature and discuss the molecular mechanisms underlying joint formation and articular development.

Environmental Oxygen Exposure Allows for the Evolution of Interdigital Cell Death in Limb Patterning

Amphibians form fingers without webbing by differential growth between digital and interdigital regions. Amniotes, however, employ interdigital cell death (ICD), an additional mechanism that contributes to a greater variation of limb shapes. Here, we investigate the role of environmental oxygen in the evolution of ICD in tetrapods. While cell death is restricted to the limb margin in amphibians with aquatic tadpoles, Eleutherodactylus coqui, a frog with terrestrial-direct-developing eggs, has cell death in the interdigital region. Chicken requires sufficient oxygen and reactive oxygen species to induce cell death, with the oxygen tension profile itself being distinct between the limbs of chicken and Xenopus laevis frogs. Notably, increasing blood vessel density in X. laevis limbs, as well as incubating tadpoles under high oxygen levels, induces ICD. We propose that the oxygen available to terrestrial eggs was an ecological feature crucial for the evolution of ICD, made possible by conserved autopod-patterning mechanisms.

Molecular signatures identify immature mesenchymal progenitors in early mouse limb buds that respond differentially to morphogen signaling

The key molecular interactions governing vertebrate limb bud development are a paradigm for studying the mechanisms controlling progenitor cell proliferation and specification during vertebrate organogenesis. However, little is known about the cellular heterogeneity of the mesenchymal progenitors in early limb buds that ultimately contribute to the chondrogenic condensations prefiguring the skeleton. We combined flow cytometric and transcriptome analyses to identify the molecular signatures of several distinct mesenchymal progenitor cell populations present in early mouse forelimb buds. In particular, jagged 1 (JAG1)-positive cells located in the posterior-distal mesenchyme were identified as the most immature limb bud mesenchymal progenitors (LMPs), which crucially depend on SHH and FGF signaling in culture. The analysis of gremlin 1 (Grem1)-deficient forelimb buds showed that JAG1-expressing LMPs are protected from apoptosis by GREM1-mediated BMP antagonism. At the same stage, the osteo-chondrogenic progenitors (OCPs) located in the core mesenchyme are already actively responding to BMP signaling. This analysis sheds light on the cellular heterogeneity of the early mouse limb bud mesenchyme and on the distinct response of LMPs and OCPs to morphogen signaling.

A single-cell transcriptomic atlas of the developing chicken limb

BACKGROUND

Through precise implementation of distinct cell type specification programs, differentially regulated in both space and time, complex patterns emerge during organogenesis. Thanks to its easy experimental accessibility, the developing chicken limb has long served as a paradigm to study vertebrate pattern formation. Through decades' worth of research, we now have a firm grasp on the molecular mechanisms driving limb formation at the tissue-level. However, to elucidate the dynamic interplay between transcriptional cell type specification programs and pattern formation at its relevant cellular scale, we lack appropriately resolved molecular data at the genome-wide level. Here, making use of droplet-based single-cell RNA-sequencing, we catalogue the developmental emergence of distinct tissue types and their transcriptome dynamics in the distal chicken limb, the so-called autopod, at cellular resolution.

RESULTS

Using single-cell RNA-sequencing technology, we sequenced a total of 17,628 cells coming from three key developmental stages of chicken autopod patterning. Overall, we identified 23 cell populations with distinct transcriptional profiles. Amongst them were small, albeit essential populations like the apical ectodermal ridge, demonstrating the ability to detect even rare cell types. Moreover, we uncovered the existence of molecularly distinct sub-populations within previously defined compartments of the developing limb, some of which have important signaling functions during autopod pattern formation. Finally, we inferred gene co-expression modules that coincide with distinct tissue types across developmental time, and used them to track patterning-relevant cell populations of the forming digits.

CONCLUSIONS

We provide a comprehensive functional genomics resource to study the molecular effectors of chicken limb patterning at cellular resolution. Our single-cell transcriptomic atlas captures all major cell populations of the developing autopod, and highlights the transcriptional complexity in many of its components. Finally, integrating our data-set with other single-cell transcriptomics resources will enable researchers to assess molecular similarities in orthologous cell types across the major tetrapod clades, and provide an extensive candidate gene list to functionally test cell-type-specific drivers of limb morphological diversification.

BMPs, TGFβ, and border security at the interzone

Synovial joints enable movement and protect the integrity of the articular cartilage. Joints form within skeletal condensations destined to undergo chondrogenesis. The suppression of this chondrogenic program in the interzone is the first morphological sign of joint formation. While we have a fairly good understanding of the essential roles of BMP and TGFβ family members in promoting chondrogenic differentiation in developing skeletal elements, we know very little about how BMP activity is suppressed specifically within the interzone, a crucial step in joint development. The function of the BMP ligand Gdf5 has been especially difficult to decipher. On the one hand, Gdf5 is required to promote chondrogenesis of articular elements. On the other hand, Gdf5 is highly expressed in the joint interzone where chondrogenesis must be suppressed for the formation of many joints. Here we review the evidence that BMP signaling must be suppressed within the joint interzone for joint morphogenesis to progress, and consider how Gdf5 exerts its divergent effects on chondrogenesis and joint formation. We also consider how TGFβ signaling impacts formation of the interzone. Finally, we propose a model whereby Gdf5 exerts distinct effects in the interzone vs. surrounding cartilage based on the repertoire of BMP receptors available in these tissues. Understanding how BMP antagonists and counteracting TGFβ signals intersect with Gdf5 to sculpt the joint interzone is essential for understanding the origin of osteoarthritis and other diseases of joint tissues.

Expression analysis of limb element markers during mouse embryonic development

Background: While data regarding expression of limb element and tissue markers during normal mouse limb development exist, few studies show expression patterns in upper and lower limbs throughout key limb development stages. A comparison to normal developmental events is essential when analyzing development of the limb in mutant mice models. Results: Expression patterns of the joint marker Gdf5, tendon and ligament marker Scleraxis, early muscle marker MyoD1, and blood vessel marker Cadherin5 (Cdh5) are presented during the most active phases of embryonic mouse limb patterning. Anti‐neurofilament staining of developing nerves in the fore‐ and hindlimbs and cartilage formation and progression also are described. Conclusions: This study demonstrates and describes a range of key morphological markers and methods that together can be used to assess normal and abnormal limb development. Developmental Dynamics 247:1217–1226, 2018. © 2018 The Authors. Developmental Dynamics published by Wiley Periodicals, Inc. on behalf of American Association of Anatomists

Expression patterns of molecular markers throughout both fore‐ and hindlimb development ‐ which can be used to assess normal and abnormal development.

Detailled description of innervation during fore‐ and hindlimb development confirming innervation first seen after limb patterning events have begun.

Description of cartilage development and progression indicates alizarin red staining not seen until E15.5 in both fore‐ and hindlimbs.

Hindlimb lags behind forelimb molecularly and morphologically until E14.5.

Detailled description of methods used to study fore‐ and hindlimb development.

Evolution of antero-posterior patterning of the limb: Insights from the chick

The developing limbs of chicken embryos have served as pioneering models for understanding pattern formation for over a century. The ease with which chick wing and leg buds can be experimentally manipulated, while the embryo is still in the egg, has resulted in the discovery of important developmental organisers, and subsequently, the signals that they produce. Sonic hedgehog (Shh) is produced by mesenchyme cells of the polarizing region at the posterior margin of the limb bud and specifies positional values across the antero-posterior axis (the axis running from the thumb to the little finger). Detailed experimental embryology has revealed the fundamental parameters required to specify antero-posterior positional values in response to Shh signaling in chick wing and leg buds. In this review, the evolution of the avian wing and leg will be discussed in the broad context of tetrapod paleontology, and more specifically, ancestral theropod dinosaur paleontology. How the parameters that dictate antero-posterior patterning could have been modulated to produce the avian wing and leg digit patterns will be considered. Finally, broader speculations will be made regarding what the antero-posterior patterning of chick limbs can tell us about the evolution of other digit patterns, including those that were found in the limbs of the earliest tetrapods.

On the Formation of Digits and Joints during Limb Development

Critical steps in forming the vertebrate limb include the positioning of digits and the positioning of joints within each digit. Recent studies have proposed that the iterative series of digits is established by a Turing-like mechanism generating stripes of chondrogenic domains. However, re-examination of available data suggest that digits are actually patterned as evenly spaced spots, not stripes, which then elongate into rod-shaped digit rays by incorporating new cells at their tips. Moreover, extension of the digit rays and the patterning of the joints occur simultaneously at the distal tip, implying that an integrated model is required to fully understand these processes.

Sonic Hedgehog Signaling in Limb Development

The gene encoding the secreted protein Sonic hedgehog (Shh) is expressed in the polarizing region (or zone of polarizing activity), a small group of mesenchyme cells at the posterior margin of the vertebrate limb bud. Detailed analyses have revealed that Shh has the properties of the long sought after polarizing region morphogen that specifies positional values across the antero-posterior axis (e.g., thumb to little finger axis) of the limb. Shh has also been shown to control the width of the limb bud by stimulating mesenchyme cell proliferation and by regulating the antero-posterior length of the apical ectodermal ridge, the signaling region required for limb bud outgrowth and the laying down of structures along the proximo-distal axis (e.g., shoulder to digits axis) of the limb. It has been shown that Shh signaling can specify antero-posterior positional values in limb buds in both a concentration- (paracrine) and time-dependent (autocrine) fashion. Currently there are several models for how Shh specifies positional values over time in the limb buds of chick and mouse embryos and how this is integrated with growth. Extensive work has elucidated downstream transcriptional targets of Shh signaling. Nevertheless, it remains unclear how antero-posterior positional values are encoded and then interpreted to give the particular structure appropriate to that position, for example, the type of digit. A distant cis-regulatory enhancer controls limb-bud-specific expression of Shh and the discovery of increasing numbers of interacting transcription factors indicate complex spatiotemporal regulation. Altered Shh signaling is implicated in clinical conditions with congenital limb defects and in the evolution of the morphological diversity of vertebrate limbs.

GLI3-related polydactyly: a review

GLI3 mutations are known to be associated with nine syndromes/conditions in which polydactyly is a feature. In this review, the embryology, pathogenesis, and animal models of GLI3-related polydactyly are discussed first. This is followed by a detailed review of the genotype-phenotype correlations. Based on our review of the literature and our clinical experiences, we recommend viewing GLI3-related syndromes/conditions as four separate entities; each characterized by a specific pattern of polydactyly. These four entities are: the preaxial polydactyly type IV-Greig-acrocallosal spectrum, postaxial polydactyly types A/B, Pallister-Hall syndrome (PHS), and oral-facial-digital overlap syndrome. We also provide illustrative clinical examples from our practice including a family with a novel GLI3 mutation causing PHS. The review also introduces the term 'Forme Fruste' preaxial polydactyly and gives several conclusions/recommendations including the recommendation to revise the current criteria for the clinical diagnosis of PHS.

John Saunders' ZPA, Sonic hedgehog and digit identity - How does it really all work?

Among John Saunders' many seminal contributions to developmental biology, his discovery of the limb 'zone of polarizing activity' (ZPA) is arguably one of the most memorable and ground-breaking. This discovery introduced the limb as a premier model for understanding developmental patterning and promoted the concept of patterning by a morphogen gradient. In the 50 years since the discovery of the ZPA, Sonic hedgehog (Shh) has been identified as the ZPA factor and the basic components of the signaling pathway and many aspects of its regulation have been elucidated. Although much has also been learned about how it regulates growth, the mechanism by which Shh patterns the limb, how it acts to instruct digit 'identity', nevertheless remains an enigma. This review focuses on what has been learned about Shh function in the limb and the outstanding puzzles that remain to be solved.