Analysis of Hox gene expression in the chick limb bud

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.

HoxA Genes and the Fin-to-Limb Transition in Vertebrates

HoxA genes encode for important DNA-binding transcription factors that act during limb development, regulating primarily gene expression and, consequently, morphogenesis and skeletal differentiation. Within these genes, HoxA11 and HoxA13 were proposed to have played an essential role in the enigmatic evolutionary transition from fish fins to tetrapod limbs. Indeed, comparative gene expression analyses led to the suggestion that changes in their regulation might have been essential for the diversification of vertebrates' appendages. In this review, we highlight three potential modifications in the regulation and function of these genes that may have boosted appendage evolution: (1) the expansion of polyalanine repeats in the HoxA11 and HoxA13 proteins; (2) the origin of +a novel long-non-coding RNA with a possible inhibitory function on HoxA11; and (3) the acquisition of cis-regulatory elements modulating 5' HoxA transcription. We discuss the relevance of these mechanisms for appendage diversification reviewing the current state of the art and performing additional comparative analyses to characterize, in a phylogenetic framework, HoxA11 and HoxA13 expression, alanine composition within the encoded proteins, long-non-coding RNAs and cis-regulatory elements.

How the embryo makes a limb: determination, polarity and identity

The vertebrate limb with its complex anatomy develops from a small bud of undifferentiated mesoderm cells encased in ectoderm. The bud has its own intrinsic polarity and can develop autonomously into a limb without reference to the rest of the embryo. In this review, recent advances are integrated with classical embryology, carried out mainly in chick embryos, to present an overview of how the embryo makes a limb bud. We will focus on how mesoderm cells in precise locations in the embryo become determined to form a limb and express the key transcription factors Tbx4 (leg/hindlimb) or Tbx5 (wing/forelimb). These Tbx transcription factors have equivalent functions in the control of bud formation by initiating a signalling cascade involving Wnts and fibroblast growth factors (FGFs) and by regulating recruitment of mesenchymal cells from the coelomic epithelium into the bud. The mesoderm that will form limb buds and the polarity of the buds is determined with respect to both antero-posterior and dorso-ventral axes of the body. The position in which a bud develops along the antero-posterior axis of the body will also determine its identity - wing/forelimb or leg/hindlimb. Hox gene activity, under the influence of retinoic acid signalling, is directly linked with the initiation of Tbx5 gene expression in the region along the antero-posterior axis of the body that will form wings/forelimbs and determines antero-posterior polarity of the buds. In contrast, Tbx4 expression in the regions that will form legs/hindlimbs is regulated by the homeoprotein Pitx1 and there is no evidence that Hox genes determine antero-posterior polarity of the buds. Bone morphogenetic protein (BMP) signalling determines the region along the dorso-ventral axis of the body in which both wings/forelimbs and legs/hindlimbs develop and dorso-ventral polarity of the buds. The polarity of the buds leads to the establishment of signalling regions - the dorsal and ventral ectoderm, producing Wnts and BMPs, respectively, the apical ectodermal ridge producing fibroblast growth factors and the polarizing region, Sonic hedgehog (Shh). These signals are the same in both wings/forelimbs and legs/hindlimbs and control growth and pattern formation by providing the mesoderm cells of the limb bud as it develops with positional information. The precise anatomy of the limb depends on the mesoderm cells in the developing bud interpreting positional information according to their identity - determined by Pitx1 in hindlimbs - and genotype. The competence to form a limb extends along the entire antero-posterior axis of the trunk - with Hox gene activity inhibiting the formation of forelimbs in the interlimb region - and also along the dorso-ventral axis.

Epigenetic mechanisms and boundaries in the regulation of mammalian Hox clusters

Hox gene expression imparts segment identity to body structures along the anterior-posterior axis and is tightly governed by higher order chromatin mechanisms. Chromatin regulatory features of the homeotic complex are best defined in Drosophila melanogaster, where multiple cis-regulatory elements have been identified that ensure collinear Hox gene expression patterns in accordance with their genomic organization. Recent studies focused on delineating the epigenetic features of the vertebrate Hox clusters have helped reveal their dynamic chromatin organization and its impact on gene expression. Enrichment for the 'activating' H3K4me3 and 'repressive' H3K27me3 histone modifications is a particularly strong read-out for transcriptional status and correlates well with the evidence for chromatin loop domain structures and stage specific topological changes at these loci. However, it is not clear how such distinct domains are imposed and regulated independent of each other. Comparative analysis of the chromatin structure and organization of the homeotic gene clusters in fly and mammals is increasingly revealing the functional conservation of chromatin mediated mechanisms. Here we discuss the case for interspersed boundary elements existing within mammalian Hox clusters along with their possible roles and mechanisms of action. Recent studies suggest a role for factors other than the well characterized vertebrate boundary factor CTCF, such as the GAGA binding factor (GAF), in maintaining chromatin domains at the Hox loci. We also present data demonstrating how such regulatory elements may be involved in organizing higher order structure and demarcating active domains of gene expression at the mammalian Hox clusters.

Data-driven modelling of a gene regulatory network for cell fate decisions in the growing limb bud

Parameter optimization coupled with model selection is a convenient approach to infer gene regulatory networks from experimental gene expression data, but so far it has been limited to single cells or static tissues where growth is not significant. Here, we present a computational study in which we determine an optimal gene regulatory network from the spatiotemporal dynamics of gene expression patterns in a complex 2D growing tissue (non-isotropic and heterogeneous growth rates). We use this method to predict the regulatory mechanisms that underlie proximodistal (PD) patterning of the developing limb bud. First, we map the expression patterns of the PD markers Meis1, Hoxa11 and Hoxa13 into a dynamic description of the tissue movements that drive limb morphogenesis. Secondly, we use reverse-engineering to test how different gene regulatory networks can interpret the opposing gradients of fibroblast growth factors (FGF) and retinoic acid (RA) to pattern the PD markers. Finally, we validate and extend the best model against various previously published manipulative experiments, including exogenous application of RA, surgical removal of the FGF source and genetic ectopic expression of Meis1. Our approach identifies the most parsimonious gene regulatory network that can correctly pattern the PD markers downstream of FGF and RA. This network reveals a new model of PD regulation which we call the “crossover model”, because the proximal morphogen (RA) controls the distal boundary of Hoxa11, while conversely the distal morphogens (FGFs) control the proximal boundary.

