Analysis of Hox gene expression in the chick limb bud

A genetic approach to the transcriptional regulation of Hox gene clusters

The evolution of vertebrate genomes was accompanied by an astounding increase in the complexity of their regulatory modalities. Genetic redundancy resulting from large-scale genome duplications at the base of the chordate tree was repeatedly exploited by the functional redeployment of paralogous genes via innovations in their regulatory circuits. As a paradigm of such regulatory evolution, we have extensively studied those control mechanisms at work in-cis over vertebrate Hox gene clusters. Here, we review the portfolio of genetic strategies that have been developed to tackle the intricate relationship between genomic topography and the transcriptional activities in this gene family, and we describe some of the mechanistic insights we gained by using the HoxD cluster as an example. We discuss the high heuristic value of this system in our general understanding of how changes in transcriptional regulation can diversify gene function and thereby fuel morphological evolution.

A regulatory archipelago controls Hox genes transcription in digits

The evolution of digits was an essential step in the success of tetrapods. Among the key players, Hoxd genes are coordinately regulated in developing digits, where they help organize growth and patterns. We identified the distal regulatory sites associated with these genes by probing the three-dimensional architecture of this regulatory unit in developing limbs. This approach, combined with in vivo deletions of distinct regulatory regions, revealed that the active part of the gene cluster contacts several enhancer-like sequences. These elements are dispersed throughout the nearby gene desert, and each contributes either quantitatively or qualitatively to Hox gene transcription in presumptive digits. We propose that this genetic system, which we call a "regulatory archipelago," provides an inherent flexibility that may partly underlie the diversity in number and morphology of digits across tetrapods, as well as their resilience to drastic variations.

HOXA13 and HOXD13 expression during development of the syndactylous digits in the marsupial Macropus eugenii

Background

Kangaroos and wallabies have specialised limbs that allow for their hopping mode of locomotion. The hindlimbs differentiate much later in development but become much larger than the forelimbs. The hindlimb autopod has only four digits, the fourth of which is greatly elongated, while digits two and three are syndactylous. We investigated the expression of two genes, HOXA13 and HOXD13, that are crucial for digit patterning in mice during formation of the limbs of the tammar wallaby.

Results

We describe the development of the tammar limbs at key stages before birth. There was marked heterochrony and the hindlimb developed more slowly than the forelimb. Both tammar HOXA13 and HOXD13 have two exons as in humans, mice and chickens. HOXA13 had an early and distal mRNA distribution in the tammar limb bud as in the mouse, but forelimb expression preceded that in the hindlimb. HOXD13 mRNA was expressed earlier in the forelimb than the hindlimb and was predominantly detected in the interdigital tissues of the forelimb. In contrast, the hindlimb had a more restricted expression pattern that appeared to be expressed at discrete points at both posterior and anterior margins of the limb bud, and was unlike expression seen in the mouse and the chicken.

Conclusions

This is the first examination of HOXA and HOXD gene expression in a marsupial. The gene structure and predicted proteins were highly conserved with their eutherian orthologues. Interestingly, despite the morphological differences in hindlimb patterning, there were no modifications to the polyalanine tract of either HOXA13 or HOXD13 when compared to those of the mouse and bat but there was a marked difference between the tammar and the other mammals in the region of the first polyserine tract of HOXD13. There were also altered expression domains for both genes in the developing tammar limbs compared to the chicken and mouse. Together these findings suggest that the timing of HOX gene expression may contribute to the heterochrony of the forelimb and hindlimb and that alteration to HOX domains may influence phenotypic differences that lead to the development of marsupial syndactylous digits.

The Devonian tetrapod Acanthostega gunnari Jarvik: postcranial anatomy, basal tetrapod interrelationships and patterns of skeletal evolution

The postcranial skeleton of Acanthostega gunnari from the Famennian of East Greenland displays a unique, transitional, mixture of features conventionally associated with fishand tetrapod-like morphologies. The rhachitomous vertebral column has a primitive, barely differentiated atlas-axis complex, encloses an unconstricted notochordal canal, and the weakly ossified neural arches have poorly developed zygapophyses. More derived axial skeletal features include caudal vertebral proliferation and, transiently, neural radials supporting unbranched and unsegmented lepidotrichia. Sacral and post-sacral ribs reiterate uncinate cervical and anterior thoracic rib morphologies: a simple distal flange supplies a broad surface for iliac attachment. The octodactylous forelimb and hindlimb each articulate with an unsutured, foraminate endoskeletal girdle. A broad-bladed femoral shaft with extreme anterior torsion and associated flattened epipodials indicates a paddle-like hindlimb function. Phylogenetic analysis places Acanthostega as the sister-group of Ichthyostega plus all more advanced tetrapods. Tulerpeton appears to be a basal stemamniote plesion, tying the amphibian-amniote split to the uppermost Devonian. Caerorhachis may represent a more derived stem-amniote plesion. Postcranial evolutionary trends spanning the taxa traditionally associated with the fish-tetrapod transition are discussed in detail. Comparison between axial skeletons of primitive tetrapods suggests that plesiomorphic fish-like morphologies were re-patterned in a cranio-caudal direction with the emergence of tetrapod vertebral regionalisation. The evolution of digited limbs lags behind the initial enlargement of endoskeletal girdles, whereas digit evolution precedes the elaboration of complex carpal and tarsal articulations. Pentadactylous limbs appear to have stabilised independently in amniote and amphibian lineages; the colosteid Greererpeton has a pentadactylous manus, indicating that basal amphibian forelimbs may not be restricted to patterns of four digits or less.

