Gene transpositions in the HoxD complex reveal a hierarchy of regulatory controls

Hox genes in development and beyond

Hox genes encode evolutionarily conserved transcription factors that are essential for the proper development of bilaterian organisms. Hox genes are unique because they are spatially and temporally regulated during development in a manner that is dictated by their tightly linked genomic organization. Although their genetic function during embryonic development has been interrogated, less is known about how these transcription factors regulate downstream genes to direct morphogenetic events. Moreover, the continued expression and function of Hox genes at postnatal and adult stages highlights crucial roles for these genes throughout the life of an organism. Here, we provide an overview of Hox genes, highlighting their evolutionary history, their unique genomic organization and how this impacts the regulation of their expression, what is known about their protein structure, and their deployment in development and beyond.

Hox dosage and morphological diversification during development and evolution

Hox genes encode for evolutionary conserved transcription factors that have long fascinated biologists since the observation of the first homeotic transformations in flies. Hox genes are developmental architects that instruct the formation of various and precise morphologies along the body axes in cnidarian and bilaterian species. In contrast to these highly specific developmental functions, Hox genes encode for proteins that display poorly selective DNA-binding properties in vitro. This "Hox paradox" has been partially solved with the discovery of the TALE-class cofactors, which interact with all Hox members and form versatile Hox/TALE protein complexes on DNA. Here, we describe the role of the Hox dosage as an additional molecular strategy contributing to further resolve the Hox paradox. We present several cases where the Hox dosage is involved in the formation of different morphologies in invertebrates and vertebrates, with a particular emphasis on flight appendages in insects. We also discuss how the Hox dosage could be interpreted in different types of target enhancers within the nuclear environment in vivo. Altogether our survey underlines the Hox dosage as a key mechanism for shaping Hox molecular function during development and evolution.

A novel microdeletion upstream of HOXD13 in a Chinese family with synpolydactyly

Synpolydactyly (SPD) is a digital malformation with the typical clinical phenotype of the webbing of 3/4 fingers and/or 4/5 toes, and combined with polydactyly. In this study, we investigated a Chinese family with SPD and genetic analysis found that all of the affected individuals in the family carry a heterozygous 11,451 bp microdeletion at chr2:176933872-176945322 (GRCh37), which is located upstream of HOXD13 gene, the known disease gene for SPD1. All the affected individuals in the family carry the heterozygous deletion variant, and the variant co-segregated with SPD in the family. Thus, we speculate that the 11,451 bp microdeletion is the disease-causing variant in the family. To date, the microdeletion associating with SPD1 which we identified is the smallest deletion upstream of the HOXD13 gene and not altering the sequence of the HOXD13 gene.

Mesomelic dysplasias associated with the HOXD locus are caused by regulatory reallocations

Human families with chromosomal rearrangements at 2q31, where the human HOXD locus maps, display mesomelic dysplasia, a severe shortening and bending of the limb. In mice, the dominant Ulnaless inversion of the HoxD cluster produces a similar phenotype suggesting the same origin for these malformations in humans and mice. Here we engineer 1 Mb inversion including the HoxD gene cluster, which positioned Hoxd13 close to proximal limb enhancers. Using this model, we show that these enhancers contact and activate Hoxd13 in proximal cells, inducing the formation of mesomelic dysplasia. We show that a secondary Hoxd13 null mutation in-cis with the inversion completely rescues the alterations, demonstrating that ectopic HOXD13 is directly responsible for this bone anomaly. Single-cell expression analysis and evaluation of HOXD13 binding sites suggests that the phenotype arises primarily by acting through genes normally controlled by HOXD13 in distal limb cells. Altogether, these results provide a conceptual and mechanistic framework to understand and unify the molecular origins of human mesomelic dysplasia associated with 2q31.

Fryns type mesomelic dysplasia of the upper limbs caused by inverted duplications of the HOXD gene cluster

The HoxD cluster is critical for vertebrate limb development. Enhancers located in both the telomeric and centromeric gene deserts flanking the cluster regulate the transcription of HoxD genes. In rare patients, duplications, balanced translocations or inversions misregulating HOXD genes are responsible for mesomelic dysplasia of the upper and lower limbs. By aCGH, whole-genome mate-pair sequencing, long-range PCR and fiber fluorescent in situ hybridization, we studied patients from two families displaying mesomelic dysplasia limited to the upper limbs. We identified microduplications including the HOXD cluster and showed that microduplications were in an inverted orientation and inserted between the HOXD cluster and the telomeric enhancers. Our results highlight the existence of an autosomal dominant condition consisting of isolated ulnar dysplasia caused by microduplications inserted between the HOXD cluster and the telomeric enhancers. The duplications likely disconnect the HOXD9 to HOXD11 genes from their regulatory sequences. This presumptive loss-of-function may have contributed to the phenotype. In both cases, however, these rearrangements brought HOXD13 closer to telomeric enhancers, suggesting that the alterations derive from the dominant-negative effect of this digit-specific protein when ectopically expressed during the early development of forearms, through the disruption of topologically associating domain structure at the HOXD locus.

The constrained architecture of mammalian Hox gene clusters

Vertebrate Hox genes are clustered. This organization has a functional relevance, as the transcription of each gene in time and space depends upon its relative position within the gene cluster. Hox clusters display a high organization, and all genes are transcribed from the same DNA strand. Here, we investigate the importance of this uniform transcriptional polarity by engineering alleles where one or several transcription units are inverted, with or without a CTCF site. We observe that inversions are likely detrimental to the proper implementation of this genetic system. We propose that the enhanced organization of Hox clusters in vertebrates evolved in conjunction with the emergence of global gene regulation to optimize a coordinated response of selected subsets of target genes.