Developmental stages for the divergence of relative limb length between a twig and a trunk-ground Anolis lizard species

The divergent evolution of niche-related traits can facilitate adaptive radiation, yet identification of the genetic or molecular mechanisms underlying such trait changes remains a major challenge in evolutionary biology. Conducting a detailed morphological comparison along growth trajectories is a powerful method for observing the formation of differences in niche-related traits. Here, we focused on hindlimb length of Anolis lizards, differences in which are related to adaptation for use of different microhabitats. We measured the length of hindlimb skeletons in different ecomorphs of anole lizards (A. sagrei, a trunk-ground ecomorph with long hindlimbs, and A. angusticeps, a twig ecomorph with short hindlimbs) from early embryonic stages to adulthood, to determine which hindlimb elements mainly differentiate the species and the timing of the formation of these differences. With respect to the digit, differences between the species mainly occurred during the embryonic stages of interdigit reduction, when the cartilage of the distal phalanges was simultaneously forming. In addition, we compared the relative length of developing autopods in early embryonic stages using whole-mount in situ hybridization before the formation of the cartilaginous bones, and the results showed that the relative growth rate of the Hoxa11-negative distal region in A. sagrei was greater than that in A. angusticeps. Our results show that there are several important developmental stages for hindlimb length differentiation between A. angusticeps and A. sagrei, depending on which hindlimb element is considered. In particular, the species differences were largely due to variations in digit length, which arose at early embryonic stages.

Structure, function and evolution of topologically associating domains (TADs) at HOX loci

Hox genes encode transcription factors necessary for patterning the major developing anterior to posterior embryonic axis. In addition, during vertebrate evolution, various subsets of this gene family were co-opted along with the emergence of novel body structures, such as the limbs or the external genitalia. The morphogenesis of these axial structures thus relies in part upon the precisely controlled transcription of specific Hox genes, a mechanism involving multiple long-range enhancers. Recently, it was reported that such regulatory mechanisms were largely shared between different developing tissues, though with some specificities, suggesting the recruitment of ancestral regulatory modalities from one tissue to another. The analysis of chromatin architectures at HoxD and HoxA loci revealed the existence of two flanking topologically associating domains (TADs), precisely encompassing the adjacent regulatory landscapes. Here, we discuss the function of these TADs in the control of Hox gene regulation and we speculate about their capacity to serve as structural frameworks for the emergence of novel enhancers. In this view, TADs may have been used as genomic niches to evolve pleiotropic regulations found at many developmental loci.

Possible association of single nucleotide polymorphisms in the 3' untranslated region of HOXB9 with acetabular overcoverage

OBJECTIVES

Excessive acetabular coverage is the most common cause of pincer-type femoroacetabular impingement. To date, an association between acetabular over-coverage and genetic variations has not been studied. In this study we investigated the association between single nucleotide polymorphisms (SNPs) of paralogous Homeobox (HOX)9 genes and acetabular coverage in Japanese individuals to identify a possible genetic variation associated with acetabular over-coverage.

METHODS

We investigated 19 total SNPs in the four HOX9 paralogs, then focused in detail on seven of those located in the 3' untranslated region of HOXB9 (rs8844, rs3826541, rs3826540, rs7405887, rs2303485, rs2303486, rs79931349) using a case-control association study. The seven HOXB9 SNPs were genotyped in 316 subjects who had all undergone radiological examination. The association study was performed by both single-locus and haplotype-based analyses.

RESULTS

The genotype and allele frequencies of the five HOXB9 SNPs showed significant association with acetabular over-coverage compared with controls (rs7405887 OR = 3.16, p = 5.29E-6, 95% CI 1.91 to 5.25). A significant difference was also detected when haplotypes were evaluated (OR = 2.59, p = 2.61E-5, 95% CI 1.65 to 4.08). The two HOXB9 SNPs (rs2303485, rs2303486) were associated with decreased acetabular coverage (rs2303485 OR = 0.524, p = 0.0091, 95% CI 0.322 to 0.855; rs2303486 OR = 0.519, p = 0.011, 95% CI 0.312 to 0.865).

CONCLUSIONS

The five HOXB9 SNPs (rs8844, rs3826541, rs3826540, rs7405887, rs79931349) were associated with acetabular over-coverage. On the other hand, the two SNPs (rs2303485 and rs2303486) were associated with the lower acetabular coverage. The association of rs2303486 would be consistent with the previous study. Therefore, the HOXB9 SNPs might be involved in the morphogenesis of acetabular coverage, and could be an independent risk factor for developing pincer-type femoroacetabular impingement. Cite this article: Bone Joint Res 2015;4:50-5.

HoxA and HoxD expression in a variety of vertebrate body plan features reveals an ancient origin for the distal Hox program

Background

Hox genes are master regulatory genes that specify positional identities during axial development in animals. Discoveries regarding their concerted expression patterns have commanded intense interest due to their complex regulation and specification of body plan features in jawed vertebrates. For example, the posterior HoxD genes switch to an inverted collinear expression pattern in the mouse autopod where HoxD13 switches from a more restricted to a less restricted domain relative to its neighboring gene on the cluster. We refer to this program as the ‘distal phase’ (DP) expression pattern because it occurs in distal regions of paired fins and limbs, and is regulated independently by elements in the 5′ region upstream of the HoxD cluster. However, few taxa have been evaluated with respect to this pattern, and most studies have focused on pectoral fin morphogenesis, which occurs relatively early in development.

Results

