A regulatory 'landscape effect' over the HoxD cluster

Adaptive tail-length evolution in deer mice is associated with differential Hoxd13 expression in early development

Variation in the size and number of axial segments underlies much of the diversity in animal body plans. Here we investigate the evolutionary, genetic and developmental mechanisms driving tail-length differences between forest and prairie ecotypes of deer mice (Peromyscus maniculatus). We first show that long-tailed forest mice perform better in an arboreal locomotion assay, consistent with tails being important for balance during climbing. We then identify six genomic regions that contribute to differences in tail length, three of which associate with caudal vertebra length and the other three with vertebra number. For all six loci, the forest allele increases tail length, indicative of the cumulative effect of natural selection. Two of the genomic regions associated with variation in vertebra number contain Hox gene clusters. Of those, we find an allele-specific decrease in Hoxd13 expression in the embryonic tail bud of long-tailed forest mice, consistent with its role in axial elongation. Additionally, we find that forest embryos have more presomitic mesoderm than prairie embryos and that this correlates with an increase in the number of neuromesodermal progenitors, which are modulated by Hox13 paralogues. Together, these results suggest a role for Hoxd13 in the development of natural variation in adaptive morphology on a microevolutionary timescale.

EVOLUTIONARY CO-OPTION OF AN ANCESTRAL CLOACAL REGULATORY LANDSCAPE DURING THE EMERGENCE OF DIGITS AND GENITALS

The transition from fins to limbs has been a rich source of discussion for more than a century. One open and important issue is understanding how the mechanisms that pattern digits arose during vertebrate evolution. In this context, the analysis of Hox gene expression and functions to infer evolutionary scenarios has been a productive approach to explain the changes in organ formation, particularly in limbs. In tetrapods, the transcription of Hoxd genes in developing digits depends on a well-characterized set of enhancers forming a large regulatory landscape1,2. This control system has a syntenic counterpart in zebrafish, even though they lack bona fide digits, suggestive of deep homology3 between distal fin and limb developmental mechanisms. We tested the global function of this landscape to assess ancestry and source of limb and fin variation. In contrast to results in mice, we show here that the deletion of the homologous control region in zebrafish has a limited effect on the transcription of hoxd genes during fin development. However, it fully abrogates hoxd expression within the developing cloaca, an ancestral structure related to the mammalian urogenital sinus. We show that similar to the limb, Hoxd gene function in the urogenital sinus of the mouse also depends on enhancers located in this same genomic domain. Thus, we conclude that the current regulation underlying Hoxd gene expression in distal limbs was co-opted in tetrapods from a preexisting cloacal program. The orthologous chromatin domain in fishes may illustrate a rudimentary or partial step in this evolutionary co-option.

Echiuran Hox genes provide new insights into the correspondence between Hox subcluster organization and collinearity pattern

In many bilaterians, Hox genes are generally clustered along the chromosomes and expressed in spatial and temporal order. In vertebrates, the expression of Hox genes follows a whole-cluster spatio-temporal collinearity (WSTC) pattern, whereas in some invertebrates the expression of Hox genes exhibits a subcluster-level spatio-temporal collinearity pattern. In bilaterians, the diversity of collinearity patterns and the cause of collinearity differences in Hox gene expression remain poorly understood. Here, we investigate genomic organization and expression pattern of Hox genes in the echiuran worm Urechis unicinctus (Annelida, Echiura). Urechis unicinctus has a split cluster with four subclusters divided by non-Hox genes: first subcluster (Hox1 and Hox2), second subcluster (Hox3), third subcluster (Hox4, Hox5, Lox5, Antp and Lox4), fourth subcluster (Lox2 and Post2). The expression of U. unicinctus Hox genes shows a subcluster-based whole-cluster spatio-temporal collinearity (S-WSTC) pattern: the anterior-most genes in each subcluster are activated in a spatially and temporally colinear manner (reminiscent of WSTC), with the subsequent genes in each subcluster then being very similar to their respective anterior-most subcluster gene. Combining genomic organization and expression profiles of Hox genes in different invertebrate lineages, we propose that the spatio-temporal collinearity of invertebrate Hox is subcluster-based.

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.

Dbx2 regulation in limbs suggests interTAD sharing of enhancers

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Geminin is required for Hox gene regulation to pattern the developing limb

