Analysis of Hox gene expression in the chick limb bud

Shaping Hox gene activity to generate morphological diversity across vertebrate phylogeny

The importance of Hox genes for the development and evolution of the vertebrate axial skeleton and paired appendages has been recognized for already several decades. The steady growth of genomic sequence data from an increasing number of vertebrate species, together with the improvement of methods to analyze genomic structure and interactions, as well as to control gene activity in various species has refined our understanding of Hox gene activity in development and evolution. Here, I will review recent data addressing the influence of Hox regulatory processes in the evolution of the fins and the emergence of the tetrapod limb. In addition, I will discuss the involvement of posterior Hox genes in the control of vertebrate axial extension, focusing on an apparently divergent activity that Hox13 paralog group genes have on the regulation of tail bud development in mouse and zebrafish embryos.

Catching a killer: Mechanisms of programmed cell death and immune activation in Amyotrophic Lateral Sclerosis

In the central nervous system (CNS), execution of programmed cell death (PCD) is crucial for proper neurodevelopment. However, aberrant activation of these pathways in adult CNS leads to neurodegenerative diseases including amyotrophic lateral sclerosis (ALS). How a cell dies is critical, as it can drive local immune activation and tissue damage. Classical apoptosis engages several mechanisms to evoke "immunologically silent" responses, whereas other forms of programmed death such as pyroptosis, necroptosis, and ferroptosis release molecules that can potentiate immune responses and inflammation. In ALS, a fatal neuromuscular disorder marked by progressive death of lower and upper motor neurons, several cell types in the CNS express machinery for multiple PCD pathways. The specific cell types engaging PCD, and ultimate mechanisms by which neuronal death occurs in ALS are not well defined. Here, we provide an overview of different PCD pathways implicated in ALS. We also examine immune activation in ALS and differentiate apoptosis from necrotic mechanisms based on downstream immunological consequences. Lastly, we highlight therapeutic strategies that target cell death pathways in the treatment of neurodegeneration and inflammation in ALS.

Untangling the relationship between developmental and evolutionary integration

Patterns of integration and modularity among organismal traits are prevalent across the tree of life, and at multiple scales of biological organization. Over the past several decades, researchers have studied these patterns at the developmental, and evolutionary levels. While their work has identified the potential drivers of these patterns at different scales, there appears to be a lack of consensus on the relationship between developmental and evolutionary integration. Here, we review and summarize key studies and build a framework to describe the conceptual relationship between these patterns across organismal scales and illustrate how, and why some of these studies may have yielded seemingly conflicting outcomes. We find that among studies that analyze patterns of integration and modularity using morphological data, the lack of consensus may stem in part from the difficulty of fully disentangling the developmental and functional causes of integration. Nonetheless, in some empirical systems, patterns of evolutionary modularity have been found to coincide with expectations based on developmental processes, suggesting that in some circumstances, developmental modularity may translate to evolutionary modularity. We also advance an extension to Hallgrímsson et al.'s palimpsest model to describe how patterns of trait modularity may shift across different evolutionary scales. Finally, we also propose some directions for future research which will hopefully be useful for investigators interested in these issues.

Retinoic acid influences the timing and scaling of avian wing development

A fundamental question in biology is how embryonic development is timed between different species. To address this problem, we compared wing development in the quail and the larger chick. We reveal that pattern formation is faster in the quail as determined by the earlier activation of 5′Hox genes, termination of developmental organizers (Shh and Fgf8), and the laying down of the skeleton (Sox9). Using interspecies tissue grafts, we show that developmental timing can be reset during a critical window of retinoic acid signaling. Accordingly, extending the duration of retinoic acid signaling switches developmental timing between the quail and the chick and the chick and the larger turkey. However, the incremental growth rate is comparable between all three species, suggesting that the pace of development primarily governs differences in the expansion of the skeletal pattern. The widespread distribution of retinoic acid could coordinate developmental timing throughout the embryo.

•

Quail wings develop faster than chick and turkey wings

•

Retinoic acid can set the species timing of wing development

•

Developmental timing is independent of growth and scales the skeletal pattern

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.

Comprehensive profiling of novel epithelial-mesenchymal transition mediators and their clinical significance in colorectal cancer

Epithelial–mesenchymal transition (EMT) is a drastic phenotypic change during cancer metastasis and is one of the most important hallmarks of aggressive cancer. Although the overexpression of some specific transcription factors explains the functional alteration of EMT-induced cells, a complete picture of this biological process is yet to be elucidated. To comprehensively profile EMT-related genes in colorectal cancer, we quantified the EMT induction ability of each gene according to its similarity to the cancer stromal gene signature and termed it “mesenchymal score.” This bioinformatic approach successfully identified 90 candidate EMT mediators, which are strongly predictive of survival in clinical samples. Among these candidates, we discovered that the neuronal gene ARC, possibly originating from the retrotransposon, unexpectedly plays a crucial role in EMT induction. Profiling of novel EMT mediators we demonstrated here may help understand the complexity of the EMT program and open up new avenues for therapeutic intervention in colorectal cancer.

Role of HOXC10 in Cancer

The HOXC10 gene, a member of the HOX genes family, plays crucial roles in mammalian physiological processes, such as limb morphological development, limb regeneration, and lumbar motor neuron differentiation. HOXC10 is also associated with angiogenesis, fat metabolism, and sex regulation. Additional evidence suggests that HOXC10 dysregulation is closely associated with various tumors. HOXC10 is an important transcription factor that can activate several oncogenic pathways by regulating various target molecules such as ERK, AKT, p65, and epithelial mesenchymal transition-related genes. HOXC10 also induces drug resistance in cancers by promoting the DNA repair pathway. In this review, we summarize HOXC10 gene structure and expression as well as the role of HOXC10 in different human cancer processes. This review will provide insight into the status of HOXC10 research and help identify novel targets for cancer therapy.