Here, we demonstrate for the first time that the DP expression pattern occurs with the posterior HoxA genes, and is therefore not solely associated with the HoxD gene cluster. Further, DP Hox expression is not confined to paired fins and limbs, but occurs in a variety of body plan features, including paddlefish barbels - sensory adornments that develop from the first mandibular arch (the former ‘Hox-free zone), and the vent (a medial structure that is analogous to a urethra). We found DP expression of HoxD13 and HoxD12 in the paddlefish barbel; and we present the first evidence for DP expression of the HoxA genes in the hindgut and vent of three ray-finned fishes. The HoxA DP expression pattern is predicted by the recent finding of a shared 5′ regulatory architecture in both the HoxA and HoxD clusters, but has not been previously observed in any body plan feature.

Conclusions

The Hox DP expression pattern appears to be an ancient module that has been co-opted in a variety of structures adorning the vertebrate bauplan. This module provides a shared genetic program that implies deep homology of a variety of distally elongated structures that has played a significant role in the evolution of morphological diversity in vertebrates

Electronic supplementary material

The online version of this article (doi:10.1186/2041-9139-5-44) contains supplementary material, which is available to authorized users.

The pisiform growth plate is lost in humans and supports a role for Hox in growth plate formation

The human pisiform is a small, nodular, although functionally significant, bone of the wrist. In most other mammals, including apes and Australopithecus afarensis, pisiforms are elongate. An underappreciated fact is that the typical mammalian pisiform forms from two ossification centers. We hypothesize that: (i) the presence of a secondary ossification center in mammalian pisiforms indicates the existence of a growth plate; and (ii) human pisiform reduction results from growth plate loss. To address these hypotheses, we surveyed African ape pisiform ossification and confirmed the presence of a late-forming secondary ossification center in chimpanzees and gorillas. Identification of the initial ossification center occurs substantially earlier in apes relative to humans, raising questions concerning the homology of the human pisiform and the two mammalian ossification centers. Second, we conducted histological and immunohistochemical analyses of pisiform ossification in mice. We confirm the presence of two ossification centers separated by organized columnar and hypertrophic chondrocyte zones. Flattened chondrocytes were highly mitotic, indicating the presence of a growth plate. Hox genes have been proposed to play a fundamental role in growth plate patterning. The existence of a pisiform growth plate presents an interesting test case for the association between Hox expression and growth plate formation, and could explain the severe effects on the pisiform observed in Hoxa11 and Hoxd11 knockout mice. Consistent with this hypothesis, we show that Hoxd11 is expressed adjacent to the pisiform in late-stage embryonic mouse limbs supporting a role for Hox genes in growth plate specification. This raises questions concerning the mechanisms regulating Hox expression in the developing carpus.

Genetic interactions between Shox2 and Hox genes during the regional growth and development of the mouse limb

The growth and development of the vertebrate limb relies on homeobox genes of the Hox and Shox families, with their independent mutation often giving dose-dependent effects. Here we investigate whether Shox2 and Hox genes function together during mouse limb development by modulating their relative dosage and examining the limb for nonadditive effects on growth. Using double mRNA fluorescence in situ hybridization (FISH) in single embryos, we first show that Shox2 and Hox genes have associated spatial expression dynamics, with Shox2 expression restricted to the proximal limb along with Hoxd9 and Hoxa11 expression, juxtaposing the distal expression of Hoxa13 and Hoxd13. By generating mice with all possible dosage combinations of mutant Shox2 alleles and HoxA/D cluster deletions, we then show that their coordinated proximal limb expression is critical to generate normally proportioned limb segments. These epistatic interactions tune limb length, where Shox2 underexpression enhances, and Shox2 overexpression suppresses, Hox-mutant phenotypes. Disruption of either Shox2 or Hox genes leads to a similar reduction in Runx2 expression in the developing humerus, suggesting their concerted action drives cartilage maturation during normal development. While we furthermore provide evidence that Hox gene function influences Shox2 expression, this regulation is limited in extent and is unlikely on its own to be a major explanation for their genetic interaction. Given the similar effect of human SHOX mutations on regional limb growth, Shox and Hox genes may generally function as genetic interaction partners during the growth and development of the proximal vertebrate limb.

Hox gene regulation and timing in embryogenesis

Hox genes are critical regulators of embryonic development in bilaterian animals. They exhibit a unique mode of transcriptional regulation where the position of the genes along the chromosome corresponds to the time and place of their expression during development. The sequential temporal activation of these genes in the primitive streak helps determining their subsequent pattern of expression along the anterior-posterior axis of the embryo, yet the precise correspondence between these two collinear processes is not fully understood. In addition, vertebrate Hox genes evolved similar modes of regulation along secondary body axes, such as the developing limbs. We review the current understanding of the mechanisms operating during activation, maintenance and silencing of Hox gene expression in these various contexts, and discuss the evolutionary significance of their genomic organization.

Direct functional consequences of ZRS enhancer mutation combine with secondary long range SHH signalling effects to cause preaxial polydactyly

Sonic hedgehog (SHH) plays a central role in patterning numerous embryonic tissues including, classically, the developing limb bud where it controls digit number and identity. This study utilises the polydactylous Silkie (Slk) chicken breed, which carries a mutation in the long range limb-specific regulatory element of SHH, the ZRS. Using allele specific SHH expression analysis combined with quantitative protein analysis, we measure allele specific changes in SHH mRNA and concentration of SHH protein over time. This confirms that the Slk ZRS enhancer mutation causes increased SHH expression in the posterior leg mesenchyme. Secondary consequences of this increased SHH signalling include increased FGF pathway signalling and growth as predicted by the SHH/GREM1/FGF feedback loop and the Growth/Morphogen models. Manipulation of Hedgehog, FGF signalling and growth demonstrate that anterior-ectopic expression of SHH and induction of preaxial polydactyly is induced secondary to increased SHH signalling and Hedgehog-dependent growth directed from the posterior limb. We predict that increased long range SHH signalling acts in combination with changes in activation of SHH transcription from the Slk ZRS allele. Through analysis of the temporal dynamics of anterior SHH induction we predict a gene regulatory network which may contribute to activation of anterior SHH expression from the Slk ZRS.

•