Development of the complex structure of the vertebrate limb requires carefully orchestrated interactions between multiple regulatory pathways and proteins. Among these, precise regulation of 5' Hox transcription factor expression is essential for proper limb bud patterning and elaboration of distinct limb skeletal elements. Here, we identified Geminin (Gmnn) as a novel regulator of this process. A conditional model of Gmnn deficiency resulted in loss or severe reduction of forelimb skeletal elements, while both the forelimb autopod and hindlimb were unaffected. 5' Hox gene expression expanded into more proximal and anterior regions of the embryonic forelimb buds in this Gmnn-deficient model. A second conditional model of Gmnn deficiency instead caused a similar but less severe reduction of hindlimb skeletal elements and hindlimb polydactyly, while not affecting the forelimb. An ectopic posterior SHH signaling center was evident in the anterior hindlimb bud of Gmnn-deficient embryos in this model. This center ectopically expressed Hoxd13, the HOXD13 target Shh, and the SHH target Ptch1, while these mutant hindlimb buds also had reduced levels of the cleaved, repressor form of GLI3, a SHH pathway antagonist. Together, this work delineates a new role for Gmnn in modulating Hox expression to pattern the vertebrate limb.

A complex regulatory landscape involved in the development of mammalian external genitals

Developmental genes are often controlled by large regulatory landscapes matching topologically associating domains (TADs). In various contexts, the associated chromatin backbone is modified by specific enhancer-enhancer and enhancer-promoter interactions. We used a TAD flanking the mouse HoxD cluster to study how these regulatory architectures are formed and deconstructed once their function achieved. We describe this TAD as a functional unit, with several regulatory sequences acting together to elicit a transcriptional response. With one exception, deletion of these sequences didn't modify the transcriptional outcome, a result at odds with a conventional view of enhancer function. The deletion and inversion of a CTCF site located near these regulatory sequences did not affect transcription of the target gene. Slight modifications were nevertheless observed, in agreement with the loop extrusion model. We discuss these unexpected results considering both conventional and alternative explanations relying on the accumulation of poorly specific factors within the TAD backbone.

The evo-devo of molluscs: Insights from a genomic perspective

Molluscs represent one of ancient and evolutionarily most successful groups of marine invertebrates, with a tremendous diversity of morphology, behavior, and lifestyle. Molluscs are excellent subjects for evo-devo studies; however, understanding of the evo-devo of molluscs has been largely hampered by incomplete fossil records and limited molecular data. Recent advancement of genomics and other technologies has greatly fueled the molluscan "evo-devo" field, and decoding of several molluscan genomes provides unprecedented insights into molluscan biology and evolution. Here, we review the recent progress of molluscan genome sequencing as well as novel insights gained from their genomes, by emphasizing how molluscan genomics enhances our understanding of the evo-devo of molluscs.

Role for first zinc finger of WT1 in DNA sequence specificity: Denys-Drash syndrome-associated WT1 mutant in ZF1 enhances affinity for a subset of WT1 binding sites

Wilms tumor protein (WT1) is a Cys2-His2 zinc-finger transcription factor vital for embryonic development of the genitourinary system. The protein contains a C-terminal DNA binding domain with four tandem zinc-fingers (ZF1-4). An alternative splicing of Wt1 can add three additional amino acids-lysine (K), threonine (T) and serine (S)-between ZF3 and ZF4. In the -KTS isoform, ZF2-4 determine the sequence-specificity of DNA binding, whereas the function of ZF1 remains elusive. Three X-ray structures are described here for wild-type -KTS isoform ZF1-4 in complex with its cognate DNA sequence. We observed four unique ZF1 conformations. First, like ZF2-4, ZF1 can be positioned continuously in the DNA major groove forming a 'near-cognate' complex. Second, while ZF2-4 make base-specific interactions with one DNA molecule, ZF1 can interact with a second DNA molecule (or, presumably, two regions of the same DNA molecule). Third, ZF1 can intercalate at the joint of two tail-to-head DNA molecules. If such intercalation occurs on a continuous DNA molecule, it would kink the DNA at the ZF1 binding site. Fourth, two ZF1 units can dimerize. Furthermore, we examined a Denys-Drash syndrome-associated ZF1 mutation (methionine at position 342 is replaced by arginine). This mutation enhances WT1 affinity for a guanine base. X-ray crystallography of the mutant in complex with its preferred sequence revealed the interactions responsible for this affinity change. These results provide insight into the mechanisms of action of WT1, and clarify the fact that ZF1 plays a role in determining sequence specificity of this critical transcription factor.