Development and evolution of the lateral plate mesoderm: comparative analysis of amphioxus and lamprey with implications for the acquisition of paired fins

Possession of paired appendages is regarded as a novelty that defines crown gnathostomes and allows sophisticated behavioral and locomotive patterns. During embryonic development, initiation of limb buds in the lateral plate mesoderm involves several steps. First, the lateral plate mesoderm is regionalized into the cardiac mesoderm (CM) and the posterior lateral plate mesoderm (PLPM). Second, in the PLPM, Hox genes are expressed in a collinear manner to establish positional values along the anterior-posterior axis. The developing PLPM splits into somatic and splanchnic layers. In the presumptive limb field of the somatic layer, expression of limb initiation genes appears. To gain insight into the evolutionary sequence leading to the emergence of paired appendages in ancestral vertebrates, we examined the embryonic development of the ventral mesoderm in the cephalochordate amphioxus Branchiostoma floridae and of the lateral plate mesoderm in the agnathan lamprey Lethenteron japonicum, and studied the expression patterns of cognates of genes known to be expressed in these mesodermal layers during amniote development. We observed that, although the amphioxus ventral mesoderm posterior to the pharynx was not regionalized into CM and posterior ventral mesoderm, the lateral plate mesoderm of lampreys was regionalized into CM and PLPM, as in gnathostomes. We also found nested expression of two Hox genes (LjHox5i and LjHox6w) in the PLPM of lamprey embryos. However, histological examination showed that the PLPM of lampreys was not separated into somatic and splanchnic layers. These findings provide insight into the sequential evolutionary changes that occurred in the ancestral lateral plate mesoderm leading to the emergence of paired appendages.

Limb regeneration: a new development?

Salamander limb regeneration is a classical model of tissue morphogenesis and patterning. Through recent advances in cell labeling and molecular analysis, a more precise, mechanistic understanding of this process has started to emerge. Long-standing questions include to what extent limb regeneration recapitulates the events observed in mammalian limb development and to what extent are adult- or salamander- specific aspects deployed. Historically, researchers studying limb development and limb regeneration have proposed different models of pattern formation. Here we discuss recent data on limb regeneration and limb development to argue that although patterning mechanisms are likely to be similar, cell plasticity and signaling from nerves play regeneration-specific roles.

Diffusible signals, not autonomous mechanisms, determine the main proximodistal limb subdivision

Vertebrate limbs develop three main proximodistal (PD) segments (upper arm, forearm, and hand) in a proximal-to-distal sequence. Despite extensive research into limb development, whether PD specification occurs through nonautonomous or autonomous mechanisms is not resolved. Heterotopic transplantation of intact and recombinant chicken limb buds identifies signals in the embryo trunk that proximalize distal limb cells to generate a complete PD axis. In these transplants, retinoic acid induces proximalization, which is counteracted by fibroblast growth factors from the distal limb bud; these related actions suggest that the first limb-bud PD regionalization results from the balance between proximal and distal signals. The plasticity of limb progenitor cell identity in response to diffusible signals provides a unifying view of PD patterning during vertebrate limb development and regeneration.

Initiation of proximal-distal patterning in the vertebrate limb by signals and growth

Two broad classes of models have been proposed to explain the patterning of the proximal-distal axis of the vertebrate limb (from the shoulder to the digit tips). Differentiating between them, we demonstrate that early limb mesenchyme in the chick is initially maintained in a state capable of generating all limb segments through exposure to a combination of proximal and distal signals. As the limb bud grows, the proximal limb is established through continued exposure to flank-derived signal(s), whereas the developmental program determining the medial and distal segments is initiated in domains that grow beyond proximal influence. In addition, the system we have developed, combining in vitro and in vivo culture, opens the door to a new level of analysis of patterning mechanisms in the limb.

Physical forces may cause Hox gene collinearity in the primary and secondary axes of the developing vertebrates

The features of spatial and temporal Hox gene collinearity along the anteroposterior and secondary axes of vertebrate development have been extensively studied. However, the understanding of these features remains problematic. Some genetic engineering experiments were performed and the consequent modifications of the Hoxd gene expressions in the vertebrate limb and trunk were presented. A two-phases model was proposed to describe the above results but still many data cannot be explained. In the present work a different mechanism is put forward in order to deal with the above experiments. This alternative mechanism (coined biophysical model), is based on the hypothesis that physical forces decondense and 'loop out' the chromatin fiber causing the observed Hox gene collinearity phenomena at the early stages of axonal development. The two models are compared in detail. The biophysical model adequately explains the data even in cases where the results are characterized as unexpected. Furthermore, the biophysical model predicts that the Hox gene expressions are entangled in space and time and this coupling is compatible with the data of the early developmental stages. Additional experiments are proposed for a direct test of this model.

Comparative analysis of 3D expression patterns of transcription factor genes and digit fate maps in the developing chick wing