Hoxd13/Bmp2-mediated mechanism involved in zebrafish finfold design

The overexpression of hoxd13a during zebrafish fin development causes distal endochondral expansion and simultaneous reduction of the finfold, mimicking the major events thought to have happened during the fin-to-limb transition in Vertebrates. We investigated the effect of hoxd13a overexpression on putative downstream targets and found it to cause downregulation of proximal fin identity markers (meis1 and emx2) and upregulation of genes involved in skeletogenesis/patterning (fbn1, dacha) and AER/Finfold maintenance (bmps). We then show that bmp2b overexpression leads to finfold reduction, recapitulating the phenotype observed in hoxd13a-overexpressing fins. In addition, we show that during the development of the long finfold in leo^t1/lof^dt1 mutants, hoxd13a and bmp2b are downregulated. Our results suggest that modulation of the transcription factor Hoxd13 during evolution may have been involved in finfold reduction through regulation of the Bmp signalling that then activated apoptotic mechanisms impairing finfold elongation.

Latent developmental potential to form limb-like skeletal structures in zebrafish

Changes in appendage structure underlie key transitions in vertebrate evolution. Addition of skeletal elements along the proximal-distal axis facilitated critical transformations, including the fin-to-limb transition that permitted generation of diverse modes of locomotion. Here, we identify zebrafish mutants that form supernumerary long bones in their pectoral fins. These new bones integrate into musculature, form joints, and articulate with neighboring elements. This phenotype is caused by activating mutations in previously unrecognized regulators of appendage patterning, vav2 and waslb, that function in a common pathway. This pathway is required for appendage development across vertebrates, and loss of Wasl in mice causes defects similar to those seen in murine Hox mutants. Concordantly, formation of supernumerary bones requires Hox11 function, and mutations in the vav2/wasl pathway drive enhanced expression of hoxa11b, indicating developmental homology with the forearm. Our findings reveal a latent, limb-like pattern ability in fins that is activated by simple genetic perturbation.

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.

Apical ectodermal ridge regulates three principal axes of the developing limb

Understanding limb development not only gives insights into the outgrowth and differentiation of the limb, but also has clinical relevance. Limb development begins with two paired limb buds (forelimb and hindlimb buds), which are initially undifferentiated mesenchymal cells tipped with a thickening of the ectoderm, termed the apical ectodermal ridge (AER). As a transitional embryonic structure, the AER undergoes four stages and contributes to multiple axes of limb development through the coordination of signalling centres, feedback loops, and other cell activities by secretory signalling and the activation of gene expression. Within the scope of proximodistal patterning, it is understood that while fibroblast growth factors (FGFs) function sequentially over time as primary components of the AER signalling process, there is still no consensus on models that would explain proximodistal patterning itself. In anteroposterior patterning, the AER has a dual-direction regulation by which it promotes the sonic hedgehog (Shh) gene expression in the zone of polarizing activity (ZPA) for proliferation, and inhibits Shh expression in the anterior mesenchyme. In dorsoventral patterning, the AER activates Engrailed-1 (En1) expression, and thus represses Wnt family member 7a (Wnt7a) expression in the ventral ectoderm by the expression of Fgfs, Sp6/8, and bone morphogenetic protein (Bmp) genes. The AER also plays a vital role in shaping the individual digits, since levels of Fgf4/8 and Bmps expressed in the AER affect digit patterning by controlling apoptosis. In summary, the knowledge of crosstalk within AER among the three main axes is essential to understand limb growth and pattern formation, as the development of its areas proceeds simultaneously.

Sarcopterygian fin ontogeny elucidates the origin of hands with digits

How the hand and digits originated from fish fins during the Devonian fin-to-limb transition remains unsolved. Controversy in this conundrum stems from the scarcity of ontogenetic data from extant lobe-finned fishes. We report the patterning of an autopod-like domain by hoxa13 during fin development of the Australian lungfish, the most closely related extant fish relative of tetrapods. Differences from tetrapod limbs include the absence of digit-specific expansion of hoxd13 and hand2 and distal limitation of alx4 and pax9, which potentially evolved through an enhanced response to shh signaling in limbs. These developmental patterns indicate that the digit program originated in postaxial fin radials and later expanded anteriorly inside of a preexisting autopod-like domain during the evolution of limbs. Our findings provide a genetic framework for the transition of fins into limbs that supports the significance of classical models proposing a bending of the tetrapod metapterygial axis.

Genetics of developmental dysplasia of the hip

In the last decade, the advances in the molecular analyses and sequencing techniques allowed researchers to study developmental dysplasia of the hip (DDH) more thoroughly. Certain chromosomes, genes, loci and polymorphisms are being associated with variable severity of this disorder. The wide range of signs and symptoms is dependent either on isolated or systemic manifestation. Phenotypes of isolated cases range from only a mild ligamental laxity, through subluxation, to a complete dislocation of the femoral head. Systemic manifestation is connected to various forms of skeletal dysplasia and other malformations characterized by significant genetic aberrations. To reveal the background of DDH heredity, multiple studies focused on large sample sizes with an emphasis on the correlation between genotype, phenotype and continuous clinical examination. Etiological risk factors that have been observed and documented in patients include genetic, environmental, and mechanical factors, which significantly contribute to the familial or nonfamilial occurrence and phenotypic variability of this disorder. Still, the multifactorial etiology and pathogenesis of DDH are not yet sufficiently clarified, explained, or understood. Formation of connective tissue, osteogenesis, chondrogenesis, and all other affected pathways and variations in the function of their individual elements contribute to the creation of the pathology in a developing human body. This review article presents an up-to-date list of known DDH associated genes, their products, and functional characteristics.

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 Logical Model of Homology for Comparative Biology