Impact of genome architecture on the functional activation and repression of Hox regulatory landscapes

BACKGROUND

The spatial organization of the mammalian genome relies upon the formation of chromatin domains of various scales. At the level of gene regulation in cis, collections of enhancer sequences define large regulatory landscapes that usually match with the presence of topologically associating domains (TADs). These domains often contain ranges of enhancers displaying similar or related tissue specificity, suggesting that in some cases, such domains may act as coherent regulatory units, with a global on or off state. By using the HoxD gene cluster, which specifies the topology of the developing limbs via highly orchestrated regulation of gene expression, as a paradigm, we investigated how the arrangement of regulatory domains determines their activity and function.

RESULTS

Proximal and distal cells in the developing limb express different levels of Hoxd genes, regulated by flanking 3' and 5' TADs, respectively. We characterized the effect of large genomic rearrangements affecting these two TADs, including their fusion into a single chromatin domain. We show that, within a single hybrid TAD, the activation of both proximal and distal limb enhancers globally occurred as when both TADs are intact. However, the activity of the 3' TAD in distal cells is generally increased in the fused TAD, when compared to wild type where it is silenced. Also, target gene activity in distal cells depends on whether or not these genes had previously responded to proximal enhancers, which determines the presence or absence of H3K27me3 marks. We also show that the polycomb repressive complex 2 is mainly recruited at the Hox gene cluster and can extend its coverage to far-cis regulatory sequences as long as confined to the neighboring TAD structure.

CONCLUSIONS

We conclude that antagonistic limb proximal and distal enhancers can exert their specific effects when positioned into the same TAD and in the absence of their genuine target genes. We also conclude that removing these target genes reduced the coverage of a regulatory landscape by chromatin marks associated with silencing, which correlates with its prolonged activity in time.

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.

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

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

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

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

Sea cucumber genome provides insights into saponin biosynthesis and aestivation regulation

Echinoderms exhibit several fascinating evolutionary innovations that are rarely seen in the animal kingdom, but how these animals attained such features is not well understood. Here we report the sequencing and analysis of the genome and extensive transcriptomes of the sea cucumber Apostichopus japonicus, a species from a special echinoderm group with extraordinary potential for saponin synthesis, aestivation and organ regeneration. The sea cucumber does not possess a reorganized Hox cluster as previously assumed for all echinoderms, and the spatial expression of Hox7 and Hox11/13b potentially guides the embryo-to-larva axial transformation. Contrary to the typical production of lanosterol in animal cholesterol synthesis, the oxidosqualene cyclase of sea cucumber produces parkeol for saponin synthesis and has "plant-like" motifs suggestive of convergent evolution. The transcriptional factors Klf2 and Egr1 are identified as key regulators of aestivation, probably exerting their effects through a clock gene-controlled process. Intestinal hypometabolism during aestivation is driven by the DNA hypermethylation of various metabolic gene pathways, whereas the transcriptional network of intestine regeneration involves diverse signaling pathways, including Wnt, Hippo and FGF. Decoding the sea cucumber genome provides a new avenue for an in-depth understanding of the extraordinary features of sea cucumbers and other echinoderms.

Investigation into the role of the germline epigenome in the transmission of glucocorticoid-programmed effects across generations

Background

Early life exposure to adverse environments affects cardiovascular and metabolic systems in the offspring. These programmed effects are transmissible to a second generation through both male and female lines, suggesting germline transmission. We have previously shown that prenatal overexposure to the synthetic glucocorticoid dexamethasone (Dex) in rats reduces birth weight in the first generation (F1), a phenotype which is transmitted to a second generation (F2), particularly through the male line. We hypothesize that Dex exposure affects developing germ cells, resulting in transmissible alterations in DNA methylation, histone marks and/or small RNA in the male germline.

Results

We profile epigenetic marks in sperm from F1 Sprague Dawley rats expressing a germ cell-specific GFP transgene following Dex or vehicle treatment of the mothers, using methylated DNA immunoprecipitation sequencing, small RNA sequencing and chromatin immunoprecipitation sequencing for H3K4me3, H3K4me1, H3K27me3 and H3K9me3. Although effects on birth weight are transmitted to the F2 generation through the male line, no differences in DNA methylation, histone modifications or small RNA were detected between germ cells and sperm from Dex-exposed animals and controls.

Conclusions