Hoxd13, Tbx2, Tbx3, Sall1 and Sall3 genes are candidates for encoding antero-posterior positional values in the developing chick wing and specifying digit identity. In order to build up a detailed profile of gene expression patterns in cell lineages that give rise to each of the digits over time, we compared 3 dimensional (3D) expression patterns of these genes during wing development and related them to digit fate maps. 3D gene expression data at stages 21, 24 and 27 spanning early bud to digital plate formation, captured from in situ hybridisation whole mounts using Optical Projection Tomography (OPT) were mapped to reference wing bud models. Grafts of wing bud tissue from GFP chicken embryos were used to fate map regions of the wing bud giving rise to each digit; 3D images of the grafts were captured using OPT and mapped on to the same models. Computational analysis of the combined computerised data revealed that Tbx2 and Tbx3 are expressed in digit 3 and 4 progenitors at all stages, consistent with encoding stable antero-posterior positional values established in the early bud; Hoxd13 and Sall1 expression is more dynamic, being associated with posterior digit 3 and 4 progenitors in the early bud but later becoming associated with anterior digit 2 progenitors in the digital plate. Sox9 expression in digit condensations lies within domains of digit progenitors defined by fate mapping; digit 3 condensations express Hoxd13 and Sall1, digit 4 condensations Hoxd13, Tbx3 and to a lesser extent Tbx2. Sall3 is only transiently expressed in digit 3 progenitors at stage 24 together with Sall1 and Hoxd13; then becomes excluded from the digital plate. These dynamic patterns of expression suggest that these genes may play different roles in digit identity either together or in combination at different stages including the digit condensation stage.

Hox expression in the American alligator and evolution of archosaurian axial patterning

The avian body plan has undergone many modifications, most associated with adaptation to flight and bipedal walking. Some of these modifications may be owing to avian-specific changes in the embryonic Hox expression code. Here, we have examined Hox expression in alligator, the closest living relative of birds, and an archosaur with a more conservative body plan. Two differences in Hox expression between chick, alligator, and other tetrapods correlate with aspects of alligator or bird-specific skeletal morphology. First, absence of a thoracic subdomain of Hoxc-8 expression in alligator correlates with morphological adaptations in crocodilian thoracic segments. Second, Hoxa-5, a gene required to pattern the cervical-thoracic transition, shows unique patterns of expression in chick, alligator, and mouse, correlating with species-specific morphological patterning of this region. Given that cervical vertebral morphologies evolved independently in the bird and mammalian lineages, the underlying developmental mechanisms, including refinement of Hox expression domains, may be distinct.

SHH propagates distal limb bud development by enhancing CYP26B1-mediated retinoic acid clearance via AER-FGF signalling

The essential roles of SHH in anteroposterior (AP) and AER-FGF signalling in proximodistal (PD) limb bud development are well understood. In addition, these morphoregulatory signals are key components of the self-regulatory SHH/GREM1/AER-FGF feedback signalling system that regulates distal progression of limb bud development. This study uncovers an additional signalling module required for coordinated progression of limb bud axis development. Transcriptome analysis using Shh-deficient mouse limb buds revealed that the expression of proximal genes was distally extended from early stages onwards, which pointed to a more prominent involvement of SHH in PD limb axis development. In particular, retinoic acid (RA) target genes were upregulated proximally, while the expression of the RA-inactivating Cyp26b1 enzyme was downregulated distally, pointing to increased RA activity in Shh-deficient mouse limb buds. Further genetic and molecular analysis established that Cyp26b1 expression is regulated by AER-FGF signalling. During initiation of limb bud outgrowth, the activation of Cyp26b1 expression creates a distal 'RA-free' domain, as indicated by complementary downregulation of a transcriptional sensor of RA activity. Subsequently, Cyp26b1 expression increases as a consequence of SHH-dependent upregulation of AER-FGF signalling. To better understand the underlying signalling interactions, computational simulations of the spatiotemporal expression patterns and interactions were generated. These simulations predicted the existence of an antagonistic AER-FGF/CYP26B1/RA signalling module, which was verified experimentally. In summary, SHH promotes distal progression of limb development by enhancing CYP26B1-mediated RA clearance as part of a signalling network linking the SHH/GREM1/AER-FGF feedback loop to the newly identified AER-FGF/CYP26B1/RA module.

Salamander limb development: integrating genes, morphology, and fossils

The development of the tetrapod limb during skeletogenesis follows a highly conservative pattern characterized by a general proximo-distal progression in the establishment of skeletal elements and a postaxial polarity in digit development. Salamanders represent the only exception to this pattern and display an early establishment of distal autopodial structures, specifically the basale commune, an amalgamation of distal carpal and tarsal 1 and 2, and a distinct preaxial polarity in digit development. This deviance from the conserved tetrapod pattern has resulted in a number of hypotheses to explain its developmental basis and evolutionary history. Here we summarize the current knowledge of salamander limb development under consideration of the fossil record to provide a deep time perspective of this evolutionary pathway and highlight what data will be needed in the future to gain a better understanding of salamander limb development specifically and tetrapod limb development and evolution more broadly.

In search of the Golden Fleece: unraveling principles of morphogenesis by studying the integrative biology of skin appendages