In many animal species with a bilateral symmetry, Hox genes are clustered either at one or at several genomic loci. This organization has a functional relevance, as the transcriptional control applied to each gene depends upon its relative position within the gene cluster. It was previously noted that vertebrate Hox clusters display a much higher level of genomic organization than their invertebrate counterparts. The former are always more compact than the latter, they are generally devoid of repeats and of interspersed genes, and all genes are transcribed by the same DNA strand, suggesting that particular factors constrained these clusters toward a tighter structure during the evolution of the vertebrate lineage. Here, we investigate the importance of uniform transcriptional orientation by engineering several alleles within the HoxD cluster, such as to invert one or several transcription units, with or without a neighboring CTCF site. We observe that the association between the tight structure of mammalian Hox clusters and their regulation makes inversions likely detrimental to the proper implementation of this complex genetic system. We propose that the consolidation of Hox clusters in vertebrates, including transcriptional polarity, evolved in conjunction with the emergence of global gene regulation via the flanking regulatory landscapes, to optimize a coordinated response of selected subsets of target genes in cis.

The HoxD cluster is a dynamic and resilient TAD boundary controlling the segregation of antagonistic regulatory landscapes

The mammalian HoxD cluster lies between two topologically associating domains (TADs) matching distinct enhancer-rich regulatory landscapes. During limb development, the telomeric TAD controls the early transcription of Hoxd genes in forearm cells, whereas the centromeric TAD subsequently regulates more posterior Hoxd genes in digit cells. Therefore, the TAD boundary prevents the terminal Hoxd13 gene from responding to forearm enhancers, thereby allowing proper limb patterning. To assess the nature and function of this CTCF-rich DNA region in embryos, we compared chromatin interaction profiles between proximal and distal limb bud cells isolated from mutant stocks where various parts of this boundary region were removed. The resulting progressive release in boundary effect triggered inter-TAD contacts, favored by the activity of the newly accessed enhancers. However, the boundary was highly resilient, and only a 400-kb deletion, including the whole-gene cluster, was eventually able to merge the neighboring TADs into a single structure. In this unified TAD, both proximal and distal limb enhancers nevertheless continued to work independently over a targeted transgenic reporter construct. We propose that the whole HoxD cluster is a dynamic TAD border and that the exact boundary position varies depending on both the transcriptional status and the developmental context.

Control of growth and gut maturation by HoxD genes and the associated lncRNA Haglr

During embryonic development, Hox genes participate in the building of a functional digestive system in metazoans, and genetic conditions involving these genes lead to important, sometimes lethal, growth retardation. Recently, this phenotype was obtained after deletion of Haglr, the Hoxd antisense growth-associated long noncoding RNA (lncRNA) located between Hoxd1 and Hoxd3 In this study, we have analyzed the function of Hoxd genes in delayed growth trajectories by looking at several nested targeted deficiencies of the mouse HoxD cluster. Mutant pups were severely stunted during the suckling period, but many recovered after weaning. After comparing seven distinct HoxD alleles, including CRISPR/Cas9 deletions involving Haglr, we identified Hoxd3 as the critical component for the gut to maintain milk-digestive competence. This essential function could be abrogated by the dominant-negative effect of HOXD10 as shown by a genetic rescue approach, thus further illustrating the importance of posterior prevalence in Hox gene function. A role for the lncRNA Haglr in the control of postnatal growth could not be corroborated.

No Sexual Dimorphism Detected in Digit Ratios of the Fire Salamander (Salamandra salamandra)

It has been proposed that digit ratio may be used as a biomarker of early developmental effects. Specifically, the second-to-fourth digit ratio (2D:4D) has been linked to the effects of sex hormones and their receptor genes, but other digit ratios have also been investigated. Across taxa, patterns of sexual dimorphism in digit ratios are ambiguous and a scarcity of studies in basal tetrapods makes it difficult to understand how ratios have evolved. Here, we focus on examining sex differences in digit ratios (2D:3D, 2D:4D, and 3D:4D) in a common amphibian, the fire salamander (Salamandra salamandra). We used graphic software to measure soft tissue digit length and digit bone length from X-rays. We found a nonsignificant tendency in males to have a lower 2D:3D than females; however, no sexual differences were detected in the other ratios. We discuss our results in the context of other studies of digit ratios, and how sex determination systems, as well as other factors, might impact patterns of sexual dimorphism, particularly in reptiles and in amphibians. Our findings suggest that caution is needed when using digit ratios as a potential indicator of prenatal hormonal effects in amphibians and highlight the need for more comparative studies to elucidate the evolutionary and genetic mechanisms implicated in sexually dimorphic patterns across taxonomic groups.

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.

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.

Long-range regulation by shared retinoic acid response elements modulates dynamic expression of posterior Hoxb genes in CNS development