There is a growing body of research on the evolution of anatomy in a wide variety of organisms. Discoveries in this field could be greatly accelerated by computational methods and resources that enable these findings to be compared across different studies and different organisms and linked with the genes responsible for anatomical modifications. Homology is a key concept in comparative anatomy; two important types are historical homology (the similarity of organisms due to common ancestry) and serial homology (the similarity of repeated structures within an organism). We explored how to most effectively represent historical and serial homology across anatomical structures to facilitate computational reasoning. We assembled a collection of homology assertions from the literature with a set of taxon phenotypes for the skeletal elements of vertebrate fins and limbs from the Phenoscape Knowledgebase. Using seven competency questions, we evaluated the reasoning ramifications of two logical models: the Reciprocal Existential Axioms (REA) homology model and the Ancestral Value Axioms (AVA) homology model. The AVA model returned all user-expected results in addition to the search term and any of its subclasses. The AVA model also returns any superclass of the query term in which a homology relationship has been asserted. The REA model returned the user-expected results for five out of seven queries. We identify some challenges of implementing complete homology queries due to limitations of OWL reasoning. This work lays the foundation for homology reasoning to be incorporated into other ontology-based tools, such as those that enable synthetic supermatrix construction and candidate gene discovery. [Homology; ontology; anatomy; morphology; evolution; knowledgebase; phenoscape.].

Hybridization Chain Reaction for Quantitative and Multiplex Imaging of Gene Expression in Amphioxus Embryos and Adult Tissues

In situ hybridization (ISH) methods remain the most popular approach for profiling the expression of a gene at high spatial resolution and have been broadly used to address many biological questions. One compelling application is in the field of evo-devo, where comparing gene expression patterns has offered insight into how vertebrate development has evolved. Gene expression profiling in the invertebrate chordate amphioxus (cephalochordate) has been particularly instrumental in this context: its key phylogenetic position as sister group to all other chordates makes it an ideal model system to compare with vertebrates and for reconstructing the ancestral condition of our phylum. However, while ISH methods have been developed extensively in vertebrate model systems to fluorescently detect the expression of multiple genes simultaneously at a cellular and subcellular resolution, amphioxus gene expression profiling is still based on single-gene nonfluorescent chromogenic methods, whose spatial resolution is often compromised by diffusion of the chromogenic product. This represents a serious limitation for reconciling gene expression dynamics between amphioxus and vertebrates and for molecularly identifying cell types, defined by their combinatorial code of gene expression, that may have played pivotal roles in evolutionary innovation. Herein we overcome these problems by describing a new protocol for application of the third-generation hybridization chain reaction (HCR) to the amphioxus, which permits fluorescent, multiplex, and quantitative detection of gene expression in situ, within the changing morphology of the developing embryo, and in adult tissues. A detailed protocol is herein provided for whole-mount preparations of embryos and vibratome sections of adult tissues.

Development, repair, and regeneration of the limb musculoskeletal system

The limb musculoskeletal system provides a primary means for locomotion, manipulation of objects and protection for most vertebrate organisms. Intricate integration of the bone, tendon and muscle tissues are required for function. These three tissues arise largely independent of one another, but the connections formed during later development are maintained throughout life and are re-established following injury. Each of these tissues also have mesenchymal stem/progenitor cells that function in maintenance and repair. Here in, we will review the major events in the development of limb skeleton, tendon, and muscle tissues, their response to injury, and discuss current knowledge regarding resident progenitor/stem cells within each tissue that participate in development, repair, and regeneration in vivo.

Similarities and differences in the regulation of HoxD genes during chick and mouse limb development

In all tetrapods examined thus far, the development and patterning of limbs require the activation of gene members of the HoxD cluster. In mammals, they are regulated by a complex bimodal process that controls first the proximal patterning and then the distal structure. During the shift from the former to the latter regulation, this bimodal regulatory mechanism allows the production of a domain with low Hoxd gene expression, at which both telomeric (T-DOM) and centromeric regulatory domains (C-DOM) are silent. These cells generate the future wrist and ankle articulations. We analyzed the implementation of this regulatory mechanism in chicken, i.e., in an animal for which large morphological differences exist between fore- and hindlimbs. We report that although this bimodal regulation is globally conserved between the mouse and the chick, some important modifications evolved at least between these two model systems, in particular regarding the activity of specific enhancers, the width of the TAD boundary separating the two regulations, and the comparison between the forelimb versus hindlimb regulatory controls. At least one aspect of these regulations seems to be more conserved between chick and bats than with mouse, which may relate to the extent to which forelimbs and hindlimbs of these various animals differ in their morphologies.

A comparison of Hox gene regulation during the development of limbs in birds and mammals reveals that whereas the characteristic bimodal regulatory system, based on large chromatin domains, is largely conserved between these morphologically distinct structures, some differences are revealed in the way this is implemented in various vertebrates.

The shapes of limbs vary greatly among tetrapod species, even between the forelimbs and hindlimbs of the same animal. Hox genes regulate the proper growth and patterning of tetrapod limbs. In order to evaluate whether variations in the complex regulation of a cluster of Hox genes—the Hoxd genes—during limb development contribute to the differences in limb shape, we compared their transcriptional control during limb bud development in the forelimbs and hindlimbs of mouse and chicken embryos. We found that the regulatory mechanism underlying Hoxd gene expression is highly conserved, but some clear differences exist. For instance, we observed a variation in the topologically associating domain (TAD; a self-interacting genomic region) boundary interval between the mouse and the chick, as well as differences in the activity of a conserved enhancer element situated within the telomeric regulatory domain. In contrast to the mouse, the chicken enhancer has a stronger activity in the forelimb buds than in the hindlimb buds, which is correlated with the striking differences in the mRNA levels of the genes. We conclude that differences in both the timing and duration of TAD activities and in the width of their boundary may parallel the important decrease in Hoxd gene transcription in chick hindlimb buds versus forelimb buds. These differences may also account for the slightly distinct regulatory strategies implemented by mammals and birds at this locus.