Although the phenotype is transmitted to a second generation, we are unable to detect specific changes in DNA methylation, common histone modifications or small RNA profiles in sperm. Dex exposure is associated with more variable 5mC levels, particularly at non-promoter loci. Although this could be one mechanism contributing to the observed phenotype, other germline epigenetic modifications or non-epigenetic mechanisms may be responsible for the transmission of programmed effects across generations in this model.

Electronic supplementary material

The online version of this article (10.1186/s13059-018-1422-4) contains supplementary material, which is available to authorized users.

Structure and function of the Nppa-Nppb cluster locus during heart development and disease

Atrial natriuretic factor and brain natriuretic peptide are two important biomarkers in clinical cardiology. These two natriuretic peptide hormones are encoded by the paralogous genes Nppa and Nppb, which are evolutionary conserved. Both genes are predominantly expressed by the heart muscle during the embryonic and fetal stages, and in particular Nppa expression is strongly reduced in the ventricles after birth. Upon cardiac stress, Nppa and Nppb are strongly upregulated in the ventricular myocardium. Much is known about the molecular and physiological ques inducing Nppa and Nppb expression; however, the transcriptional regulatory mechanisms of the Nppa-Nppb cluster in vivo has proven to be quite complex and is not well understood. In this review, we will provide recent insights into the dynamic and complex regulation of Nppa and Nppb during heart development and hypertrophy, and the association of this gene cluster with the cardiomyocyte-intrinsic program of heart regeneration.

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.

Scallop genome provides insights into evolution of bilaterian karyotype and development

Reconstructing the genomes of bilaterian ancestors is central to our understanding of animal evolution, where knowledge from ancient and/or slow-evolving bilaterian lineages is critical. Here we report a high-quality, chromosome-anchored reference genome for the scallop Patinopecten yessoensis, a bivalve mollusc that has a slow-evolving genome with many ancestral features. Chromosome-based macrosynteny analysis reveals a striking correspondence between the 19 scallop chromosomes and the 17 presumed ancestral bilaterian linkage groups at a level of conservation previously unseen, suggesting that the scallop may have a karyotype close to that of the bilaterian ancestor. Scallop Hox gene expression follows a new mode of subcluster temporal co-linearity that is possibly ancestral and may provide great potential in supporting diverse bilaterian body plans. Transcriptome analysis of scallop mantle eyes finds unexpected diversity in phototransduction cascades and a potentially ancient Pax2/5/8-dependent pathway for noncephalic eyes. The outstanding preservation of ancestral karyotype and developmental control makes the scallop genome a valuable resource for understanding early bilaterian evolution and biology.

The origins, scaling and loss of tetrapod digits

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

Reorganisation of Hoxd regulatory landscapes during the evolution of a snake-like body plan

Within land vertebrate species, snakes display extreme variations in their body plan, characterized by the absence of limbs and an elongated morphology. Such a particular interpretation of the basic vertebrate body architecture has often been associated with changes in the function or regulation of Hox genes. Here, we use an interspecies comparative approach to investigate different regulatory aspects at the snake HoxD locus. We report that, unlike in other vertebrates, snake mesoderm-specific enhancers are mostly located within the HoxD cluster itself rather than outside. In addition, despite both the absence of limbs and an altered Hoxd gene regulation in external genitalia, the limb-associated bimodal HoxD chromatin structure is maintained at the snake locus. Finally, we show that snake and mouse orthologous enhancer sequences can display distinct expression specificities. These results show that vertebrate morphological evolution likely involved extensive reorganisation at Hox loci, yet within a generally conserved regulatory framework.

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

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

Evidence-Based Theory for Integrated Genome Regulation of Ontogeny--An Unprecedented Role of Nuclear FGFR1 Signaling

Genetic experiments have positioned the fgfr1 gene at the top of the gene hierarchy that governs gastrulation, as well as the subsequent development of the major body axes, nervous system, muscles, and bones, by affecting downstream genes that control the cell cycle, pluripotency, and differentiation, as well as microRNAs. Studies show that this regulation is executed by a single protein, the nuclear isoform of FGFR1 (nFGFR1), which integrates signals from development‐initiating factors, such as retinoic acid (RA), and operates at the interface of genomic and epigenomic information. nFGFR1 cooperates with a multitude of transcriptional factors (TFs), and targets thousands of genes encoding for mRNAs, as well as miRNAs in top ontogenic networks. nFGFR1 binds to the promoters of ancient proto‐oncogenes and tumor suppressor genes, in addition to binding to metazoan morphogens that delineate body axes, and construct the nervous system, as well as mesodermal and endodermal tissues. The discovery of pan‐ontogenic gene programming by integrative nuclear FGFR1 signaling (INFS) impacts our understanding of ontogeny, as well as developmental pathologies, and holds new promise for reconstructive medicine, and cancer therapy. J. Cell. Physiol. 231: 1199–1218, 2016. © 2016 The Authors. Journal of Cellular Physiology published by Wiley Periodicals, Inc.