Retinoic acid (RA) signaling plays an important role in determining the anterior boundary of Hox gene expression in the neural tube during embryogenesis. In particular, RA signaling is implicated in a rostral expansion of the neural expression domain of 5׳ Hoxb genes (Hoxb9-Hoxb5) in mice. However, underlying mechanisms for this gene regulation have remained elusive due to the lack of RA responsive element (RARE) in the 5׳ half of the HoxB cluster. To identify cis-regulatory elements required for the rostral expansion, we developed a recombineering technology to serially label multiple genes with different reporters in a single bacterial artificial chromosome (BAC) vector containing the mouse HoxB cluster. This allowed us to simultaneously monitor the expression of multiple genes. In contrast to plasmid-based reporters, transgenic BAC reporters faithfully recapitulated endogenous gene expression patterns of the Hoxb genes including the rostral expansion. Combined inactivation of two RAREs, DE-RARE and ENE-RARE, in the BAC completely abolished the rostral expansion of the 5׳ Hoxb genes. Knock-out of endogenous DE-RARE lead to significantly reduced expression of multiple Hoxb genes and attenuated Hox gene response to exogenous RA treatment in utero. Regulatory potential of DE-RARE was further demonstrated by its ability to anteriorize 5׳ Hoxa gene expression in the neural tube when inserted into a HoxA BAC reporter. Our data demonstrate that multiple RAREs cooperate to remotely regulate 5׳ Hoxb genes during CNS development, providing a new insight into the mechanisms for gene regulation within the Hox clusters.

The genetics of murine Hox loci: TAMERE, STRING, and PANTHERE to engineer chromosome variants

Following their duplications at the base of the vertebrate clade, Hox gene clusters underwent remarkable sub- and neo-functionalization events. Many of these evolutionary innovations can be associated with changes in the transcriptional regulation of their genes, where an intricate relationship between the structure of the gene cluster and the architecture of the surrounding genomic landscape is at play. Here, we report on a portfolio of in vivo genome engineering strategies in mice, which have been used to probe and decipher the genetic and molecular underpinnings of the complex regulatory mechanisms implemented at these loci.

Transgene- and locus-dependent imprinting reveals allele-specific chromosome conformations

When positioned into the integrin α-6 gene, an Hoxd9lacZ reporter transgene displayed parental imprinting in mouse embryos. While the expression from the paternal allele was comparable with patterns seen for the same transgene when present at the neighboring HoxD locus, almost no signal was scored at this integration site when the transgene was inherited from the mother, although the Itga6 locus itself is not imprinted. The transgene exhibited maternal allele-specific DNA hypermethylation acquired during oogenesis, and its expression silencing was reversible on passage through the male germ line. Histone modifications also corresponded to profiles described at known imprinted loci. Chromosome conformation analyses revealed distinct chromatin microarchitectures, with a more compact structure characterizing the maternally inherited repressed allele. Such genetic analyses of well-characterized transgene insertions associated with a de novo-induced parental imprint may help us understand the molecular determinants of imprinting.

HOX and TALE signatures specify human stromal stem cell populations from different sources

Human stromal stem cell populations reside in different tissues and anatomical sites, however a critical question related to their efficient use in regenerative medicine is whether they exhibit equivalent biological properties. Here, we compared cellular and molecular characteristics of stromal stem cells derived from the bone marrow, at different body sites (iliac crest, sternum, and vertebrae) and other tissues (dental pulp and colon). In particular, we investigated whether homeobox genes of the HOX and TALE subfamilies might provide suitable markers to identify distinct stromal cell populations, as HOX proteins control cell positional identity and, together with their co-factors TALE, are involved in orchestrating differentiation of adult tissues. Our results show that stromal populations from different sources, although immunophenotypically similar, display distinct HOX and TALE signatures, as well as different growth and differentiation abilities. Stromal stem cells from different tissues are characterized by specific HOX profiles, differing in the number and type of active genes, as well as in their level of expression. Conversely, bone marrow-derived cell populations can be essentially distinguished for the expression levels of specific HOX members, strongly suggesting that quantitative differences in HOX activity may be crucial. Taken together, our data indicate that the HOX and TALE profiles provide positional, embryological and hierarchical identity of human stromal stem cells. Furthermore, our data suggest that cell populations derived from different body sites may not represent equivalent cell sources for cell-based therapeutical strategies for regeneration and repair of specific tissues.

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.

Vertebrate GAGA factor associated insulator elements demarcate homeotic genes in the HOX clusters

Background

Hox genes impart segment identity to body structures along the anterior-posterior axis and are crucial for the proper development of all organisms. Multiple regulatory elements, best defined in Drosophila melanogaster, ensure that Hox expression patterns follow the spatial and temporal colinearity reflected in their tight genomic organization. However, the precise mechanisms that regulate colinear patterns of Hox gene expression remain unclear, especially in higher vertebrates where it is not fully determined how the distinct activation domains of the tightly clustered Hox genes are defined independently of each other. Here, we report the identification of a large number of novel cis-elements at mammalian Hox clusters that can help in regulating their precise expression pattern.

Results

We have identified DNA elements at all four murine Hox clusters that show poor association with histone H3 in chromatin immunoprecipitation (ChIP)-chip tiling arrays. The majority of these elements lie in the intergenic regions segregating adjacent Hox genes; we demonstrate that they possess efficient enhancer-blocking activity in mammalian cells. Further, we find that these histone-free intergenic regions bear GA repeat motifs and associate with the vertebrate homolog of the GAGA binding boundary factor. This suggests that they can act as GAGA factor-dependent chromatin boundaries that create independent domains, insulating each Hox gene from the influence of neighboring regulatory elements.

Conclusions

Our results reveal a large number of potential regulatory elements throughout the murine Hox clusters. We further demarcate the precise location of several novel cis-elements bearing chromatin boundary activity that appear to segregate successive Hox genes. This reflects a pattern reminiscent of the organization of homeotic genes in Drosophila, where such regulatory elements have been characterized. Our findings thus provide new insights into the regulatory processes and evolutionarily conserved epigenetic mechanisms that control homeotic gene expression.

Impact of copy number variations (CNVs) on long-range gene regulation at the HoxD locus