The microRNA-23a cluster regulates the developmental HoxA cluster function during osteoblast differentiation

MicroRNAs (miRs) and Hox transcription factors have decisive roles in postnatal bone formation and homeostasis. In silico analysis identified extensive interaction between HOXA cluster mRNA and microRNAs from the miR-23a cluster. However, Hox regulation by the miR-23a cluster during osteoblast differentiation remains undefined. We examined this regulation in preosteoblasts and in a novel miR-23a cluster knockdown mouse model. Overexpression and knockdown of the miR-23a cluster in preosteoblasts decreased and increased, respectively, the expression of the proteins HOXA5, HOXA10, and HOXA11; these proteins' mRNAs exhibited significant binding with the miR-23a cluster miRNAs, and miRNA 3'-UTR reporter assays confirmed repression. Importantly, during periods correlating with development and differentiation of bone cells, we found an inverse pattern of expression between HoxA factors and members of the miR-23a cluster. HOXA5 and HOXA11 bound to bone-specific promoters, physically interacted with transcription factor RUNX2, and regulated bone-specific genes. Depletion of HOXA5 or HOXA11 in preosteoblasts also decreased cellular differentiation. Additionally, stable overexpression of the miR-23a cluster in osteoblasts decreased the recruitment of HOXA5 and HOXA11 to osteoblast gene promoters, significantly inhibiting histone H3 acetylation. Heterozygous miR-23a cluster knockdown female mice (miR-23a Cl^WT/ZIP) had significantly increased trabecular bone mass when compared with WT mice. Furthermore, miR-23a cluster knockdown in calvarial osteoblasts of these mice increased the recruitment of HOXA5 and HOXA11, with a substantial enrichment of promoter histone H3 acetylation. Taken together, these findings demonstrate that the miR-23a cluster is required for maintaining stage-specific HoxA factor expression during osteogenesis.

Heterogeneous combinatorial expression of Hoxd genes in single cells during limb development

BACKGROUND

Global analyses of gene expression during development reveal specific transcription patterns associated with the emergence of various cell types, tissues, and organs. These heterogeneous patterns are instrumental to ensure the proper formation of the different parts of our body, as shown by the phenotypic effects generated by functional genetic approaches. However, variations at the cellular level can be observed within each structure or organ. In the developing mammalian limbs, expression of Hox genes from the HoxD cluster is differentially controlled in space and time, in cells that will pattern the digits and the forearms. While the Hoxd genes broadly share a common regulatory landscape and large-scale analyses have suggested a homogenous Hox gene transcriptional program, it has not previously been clear whether Hoxd genes are expressed together at the same levels in the same cells.

RESULTS

We report a high degree of heterogeneity in the expression of the Hoxd11 and Hoxd13 genes. We analyzed single-limb bud cell transcriptomes and show that Hox genes are expressed in specific combinations that appear to match particular cell types. In cells giving rise to digits, we find that the expression of the five relevant Hoxd genes (Hoxd9 to Hoxd13) is unbalanced, despite their control by known global enhancers. We also report that specific combinatorial expression follows a pseudo-time sequence, which is established based on the transcriptional diversity of limb progenitors.

CONCLUSIONS

Our observations reveal the existence of distinct combinations of Hoxd genes at the single-cell level during limb development. In addition, we document that the increasing combinatorial expression of Hoxd genes in this developing structure is associated with specific transcriptional signatures and that these signatures illustrate a temporal progression in the differentiation of these cells.

Expression of meis and hoxa11 in dipnoan and teleost fins provides new insights into the evolution of vertebrate appendages

Background

The concerted activity of Meis and Hoxa11 transcription factors is essential for the subdivision of tetrapod limbs into proximo-distal (PD) domains; however, little is know about the evolution of this patterning mechanism. Here, we aim to study the expression of meis and hoxa11 orthologues in the median and paired rayed fins of zebrafish and in the lobed fins of the Australian lungfish.

Results

First, a late phase of expression of meis1.1 and hoxa11b in zebrafish dorsal and anal fins relates with segmentation of endochondral elements in proximal and distal radials. Second, our zebrafish in situ hybridization results reveal spatial and temporal changes between pectoral and pelvic fins. Third, in situ analysis of meis1, meis3 and hoxa11 genes in Neoceratodus pectoral fins identifies decoupled domains of expression along the PD axis.

Conclusions

Our data raise the possibility that the origin of stylopod and zeugopod lies much deeper in gnathostome evolution and that variation in meis and hoxa11 expression has played a substantial role in the transformation of appendage anatomy. Moreover, these observations provide evidence that the Meis/Hoxa11 profile considered a hallmark of stylopod/zeugopod patterning is present in Neoceratodus.

Electronic supplementary material

The online version of this article (10.1186/s13227-018-0099-9) contains supplementary material, which is available to authorized users.

Flight feather development: its early specialization during embryogenesis

BACKGROUND

Flight feathers, a type of feather that is unique to extant/extinct birds and some non-avian dinosaurs, are the most evolutionally advanced type of feather. In general, feather types are formed in the second or later generation of feathers at the first and following molting, and the first molting begins at around two weeks post hatching in chicken. However, it has been stated in some previous reports that the first molting from the natal down feathers to the flight feathers is much earlier than that for other feather types, suggesting that flight feather formation starts as an embryonic event. The aim of this study was to determine the inception of flight feather morphogenesis and to identify embryological processes specific to flight feathers in contrast to those of down feathers.

RESULTS