Nanoscale spatial organization of the HoxD gene cluster in distinct transcriptional states

Chromatin condensation plays an important role in the regulation of gene expression. Recently, it was shown that the transcriptional activation of Hoxd genes during vertebrate digit development involves modifications in 3D interactions within and around the HoxD gene cluster. This reorganization follows a global transition from one set of regulatory contacts to another, between two topologically associating domains (TADs) located on either side of the HoxD locus. Here, we use 3D DNA FISH to assess the spatial organization of chromatin at and around the HoxD gene cluster and report that although the two TADs are tightly associated, they appear as spatially distinct units. We measured the relative position of genes within the cluster and found that they segregate over long distances, suggesting that a physical elongation of the HoxD cluster can occur. We analyzed this possibility by super-resolution imaging (STORM) and found that tissues with distinct transcriptional activity exhibit differing degrees of elongation. We also observed that the morphological change of the HoxD cluster in developing digits is associated with its position at the boundary between the two TADs. Such variations in the fine-scale architecture of the gene cluster suggest causal links among its spatial configuration, transcriptional activation, and the flanking chromatin context.

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.

Overexpression of HOX genes is prevalent in Ewing sarcoma and is associated with altered epigenetic regulation of developmental transcription programs

The polycomb proteins BMI-1 and EZH2 are highly overexpressed by Ewing sarcoma (ES), a tumor of stem cell origin that is driven by EWS-ETS fusion oncogenes, most commonly EWS-FLI1. In the current study we analyzed expression of transcription programs that are controlled by polycomb proteins during embryonic development to determine if they are abnormal in ES. Our results show that polycomb target gene expression in ES deviates from normal tissues and stem cells and that, as expected, most targets are relatively repressed. However, we also discovered a paradoxical up regulation of numerous polycomb targets and these were highly enriched for homeobox (HOX) genes. Comparison of HOX profiles between malignant and non-malignant tissues revealed a distinctive HOX profile in ES, which was characterized by overexpression of posterior HOXD genes. In addition, ectopic expression of EWS-FLI1 during stem cell differentiation led to aberrant up regulation of posterior HOXD genes. Mechanistically, this up regulation was associated with altered epigenetic regulation. Specifically, ES and EWS-FLI1+ stem cells displayed a relative loss of polycomb-dependent H3K27me3 and gain of trithorax-dependent H3K4me3 at the promoters of posterior HOXD genes and also at the HOXD11.12 polycomb response element. In addition, a striking correlation was evident between HOXD13 and other genes whose regulation is coordinately regulated during embryonic development by distal enhancer elements. Together, these studies demonstrate that epigenetic regulation of polycomb target genes, in particular HOXD genes, is altered in ES and that these changes are mediated downstream of EWS-FLI1.

MicroRNA-196b is transcribed from an autonomous promoter and is directly regulated by Cdx2 and by posterior Hox proteins during embryogenesis

The miR-196 miRNA gene family located within the Hox gene clusters has been shown to function during embryogenesis and to be aberrantly expressed in various malignancies, including leukaemia, melanoma, and colorectal cancer. Despite its involvement in numerous biological processes, the control of miR-196 expression is still poorly defined. We identified the miR-196b promoter and found that the mature miR-196b originates from a large, non-coding primary transcript, which starts within an autonomous TATA box promoter and is not in physical continuity with either the Hoxa10 or Hoxa9 main primary transcripts. A ~680bp genomic fragment, spanning the pri-miR-196b transcription start site, is sufficient to recapitulate the neural tube expression pattern of miR-196 during embryogenesis. This region contains potential binding sites for Cdx and 5'Hox transcription factors. Two of these sites revealed to be necessary for neural tube expression and were bound in vivo by Cdx2 and Hoxd13. We show that Cdx2 is required for miR-196 expression and that both Cdx2 and 5'Hox, but not 3'Hox, are able to activate the miR-196b promoter. The possible role of Cdx2- and 5'Hox-mediated regulation of miR-196 expression in vertebrate anterior-posterior (AP) axis formation during embryogenesis is discussed.

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

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

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.