Copy number variations are genomic structural variants that are frequently associated with human diseases. Among these copy number variations, duplications of DNA segments are often assumed to lead to dosage effects by increasing the copy number of either genes or their regulatory elements. We produced a series of large targeted duplications within a conserved gene desert upstream of the murine HoxD locus. This DNA region, syntenic to human 2q31-32, contains a range of regulatory elements required for Hoxd gene transcription, and it is often disrupted and/or reorganized in human genetic conditions collectively known as the 2q31 syndrome. Unexpectedly, one such duplication led to a transcriptional down-regulation in developing digits by impairing physical interactions between the target genes and their upstream regulatory elements, thus phenocopying the effect obtained when these enhancer sequences are deleted. These results illustrate the detrimental consequences of interrupting highly conserved regulatory landscapes and reveal a mechanism where genomic duplications lead to partial loss of function of nearby located genes.

A function for all posterior Hoxd genes during digit development?

BACKGROUND

Four posterior Hoxd genes, from Hoxd13 to Hoxd10, are collectively regulated during the development of tetrapod digits. Besides the well-documented role of Hoxd13, the function of the neighboring genes has been difficult to evaluate due to the close genetic linkage and potential regulatory interferences. We used a combination of five small nested deletions in cis, involving from two to four consecutive genes of the Hoxd13 to Hoxd9 loci, in mice, to evaluate their combined functional importance.

RESULTS

We show that deletions leading to a gain of function of Hoxd13, via regulatory re-allocation, generate abnormal phenotypes, in agreement with the dominant negative role of this gene. We also show that Hoxd10, Hoxd11, and Hoxd12 all seem to play a genuine role in digit development, though less compelling than that of Hoxd13. In contrast, the nearby Hoxd9 contributed no measurable function in digits.

CONCLUSIONS

We conclude that a slight and transient deregulation of Hoxd13 expression can readily affect the relative lengths of limb segments and that all posterior Hoxd genes likely contribute to the final limb morphology. We discuss the difficulty to clearly assess the functional share of individual genes within such a gene family, where closely located neighbors, coding for homologous proteins, are regulated by a unique circuitry and all contribute to shape the distal parts of our appendages.

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.

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.

Comparative analyses of vertebrate posterior HoxD clusters reveal atypical cluster architecture in the caecilian Typhlonectes natans

BACKGROUND

The posterior genes of the HoxD cluster play a crucial role in the patterning of the tetrapod limb. This region is under the control of a global, long-range enhancer that is present in all vertebrates. Variation in limb types, as is the case in amphibians, can probably not only be attributed to variation in Hox genes, but is likely to be the product of differences in gene regulation. With a collection of vertebrate genome sequences available today, we used a comparative genomics approach to study the posterior HoxD cluster of amphibians. A frog and a caecilian were included in the study to compare coding sequences as well as to determine the gain and loss of putative regulatory sequences.

RESULTS

We sequenced the posterior end of the HoxD cluster of a caecilian and performed comparative analyses of this region using HoxD clusters of other vertebrates. We determined the presence of conserved non-coding sequences and traced gains and losses of these footprints during vertebrate evolution, with particular focus on amphibians. We found that the caecilian HoxD cluster is almost three times larger than its mammalian counterpart. This enlargement is accompanied with the loss of one gene and the accumulation of repeats in that area. A similar phenomenon was observed in the coelacanth, where a different gene was lost and expansion of the area where the gene was lost has occurred. At least one phylogenetic footprint present in all vertebrates was lost in amphibians. This conserved region is a known regulatory element and functions as a boundary element in neural tissue to prevent expression of Hoxd genes.

CONCLUSION

The posterior part of the HoxD cluster of Typhlonectes natans is among the largest known today. The loss of Hoxd-12 and the expansion of the intergenic region may exert an influence on the limb enhancer, by having to bypass a distance seven times that of regular HoxD clusters. Whether or not there is a correlation with the loss of limbs remains to be investigated. These results, together with data on other vertebrates show that the tetrapod Hox clusters are more variable than previously thought.

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.

Uncoupling time and space in the collinear regulation of Hox genes

During development of the vertebrate body axis, Hox genes are transcribed sequentially, in both time and space, following their relative positions within their genomic clusters. Analyses of animal genomes support the idea that Hox gene clustering is essential for coordinating the various times of gene activations. However, the eventual collinear ordering of the gene specific transcript domains in space does not always require genomic clustering. We analyzed these complex regulatory relationships by using mutant alleles at the mouse HoxD locus, including one that splits the cluster into two pieces. We show that both positive and negative regulatory influences, located on either side of the cluster, control an early phase of collinear expression in the trunk. Interestingly, this early phase does not systematically impact upon the subsequent expression patterns along the main body axis, indicating that the mechanism underlying temporal collinearity is distinct from those acting during the second phase. We discuss the potential functions and evolutionary origins of these mechanisms, as well as their relationship with similar processes at work during limb development.

Hox genes encode proteins that control embryonic development along the head-to-tail axis. These genes are clustered in one site on the chromosome and their respective positions within the cluster determine their time and place of activation. Here, by using a large set of targeted mutations disturbing the integrity of the gene cluster, we show that the spatial organization of expression domains does not directly depend upon the timing of activation as was previously suggested. This uncoupling between space and time in the regulation of these Hox genes coincides with the existence of two major phases of regulation. The first is time-dependent and involves global regulatory influences, located outside the gene cluster, whereas the second relies upon more local regulatory elements, likely interspersed between the genes, inside the cluster. These results provide the bases for future analyses of collinear mechanisms and indicate that different types of collinearities are not necessarily related, neither in function, nor in their evolutionary histories.