The mythological story of the Golden Fleece symbolizes the magical regenerative power of skin appendages. Similar to the adventurous pursuit of the Golden Fleece by the multi-talented Argonauts, today we also need an integrated multi-disciplined approach to understand the cellular and molecular processes during development, regeneration and evolution of skin appendages. To this end, we have explored several aspects of skin appendage biology that contribute to the Turing activator/inhibitor model in feather pattern formation, the topo-biological arrangement of stem cells in organ shape determination, the macro-environmental regulation of stem cells in regenerative hair waves, and potential novel molecular pathways in the morphological evolution of feathers. Here we show our current integrative biology efforts to unravel the complex cellular behavior in patterning stem cells and the control of regional specificity in skin appendages. We use feather/scale tissue recombination to demonstrate the timing control of competence and inducibility. Feathers from different body regions are used to study skin regional specificity. Bioinformatic analyses of transcriptome microarrays show the potential involvement of candidate molecular pathways. We further show Hox genes exhibit some region specific expression patterns. To visualize real time events, we applied time-lapse movies, confocal microscopy and multiphoton microscopy to analyze the morphogenesis of cultured embryonic chicken skin explants. These modern imaging technologies reveal unexpectedly complex cellular flow and organization of extracellular matrix molecules in three dimensions. While these approaches are in preliminary stages, this perspective highlights the challenges we face and new integrative tools we will use. Future work will follow these leads to develop a systems biology view and understanding in the morphogenetic principles that govern the development and regeneration of ectodermal organs.

A long noncoding RNA maintains active chromatin to coordinate homeotic gene expression

The genome is extensively transcribed into long intergenic noncoding RNAs (lincRNAs), many of which are implicated in gene silencing. Potential roles of lincRNAs in gene activation are much less understood. Development and homeostasis require coordinate regulation of neighbouring genes through a process termed locus control. Some locus control elements and enhancers transcribe lincRNAs, hinting at possible roles in long-range control. In vertebrates, 39 Hox genes, encoding homeodomain transcription factors critical for positional identity, are clustered in four chromosomal loci; the Hox genes are expressed in nested anterior-posterior and proximal-distal patterns colinear with their genomic position from 3' to 5'of the cluster. Here we identify HOTTIP, a lincRNA transcribed from the 5' tip of the HOXA locus that coordinates the activation of several 5' HOXA genes in vivo. Chromosomal looping brings HOTTIP into close proximity to its target genes. HOTTIP RNA binds the adaptor protein WDR5 directly and targets WDR5/MLL complexes across HOXA, driving histone H3 lysine 4 trimethylation and gene transcription. Induced proximity is necessary and sufficient for HOTTIP RNA activation of its target genes. Thus, by serving as key intermediates that transmit information from higher order chromosomal looping into chromatin modifications, lincRNAs may organize chromatin domains to coordinate long-range gene activation.

Evolution of vertebrate appendicular structures: Insight from genetic and palaeontological data

The new body of evidence from fossils and comparative-developmental analysis of subset of appendicular patterning genes has revealed that limb elements seen in tetrapods are assembled in fish fin over evolutionary time. However, despite of deep homology in basic structure and underlying developmental system, there remains a large morphological gap between distal elements of tetrapod limb and distal fin skeleton of tetrapodomorph fish. Understanding the genetic basis of major transformations in distal-limb morphology is the next challenge for evolutionary developmental biologists. Here by integrating data from fossils, comparative-developmental and genetic studies, models are proposed describing the evolution of cis-regulatory elements as a basis for diversification of appendicular architecture. Instead of emphasizing the subset of developmental genes, for instance Hoxd genes, the focus here is on the significance of elucidating cis-regulatory elements for multiple other key molecular players of limb/fin development and genetic/molecular interactions among them, for a better understanding of the developmental and genetic basis of limb evolution.

Specific age-associated DNA methylation changes in human dermal fibroblasts

Epigenetic modifications of cytosine residues in the DNA play a critical role for cellular differentiation and potentially also for aging. In mesenchymal stromal cells (MSC) from human bone marrow we have previously demonstrated age-associated methylation changes at specific CpG-sites of developmental genes. In continuation of this work, we have now isolated human dermal fibroblasts from young (<23 years) and elderly donors (>60 years) for comparison of their DNA methylation profiles using the Infinium HumanMethylation27 assay. In contrast to MSC, fibroblasts could not be induced towards adipogenic, osteogenic and chondrogenic lineage and this is reflected by highly significant differences between the two cell types: 766 CpG sites were hyper-methylated and 752 CpG sites were hypo-methylated in fibroblasts in comparison to MSC. Strikingly, global DNA methylation profiles of fibroblasts from the same dermal region clustered closely together indicating that fibroblasts maintain positional memory even after in vitro culture. 75 CpG sites were more than 15% differentially methylated in fibroblasts upon aging. Very high hyper-methylation was observed in the aged group within the INK4A/ARF/INK4b locus and this was validated by pyrosequencing. Age-associated DNA methylation changes were related in fibroblasts and MSC but they were often regulated in opposite directions between the two cell types. In contrast, long-term culture associated changes were very consistent in fibroblasts and MSC. Epigenetic modifications at specific CpG sites support the notion that aging represents a coordinated developmental mechanism that seems to be regulated in a cell type specific manner.

Spatiotemporal changes in cell adhesiveness during vertebrate limb morphogenesis