Temporal dynamics and developmental memory of 3D chromatin architecture at Hox gene loci

Hox genes are essential regulators of embryonic development. Their step-wise transcriptional activation follows their genomic topology and the various states of activation are subsequently memorized into domains of progressively overlapping gene products. We have analyzed the 3D chromatin organization of Hox clusters during their early activation in vivo, using high-resolution circular chromosome conformation capture. Initially, Hox clusters are organized as single chromatin compartments containing all genes and bivalent chromatin marks. Transcriptional activation is associated with a dynamic bi-modal 3D organization, whereby the genes switch autonomously from an inactive to an active compartment. These local 3D dynamics occur within a framework of constitutive interactions within the surrounding Topological Associated Domains, indicating that this regulation process is mostly cluster intrinsic. The step-wise progression in time is fixed at various body levels and thus can account for the chromatin architectures previously described at a later stage for different anterior to posterior levels.DOI: http://dx.doi.org/10.7554/eLife.02557.001.

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.

Epigenome-wide DNA methylation landscape of melanoma progression to brain metastasis reveals aberrations on homeobox D cluster associated with prognosis

Melanoma brain metastasis (MBM) represents a frequent complication of cutaneous melanoma. Despite aggressive multi-modality therapy, patients with MBM often have a survival rate of <1 year. Alteration in DNA methylation is a major hallmark of tumor progression and metastasis; however, it remains largely unexplored in MBM. In this study, we generated a comprehensive DNA methylation landscape through the use of genome-wide copy number, DNA methylation and gene expression data integrative analysis of melanoma progression to MBM. A progressive genome-wide demethylation in low CpG density and an increase in methylation level of CpG islands according to melanoma progression were observed. MBM-specific partially methylated domains (PMDs) affecting key brain developmental processes were identified. Differentially methylated CpG sites between MBM and lymph node metastasis (LNM) from patients with good prognosis were identified. Among the most significantly affected genes were the HOX family members. DNA methylation of HOXD9 gene promoter affected transcript and protein expression and was significantly higher in MBM than that in early stages. A MBM-specific PMD was identified in this region. Low methylation level of this region was associated with active HOXD9 expression, open chromatin and histone modifications associated with active transcription. Demethylating agent induced HOXD9 expression in melanoma cell lines. The clinical relevance of this finding was verified in an independent large cohort of melanomas (n = 145). Patients with HOXD9 hypermethylation in LNM had poorer disease-free and overall survival. This epigenome-wide study identified novel methylated genes with functional and clinical implications for MBM patients.

Evolutionary constraints favor a biophysical model explaining hox gene collinearity

The Hox gene collinearity enigma has often been approached using models based on biomolecular mechanisms. The biophysical model is an alternative approach based on the hypothesis that collinearity is caused by physical forces pulling the Hox genes from a territory where they are inactive to a distinct spatial domain where they are activated in a step by step manner. Such Hox gene translocations have recently been observed in support of the biophysical model. Genetic engineering experiments, performed on embryonic mice, gave rise to several unexpected mutant expressions that the biomolecular models cannot predict. On the contrary, the biophysical model offers convincing explanation. Evolutionary constraints consolidate the Hox clusters and as a result, denser and well organized clusters may create more efficient physical forces and a more emphatic manifestation of gene collinearity. This is demonstrated by stochastic modeling with white noise perturbing the expression of Hox genes. As study cases the genomes of mouse and amphioxus are used. The results support the working hypothesis that vertebrates have adopted their comparably more compact Hox clustering as a tool needed to develop more complex body structures. Several experiments are proposed in order to test further the physical forces hypothesis.

A genetic approach to the recruitment of PRC2 at the HoxD locus

Polycomb group (PcG) proteins are essential for the repression of key factors during early development. In Drosophila, the polycomb repressive complexes (PRC) associate with defined polycomb response DNA elements (PREs). In mammals, however, the mechanisms underlying polycomb recruitment at targeted loci are poorly understood. We have used an in vivo approach to identify DNA sequences of importance for the proper recruitment of polycomb proteins at the HoxD locus. We report that various genomic re-arrangements of the gene cluster do not strongly affect PRC2 recruitment and that relatively small polycomb interacting sequences appear necessary and sufficient to confer polycomb recognition and targeting to ectopic loci. In addition, a high GC content, while not sufficient to recruit PRC2, may help its local spreading. We discuss the importance of PRC2 recruitment over Hox gene clusters in embryonic stem cells, for their subsequent coordinated transcriptional activation during development.