Establishment of Hox vertebral identities in the embryonic spine precursors

The vertebrate spine exhibits two striking characteristics. The first one is the periodic arrangement of its elements-the vertebrae-along the anteroposterior axis. This segmented organization is the result of somitogenesis, which takes place during organogenesis. The segmentation machinery involves a molecular oscillator-the segmentation clock-which delivers a periodic signal controlling somite production. During embryonic axis elongation, this signal is displaced posteriorly by a system of traveling signaling gradients-the wavefront-which depends on the Wnt, FGF, and retinoic acid pathways. The other characteristic feature of the spine is the subdivision of groups of vertebrae into anatomical domains, such as the cervical, thoracic, lumbar, sacral, and caudal regions. This axial regionalization is controlled by a set of transcription factors called Hox genes. Hox genes exhibit nested expression domains in the somites which reflect their linear arrangement along the chromosomes-a property termed colinearity. The colinear disposition of Hox genes expression domains provides a blueprint for the regionalization of the future vertebral territories of the spine. In amniotes, Hox genes are activated in the somite precursors of the epiblast in a temporal colinear sequence and they were proposed to control their progressive ingression into the nascent paraxial mesoderm. Consequently, the positioning of the expression domains of Hox genes along the anteroposterior axis is largely controlled by the timing of Hox activation during gastrulation. Positioning of the somitic Hox domains is subsequently refined through a crosstalk with the segmentation machinery in the presomitic mesoderm. In this review, we focus on our current understanding of the embryonic mechanisms that establish vertebral identities during vertebrate development.

Distinct roles and regulations for HoxD genes in metanephric kidney development

Hox genes encode homeodomain-containing proteins that control embryonic development in multiple contexts. Up to 30 Hox genes, distributed among all four clusters, are expressed during mammalian kidney morphogenesis, but functional redundancy between them has made a detailed functional account difficult to achieve. We have investigated the role of the HoxD cluster through comparative molecular embryological analysis of a set of mouse strains carrying targeted genomic rearrangements such as deletions, duplications, and inversions. This analysis allowed us to uncover and genetically dissect the complex role of the HoxD cluster. Regulation of metanephric mesenchyme-ureteric bud interactions and maintenance of structural integrity of tubular epithelia are differentially controlled by some Hoxd genes during renal development, consistent with their specific expression profiles. We also provide evidence for a kidney-specific form of colinearity that underlies the differential expression of two distinct sets of genes located on both sides and overlapping at the Hoxd9 locus. These insights further our knowledge of the genetic control of kidney morphogenesis and may contribute to understanding certain congenital kidney malformations, including polycystic kidney disease and renal hypoplasia.

Hox genes encode proteins that control embryonic development along the head-to-tail axis and in multiple organs. Here, we show that several members of this gene family are necessary for the normal development of the mammalian kidneys. These genes are clustered in one site on the chromosome and their respective positions within the group determine which component of the kidneys they will contribute to. Using a large collection of engineered mutations in this system, we show that these genes are required both for the growth of the kidneys and for their proper organization, such that mutations in some genes reduce the size of the organs, whereas mutations in others induce polycystic kidneys. Our set of genetic rearrangements also allowed us to localize the position of regulatory sequences, which control the expression of these genes during kidney development.

Ectopic nuclear reorganisation driven by a Hoxb1 transgene transposed into Hoxd

The extent to which the nuclear organisation of a gene impacts on its ability to be expressed, or whether nuclear organisation merely reflects gene expression states, remains an important but unresolved issue. A model system that has been instrumental in investigating this question utilises the murine Hox gene clusters encoding homeobox-containing proteins. Nuclear reorganisation and chromatin decondensation, initiated towards the 3' end of the clusters, accompanies activation of Hox genes in both differentiation and development, and might be linked to mechanisms underlying colinearity. To investigate this, and to delineate the cis-acting elements involved, here we analyse the nuclear behaviour of a 3' Hoxb1 transgene transposed to the 5' end of the Hoxd cluster. We show that this transgene contains the cis-acting elements sufficient to initiate ectopic local nuclear reorganisation and chromatin decondensation and to break Hoxd colinearity in the primitive streak region of the early embryo. Significantly, in rhombomere 4, the transgene is able to induce attenuated nuclear reorganisation and decondensation of Hoxd even though there is no detectable expression of the transgene at this site. This shows that reorganisation of chromosome territories and chromatin decondensation can be uncoupled from transcription itself and suggests that they can therefore operate upstream of gene expression.

Long-range gene control and genetic disease

The past two decades have seen great progress in the elucidation of the genetic basis of human genetic disease. Many clinical phenotypes have been linked with mutations or deletions in specific causative genes. However, it is often less recognized that in addition to the integrity of the protein-coding sequences, human health critically also depends on the spatially, temporally, and quantitatively correct expression of those genes. Genetic disease can therefore equally be caused by disruption of the regulatory mechanisms that ensure proper gene expression. The term "position effect" is used in those situations where the expression level of a gene is deleteriously affected by an alteration in its chromosomal environment, while maintaining an intact transcription unit. Here, we review recent advances in our understanding of the possible mechanisms of a number of "position effect" disease cases and discuss the findings with respect to current models for genome organization and long-range control of gene expression.

Global control regions and regulatory landscapes in vertebrate development and evolution

During the course of evolution, many genes that control the development of metazoan body plans were co-opted to exert novel functions, along with the emergence or modification of structures. Gene amplification and/or changes in the cis-regulatory modules responsible for the transcriptional activity of these genes have certainly contributed in a major way to evolution of gene functions. In some cases, these processes led to the formation of groups of adjacent genes that appear to be controlled by both global and shared mechanisms.