During vertebrate limb development, various molecules are expressed in the presumptive limb field or the limb bud in a spatiotemporal-specific manner. The combination of these molecules regulates cellular properties that affect limb initiation and its morphogenesis, especially cartilage formation. Cell adhesiveness of the limb mesenchyme is a key factor in the regulation of cell distribution. Differential adhesiveness of mesenchymal cells is first observed between cells in the presumptive limb field and flank region, and the adhesiveness of the cells in the limb field is higher than that of cells in the flank region. In the limb bud, the adhesiveness of mesenchymal cells shows spatiotemporal difference, which reflects the positional identity of the cells. Position-dependent cell adhesiveness is also observed in blastema cells of the regenerating limb. Therefore, local changes in cell adhesiveness are observed during limb development and regeneration, suggesting significant roles for cell adhesiveness in limb morphogenesis.

A regulatory 'landscape effect' over the HoxD cluster

Faithful expression of Hox genes in both time and space is essential for proper patterning of the primary body axis. Transgenic approaches in vertebrates have suggested that this collinear activation process is regulated in a largely gene cluster-autonomous manner. In contrast, more recently co-opted expression specificities, required in other embryonic structures, depend upon long-range enhancer sequences acting from outside the gene clusters. This regulatory dichotomy was recently questioned, since gene activation along the trunk seems to be partially regulated by signals located outside of the cluster. We investigated these alternative regulatory strategies by engineering a large inversion that precisely separates the murine HoxD complex from its centromeric neighborhood. Mutant animals displayed posterior transformations along with subtle deregulations of Hoxd genes, indicating an impact of the centromeric landscape on the fine-tuning of Hoxd gene expression. Proximal limbs were also affected, suggesting that this 'landscape effect' is generic and impacts upon regulatory mechanisms of various qualities and evolutionary origins.

Wdr5 is required for chick skeletal development

Wdr5, a bone morphogenetic protein 2 (BMP-2)-induced protein belonging to the family of the WD repeat proteins, is expressed in proliferating and hypertrophic chondrocytes of the growth plate and in osteoblasts. Although previous studies have provided insight into the mechanisms by which Wdr5 affects chondrocyte and osteoblast differentiation, whether Wdr5 is required in vivo for endochondral bone development has not been addressed. In this study, using an avian replication competent retrovirus (RCAS) system delivering Wdr5 short hairpin (sh) RNA to silence Wdr5 in the developing limb, we report that reduction of Wdr5 levels delays endochondral bone development and consequently results in shortening of the skeletal elements. Shortening of the skeletal elements was due to impaired chondrocyte maturation, evidenced by a significant reduction of Runx2, type X collagen, and osteopontin expression. A decrease in Runx2, type collagen I, and ostepontin expression in osteoblasts and a subsequent defect in mineralized bone was observed as well when Wdr5 levels were reduced. Most important, retroviral misexpression of Runx2 rescued the phenotype induced by Wdr5 shRNA. These findings suggest that during limb development, Wdr5 is required for endochondral bone formation and that Wdr5 influences this process, at least in part, by regulating Runx2 expression.

Identification of genes downstream of the Shh signalling in the developing chick wing and syn-expressed with Hoxd13 using microarray and 3D computational analysis

Sonic hedgehog (Shh) signalling by the polarizing region at the posterior margin of the chick wing bud is pivotal in patterning the digits but apart from a few key downstream genes, such as Hoxd13, which is expressed in the posterior region of the wing that gives rise to the digits, the genes that mediate the response to Shh signalling are not known. To find genes that are co-expressed with Hoxd13 in the posterior of chick wing buds and regulated in the same way, we used microarrays to compare gene expression between anterior and posterior thirds of wing buds from normal chick embryos and from polydactylous talpid³ mutant chick embryos, which have defective Shh signalling due to lack of primary cilia. We identified 1070 differentially expressed gene transcripts, which were then clustered. Two clusters contained genes predominantly expressed in posterior thirds of normal wing buds; in one cluster, genes including Hoxd13, were expressed at high levels in anterior and posterior thirds in talpid³ wing buds, in the other cluster, genes including Ptc1, were expressed at low levels in anterior and posterior thirds in talpid³ wing buds. Expression patterns of genes in these two clusters were validated in normal and talpid³ mutant wing buds by in situ hybridisation and demonstrated to be responsive to application of Shh. Expression of several genes in the Hoxd13 cluster was also shown to be responsive to manipulation of protein kinase A (PKA) activity, thus demonstrating regulation by Gli repression. Genes in the Hoxd13 cluster were then sub-clustered by computational comparison of 3D expression patterns in normal wing buds to produce syn-expression groups. Hoxd13 and Sall1 are syn-expressed in the posterior region of early chick wing buds together with 6 novel genes which are likely to be functionally related and represent secondary targets of Shh signalling. Other groups of syn-expressed genes were also identified, including a group of genes involved in vascularisation.

Predicting embryonic patterning using mutual entropy fitness and in silico evolution

During vertebrate embryogenesis, the expression of Hox genes that define anterior-posterior identity follows general rules: temporal colinearity and posterior prevalence. A mathematical measure for the quality or fitness of the embryonic pattern produced by a gene regulatory network is derived. Using this measure and in silico evolution we derive gene interaction networks for anterior-posterior (AP) patterning under two developmental paradigms. For patterning during growth (paradigm I), which is appropriate for vertebrates and short germ-band insects, the algorithm creates gene expression patterns reminiscent of Hox gene expression. The networks operate through a timer gene, the level of which measures developmental progression (a candidate is the widely conserved posterior morphogen Caudal). The timer gene provides a simple mechanism to coordinate patterning with growth rate. The timer, when expressed as a static spatial gradient, functions as a classical morphogen (paradigm II), providing a natural way to derive the AP patterning, as seen in long germ-band insects that express their Hox genes simultaneously, from the ancestral short germ-band system. Although the biochemistry of Hox regulation in higher vertebrates is complex, the actual spatiotemporal expression phenotype is not, and simple activation and repression by Hill functions suffices in our model. In silico evolution provides a quantitative demonstration that continuous positive selection can generate complex phenotypes from simple components by incremental evolution, as Darwin proposed.