We found that the second generation of feather that shows a flight feather-type arrangement has already started developing by chick embryonic day 18, deep in the skin of the flight feather-forming region. This was confirmed by shh gene expression that shows barb pattern, and the expression pattern revealed that the second generation of feather development in the flight feather-forming region seems to start by embryonic day 14. The first stage at which we detected a specific morphology of the feather bud in the flight feather-forming region was embryonic day 11, when internal invagination of the feather bud starts, while the external morphology of the feather bud is radial down-type.

CONCLUSION

The morphogenesis for the flight feather, the most advanced type of feather, has been drastically modified from the beginning of feather morphogenesis, suggesting that early modification of the embryonic morphogenetic process may have played a crucial role in the morphological evolution of this key innovation. Co-optation of molecular cues for axial morphogenesis in limb skeletal development may be able to modify morphogenesis of the feather bud, giving rise to flight feather-specific morphogenesis of traits.

The role of HoxA11 and HoxA13 in the evolution of novel fin morphologies in a representative batoid (Leucoraja erinacea)

Background

Batoids exhibit unique body plans with derived fin morphologies, such as the anteriorly expanded pectoral fins that fuse to the head, or distally extended anterior pelvic fin lobes used for a modified swimming technique utilized by skates (Rajidae). The little skate (Leucoraja erinacea), exhibits both of these unique fin morphologies. These fin modifications are not present in a typical shark body plan, and little is known regarding the mechanisms underlying their development. A recent study identified a novel apical ectodermal ridge (AER) associated with the development of the anterior pectoral fin in the little skate, but the role of the posterior HoxA genes was not featured during skate fin development.

Results

We present the first evidence for HoxA expression (HoxA11 and HoxA13) in novel AER domains associated with the development of three novel fin morphologies in a representative batoid, L. erinacea. We found HoxA13 expression associated with the recently described novel AER in the anterior pectoral fin, and HoxA11 expression in a novel AER domain in the anterior pelvic fin that we describe here. We find that both HoxA11 and HoxA13 are expressed in claspers, and while HoxA11 is expressed in pelvic fins and claspers, HoxA13 is expressed exclusively in developing claspers of males. Finally, HoxA11 expression is associated with the developing fin rays in paired fins.

Conclusion

Overall, these results indicate that the posterior HoxA genes play an important role in the morphological evolution of paired fins in a representative batoid. These data suggest that the batoids utilize a unique Hox code, where the posterior HoxA genes exhibit distinct expression patterns that are likely associated with specification of novel fin morphologies.

Electronic supplementary material

The online version of this article (10.1186/s13227-017-0088-4) contains supplementary material, which is available to authorized users.

Embryonic timing, axial stem cells, chromatin dynamics, and the Hox clock

Collinear regulation of Hox genes in space and time has been an outstanding question ever since the initial work of Ed Lewis in 1978. Here we discuss recent advances in our understanding of this phenomenon in relation to novel concepts associated with large-scale regulation and chromatin structure during the development of both axial and limb patterns. We further discuss how this sequential transcriptional activation marks embryonic stem cell-like axial progenitors in mammals and, consequently, how a temporal genetic system is further translated into spatial coordinates via the fate of these progenitors. In this context, we argue the benefit and necessity of implementing this unique mechanism as well as the difficulty in evolving an alternative strategy to deliver this critical positional information.

Alterations in anterior-posterior patterning and its accompanying changes along the proximal-distal axis during the fin-to-limb transition

The evolution from fins to limbs was one of the most successful innovations for vertebrates, allowing them to vastly expand their behaviors and habitats. Fossil records suggest that morphological changes occurred not only along the proximal-distal axis included appearance of the autopod, but also occurred along the anterior-posterior axis included reductions in the size and number of basal bones and digits. This review focuses on recent progress in developmental and genetic studies aimed at elucidating the mechanisms underlying alteration of anterior-posterior patterning and its accompanying changes along the proximal-distal axis during the fin-to-limb transition.

Preaxial polydactyly in an individual with Wiedemann-Steiner syndrome caused by a novel nonsense mutation in KMT2A

Wiedemann-Steiner syndrome (WDSTS) is an autosomal dominant disorder characterized by hypertrichosis, intellectual disability, and dysmorphic facial appearances (down-slanted vertically narrow palpebral fissures, wide nasal bridge, broad nasal tip, and thick eyebrows). In 2012, Jones and co-workers identified heterozygous mutations in KMT2A (lysine methyltransferase 2A) as the molecular cause of WDSTS. Although the phenotype of this syndrome continues to expand, the associated features are not fully understood. Here, we report WDSTS in a 12-year-old Japanese boy with a novel nonsense mutation in KMT2A. He had right preaxial polydactyly, which has not been previously reported in WDSTS. We could not identify a causal relationship between the KMT2A mutation and preaxial polydactyly, and cannot exclude the preaxial polydactyly is a simple coincidence. We summarized the clinical features of WDSTS associated with KMT2A mutation and discussed the cardinal symptoms in detail.

Pitx1 directly modulates the core limb development program to implement hindlimb identity

Forelimbs (FLs) and hindlimbs (HLs) develop complex musculoskeletal structures that rely on the deployment of a conserved developmental program. Pitx1, a transcription factor gene with expression restricted to HL and absent from FL, plays an important role in generating HL features. The genomic mechanisms by which Pitx1 effects HL identity remain poorly understood. Here, we use expression profiling and analysis of direct Pitx1 targets to characterize the HL- and FL-restricted genetic programs in mouse and situate the Pitx1-dependent gene network within the context of limb-specific gene regulation. We show that Pitx1 is a crucial component of a narrow network of HL-restricted regulators, acting on a developmental program that is shared between FL and HL. Pitx1 targets sites that are in a similar chromatin state in FL and HL and controls expression of patterning genes as well as the chondrogenic program, consistent with impaired chondrogenesis in Pitx1^-/- HL. These findings support a model in which multifactorial actions of a limited number of HL regulators redirect the generic limb development program in order to generate the unique structural features of the limb.

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.