Hoxd and Gli3 interactions modulate digit number in the amniote limb

During limb development, Sonic hedgehog (SHH) and HOX proteins are considered among the most important factors regulating digit number and identity. SHH signaling prevents the processing of GLI3 into a short form that functions as a strong transcriptional repressor. Gli3 mutant limbs are characterized by a severe polydactyly and associated ectopic anterior expression of 5'Hoxd genes. To genetically determine the involvement of 5'Hoxd genes in the polydactyly of Gli3 mutants, we have generated a compound mutant that simultaneously removes the three most 5'-located Hoxd genes and Gli3. Remarkably, the limbs that form in the absence of all four of these genes show the most severe polydactyly so far reported in the mouse. The analysis of gene expression performed in compound mutants allows us to propose that the increase in the number of digits is mediated by the gain in function of Hoxd10 and Hoxd9. Our results also support the notion that an adequate balance between positive and negative effects of different Hoxd genes is required for pentadactyly.

Mining of putative cis-acting elements for chromatin mediated regulation of Hox genes in mammals by in-silico analysis

The remarkable conservation in developmental strategies across phyla is well reflected in the conservation of the homeotic gene complexes responsible for establishing the body plan in embryonic development. On the other hand, changes in the strategy of transcription regulation are believed to form one of the major factors in the evolution of developmental mechanisms and phenotypic evolution of species. Apart from transcription regulation by gene specific transcription factors, the role of regulators mediating modifications of chromatin proteins, especially of HOX gene clusters in Drosophila is well documented. By comparative genomics we have identified novel motifs conserved in mouse, chimpanzee and human in the noncoding upstream/intronic sequences of Hox genes by in silico analysis. These motifs lack the binding sites for known transcription factors and are significantly over represented in the target genes of one of the core components of Polycomb Repressive Complex namely Supressor of zeste 12 (SUZ12) in human embryonic cells reported by Lee et al. [2006a. Cell 125:301-313]. Therefore, we predict that they could be the sites of interaction of chromatin modifying complexes for epigenetic regulation.

Interactions between HOXD and Gli3 genes control the limb apical ectodermal ridge via Fgf10

The development of the vertebrate limb is dependent upon two signaling centers, the apical ectodermal ridge (AER), which provides the underlying mesenchyme with essential growth factors, and the zone of polarizing activity (ZPA), the source of the Sonic hedgehog (SHH) product. Recent work involving gain and loss of function of Hox genes has emphasized their impact both on AER maintenance and Shh transcriptional activation. Here, we describe antagonistic interactions between posterior Hoxd genes and Gli3, suggesting that the latter product protects the AER from the deleterious effect of the formers, and we present evidence that Fgf10 is the mediator of HOX-dependent AER expansion. Furthermore, the striking similarity between some of the hereby observed Hox/Gli3-dependent morphogenetic defects and those displayed by fetuses with severely altered retinoic acid metabolism suggests a tight connection between these various pathways. The nature of these potential interactions is discussed in the context of proximal-distal growth and patterning.

Hox genes in time and space during vertebrate body formation

Vertebrae display distinct morphological features at different levels of the body axis. Links between collinear Hox gene activation and the progressive mode of body axis elongation have provided a fascinating blueprint of the mechanisms for establishing these morphological identities. In this review, we first discuss the regulation and possible role of collinear Hox gene activation during body formation and then highlight the direct role of Hox genes in controlling cellular movements during gastrulation, therefore contributing to body formation. Additional related research aspects, such as imaging of chromatin regulation, roles of micro RNAs and evolutional findings are also discussed.

Transgenic analysis of Hoxd gene regulation during digit development

In tetrapods, posterior Hoxd genes (from groups 10 to 13) are necessary to properly pattern the developing autopods, including the number and identities of digits. Their coordinated expression is achieved by sharing a global control region (GCR), which was isolated and localized 200 kb 5' (centromeric) of the gene cluster. However, in transgenic assays, the GCR was unable to fully recapitulate all aspects of the endogenous Hoxd expression patterns during distal limb development. In this paper, we further analyze the regulatory potential of this locus and report the characterization of Prox, a second enhancer element that contributes to the transcriptional activity of posterior Hoxd genes in developing distal limb buds. We show that the GCR and Prox elements complement each other and work in combination to correctly establish the late phase of Hoxd genes expression. Based on DNA sequence conservation and transgenic assays, we discuss the functions of these regulatory regions as well as a potential evolutionary scheme accounting for their emergence along with the evolution of tetrapod limbs.

Nuclear reorganisation and chromatin decondensation are conserved, but distinct, mechanisms linked to Hox gene activation

The relocalisation of some genes to positions outside chromosome territories, and the visible decondensation or unfolding of interphase chromatin, are two striking facets of nuclear reorganisation linked to gene activation that have been assumed to be related to each other. Here, in a study of nuclear reorganisation around the Hoxd cluster, we suggest that this may not be the case. Despite its very different genomic environment from Hoxb, Hoxd also loops out from its chromosome territory, and unfolds, upon activation in differentiating embryonic stem (ES) cells and in the tailbud of the embryo. However, looping out and decondensation are not simply two different manifestations of the same underlying change in chromatin structure. We show that, in the limb bud of the embryonic day 9.5 embryo,where Hoxd is also activated, there is visible decondensation of chromatin but no detectable movement of the region out from the chromosome territory. During ES cell differentiation, decondensed alleles can also be found inside of chromosome territories, and loci that have looped out of the territories can appear to still be condensed. We conclude that evolutionarily conserved chromosome remodelling mechanisms, predating the duplication of mammalian Hox loci, underlie Hox regulation along the rostrocaudal embryonic axis. However, we suggest that separate modes of regulation can modify Hoxd chromatin in different ways in different developmental contexts.