Comparison of Iroquois gene expression in limbs/fins of vertebrate embryos

In Drosophila, Iroquois (Irx) genes have various functions including the specification of the identity of wing veins. Vertebrate Iroquois (Irx) genes have been reported to be expressed in the developing digits of mouse limbs. Here we carry out a phylogenetic analysis of vertebrate Irx genes and compare expression in developing limbs of mouse, chick and human embryos and in zebrafish pectoral fin buds. We confirm that the six Irx gene families in vertebrates are well defined and that Clusters A and B are duplicates; in contrast, Irx1 and 3, Irx2 and 5, and Irx4 and 6 are paralogs. All Irx genes in mouse and chick are expressed in developing limbs. Detailed comparison of the expression patterns in mouse and chick shows that expression patterns of genes in the same cluster are generally similar but paralogous genes have different expression patterns. Mouse and chick Irx1 are expressed in digit condensations, whereas mouse and chick Irx6 are expressed interdigitally. The timing of Irx1 expression in individual digits in mouse and chick is different. Irx1 is also expressed in digit condensations in developing human limbs, thus showing conservation of expression of this gene in higher vertebrates. In zebrafish, Irx genes of all but six of the families are expressed in early stage pectoral fin buds but not at later stages, suggesting that these genes are not involved in patterning distal structures in zebrafish fins.

Bare bones pattern formation: a core regulatory network in varying geometries reproduces major features of vertebrate limb development and evolution

BACKGROUND

Major unresolved questions regarding vertebrate limb development concern how the numbers of skeletal elements along the proximodistal (P-D) and anteroposterior (A-P) axes are determined and how the shape of a growing limb affects skeletal element formation. There is currently no generally accepted model for these patterning processes, but recent work on cartilage development (chondrogenesis) indicates that precartilage tissue self-organizes into nodular patterns by cell-molecular circuitry with local auto-activating and lateral inhibitory (LALI) properties. This process is played out in the developing limb in the context of a gradient of fibroblast growth factor (FGF) emanating from the apical ectodermal ridge (AER).

RESULTS

We have simulated the behavior of the core chondrogenic mechanism of the developing limb in the presence of an FGF gradient using a novel computational environment that permits simulation of LALI systems in domains of varying shape and size. The model predicts the normal proximodistal pattern of skeletogenesis as well as distal truncations resulting from AER removal. Modifications of the model's parameters corresponding to plausible effects of Hox proteins and formins, and of the reshaping of the model limb, bud yielded simulated phenotypes resembling mutational and experimental variants of the limb. Hypothetical developmental scenarios reproduce skeletal morphologies with features of fossil limbs.

CONCLUSIONS

The limb chondrogenic regulatory system operating in the presence of a gradient has an inherent, robust propensity to form limb-like skeletal structures. The bare bones framework can accommodate ancillary gene regulatory networks controlling limb bud shaping and establishment of Hox expression domains. This mechanism accounts for major features of the normal limb pattern and, under variant geometries and different parameter values, those of experimentally manipulated, genetically aberrant and evolutionary early forms, with no requirement for an independent system of positional information.

Control of vertebrate Hox clusters by remote and global cis-acting regulatory sequences

Despite apparently shared structural organisation and functional roles, vertebrate Hox genes are controlled by regulatory mechanisms rather distinct from those of the prototypic Drosophila Antennapedia (ANT-C) and Bithorax (BX-C) Complexes. If individual regulatory modules have been shown to recapitulate specific Hox expression patterns, other experimental studies underscore that vertebrate Hox clusters are controlled in many of their functions in a global manner, through distinct mechanisms. We will discuss the different models that have been proposed to account for these global regulatory modes. In this context, the studies of the regulation of the HoxD complex during limb development highlighted the role of global regulatory elements and the different mechanisms associated to transform a structural organisation into distinct temporal and spatial expression domains. We will further discuss how these mechanisms may have benefited from the structure of the vertebrate homeotic clusters and reciprocally contribute to shape their evolution towards an increased level of organisation and compaction.

Analysis of hoxa11 and hoxa13 expression during patternless limb regeneration in Xenopus

During limb regeneration, anuran tadpoles and urodele amphibians generate pattern-organizing, multipotent, mesenchymal blastema cells, which give rise to a replica of the lost limb including patterning in three dimensions. To facilitate the regeneration of nonregenerative limbs in other vertebrates, it is important to elucidate the molecular differences between blastema cells that can regenerate the pattern of limbs and those that cannot. In Xenopus froglet (soon after metamorphosis), an amputated limb generates blastema cells that do not produce proper patterning, resulting in a patternless regenerate, a spike, regardless of the amputation level. We found that re-expression of hoxa11 and hoxa13 in the froglet blastema is initiated although the subsequent proximal-distal patterning, including separation of the hoxa11 and hoxa13 expression domains, is disrupted. We also observed an absence of EphA4 gene expression in the froglet blastema and a failure of position-dependent cell sorting, which correlated with the altered hoxa11 and hoxa13 expression. Quantitative analysis of hoxa11 and hoxa13 expression revealed that hoxa13 transcript levels were reduced in the froglet blastema compared with the tadpole blastema. Moreover, the expression of sox9, an important regulator of chondrogenic differentiation, was detected earlier in patternless blastemas than in tadpole blastemas. These results suggest that appropriate spatial, temporal, and quantitative gene expression is necessary for pattern regeneration by blastema cells.