Hox genes are essential for the proper organization of structures along the developing vertebrate body axis. These genes must be activated at a precise time and their premature transcription is deleterious to the organism. Early on, Hox gene clusters are covered by Polycomb Repressive protein Complexes (PRCs), which help keep these genes silent. However, the mechanism(s) that selectively recruit PRCs to these particular genomic loci remains elusive. We have used a collection of mutant mice carrying a set of deletions inside and outside the HoxD cluster to try and detect the presence of any DNA sequence of particular importance in this mechanism. We conclude that a range of low affinity sequences synergize to recruit PRCs over the gene cluster, which makes this process very robust and resistant to genetic perturbations.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

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.

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

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

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.

The origin of the tetrapod limb: from expeditions to enhancers

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

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.

Comparison of models for the collinearity of hox genes in the developmental axes of vertebrates

Hox gene clusters are very frequent in many animal genomes and their role in development is pivotal. Particularly in vertebrates, intensive efforts have established several properties of Hox clusters. The collinearity of Hox gene expressions (spatial, temporal and quantitative) is a common feature of the vertebrates. During the last decade, genetic engineering experiments have revealed some important facets of collinearity during limb and trunk development in mice. Two models have been proposed to explain all these properties. On one hand the ‘two-phases model’ makes use of the molecular regulatory mechanisms acting on the Hox genes. On the other hand, the’biophysical model’ is based on the signals transduced inside the cell nucleus and the generation of forces which apply on the cluster and lead to a coordinated activation of Hox genes. The two models differ fundamentally and a critical and detailed comparison is presented. Furthermore, experiments are proposed for which the two models provide divergent predictions. The outcome of these experiments will help to decide which of the two models is valid (if any).

Anterior-posterior differences in HoxD chromatin topology in limb development

A late phase of HoxD activation is crucial for the patterning and growth of distal structures across the anterior-posterior (A-P) limb axis of mammals. Polycomb complexes and chromatin compaction have been shown to regulate Hox loci along the main body axis in embryonic development, but the extent to which they have a role in limb-specific HoxD expression, an evolutionary adaptation defined by the activity of distal enhancer elements that drive expression of 5' Hoxd genes, has yet to be fully elucidated. We reveal two levels of chromatin topology that differentiate distal limb A-P HoxD activity. Using both immortalised cell lines derived from posterior and anterior regions of distal E10.5 mouse limb buds, and analysis in E10.5 dissected limb buds themselves, we show that there is a loss of polycomb-catalysed H3K27me3 histone modification and a chromatin decompaction over HoxD in the distal posterior limb compared with anterior. Moreover, we show that the global control region (GCR) long-range enhancer spatially colocalises with the 5' HoxD genomic region specifically in the distal posterior limb. This is consistent with the formation of a chromatin loop between 5' HoxD and the GCR regulatory module at the time and place of distal limb bud development when the GCR participates in initiating Hoxd gene quantitative collinearity and Hoxd13 expression. This is the first example of A-P differences in chromatin compaction and chromatin looping in the development of the mammalian secondary body axis (limb).

A fluorescence spotlight on the clockwork development and metabolism of bone

Biological phenomena that exhibit periodic activity are often referred as biorhythms or biological clocks. Among these, circadian rhythms, cyclic patterns reflecting a 24-h cycle, are the most obvious in many physiological activities including bone growth and metabolism. In the late 1990s, several clock genes were isolated and their primary structures and functions were identified. The feedback loop model of transcriptional factors was proposed to work as a circadian core oscillator not only in the suprachiasmatic nuclei of the anterior hypothalamus, which is recognized as the mammalian central clock, but also in various peripheral tissues including cartilage and bone. Looking back to embryonic development, the fundamental architecture of skeletal patterning is regulated by ultradian clocks that are defined as biorhythms that cycle more than once every 24 h. As post-genomic approaches, transcriptome analysis by micro-array and bioimaging assays to detect luminescent and fluorescent signals have been exploited to uncover a more comprehensive set of genes and spatio-temporal regulation of the clockwork machinery in animal models. In this review paper, we provide an overview of topics related to these molecular clocks in skeletal biology and medicine, and discuss how fluorescence imaging approaches can contribute to widening our views of this realm of biomedical science.

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.