The 2017 ABJS Nicolas Andry Award: Advancing Personalized Medicine for Clubfoot Through Translational Research

BACKGROUND

Clubfoot is one of the most common pediatric orthopaedic disorders. While the Ponseti method has revolutionized clubfoot treatment, it is not effective for all patients. When the Ponseti method does not correct the foot, patients are at risk for lifelong disability and may require more-extensive surgery.

QUESTIONS/PURPOSES

(1) What genetic and morphologic abnormalities contribute to the development of clubfoot? (2) How can this information be used to devise personalized treatment paradigms for patients with clubfoot?

METHODS

Human gene sequencing, molecular genetic engineering of mouse models of clubfoot, MRI of clubfoot, and development of new treatment methods all have been used by our group to understand the biological basis and improve therapy for this group of disorders.

RESULTS

We gained new insight into clubfoot pathogenesis from our discovery that mutations in the PITX1-TBX4-HOXC transcriptional pathway cause familial clubfoot and vertical talus in a small number of families, with the unique lower limb expression of these genes providing an explanation for the lack of upper extremity involvement in these disorders. MRI studies revealed corresponding morphologic abnormalities, including hypomorphic muscle, bone, and vasculature, that are not only associated with these gene mutations, but also are biomarkers for treatment-resistant clubfoot.

CONCLUSIONS

Based on an understanding of the underlying biology, we improved treatment methods for neglected and syndromic clubfoot, developed new treatment for congenital vertical talus based on the principles of the Ponseti method, and designed a new dynamic clubfoot brace to improve strength and compliance.

Dorsal fin development in flounder, Paralichthys olivaceus: Bud formation and its cellular origin

The development of the median fin has not been investigated extensively in teleosts, although in other fishes it has been proposed that it involves the same genetic programs operating in the paired appendages. Adult median fins develop from the larval bud; therefore an investigation of fin bud formation and its cellular origin is essential to understanding the maturation mechanisms. In Paralichthys olivaceus, skeletogenesis proceeds from an anterior to posterior direction providing a good opportunity to study the formation of dorsal fin bud. An apical ectodermal ridge appeared at the basal stratum of the presumptive dorsal fin was first observed at 3 days post hatching. Then the apical ectodermal fold formed as the bud outgrew in 6 days post-hatch larvae. The bud continued to grow, breaking through the dorsal fin fold in 9 days post-hatch larvae. At 13 days post-hatch, the bud grew beyond the edge of the fin fold and formed into the four future rays. Molecular markers of cell type showed the existence of neural crest cells, scleroblasts and sclerotomes in the dorsal fin bud. The earliest gene expression in the dorsal fin bud was Hoxd10 at 3 days post-hatch larvae, then Hoxd9, Hoxd11 and Hoxd12. This indicates Hoxd10 might be a candidate molecular marker of the bud formation site. Some key molecular markers for paired appendage development, such as FGF8, Wnt7, and Shh were expressed at the apical ectodermal ridge and later the apical ectodermal fold. Moreover, the form of the dorsal fin bud could be inhibited by Hh pathway inhibitor, further indicating that common basic molecular mechanisms might be utilized by median fins.

Upstream regulation for initiation of restricted Shh expression in the chick limb bud

BACKGROUND

The organizing center, which serves as a morphogen source, has crucial functions in morphogenesis in animal development. The center is necessarily located in a certain restricted area in the morphogenetic field, and there are several ways in which an organizing center can be restricted. The organizing center for limb morphogenesis, the ZPA (zone of polarizing activity), specifically expresses the Shh gene and is restricted to the posterior region of the developing limb bud.

RESULTS

The pre-pattern along the limb anteroposterior axis, provided by anterior Gli3 expression and posterior Hand2 expression, seems insufficient for the initiation of Shh expression restricted to a narrow, small spot in the posterior limb field. Comparison of the spatiotemporal patterns of gene expression between Shh and some candidate genes (Fgf8, Hoxd10, Hoxd11, Tbx2, and Alx4) upstream of Shh expression suggested that a combination of these genes' expression provides the restricted initiation of Shh expression.

CONCLUSIONS

Taken together with results of functional assays, we propose a model in which positive and negative transcriptional regulatory networks accumulate their functions in the intersection area of their expression regions to provide a restricted spot for the ZPA, the source of morphogen, Shh. Developmental Dynamics 246:417-430, 2017. © 2017 Wiley Periodicals, Inc.

Integration of Shh and Fgf signaling in controlling Hox gene expression in cultured limb cells

During embryonic development, fields of progenitor cells form complex structures through dynamic interactions with external signaling molecules. How complex signaling inputs are integrated to yield appropriate gene expression responses is poorly understood. In the early limb bud, for instance, Sonic hedgehog (Shh) is expressed in the distal posterior mesenchyme, where it acts as a mediator of anterior to posterior (AP) patterning, whereas fibroblast growth factor 8 (Fgf8) is produced by the apical ectodermal ridge (AER) at the distal tip of the limb bud to direct outgrowth along the proximal to distal (PD) axis. Here we use cultured limb mesenchyme cells to assess the response of the target Hoxd genes to these two factors. We find that they act synergistically and that both factors are required to activate Hoxd13 in limb mesenchymal cells. However, the analysis of the enhancer landscapes flanking the HoxD cluster reveals that the bimodal regulatory switch observed in vivo is only partially achieved under these in vitro conditions, suggesting an additional requirement for other factors.

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.

The origins, scaling and loss of tetrapod digits