Overexpression of posterior SHH in the limb bud can cause preaxial polydactyly.

•

Increased activation of SHH/GREM/FGF feedback and growth induces Slk preaxial polydactyly.

•

Autoregulated expression of SHH can occur within 1.5–2 h in the limb bud.

Diffusible signals and epigenetic timing cooperate in late proximo-distal limb patterning

Developing vertebrate limbs initiate proximo-distal patterning by interpreting opposing gradients of diffusible signaling molecules. We report two thresholds of proximo-distal signals in the limb bud: a higher threshold that establishes the upper-arm to forearm transition; and a lower one that positions a later transition from forearm to hand. For this last transition to happen, however, the signal environment seems to be insufficient, and we show that a timing mechanism dependent on histone acetylation status is also necessary. Therefore, as a consequence of the time dependence, the lower signaling threshold remains cryptic until the timing mechanism reveals it. We propose that this timing mechanism prevents the distal transition from happening too early, so that the prospective forearm has enough time to expand and form a properly sized segment. Importantly, the gene expression changes provoked by the first transition further regulate proximo-distal signal distribution, thereby coordinating the positioning of the two thresholds, which ensures robustness. This model is compatible with the most recent genetic analyses and underscores the importance of growth during the time-dependent patterning phase, providing a new mechanistic framework for understanding congenital limb defects.

A combination of activation and repression by a colinear Hox code controls forelimb-restricted expression of Tbx5 and reveals Hox protein specificity

Tight control over gene expression is essential for precision in embryonic development and acquisition of the regulatory elements responsible is the predominant driver for evolution of new structures. Tbx5 and Tbx4, two genes expressed in forelimb and hindlimb-forming regions respectively, play crucial roles in the initiation of limb outgrowth. Evolution of regulatory elements that activate Tbx5 in rostral LPM was essential for the acquisition of forelimbs in vertebrates. We identified such a regulatory element for Tbx5 and demonstrated Hox genes are essential, direct regulators. While the importance of Hox genes in regulating embryonic development is clear, Hox targets and the ways in which each protein executes its specific function are not known. We reveal how nested Hox expression along the rostro-caudal axis restricts Tbx5 expression to forelimb. We demonstrate that Hoxc9, which is expressed in caudal LPM where Tbx5 is not expressed, can form a repressive complex on the Tbx5 forelimb regulatory element. This repressive capacity is limited to Hox proteins expressed in caudal LPM and carried out by two separate protein domains in Hoxc9. Forelimb-restricted expression of Tbx5 and ultimately forelimb formation is therefore achieved through co-option of two characteristics of Hox genes; their colinear expression along the body axis and the functional specificity of different paralogs. Active complexes can be formed by Hox PG proteins present throughout the rostral-caudal LPM while restriction of Tbx5 expression is achieved by superimposing a dominant repressive (Hoxc9) complex that determines the caudal boundary of Tbx5 expression. Our results reveal the regulatory mechanism that ensures emergence of the forelimbs at the correct position along the body. Acquisition of this regulatory element would have been critical for the evolution of limbs in vertebrates and modulation of the factors we have identified can be molecular drivers of the diversity in limb morphology.

The acquisition of limbs during vertebrate evolution was a very successful innovation that enabled this group of species to diversify and colonise land. It has become clear recently that the primary driver behind the evolution of new structures, such as limbs, is the acquisition of novel regulatory elements that control when and where genes are activated rather than the proteins encoded by the genes themselves acquiring novel functions. We have identified the regulatory element from a gene, Tbx5. Activation of Tbx5 in the forelimb-forming region of the developing embryos is essential for forelimbs to form and disruption of human TBX5 causes limb abnormalities. We show that activation of Tbx5 in a restricted territory is achieved through a combination of activation inputs that are present broadly throughout the embryo flank and dominant, repressive inputs present only in more caudal regions of the flank. The sum of these inputs yields restricted activation in the rostral, forelimb-forming flank. Our results explain how the regulatory switches that were harnessed for the acquisition of limbs during evolution operate and how they can be turned off during the evolution of limblessness in species such as the snake.

Association of a single nucleotide polymorphism in HOXB9 with developmental dysplasia of the hip: a case-control study

Developmental dysplasia of the hip (DDH) is one of the most common skeletal disorders. It comprises a spectrum of abnormalities, including shallow acetabulum and decreased coverage of the femoral head. Genetic component plays a considerable role in the aetiology of DDH. HOXB9 may be involved in the aetiology and pathogenesis of DDH, as it plays an important role in the development of the limbs. Our objective was to evaluate whether single nucleotide polymorphisms (SNPs) of HOXB9 (rs2303486 and rs8844) were associated with DDH in Chinese population. The HOXB9 tag SNPs were genotyped in 460 DDH cases and 562 control subjects by Taqman assay, and their association was examined. rs8844 was not associated with DDH. rs2303486 was associated with DDH in the dominant genetic model (p = 0.037; odds ratio (OR) = 1.32; 95% confidence interval (CI) = 1.02–1.71). After stratification by sex, significant association of the dominant genetic model still existed in the female subjects (p = 0.015; OR = 1.46; 95% CI = 1.08–1.98),but not in the male subjects. After stratification by severity, we discovered an association with hip dislocation in the dominant model(p=0.042; OR = 1.35; 95% CI = 1.01–1.80), but not with subluxation or instability. HOXB9 is associated with DDH in Chinese.

Conservation and divergence of regulatory strategies at Hox Loci and the origin of tetrapod digits

During development, expression of the Hoxa and Hoxd genes in zebrafish fins and mouse limbs are regulated via a conserved chromatin structure. However, zebrafish lack certain regulatory elements required to produce digits, revealing that radials—the fin's bony elements—are likely not homologous to tetrapod digits.