Duplicated Abd-B class genes in medaka hoxAa and hoxAb clusters exhibit differential expression patterns in pectoral fin buds

Hox genes form clusters. Invertebrates and Amphioxus have only one hox cluster, but in vertebrates, they are multiple, i.e., four in the basal teleost fish Polyodon and tetrapods (HoxA, B, C, D), but seven or eight in common teleosts. We earlier completely sequenced the entire hox gene loci in medaka fish, showing a total of 46 hox genes to be encoded in seven clusters (hoxAa, Ab, Ba, Bb, Ca, Da, Db). Among them, hoxAa, hoxAb and hoxDa clusters are presumed to be important for fin-to-limb evolution because of their key role in forelimb and pectoral fin development. In the present study, we compared genome organization and nucleotide sequences of the hoxAa and hoxAb clusters to these of tetrapod HoxA clusters, and found greater similarity in hoxAa case. We then analyzed expression of Abd-B family genes in the clusters. In the trunk, those from the hoxAa cluster, i.e., hoxA9a, hoxA10a, hoxA11a and hoxA13a, were expressed in a manner keeping the colinearity rule of the hox expression as those of tetrapods, while those from the hoxAb cluster, i.e., hoxA9b, hoxA10b, hoxA11b and hoxA13b, were not. In the pectoral fins, the hoxAa cluster was expressed in split domains and did not obey the rule. By contrast, those from the hoxAb and hoxDa clusters were expressed in a manner keeping the rule, i.e., an ancestral pattern similar to those of tetrapods. It is plausible that this differential expression of the two clusters is caused by changes occurred in global control regions after cluster duplications.

Regulatory constraints in the evolution of the tetrapod limb anterior-posterior polarity

The anterior to posterior (A-P) polarity of the tetrapod limb is determined by the confined expression of Sonic hedgehog (Shh) at the posterior margin of developing early limb buds, under the control of HOX proteins encoded by gene members of both the HoxA and HoxD clusters. Here, we use a set of partial deletions to show that only the last four Hox paralogy groups can elicit this response: that is, precisely those genes whose expression is excluded from most anterior limb bud cells owing to their collinear transcriptional activation. We propose that the limb A-P polarity is produced as a collateral effect of Hox gene collinearity, a process highly constrained by its crucial importance during trunk development. In this view, the co-option of the trunk collinear mechanism, along with the emergence of limbs, imposed an A-P polarity to these structures as the most parsimonious solution. This in turn further contributed to stabilize the architecture and operational mode of this genetic system.

Anteroposterior patterning in the limb and digit specification: contribution of mouse genetics

The limb has been a privileged object of investigation and reflection for scientists over the past two centuries and continues to provide a heuristic framework to analyze vertebrate development. Recently, accumulation of new data has significantly changed our view on the mechanisms of limb patterning, in particular along the anterior-posterior axis. These data have led us to revisit the mode of action of the zone of polarizing activity. They shed light on the molecular and cellular mechanisms of patterning linked to the Shh-Gli3 signaling pathway and give insights into the mechanism of activation of these cardinal factors, as well as the consequences of their activity. These new data are in good part the result of systematic Application of tools used in contemporary mouse molecular genetics. These have extended the power of mouse genetics by introducing mutational strategies that allow fine-tuned modulation of gene expression, interchromosomal deletions and duplication. They have even made the mouse embryo amenable to cell lineage analysis that used to be the realm of chick embryos. In this review, we focus on the data acquired over the last five years from the analysis of mouse limb development and discuss new perspectives opened by these results.

A group 13 homeodomain is neither necessary nor sufficient for posterior prevalence in the mouse limb

Posterior prevalence is the general property attributed to HOX proteins describing the dominant effect of more posterior HOX proteins over the function of anterior orthologs in common areas of expression. To explore the HOX group 13 protein domains required for this property, we used the mouse Prx-1 promoter to drive transgenic expression of Hox constructs throughout the entire limb bud during development. This system allowed us to conclusively demonstrate a hierarchy of Hox function in developing limbs. Furthermore, by substituting the HOXD11 or HOXA9 homeodomain for that of HOXD13, we show that a HOXD13 homeodomain is not necessary for posterior prevalence. Proximal expression of these chimeric proteins unexpectedly caused defects consistent with wild-type HOXD13 mediated posterior prevalence. Moreover, group 13 non-homeodomain residues appear to confer the property as proximal expression of HOXA9 containing the HOXD13 homeodomain did not result in limb reductions characteristic of HOXD13. These data are most compatible with models of posterior prevalence based on protein-protein interactions and support examination of the N-terminal non-homeodomain regions of Hox group 13 proteins as necessary agents for posterior prevalence.

Differentiation of trophoblast giant cells and their metabolic functions are dependent on peroxisome proliferator-activated receptor beta/delta