Many of the great morphologists of the nineteenth century marvelled at similarities between the limbs of diverse species, and Charles Darwin noted these homologies as significant supporting evidence for descent with modification from a common ancestor. Sir Richard Owen also took great care to highlight each of the elements of the forelimb and hindlimb in a multitude of species with focused attention on the homology between the hoof of the horse and the middle digit of man. The ensuing decades brought about a convergence of palaeontology, experimental embryology and molecular biology to lend further support to the homologies of tetrapod limbs and their developmental origins. However, for all that we now understand about the conserved mechanisms of limb development and the development of gross morphological disturbances, little of what is presented in the experimental or medical literature reflects the remarkable diversity resulting from the 450 million year experiment of natural selection. An understanding of conserved and divergent limb morphologies in this new age of genomics and genome engineering promises to reveal more of the developmental potential residing in all limbs and to unravel the mechanisms of evolutionary variation in limb size and shape. In this review, we present the current state of our rapidly advancing understanding of the evolutionary origin of hands and feet and highlight what is known about the mechanisms that shape diverse limbs.This article is part of the themed issue 'Evo-devo in the genomics era, and the origins of morphological diversity'.

Intrinsic properties of limb bud cells can be differentially reset

An intrinsic timing mechanism specifies the positional values of the zeugopod (i.e. radius/ulna) and then autopod (i.e. wrist/digits) segments during limb development. Here, we have addressed whether this timing mechanism ensures that patterning events occur only once by grafting GFP-expressing autopod progenitor cells to the earlier host signalling environment of zeugopod progenitor cells. We show by detecting Hoxa13 expression that early and late autopod progenitors fated for the wrist and phalanges, respectively, both contribute to the entire host autopod, indicating that the autopod positional value is irreversibly determined. We provide evidence that Hoxa13 provides an autopod-specific positional value that correctly allocates cells into the autopod, most likely through the control of cell-surface properties as shown by cell-cell sorting analyses. However, we demonstrate that only the earlier autopod cells can adopt the host proliferation rate to permit normal morphogenesis. Therefore, our findings reveal that the ability of embryonic cells to differentially reset their intrinsic behaviours confers robustness to limb morphogenesis. We speculate that this plasticity could be maintained beyond embryogenesis in limbs with regenerative capacity.

Genome-Wide Binding of Posterior HOXA/D Transcription Factors Reveals Subgrouping and Association with CTCF

Homeotic genes code for key transcription factors (HOX-TFs) that pattern the animal body plan. During embryonic development, Hox genes are expressed in overlapping patterns and function in a partially redundant manner. In vitro biochemical screens probing the HOX-TF sequence specificity revealed largely overlapping sequence preferences, indicating that co-factors might modulate the biological function of HOX-TFs. However, due to their overlapping expression pattern, high protein homology, and insufficiently specific antibodies, little is known about their genome-wide binding preferences. In order to overcome this problem, we virally expressed tagged versions of limb-expressed posterior HOX genes (HOXA9-13, and HOXD9-13) in primary chicken mesenchymal limb progenitor cells (micromass). We determined the effect of each HOX-TF on cellular differentiation (chondrogenesis) and gene expression and found that groups of HOX-TFs induce distinct regulatory programs. We used ChIP-seq to determine their individual genome-wide binding profiles and identified between 12,721 and 28,572 binding sites for each of the nine HOX-TFs. Principal Component Analysis (PCA) of binding profiles revealed that the HOX-TFs are clustered in two subgroups (Group 1: HOXA/D9, HOXA/D10, HOXD12, and HOXA13 and Group 2: HOXA/D11 and HOXD13), which are characterized by differences in their sequence specificity and by the presence of cofactor motifs. Specifically, we identified CTCF binding sites in Group 1, indicating that this subgroup of HOX-proteins cooperates with CTCF. We confirmed this interaction by an independent biological assay (Proximity Ligation Assay) and demonstrated that CTCF is a novel HOX cofactor that specifically associates with Group 1 HOX-TFs, pointing towards a possible interplay between HOX-TFs and chromatin architecture.

Hox genes encode transcription factors that determine the vertebrate body plan and pattern structures and organs in the developing embryo. Despite decades of effort and research on Hox genes, little is known about the HOX-DNA binding properties in vivo. This lack of knowledge is mainly due to the absence of appropriate antibodies to distinguish between different HOX transcription factors. Here, we adapt a cell culture system that allows us to investigate HOX-DNA binding on a genome-wide scale. With this approach, we define and directly compare the genome-wide binding sites of nine posterior HOXA and HOXD transcription factors. We report that the in vivo HOX binding specificity differs from the in vitro specificity and find that HOX-TFs largely rely on co-factor binding and not only on direct HOX-DNA binding. Finally, we identify a novel HOX co-factor, a genome architecture protein, CTCF suggesting a possible interplay between HOX-TF function and chromatin architecture.

Inactivation of Sonic Hedgehog Signaling and Polydactyly in Limbs of Hereditary Multiple Malformation, a Novel Type of Talpid Mutant

Hereditary Multiple Malformation (HMM) is a naturally occurring, autosomal recessive, homozygous lethal mutation found in Japanese quail. Homozygote embryos (hmm^−/−) show polydactyly similar to talpid² and talpid³ mutants. Here we characterize the molecular profile of the hmm^−/− limb bud and identify the cellular mechanisms that cause its polydactyly. The hmm^−/− limb bud shows a severe lack of sonic hedgehog (SHH) signaling, and the autopod has 4 to 11 unidentifiable digits with syn-, poly-, and brachydactyly. The Zone of Polarizing Activity (ZPA) of the hmm^−/− limb bud does not show polarizing activity regardless of the presence of SHH protein, indicating that either the secretion pathway of SHH is defective or the SHH protein is dysfunctional. Furthermore, mesenchymal cells in the hmm^−/− limb bud do not respond to ZPA transplanted from the normal limb bud, suggesting that signal transduction downstream of SHH is also defective. Since primary cilia are present in the hmm^−/− limb bud, the causal gene must be different from talpid² and talpid³. In the hmm^−/− limb bud, a high amount of GLI3A protein is expressed and GLI3 protein is localized to the nucleus. Our results suggest that the regulatory mechanism of GLI3 is disorganized in the hmm^−/− limb bud.