A reevaluation of X-irradiation-induced phocomelia and proximodistal limb patterning

Phocomelia is a devastating, rare congenital limb malformation in which the long bones are shorter than normal, with the upper portion of the limb being most severely affected. In extreme cases, the hands or fingers are attached directly to the shoulder and the most proximal elements (those closest to the shoulder) are entirely missing. This disorder, previously known in both autosomal recessive and sporadic forms, showed a marked increase in incidence in the early 1960s due to the tragic toxicological effects of the drug thalidomide, which had been prescribed as a mild sedative. This human birth defect is mimicked in developing chick limb buds exposed to X-irradiation. Both X-irradiation and thalidomide-induced phocomelia have been interpreted as patterning defects in the context of the progress zone model, which states that a cell's proximodistal identity is determined by the length of time spent in a distal limb region termed the 'progress zone'. Indeed, studies of X-irradiation-induced phocomelia have served as one of the two major experimental lines of evidence supporting the validity of the progress zone model. Here, using a combination of molecular analysis and lineage tracing in chick, we show that X-irradiation-induced phocomelia is fundamentally not a patterning defect, but rather results from a time-dependent loss of skeletal progenitors. Because skeletal condensation proceeds from the shoulder to fingers (in a proximal to distal direction), the proximal elements are differentially affected in limb buds exposed to radiation at early stages. This conclusion changes the framework for considering the effect of thalidomide and other forms of phocomelia, suggesting the possibility that the aetiology lies not in a defect in the patterning process, but rather in progenitor cell survival and differentiation. Moreover, molecular evidence that proximodistal patterning is unaffected after X-irradiation does not support the predictions of the progress zone model.

Duplication

Congenital limb duplications include pre- and post-axial polydactyly, central polydactyly, and the mirror-hand spectrum. Treatment of these duplications constitutes a significant functional and aesthetic challenge for the reconstructive hand surgeon. This article provides an inclusive review of the embryologic and molecular mechanisms underlying these deformities and focuses on their clinical treatment. The anatomic variation, classification, surgical treatment, and outcomes of surgical intervention are reviewed for each of the disorders of duplication.

MicroRNAs in the Hox network: an apparent link to posterior prevalence

Homeobox (Hox) transcription factors confer anterior-posterior (AP) axial coordinates to vertebrate embryos. Hox genes are found in clusters that also contain genes for microRNAs (miRNAs). Our analysis of predicted miRNA targets indicates that Hox cluster-embedded miRNAs preferentially target Hox mRNAs. Moreover, the presumed Hox target genes are predominantly situated on the 3' side of each Hox miRNA locus. These results suggest that Hox miRNAs help repress more anterior programmes, thereby reinforcing posterior prevalence, which is the hierarchical dominance of posterior over anterior Hox gene function that is observed in bilaterians. In this way, miRNA-mediated regulation seems to recapitulate interactions at other levels of gene expression, some more ancestral, within a network under stabilizing selection.

The evolution of HoxD-11 expression in the bird wing: insights from Alligator mississippiensis

Background

Comparative morphology identifies the digits of the wing of birds as 1,2 and 3, but they develop at embryological positions that become digits 2, 3 and 4 in other amniotes. A hypothesis to explain this is that a homeotic frame shift of digital identity occurred in the evolution of the bird wing, such that digits 1,2 and 3 are developing from embryological positions 2, 3 and 4. Digit 1 of the mouse is the only digit that shows no late expression of HoxD-11. This is also true for the anterior digit of the bird wing, suggesting this digit is actually a digit 1. If this is the case, we can expect closer relatives of birds to show no HoxD-11 expression only in digit 1. To test this prediction we investigate HoxD-11 expression in crocodilians, the closest living relatives of birds.

Methodology/Principal Findings

Using degenerate primers we cloned a 606 nucleotide fragment of exon 1 of the alligator HoxD-11 gene and used it for whole-mount in-situ detection in alligator embryos. We found that in the pentadactyl forelimbs of alligator, as in the mouse, late expression of HoxD-11 is absent only in digit 1.

Conclusions/Significance

The ancestral condition for amniotes is that late-phase HoxD-11 expression is absent only in digit 1. The biphalangeal morphology and lack of HoxD-11 expression of the anterior digit of the wing is like digit 1 of alligator and mouse, but its embryological position as digit 2 is derived. HoxD-11 expression in alligator is consistent with the hypothesis that both digit morphology as well as HoxD-11 expression are shifted towards posterior in the bird wing.

Tri-phasic expression of posterior Hox genes during development of pectoral fins in zebrafish: implications for the evolution of vertebrate paired appendages