The evolution of tetrapod limbs from fish fins enabled the conquest of land by vertebrates and thus represents a key step in evolution. Despite the use of comparative gene expression analyses, critical aspects of this transformation remain controversial, in particular the origin of digits. Hoxa and Hoxd genes are essential for the specification of the different limb segments and their functional abrogation leads to large truncations of the appendages. Here we show that the selective transcription of mouse Hoxa genes in proximal and distal limbs is related to a bimodal higher order chromatin structure, similar to that reported for Hoxd genes, thus revealing a generic regulatory strategy implemented by both gene clusters during limb development. We found the same bimodal chromatin architecture in fish embryos, indicating that the regulatory mechanism used to pattern tetrapod limbs may predate the divergence between fish and tetrapods. However, when assessed in mice, both fish regulatory landscapes triggered transcription in proximal rather than distal limb territories, supporting an evolutionary scenario whereby digits arose as tetrapod novelties through genetic retrofitting of preexisting regulatory landscapes. We discuss the possibility to consider regulatory circuitries, rather than expression patterns, as essential parameters to define evolutionary synapomorphies.

Our upper limbs differ from fish fins, notably by their subdivision into arm and hand regions, which are separated by a complex articulation, the wrist. The development of this anatomy is associated with two distinct waves of expression of the Hoxa and Hoxd genes during development. Would such a shared expression pattern be sufficient to infer homology between fish fins and mouse limbs? We investigated this question here, looking at whether the two phases of Hox gene transcription that are observed during tetrapod limb development also occur during zebrafish fin development. We find the answer is “not quite.” For although the mechanisms that regulate the expression of Hoxa and Hoxd are comparable between zebrafish fins and mouse limbs, when the genomic regions that regulate Hox gene expression in fish fins are introduced into transgenic mice, they trigger Hox gene expression in only the proximal limb segment (the segment nearest the body) and not in the presumptive digits. We conclude that although fish have the Hox regulatory toolkit to produce digits, this potential is not utilized as it is in tetrapods, and as a result we propose that fin radials—the bony elements of fins—are not homologous to tetrapod digits.

Progressive specification rather than intercalation of segments during limb regeneration

An amputated salamander limb regenerates the correct number of segments. Models explaining limb regeneration were largely distinct from those for limb development, despite the presence of common patterning molecules. Intercalation has been an important concept to explain salamander limb regeneration, but clear evidence supporting or refuting this model was lacking. In the intercalation model, the first blastema cells acquire fingertip identity, creating a gap in positional identity that triggers regeneration of the intervening region from the stump. We used HOXA protein analysis and transplantation assays to show that axolotl limb blastema cells acquire positional identity in a proximal-to-distal sequence. Therefore, intercalation is not the primary mechanism for segment formation during limb regeneration in this animal. Patterning in development and regeneration uses similar mechanisms.

"Self-regulation," a new facet of Hox genes' function

BACKGROUND

Precise temporal and spatial expression of the clustered Hox genes is essential for patterning the developing embryo. Temporal activation of Hox genes was shown to be cluster-autonomous. However, gene clustering appears dispensable for spatial colinear expression. Generally, a set of Hox genes expressed in a group of cells instructs these cells about their fate such that the differential expression of Hox genes results in morphological diversity. The spatial colinearity is considered to rely both on local and long-range cis regulation.

RESULTS

Here, we report on the global deregulation of HoxA and HoxD expression patterns upon inactivation of a subset of HOXA and HOXD proteins.

CONCLUSIONS

Our data suggest the existence of a "self-regulation" mechanism, a process by which HOX proteins establish and/or maintain the spatial domains of the Hox gene family and we propose that the functionally dominant HOX proteins could contribute to generating the spatial parameters of Hox expression in a given tissue, i.e., HOX controlling the establishment of the ultimate HOX code.

Hox genes regulation in vertebrates

Hox genes encode transcription factors defining cellular identities along the major and secondary body axes. Their coordinated expression in both space and time is critical for embryonic patterning. Accordingly, Hox genes transcription is tightly controlled at multiple levels, and involves an intricate combination of local and long-range cis-regulatory elements. Recent studies revealed that in addition to transcription factors, dynamic patterns of histone marks and higher-order chromatin structure are important determinants of Hox gene regulation. Furthermore, the emerging picture suggests an involvement of various species of non-coding RNA in targeting activating and repressive complexes to Hox clusters. I review these recent developments and discuss their relevance to the control of Hox gene expression in vivo, as well as to our understanding of transcriptional regulatory mechanisms.

Transcriptional components of anteroposterior positional information during zebrafish fin regeneration

Many fish and salamander species regenerate amputated fins or limbs, restoring the size and shape of the original appendage. Regeneration requires that spared cells retain or recall information encoding pattern, a phenomenon termed positional memory. Few factors have been implicated in positional memory during vertebrate appendage regeneration. Here, we investigated potential regulators of anteroposterior (AP) pattern during fin regeneration in adult zebrafish. Sequence-based profiling from tissues along the AP axis of uninjured pectoral fins identified many genes with region-specific expression, several of which encoded transcription factors with known AP-specific expression or function in developing embryonic pectoral appendages. Transgenic reporter strains revealed that regulatory sequences of the transcription factor gene alx4a activated expression in fibroblasts and osteoblasts within anterior fin rays, whereas hand2 regulatory sequences activated expression in these same cell types within posterior rays. Transgenic overexpression of hand2 in all pectoral fin rays did not affect formation of the proliferative regeneration blastema, yet modified the lengths and widths of regenerating bones. Hand2 influenced the character of regenerated rays in part by elevation of the vitamin D-inactivating enzyme encoded by cyp24a1, contributing to region-specific regulation of bone metabolism. Systemic administration of vitamin D during regeneration partially rescued bone defects resulting from hand2 overexpression. Thus, bone-forming cells in a regenerating appendage maintain expression throughout life of transcription factor genes that can influence AP pattern, and differ across the AP axis in their expression signatures of these and other genes. These findings have implications for mechanisms of positional memory in vertebrate tissues.

A switch between topological domains underlies HoxD genes collinearity in mouse limbs