Mutation of the nuclear receptor peroxisome proliferator-activated receptor beta/delta (PPARbeta/delta) severely affects placenta development, leading to embryonic death at embryonic day 9.5 (E9.5) to E10.5 of most, but not all, PPARbeta/delta-null mutant embryos. While very little is known at present about the pathway governed by PPARbeta/delta in the developing placenta, this paper demonstrates that the main alteration of the placenta of PPARbeta/delta-null embryos is found in the giant cell layer. PPARbeta/delta activity is in fact essential for the differentiation of the Rcho-1 cells in giant cells, as shown by the severe inhibition of differentiation once PPARbeta/delta is silenced. Conversely, exposure of Rcho-1 cells to a PPARbeta/delta agonist triggers a massive differentiation via increased expression of 3-phosphoinositide-dependent kinase 1 and integrin-linked kinase and subsequent phosphorylation of Akt. The links between PPARbeta/delta activity in giant cells and its role on Akt activity are further strengthened by the remarkable pattern of phospho-Akt expression in vivo at E9.5, specifically in the nucleus of the giant cells. In addition to this phosphatidylinositol 3-kinase/Akt main pathway, PPARbeta/delta also induced giant cell differentiation via increased expression of I-mfa, an inhibitor of Mash-2 activity. Finally, giant cell differentiation at E9.5 is accompanied by a PPARbeta/delta-dependent accumulation of lipid droplets and an increased expression of the adipose differentiation-related protein (also called adipophilin), which may participate to lipid metabolism and/or steroidogenesis. Altogether, this important role of PPARbeta/delta in placenta development and giant cell differentiation should be considered when contemplating the potency of PPARbeta/delta agonist as therapeutic agents of broad application.

Genomic annotation and transcriptome analysis of the zebrafish (Danio rerio) hox complex with description of a novel member, hoxb13a

The zebrafish (Danio rerio) is an important model in evolutionary developmental biology, and its study is being revolutionized by the zebrafish genome project. Sequencing is at an advanced stage, but annotation is largely the result of in silico analyses. We have performed genomic annotation, comparative genomics, and transcriptional analysis using microarrays of the hox homeobox-containing transcription factors. These genes have important roles in specifying the body plan. Candidate sequences were located in version Z v 4 of the Ensembl genome database by TBLASTN searching with Danio and other vertebrate published Hox protein sequences. Homologies were confirmed by alignment with reference sequences, and by the relative position of genes along each cluster. RT-PCR using adult Tübingen cDNA was used to confirm annotations, to check the genomic sequence and to confirm expression in vivo. Our RT-PCR and microarray data show that all 49 hox genes are expressed in adult zebrafish. Significant expression for all known hox genes could be detected in our microarray analysis. We also find significant expression of hox 8 paralogs and hox b 7 a in the anti-sense direction. A novel gene, D. rerio hox b 13 a, was identified, and a preliminary characterization by in situ hybridization showed expression at 24 hpf at the tip of the developing tail. We are currently characterizing this gene at the functional level. We argue that the oligo design for microarrays can be greatly enhanced by the availability of genomic sequences.

Globalisation reaches gene regulation: the case for vertebrate limb development

Analysis of key regulators of vertebrate limb development has revealed that the cis-regulatory regions controlling their expression are often located several hundred kilobases upstream of the transcription units. These far up- or down-stream cis-regulatory regions tend to reside within rather large, functionally and structurally unrelated genes. Molecular analysis is beginning to reveal the complexity of these large genomic landscapes, which control the co-expression of clusters of diverse genes by this novel type of long-range and globally acting cis-regulatory region. An increasing number of spontaneous mutations in vertebrates, including humans, are being discovered inactivating or altering such global control regions. Thereby, the functions of a seemingly distant but essential gene are disrupted rather than the closest.

Comparative analysis of genes downstream of the Hoxd cluster in developing digits and external genitalia

Mammalian Hox genes encode transcription factors that are crucial for proper morphogenesis along the various body axes. Despite their extensive structural and functional characterization, the nature of their target genes remains elusive. We have addressed this question by using DNA microarrays to screen for genes whose expression in developing distal forelimbs and genital eminences was significantly modified in the absence of the full Hoxd gene complement. This comparative approach not only identified specific candidate genes, but also allowed the examination of whether a similar Hox expression pattern in distinct tissues leads to the modulation of the same or different downstream genes. We report here a set of potential target genes, most of which were not previously known to play a role in the early stages of either limb or genital bud development. Interestingly, we find that the majority of these candidate genes are differentially expressed in both structures,although often at different times. This supports the idea that both appendices involve similar genetic controls, both upstream and downstream of the Hox gene family. These results highlight the surprising mechanistic relationship between these rather different body parts and suggest a common developmental strategy to build up the most distal appendicular structures of the body, i.e. the digits and the penis/clitoris.

Long-range control of gene expression: emerging mechanisms and disruption in disease

Transcriptional control is a major mechanism for regulating gene expression. The complex machinery required to effect this control is still emerging from functional and evolutionary analysis of genomic architecture. In addition to the promoter, many other regulatory elements are required for spatiotemporally and quantitatively correct gene expression. Enhancer and repressor elements may reside in introns or up- and downstream of the transcription unit. For some genes with highly complex expression patterns--often those that function as key developmental control genes--the cis-regulatory domain can extend long distances outside the transcription unit. Some of the earliest hints of this came from disease-associated chromosomal breaks positioned well outside the relevant gene. With the availability of wide-ranging genome sequence comparisons, strong conservation of many noncoding regions became obvious. Functional studies have shown many of these conserved sites to be transcriptional regulatory elements that sometimes reside inside unrelated neighboring genes. Such sequence-conserved elements generally harbor sites for tissue-specific DNA-binding proteins. Developmentally variable chromatin conformation can control protein access to these sites and can regulate transcription. Disruption of these finely tuned mechanisms can cause disease. Some regulatory element mutations will be associated with phenotypes distinct from any identified for coding-region mutations.