Saunders's framework for understanding limb development as a platform for investigating limb evolution

John W. Saunders, Jr. made seminal discoveries unveiling how chick embryos develop their limbs. He discovered the apical ectodermal ridge (AER), the zone of polarizing activity (ZPA), and the domains of interdigital cell death within the developing limb and determined their function through experimental analysis. These discoveries provided the basis for subsequent molecular understanding of how vertebrate limbs are induced, patterned, and differentiated. These mechanisms are strongly conserved among the vast diversity of tetrapod limbs suggesting that relatively minor changes and tweaks to the molecular cascades are responsible for the diversity observed in nature. Analysis of the pathway systems first identified by Saunders in the context of animals displaying limb reduction show how alterations in these pathways have resulted in multiple mechanisms of limb and digit loss. Other classes of modification to these same patterning systems are seen at the root of other, novel limb morphological alterations and elaborations.

Critical Timing without a Timer for Embryonic Development

Timing of embryonic development is precisely controlled, but the mechanisms underlying biological timers are still unclear. Here, a validated model for timing under control of Sonic Hedgehog is revisited and generalized to an arbitrary number of genes. The developmental dynamics where a temporal sequence of gene expression recapitulates a steady-state spatial pattern can be realized through a simple network close to criticality, controlled by the duration of exposure to a morphogen. Criticality simultaneously accounts for many observed biological properties, such as timing, multistability, and canalization of genetic expression. This process can be parsimoniously generalized in many dimensions with a minimum number of genes, all repressing each other with asymmetrical strengths, which also explains sequential activation of different fates. Separation of timescales allows for a simple analytical interpretation. Finally, it is shown that even in the presence of noise, coupling between cells preserves criticality and robust patterning. The model offers a simple theoretical framework for the study of emergent developmental timers.

Positional information from oscillatory phase shifts : insights from in silico evolution

Complex cellular decisions are based on temporal dynamics of pathways, including genetic oscillators. In development, recent works on vertebrae formation have suggested that relative phase of genetic oscillators encode positional information, including differentiation front defining vertebrae positions. Precise mechanisms for this are still unknown. Here, we use computational evolution to find gene network topologies that can compute the phase difference between oscillators and convert it into a decoder morphogen concentration. Two types of networks are discovered, based on symmetry properties of the decoder gene. So called asymmetric networks are studied, and two submodules are identified converting phase information into an amplitude variable. Those networks naturally display a 'shock' for a well defined phase difference, that can be used to define a wavefront of differentiation. We show how implementation of these ideas reproduce experimental features of vertebrate segmentation.

A role for HOX13 proteins in the regulatory switch between TADs at the HoxD locus

During vertebrate limb development, Hoxd genes are regulated following a bimodal strategy involving two topologically associating domains (TADs) located on either side of the gene cluster. These regulatory landscapes alternatively control different subsets of Hoxd targets, first into the arm and subsequently into the digits. We studied the transition between these two global regulations, a switch that correlates with the positioning of the wrist, which articulates these two main limb segments. We show that the HOX13 proteins themselves help switch off the telomeric TAD, likely through a global repressive mechanism. At the same time, they directly interact with distal enhancers to sustain the activity of the centromeric TAD, thus explaining both the sequential and exclusive operating processes of these two regulatory domains. We propose a model in which the activation of Hox13 gene expression in distal limb cells both interrupts the proximal Hox gene regulation and re-enforces the distal regulation. In the absence of HOX13 proteins, a proximal limb structure grows without any sign of wrist articulation, likely related to an ancestral fish-like condition.

Developmental Mechanism of Limb Field Specification along the Anterior-Posterior Axis during Vertebrate Evolution

In gnathostomes, limb buds arise from the lateral plate mesoderm at discrete positions along the body axis. Specification of these limb-forming fields can be subdivided into several steps. The lateral plate mesoderm is regionalized into the anterior lateral plate mesoderm (ALPM; cardiac mesoderm) and the posterior lateral plate mesoderm (PLPM). Subsequently, Hox genes appear in a nested fashion in the PLPM and provide positional information along the body axis. The lateral plate mesoderm then splits into the somatic and splanchnic layers. In the somatic layer of the PLPM, the expression of limb initiation genes appears in the limb-forming region, leading to limb bud initiation. Furthermore, past and current work in limbless amphioxus and lampreys suggests that evolutionary changes in developmental programs occurred during the acquisition of paired fins during vertebrate evolution. This review presents these recent advances and discusses the mechanisms of limb field specification during development and evolution, with a focus on the role of Hox genes in this process.

HoxD expression in the fin-fold compartment of basal gnathostomes and implications for paired appendage evolution

The role of Homeobox transcription factors during fin and limb development have been the focus of recent work investigating the evolutionary origin of limb-specific morphologies. Here we characterize the expression of HoxD genes, as well as the cluster-associated genes Evx2 and LNP, in the paddlefish Polyodon spathula, a basal ray-finned fish. Our results demonstrate a collinear pattern of nesting in early fin buds that includes HoxD14, a gene previously thought to be isolated from global Hox regulation. We also show that in both Polyodon and the catshark Scyliorhinus canicula (a representative chondrichthyan) late phase HoxD transcripts are present in cells of the fin-fold and co-localize with And1, a component of the dermal skeleton. These new data support an ancestral role for HoxD genes in patterning the fin-folds of jawed vertebrates, and fuel new hypotheses about the evolution of cluster regulation and the potential downstream differentiation outcomes of distinct HoxD-regulated compartments.