Hox genes are major determinants of the animal body plan, where they organize structures along both the trunk and appendicular axes. During mouse limb development, Hoxd genes are transcribed in two waves: early on, when the arm and forearm are specified, and later, when digits form. The transition between early and late regulations involves a functional switch between two opposite topological domains. This switch is reflected by a subset of Hoxd genes mapping centrally into the cluster, which initially interact with the telomeric domain and subsequently swing toward the centromeric domain, where they establish new contacts. This transition between independent regulatory landscapes illustrates both the modularity of the limbs and the distinct evolutionary histories of its various pieces. It also allows the formation of an intermediate area of low HOX proteins content, which develops into the wrist, the transition between our arms and our hands. This regulatory strategy accounts for collinear Hox gene regulation in land vertebrate appendages.

Misexpression of Pknox2 in mouse limb bud mesenchyme perturbs zeugopod development and deltoid crest formation

The TALE (Three Amino acid Loop Extension) family consisting of Meis, Pbx and Pknox proteins is a group of transcriptional co-factors with atypical homeodomains that play pivotal roles in limb development. Compared to the in-depth investigations of Meis and Pbx protein functions, the role of Pknox2 in limb development remains unclear. Here, we showed that Pknox2 was mainly expressed in the zeugopod domain of the murine limb at E10.5 and E11.5. Misexpression of Pknox2 in the limb bud mesenchyme of transgenic mice led to deformities in the zeugopod and forelimb stylopod deltoid crest, but left the autopod and other stylopod skeletons largely intact. These malformations in zeugopod skeletons were recapitulated in mice overexpressing Pknox2 in osteochondroprogenitor cells. Molecular and cellular analyses indicated that the misexpression of Pknox2 in limb bud mesenchyme perturbed the Hox10-11 gene expression profiles, decreased Col2 expression and Bmp/Smad signaling activity in the limb. These results indicated that Pknox2 misexpression affected mesenchymal condensation and early chondrogenic differentiation in the zeugopod skeletons of transgenic embryos, suggesting Pknox2 as a potential regulator of zeugopod and deltoid crest formation.

Chromatin organization and global regulation of Hox gene clusters

During development, a properly coordinated expression of Hox genes, within their different genomic clusters is critical for patterning the body plans of many animals with a bilateral symmetry. The fascinating correspondence between the topological organization of Hox clusters and their transcriptional activation in space and time has served as a paradigm for understanding the relationships between genome structure and function. Here, we review some recent observations, which revealed highly dynamic changes in the structure of chromatin at Hox clusters, in parallel with their activation during embryonic development. We discuss the relevance of these findings for our understanding of large-scale gene regulation.

Metapodial or phalanx? An evolutionary and developmental perspective on the homology of the first ray's proximal segment

The first mammalian metapodial (MP1) has periodically been argued to actually be a phalanx, because the first ray has one less element than the four posterior rays, and because the MP1 growth plate is proximal like those of all phalanges, rather than distal as in metapodials 2-5. However, growth plates are formed at both ends in non-therian tetrapod metapodials, and phylogenetic analysis demonstrates that growth plate loss is a therian synapomorphy that postdates the establishment of the mammalian phalangeal formula. These data, along with results of developmental and morphological studies, suggest that the MP1 is not a phalanx. The singular, proximal growth plates in MPs 2-5 are likely to be an adaptation to dynamic erect quadrupedal gait which was characterized by conversion of the posterior metapodials into rigid struts with the carpus/tarsus. While the adaptive significance of the reversed ossification of MP1 is less clear, we present three functional/developmental hypotheses.

Wnt and BMP signaling cooperate with Hox in the control of Six2 expression in limb tendon precursor

The number and shape of limb tendons vary along the proximodistal axis, and the autopod contains more tendons than the zeugopod. The transcription factor Six2 is expressed in the developing tendons, and its expression can be traced back to a group of limb mesenchymal cells that are thought to be tendon precursor cells. We tried to elucidate the mechanism controlling position-specific tendon pattern formation using Six2 as a tendon marker. Six2 expression was always found in cells between the limb cartilage and ectoderm. Administration of BMP-2 or BMP antagonist Noggin to the limb bud, respectively repressed or facilitated Six2 expression. Removal of the ectoderm or administration of the Wnt antagonist sFRP-2 abolished Six2 expression and ectopic Wnt expression induced ectopic Six2 expression. Taken together, Six2 expression is induced in the cells located at the point where cartilage-derived Noggin and ectoderm-derived Wnt signals meet. Misexpression of the autopod-specific Hoxa-13 or Hoxd-13 induced ectopic expression of Six2 in the zeugopodal mesenchymal cells of the chick limb bud. Six2 expression in the dorsal autopodal mesenchyme was not detected in Hoxa-13(-/-);HoxD(del/del) mice, indicating that autopod-specific Hox is required for the regulation of Six2 expression. Misexpression of Wnt in the autopod induced ectopic Six2 expression in the autopod. On the other hand, Wnt misexpression alone never induced Six2 expression in the zeugopod, yet co-misexpression of Hoxa-13 and Wnt in the zeugopod enhanced ectopic Six2 expression. Our results indicate that autopodal Hox genes regulate Six2 expression in the autopodal tendon precursor in cooperation with the factors from cartilage and ectoderm.

The origin of the tetrapod limb: from expeditions to enhancers

More than three centuries ago natural philosophers, and later anatomists, recognized a fundamental organization to the skeleton of tetrapod limbs. Composed of three segments, stylopod, zeugopod, and autopod, this pattern has served as the basis for a remarkably broad adaptive radiation from wings and flippers to hands and digging organs. A central area of inquiry has been tracing the origins of the elements of this Bauplan in the fins of diverse fish. Can equivalents of the three segments, and the developmental processes that pattern them, be seen in fish fins? In addition, if so, how do these data inform theories of the transformation of fins into limbs? Answers to these questions come from linking discoveries in paleontology with those of developmental biology and genetics. Burgeoning discoveries in the regulatory biology of developmental genes and in the genomics of diverse species offer novel data to investigate these classical questions.

The evolutionary history of the development of the pelvic fin/hindlimb