Deregulated HOX genes in ameloblastomas are located in physical contiguity to keratin genes

The expression of the HOX gene network in mid-stage human tooth development mostly concerns the epithelial tooth germ compartment and involves the C and D HOX loci. To further dissect the HOX gene implication with tooth epithelium differentiation we compared the expression of the whole HOX network in human ameloblastomas, as paradigm of epithelial odontogenic tumors, with tooth germs. We identified two ameloblastoma molecular types with respectively low and high number of active HOX C genes. The highly expressing HOX C gene ameloblastomas were characterized by a strong keratinized phenotype. Locus C HOX genes are located on chromosome 12q13-15 in physical contiguity with one of the two keratin gene clusters included in the human genome. The most posterior HOX C gene, HOX C13, is capable to interact with hair keratin genes located on the other keratin gene cluster in physical contiguity with the HOX B locus on chromosome 17q21-22. Inside the HOX C locus, a 2.2 kb ncRNA (HOTAIR) able to repress transcription, in cis, along the entire HOX C locus and, in trans, at the posterior region of the HOX D locus has recently been identified. Interestingly both loci are deregulated in ameloblastomas. Our finding support an important role of the HOX network in characterizing the epithelial tooth compartment. Furthermore, the physical contiguity between locus C HOX and keratin genes in normal tooth epithelium and their deregulation in the neoplastic counterparts suggest they may act on the same mechanism potentially involved with epithelial tumorigenesis.

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.

Concerted involvement of Cdx/Hox genes and Wnt signaling in morphogenesis of the caudal neural tube and cloacal derivatives from the posterior growth zone

Decrease in Cdx dosage in an allelic series of mouse Cdx mutants leads to progressively more severe posterior vertebral defects. These defects are corrected by posterior gain of function of the Wnt effector Lef1. Precocious expression of Hox paralogous 13 genes also induces vertebral axis truncation by antagonizing Cdx function. We report here that the phenotypic similarity also applies to patterning of the caudal neural tube and uro-rectal tracts in Cdx and Wnt3a mutants, and in embryos precociously expressing Hox13 genes. Cdx2 inactivation after placentation leads to posterior defects, including incomplete uro-rectal septation. Compound mutants carrying one active Cdx2 allele in the Cdx4-null background (Cdx2/4), transgenic embryos precociously expressing Hox13 genes and a novel Wnt3a hypomorph mutant all manifest a comparable phenotype with similar uro-rectal defects. Phenotype and transcriptome analysis in early Cdx mutants, genetic rescue experiments and gene expression studies lead us to propose that Cdx transcription factors act via Wnt signaling during the laying down of uro-rectal mesoderm, and that they are operative in an early phase of these events, at the site of tissue progenitors in the posterior growth zone of the embryo. Cdx and Wnt mutations and premature Hox13 expression also cause similar neural dysmorphology, including ectopic neural structures that sometimes lead to neural tube splitting at caudal axial levels. These findings involve the Cdx genes, canonical Wnt signaling and the temporal control of posterior Hox gene expression in posterior morphogenesis in the different embryonic germ layers. They shed a new light on the etiology of the caudal dysplasia or caudal regression range of human congenital defects.

Reshuffling genomic landscapes to study the regulatory evolution of Hox gene clusters

The emergence of Vertebrata was accompanied by two rounds of whole-genome duplications. This enabled paralogous genes to acquire novel functions with high evolutionary potential, a process suggested to occur mostly by changes in gene regulation, rather than in protein sequences. In the case of Hox gene clusters, such duplications favored the appearance of distinct global regulations. To assess the impact of such "regulatory evolution" upon neo-functionalization, we developed PANTHERE (PAN-genomic Translocation for Heterologous Enhancer RE-shuffling) to bring the entire megabase-scale HoxD regulatory landscape in front of the HoxC gene cluster via a targeted translocation in vivo. At this chimeric locus, Hoxc genes could both interpret this foreign regulation and functionally substitute for their Hoxd counterparts. Our results emphasize the importance of evolving regulatory modules rather than their target genes in the process of neo-functionalization and offer a genetic tool to study the complexity of the vertebrate regulatory genome.