During development of the limbs, Hox genes belonging to the paralogous groups 9-13 are expressed in three distinct phases, which play key roles in the segmental patterning of limb skeletons. In teleost fishes, which have a very different organization in their fin skeletons, it is not clear whether a similar patterning mechanism is at work. To determine whether Hox genes are also expressed in several distinct phases during teleost paired fin development, we re-analyzed the expression patterns of hox9-13 genes during development of pectoral fins in zebrafish. We found that, similar to tetrapod Hox genes, expression of hoxa/d genes in zebrafish pectoral fins occurs in three distinct phases, in which the most distal/third phase is correlated with the development of the most distal structure of the fin, the fin blade. Like in tetrapods, hox gene expression in zebrafish pectoral fins during the distal/third phase is dependent upon sonic hedgehog signaling (hoxa and hoxd genes) and the presence of a long-range enhancer (hoxa genes), which indicates that the regulatory mechanisms underlying tri-phasic expression of Hox genes have remained relatively unchanged during evolution. Our results suggest that, although simpler in organization, teleost fins do have a distal structure that might be considered comparable to the autopod region of limbs.

Modeling Hox gene regulation in digits: reverse collinearity and the molecular origin of thumbness

During the development of mammalian digits, clustered Hoxd genes are expressed following a collinear regulatory strategy, leading to both the growth of digits and their morphological identities. Because gene dosage is a key parameter in this system, we used a quantitative approach, associated with a collection of mutant stocks, to investigate the nature of the underlying regulatory mechanism(s). In parallel, we elaborated a mathematical model of quantitative collinearity, which was progressively challenged and validated by the experimental approach. This combined effort suggested a two-step mechanism, which involves initially the looping and recognition of the cluster by a complex including two enhancer sequences, followed by a second step of microscanning of genes located nearby. In this scenario, the respective rank of the genes, with respect to the 5' extremity of the cluster, is primordial, as well as different gene-specific affinities. This model accounts for the quantitative variations observed in our many mutant strains, and reveals the molecular constraint leading to thumbness; i.e., why a morphological difference must occur between the most anterior digit and the others. We also show that the same model applies to the collinear regulation of Hox genes during the emergence of external genitalia, though with some differences likely illustrating the distinct functionalities of these structures in adults.

Hoxd13 binds in vivo and regulates the expression of genes acting in key pathways for early limb and skeletal patterning

5' HoxD genes are required for the correct formation of limb skeletal elements. Hoxd13, the most 5'-located HoxD gene, is important for patterning the most distal limb region, and its mutation causes human limb malformation syndromes. The mechanisms underlying the control of developmental processes by Hoxd13, and by Hox genes in general, are still elusive, due to the limited knowledge on their direct downstream target genes. We identified by ChIP-on-chip 248 known gene loci bound invivo by Hoxd13. Genes relevant to limb patterning and skeletogenesis were further analysed. We found that Hoxd13 binds invivo, in developing limbs, the loci of Hand2, a gene crucial to limb AP axis patterning, of Meis1 and Meis2, involved in PD patterning, of the Sfrp1, Barx1, and Fbn1 genes, involved in skeletogenesis, and of the Dach1, Bmp2, Bmp4, andEmx2 genes. We show that Hoxd13 misexpression in developing chick limbs alters the expression of the majority of these genes, supporting the conclusion that Hoxd13 directly regulates their transcription. Our results indicate that 5' Hox proteins regulate directly both key genes for early limb AP and PD axis patterning and genes involved, at later stages, in skeletal patterning.

The migrating gubernaculum grows like a "limb bud"

BACKGROUND

The gubernaculum is crucial for testicular descent, and in the second, or inguinoscrotal, phase of descent it has no caudal attachments. Cranially, it is attached to the testis, but its caudal free tip migrates to the scrotum controlled by the genitofemoral nerve. Recent studies show active proliferation in the tip. We hypothesized that the gubernacular tip may grow like a limb bud.

METHODS

We performed whole-mount in situ hybridization studies on male and female fetal mice (ages, E14.5-E18.5; n = 162) looking for limb bud regulatory factors.

RESULTS

Our results showed that a member of the fibroblast growth factor (Fgf) family, Fgf10, and Hoxa10 were both expressed in the male gubernaculum at E14.5, and Hoxa10 was also expressed in the E16.5 mice. Weak staining was seen in the female gubernaculum for Hoxa10 on days E14.5 and E16.5, whereas no staining for Fgf10 was seen in the female gubernaculums.

CONCLUSIONS

These studies, although preliminary, suggest limb bud regulators are essential for gubernacular growth. Hox genes and Fgfs may be fruitful areas of research to unravel the molecular control of gubernacular migration during testicular descent.

Hox patterning of the vertebrate axial skeleton

The axial skeleton in all vertebrates is composed of similar components that extend from anterior to posterior along the body axis: the occipital skull bones and cervical, thoracic, lumbar, sacral, and caudal vertebrae. Despite significant changes in the number and size of these elements during evolution, the basic character of these anatomical elements, as well as the order in which they appear in vertebrate skeletons, have remained remarkably similar. Through extensive expression analyses, classic morphological perturbation experiments in chicken and targeted loss-of-function analyses in mice, Hox genes have proven to be critical regulators in the establishment of axial skeleton morphology. The convergence of these studies to date allows an emerging understanding of Hox gene function in patterning the vertebrate axial skeleton. This review summarizes genetic and embryologic findings regarding the role of Hox genes in establishing axial morphology and how these combined results impact our current understanding of the vertebrate Hox code.