The arms and legs of man are evolutionarily derived from the paired fins of primitive jawed fish. Few evolutionary changes have attracted as much attention as the origin of tetrapod limbs from the paired fins of ancestral fish. The hindlimbs of tetrapods are derived from the pelvic fins of ancestral fish. These evolutionary origins can be seen in the examination of shared gene and protein expression patterns during the development of pelvic fins and tetrapod hindlimbs. The pelvic fins of fish express key limb positioning, limb bud induction and limb outgrowth genes in a similar manner to that seen in hindlimb development of higher vertebrates. We are now at a point where many of the key players in the development of pelvic fins and vertebrate hindlimbs have been identified and we can now readily examine and compare mechanisms between species. This is yielding fascinating insights into how the developmental programme has altered during evolution and how that relates to anatomical change. The role of pelvic fins has also drastically changed over evolutionary history, from playing a minor role during swimming to developing into robust weight-bearing limbs. In addition, the pelvic fins/hindlimbs have been lost repeatedly in diverse species over evolutionary time. Here we review the evolution of pelvic fins and hindlimbs within the context of the changes in anatomical structure and the molecular mechanisms involved.

Molecular and evolutionary basis of limb field specification and limb initiation

Specification of limb field and initiation of limb development involve multiple steps, each of which is tightly regulated both spatially and temporally. Recent developmental analyses on various vertebrates have provided insights into the molecular mechanisms that specify limb field and have revealed several genetic interactions of signals involved in limb initiation processes. Furthermore, new approaches to the study of the developmental mechanisms of the lateral plate mesoderm of amphioxus and lamprey embryos have given us clues to understand the evolutionary scenarios that led to the acquisition of paired appendages during evolution. This review highlights such recent findings and discusses the mechanisms of limb field specification and limb bud initiation during development and evolution.

Mathematical modeling of vertebrate limb development

In this paper, we review the major mathematical and computational models of vertebrate limb development and their roles in accounting for different aspects of this process. The main aspects of limb development that have been modeled include outgrowth and shaping of the limb bud, establishment of molecular gradients within the bud, and formation of the skeleton. These processes occur interdependently during development, although (as described in this review), there are various interpretations of the biological relationships among them. A wide range of mathematical and computational methods have been used to study these processes, including ordinary and partial differential equation systems, cellular automata and discrete, stochastic models, finite difference methods, finite element methods, the immersed boundary method, and various combinations of the above. Multiscale mathematical modeling and associated computational simulation have become integrated into the study of limb morphogenesis and pattern formation to an extent with few parallels in the field of developmental biology. These methods have contributed to the design and analysis of experiments employing microsurgical and genetic manipulations, evaluation of hypotheses for limb bud outgrowth, interpretation of the effects of natural mutations, and the formulation of scenarios for the origination and evolution of the limb skeleton.

A misplaced lncRNA causes brachydactyly in humans

Translocations are chromosomal rearrangements that are frequently associated with a variety of disease states and developmental disorders. We identified 2 families with brachydactyly type E (BDE) resulting from different translocations affecting chromosome 12p. Both translocations caused downregulation of the parathyroid hormone-like hormone (PTHLH) gene by disrupting the cis-regulatory landscape. Using chromosome conformation capturing, we identified a regulator on chromosome 12q that interacts in cis with PTHLH over a 24.4-megabase distance and in trans with the sex-determining region Y-box 9 (SOX9) gene on chromosome 17q. The element also harbored a long noncoding RNA (lncRNA). Silencing of the lncRNA, PTHLH, or SOX9 revealed a feedback mechanism involving an expression-dependent network in humans. In the BDE patients, the human lncRNA was upregulated by the disrupted chromosomal association. Moreover, the lncRNA occupancy at the PTHLH locus was reduced. Our results document what we believe to be a novel in cis- and in trans-acting DNA and lncRNA regulatory feedback element that is reciprocally regulated by coding genes. Furthermore, our findings provide a systematic and combinatorial view of how enhancers encoding lncRNAs may affect gene expression in normal development.

Molecular profiling of synovial joints: use of microarray analysis to identify factors that direct the development of the knee and elbow

BACKGROUND

Synovial joints develop from the interzone, a dense layer of mesenchymal progenitor cells that marks the site of the future joint. During the morphogenic events that follow, joints attain their distinct shape and organization. The molecular mechanisms controlling the initial specification of synovial joints has been studied, but the question of how individual joints attain the specific structure required for their unique functions remains largely unresolved. Here, we use microarray analysis to compare knee and elbow formation to identify factors involved in the development of specific joints.

RESULTS

The knee is enriched for the hindlimb patterning genes Hoxc9, Hoxc10, and Tbx4 and for Tgfbi, Rspo2, and Sfrp2, factors involved in transforming growth factor-beta/bone morphogenetic protein (TGFβ/BMP) and Wnt signaling. Consistent with these findings, we show that TGFβ signaling directs knee morphogenesis, and is necessary for meniscus development. The tissue surrounding the elbow is highly enriched for genes involved in muscle specification and differentiation, and in splotch-delayed muscleless mutants, elbow, but not knee morphogenesis is disrupted.

CONCLUSIONS

Our results suggest there are fundamental differences in how individual joints develop after interzone formation. Our microarray analyses provides a new resource for further investigation of the pathways involved in the morphogenesis of specific synovial joints.

Copy number analysis of 413 isolated talipes equinovarus patients suggests role for transcriptional regulators of early limb development

Talipes equinovarus is one of the most common congenital musculoskeletal anomalies and has a worldwide incidence of 1 in 1000 births. A genetic predisposition to talipes equinovarus is evidenced by the high concordance rate in twin studies and the increased risk to first-degree relatives. Despite the frequency of isolated talipes equinovarus and the strong evidence of a genetic basis for the disorder, few causative genes have been identified. To identify rare and/or recurrent copy number variants, we performed a genome-wide screen for deletions and duplications in 413 isolated talipes equinovarus patients using the Affymetrix 6.0 array. Segregation analysis within families and gene expression in mouse E12.5 limb buds were used to determine the significance of copy number variants. We identified 74 rare, gene-containing copy number variants that were present in talipes equinovarus probands and not present in 759 controls or in the Database of Genomic Variants. The overall frequency of copy number variants was similar between talipes equinovarus patients compared with controls. Twelve rare copy number variants segregate with talipes equinovarus in multiplex pedigrees, and contain the developmentally expressed transcription factors and transcriptional regulators PITX1, TBX4, HOXC13, UTX, CHD (chromodomain protein)1, and RIPPLY2. Although our results do not support a major role for recurrent copy number variations in the etiology of isolated talipes equinovarus, they do suggest a role for genes involved in early embryonic patterning in some families that can now be tested with large-scale sequencing methods.

Heterochronically early decline of Hox expression prior to cartilage formation in the avian hindlimb zeugopod

The fibula, a zeugopod bone in the hindlimb, exhibits various morphologies in tetrapod species. The fibula in some species has a similar length with the other zeugopod element, the tibia, while other species have obvious differences in the sizes of the two elements. In the avian hindlimb, for example, the fibula is extremely short, thin, and truncated. Basic morphology of the fibula is established during development, and cartilage primordium of the bone emerges in a certain region defined by a distinct combination of expression of Hox genes (Hox code). In order to elucidate how the different morphologies are produced from a region that is defined as the fixed Hox code, we examined spatial and temporal patterns of Hoxd11/Hoxd12 expression in the developing limb bud, which defines the region from which the fibula emerges, in comparison with the sites of precartilaginous mesenchymal condensations representing regions for cartilage formation among chick, mouse, and gecko embryos. We found that in the chick hindlimb, expression of Hoxd11/Hoxd12 decreased and disappeared from the presumptive zeugopod region before cartilage formation. This heterochronically early decline of expression of Hox genes is strongly correlated with the peculiar trait of the fibula in the avian hindlimb, since in the other species examined, expression of those genes continued after the onset of cartilage formation. This is morphological phenotype-related because the early disappearance was not seen in the chick forelimb. Our results suggest that temporal change of the Hox code governs diversification in morphology of homologous structures among related species.

A novel 13 base pair insertion in the sonic hedgehog ZRS limb enhancer (ZRS/LMBR1) causes preaxial polydactyly with triphalangeal thumb

Mutations in the Sonic hedgehog limb enhancer, the zone of polarizing activity regulatory sequence (ZRS, located within the gene LMBR1), commonly called the ZRS), cause limb malformations. In humans, three classes of mutations have been proposed based on the limb phenotype; single base changes throughout the region cause preaxial polydactyly (PPD), single base changes at one specific site cause Werner mesomelic syndrome, and large duplications cause polysyndactyly. This study presents a novel mutation-a small insertion. In a Swedish family with autosomal-dominant PPD, we found a 13 base pair insertion within the ZRS, NG_009240.1:g.106934_106935insTAAGGAAGTGATT (traditional nomenclature: ZRS603ins13). Computational transcription factor-binding site predictions suggest that this insertion creates new binding sites and a mouse enhancer assay shows that this insertion causes ectopic gene expression. This study is the first to discover a small insertion in an enhancer that causes a human limb malformation and suggests a potential mechanism that could explain the ectopic expression caused by this mutation.

Association analysis between HOXD9 genes and the development of developmental dysplasia of the hip in Chinese female Han population

Background

Developmental dysplasia of the hip (DDH) is a congenital or acquired deformation or misalignment of the hip joint which affects mainly females. We hypothesized that HOXD9 gene could be regulated in acetabular size or shape and related in DDH developing.

Methods

Two hundred and nine Chinese Han female DDH patients and 173 ethnic, age matched healthy female controls were genotyped for HOXD9 two tag SNPs using sequenom method.

Results

One of the two tag SNPs, rs711822, was not shown significantly differences in genotypic or allelic distribution between case and control group. Comparing the genotypic distribution of rs711819, there was significant differences between DDH patients group and control group (χ2 = 7.54, df =2, P =0.023), and the association to DDH developing reached significance (P =0.045, OR =1.79, 95 % CI: 1.01-3.17 by dominant mode).

Conclusion

In conclusion, the association between one tag SNP of HOXD9 gene and the development of DDH reach significant in our study population, this result indicate the positive correlation between HOXD9 gene and DDH developing. Further study in larger sample size and different population as well as functional studies will help to understand the pathogenesis of DDH.

The making of differences between fins and limbs

'Evo-devo', an interdisciplinary field based on developmental biology, includes studies on the evolutionary processes leading to organ morphologies and functions. One fascinating theme in evo-devo is how fish fins evolved into tetrapod limbs. Studies by many scientists, including geneticists, mathematical biologists, and paleontologists, have led to the idea that fins and limbs are homologous organs; now it is the job of developmental biologists to integrate these data into a reliable scenario for the mechanism of fin-to-limb evolution. Here, we describe the fin-to-limb transition based on key recent developmental studies from various research fields that describe mechanisms that may underlie the development of fins, limb-like fins, and limbs.

Three Drosophila Hox complex microRNAs do not have major effects on expression of evolutionarily conserved Hox gene targets during embryogenesis

The discovery of microRNAs has resulted in a major expansion of the number of molecules known to be involved in gene regulation. Elucidating the functions of animal microRNAs has posed a significant challenge as their target interactions with messenger RNAs do not adhere to simple rules. Of the thousands of known animal microRNAs, relatively few microRNA:messenger RNA regulatory interactions have been biologically validated in an normal organismal context. Here we present evidence that three microRNAs from the Hox complex in Drosophila (miR-10-5p, miR-10-3p, miR-iab-4-5p) do not have significant effects during embryogenesis on the expression of Hox genes that contain high confidence microRNAs target sites in the 3′ untranslated regions of their messenger RNAs. This is significant, in that it suggests that many predicted microRNA-target interactions may not be biologically relevant, or that the outcomes of these interactions may be so subtle that mutants may only show phenotypes in specific contexts, such as in environmental stress conditions, or in combinations with other microRNA mutations.