A CTCF-dependent mechanism underlies the Hox timer: relation to a segmented body plan

During gastrulation, Hox genes are activated in a time-sequence that follows the order of the genes along their clusters. This property, which is observed in all animals that develop following a progressive rostral-to-caudal morphogenesis, is associated with changes in the chromatin structure and epigenetic profiles of Hox clusters, suggesting a process at least partly based on sequential gene accessibility. Here, we discuss recent work on this issue, as well as a possible mechanism based on the surprising conservation in both the distribution and orientation of CTCF sites inside vertebrate Hox clusters.

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.

Sequential and directional insulation by conserved CTCF sites underlies the Hox timer in stembryos

During development, Hox genes are temporally activated according to their relative positions on their clusters, contributing to the proper identities of structures along the rostrocaudal axis. To understand the mechanism underlying this Hox timer, we used mouse embryonic stem cell-derived stembryos. Following Wnt signaling, the process involves transcriptional initiation at the anterior part of the cluster and a concomitant loading of cohesin complexes enriched on the transcribed DNA segments, that is, with an asymmetric distribution favoring the anterior part of the cluster. Chromatin extrusion then occurs with successively more posterior CTCF sites acting as transient insulators, thus generating a progressive time delay in the activation of more posterior-located genes due to long-range contacts with a flanking topologically associating domain. Mutant stembryos support this model and reveal that the presence of evolutionary conserved and regularly spaced intergenic CTCF sites controls the precision and the pace of this temporal mechanism.

Genetic cold cases; lessons from solving complex congenital limb disorders

Congenital genetic disorders affecting limb morphology in humans and other mammals are particularly well described, due to both their rather high frequencies of occurrence and the ease of their detection when expressed as severe forms. In most cases, their molecular and cellular etiology remained unknown long after their initial description, often for several decades, and sometimes close to a century. Over the past 20 yr, however, experimental and conceptual advances in our understanding of gene regulation, in particular over large genomic distances, have allowed these cold cases to be reopened and, eventually, for some of them to be solved. These investigations led not only to the isolation of the culprit genes and mechanisms, but also to the understanding of the often complex regulatory processes that are disturbed in such mutant genetic configurations. Here, we present several cases in which dormant regulatory mutations have been retrieved from the archives, starting from a historical perspective up to their molecular explanations. While some cases remain open, waiting for new tools and/or concepts to bring their investigations to an end, the solutions to others have contributed to our understanding of particular features often found in the regulation of developmental genes and hence can be used as benchmarks to address the impact of noncoding variants in the future.

Developmental and evolutionary comparative analysis of a regulatory landscape in mouse and chicken

Modifications in gene regulation are driving forces in the evolution of organisms. Part of these changes involve cis-regulatory elements (CREs), which contact their target genes through higher-order chromatin structures. However, how such architectures and variations in CREs contribute to transcriptional evolvability remains elusive. We use Hoxd genes as a paradigm for the emergence of regulatory innovations, as many relevant enhancers are located in a regulatory landscape highly conserved in amniotes. Here, we analysed their regulation in murine vibrissae and chicken feather primordia, two skin appendages expressing different Hoxd gene subsets, and compared the regulation of these genes in these appendages with that in the elongation of the posterior trunk. In the two former structures, distinct subsets of Hoxd genes are contacted by different lineage-specific enhancers, probably as a result of using an ancestral chromatin topology as an evolutionary playground, whereas the gene regulation that occurs in the mouse and chicken embryonic trunk partially relies on conserved CREs. A high proportion of these non-coding sequences active in the trunk have functionally diverged between species, suggesting that transcriptional robustness is maintained, despite considerable divergence in enhancer sequences.

Summary: Analyses of the relationships between chromatin architecture and regulatory activities at the HoxD locus show that ancestral transcription patterns can be maintained while new regulations evolve.

Sequential in cis mutagenesis in vivo reveals various functions for CTCF sites at the mouse HoxD cluster

Mammalian Hox gene clusters contain a range of CTCF binding sites. In addition to their importance in organizing a TAD border, which isolates the most posterior genes from the rest of the cluster, the positions and orientations of these sites suggest that CTCF may be instrumental in the selection of various subsets of contiguous genes, which are targets of distinct remote enhancers located in the flanking regulatory landscapes. We examined this possibility by producing an allelic series of cumulative in cis mutations in these sites, up to the abrogation of CTCF binding in the five sites located on one side of the TAD border. In the most impactful alleles, the global chromatin architecture of the locus was modified, yet not drastically, illustrating that CTCF sites located on one side of a strong TAD border are sufficient to organize at least part of this insulation. Spatial colinearity in the expression of these genes along the major body axis was nevertheless maintained, despite abnormal expression boundaries. In contrast, strong effects were scored in the selection of target genes responding to particular enhancers, leading to the misregulation of Hoxd genes in specific structures. Altogether, while most enhancer-promoter interactions can occur in the absence of this series of CTCF sites, the binding of CTCF in the Hox cluster is required to properly transform a rather unprecise process into a highly discriminative mechanism of interactions, which is translated into various patterns of transcription accompanied by the distinctive chromatin topology found at this locus. Our allelic series also allowed us to reveal the distinct functional contributions for CTCF sites within this Hox cluster, some acting as insulator elements, others being necessary to anchor or stabilize enhancer-promoter interactions, and some doing both, whereas they all together contribute to the formation of a TAD border. This variety of tasks may explain the amazing evolutionary conservation in the distribution of these sites among paralogous Hox clusters or between various vertebrates.

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.

Analysis of Polycerate Mutants Reveals the Evolutionary Co-option of HOXD1 for Horn Patterning in Bovidae

In the course of evolution, pecorans (i.e., higher ruminants) developed a remarkable diversity of osseous cranial appendages, collectively referred to as “headgear,” which likely share the same origin and genetic basis. However, the nature and function of the genetic determinants underlying their number and position remain elusive. Jacob and other rare populations of sheep and goats are characterized by polyceraty, the presence of more than two horns. Here, we characterize distinct POLYCERATE alleles in each species, both associated with defective HOXD1 function. We show that haploinsufficiency at this locus results in the splitting of horn bud primordia, likely following the abnormal extension of an initial morphogenetic field. These results highlight the key role played by this gene in headgear patterning and illustrate the evolutionary co-option of a gene involved in the early development of bilateria to properly fix the position and number of these distinctive organs of Bovidae.

Induction of a chromatin boundary in vivo upon insertion of a TAD border

Mammalian genomes are partitioned into sub-megabase to megabase-sized units of preferential interactions called topologically associating domains or TADs, which are likely important for the proper implementation of gene regulatory processes. These domains provide structural scaffolds for distant cis regulatory elements to interact with their target genes within the three-dimensional nuclear space and architectural proteins such as CTCF as well as the cohesin complex participate in the formation of the boundaries between them. However, the importance of the genomic context in providing a given DNA sequence the capacity to act as a boundary element remains to be fully investigated. To address this question, we randomly relocated a topological boundary functionally associated with the mouse HoxD gene cluster and show that it can indeed act similarly outside its initial genomic context. In particular, the relocated DNA segment recruited the required architectural proteins and induced a significant depletion of contacts between genomic regions located across the integration site. The host chromatin landscape was re-organized, with the splitting of the TAD wherein the boundary had integrated. These results provide evidence that topological boundaries can function independently of their site of origin, under physiological conditions during mouse development.

During development, enhancer sequences tightly regulate the spatio-temporal expression of target genes often located hundreds of kilobases away. This complex process is made possible by the folding of chromatin into domains, which are separated from one another by specific genomic regions referred to as boundaries. In order to understand whether such boundary sequences require their particular genomic contexts to achieve their isolating effect, we analyzed the impact of introducing one such boundary, taken from the HoxD locus, into a distinct topological domain. We show that this ectopic boundary splits the host domain into two sub-domains and affects the expression levels of a neighboring gene. We conclude that this sequence can work independently from its genomic context and thus carries all the information necessary to act as a boundary element.

Time-sequenced transcriptomes of developing distal mouse limb buds: A comparative tissue layer analysis

BACKGROUND

The development of the amniote limb has been an important model system to study patterning mechanisms and morphogenesis. For proper growth and patterning, it requires the interaction between the distal sub-apical mesenchyme and the apical ectodermal ridge (AER) that involve the separate implementation of coordinated and tissue-specific genetic programs.

RESULTS

Here, we produce and analyze the transcriptomes of both distal limb mesenchymal progenitors and the overlying ectodermal cells, following time-coursed dissections that cover from limb bud initiation to fully patterned limbs. The comparison of transcriptomes within each layer as well as between layers over time, allowed the identification of specific transcriptional signatures for each of the developmental stages. Special attention was given to the identification of genes whose transcription dynamics suggest a previously unnoticed role in the context of limb development and also to signaling pathways enriched between layers.

CONCLUSION

We interpret the transcriptomic data in light of the known development pattern and we conclude that a major transcriptional transition occurs in distal limb buds between E9.5 and E10.5, coincident with the switch from an early phase continuation of the signature of trunk progenitors, related to the initial proximo distal specification, to a late intrinsic phase of development.

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.

Chromatin topology and the timing of enhancer function at the HoxD locus

The HoxD gene cluster is critical for proper limb formation in tetrapods. In the emerging limb buds, different subgroups of Hoxd genes respond first to a proximal regulatory signal, then to a distal signal that organizes digits. These two regulations are exclusive from one another and emanate from two distinct topologically associating domains (TADs) flanking HoxD, both containing a range of appropriate enhancer sequences. The telomeric TAD (T-DOM) contains several enhancers active in presumptive forearm cells and is divided into two sub-TADs separated by a CTCF-rich boundary, which defines two regulatory submodules. To understand the importance of this particular regulatory topology to control Hoxd gene transcription in time and space, we either deleted or inverted this sub-TAD boundary, eliminated the CTCF binding sites, or inverted the entire T-DOM to exchange the respective positions of the two sub-TADs. The effects of such perturbations on the transcriptional regulation of Hoxd genes illustrate the requirement of this regulatory topology for the precise timing of gene activation. However, the spatial distribution of transcripts was eventually resumed, showing that the presence of enhancer sequences, rather than either their exact topology or a particular chromatin architecture, is the key factor. We also show that the affinity of enhancers to find their natural target genes can overcome the presence of both a strong TAD border and an unfavorable orientation of CTCF sites.

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 regulatory landscapes of developmental genes

Regulatory landscapes have been defined in vertebrates as large DNA segments containing diverse enhancer sequences that produce coherent gene transcription. These genomic platforms integrate multiple cellular signals and hence can trigger pleiotropic expression of developmental genes. Identifying and evaluating how these chromatin regions operate may be difficult as the underlying regulatory mechanisms can be as unique as the genes they control. In this brief article and accompanying poster, we discuss some of the ways in which regulatory landscapes operate, illustrating these mechanisms using genes important for vertebrate development as examples. We also highlight some of the techniques available to researchers for analysing regulatory landscapes.

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.

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.

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.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Rescue of an aggressive female sexual courtship in mice by CRISPR/Cas9 secondary mutation in vivo

OBJECTIVE

We had previously reported a mouse line carrying the Atypical female courtship (HoxD ^Afc ) allele, where an ectopic accumulation of Hoxd10 transcripts was observed in a sparse population of cells in the adult isocortex, as a result of a partial deletion of the HoxD gene cluster. Female mice carrying this allele displayed an exacerbated paracopulatory behavior, culminating in a severe mutilation of the studs' external genitals. To unequivocally demonstrate that this intriguing phenotype was indeed caused by an illegitimate function of the HOXD10 protein, we use CRISPR/Cas9 technology to induce a microdeletion into the homeobox of the Hoxd10 gene in cis with the HoxD ^Afc allele.

RESULTS

Females carrying this novel HoxD^{Del(1-9)d10hd} allele no longer mutilate males. We conclude that a brain malfunction leading to a severe pathological behavior can be caused by the mere binding to DNA of a transcription factor expressed ectopically. We also show that in HoxD ^Afc mice, Hoxd10 was expressed in cells containing glutamate decarboxylase (Gad1) and Cholecystokinin (Cck) transcripts, corroborating our proposal that a small fraction of GABAergic neurons in adult hippocampus may participate to some aspects of female courtship.

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

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

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.

Topological Domains, Metagenes, and the Emergence of Pleiotropic Regulations at Hox Loci

Hox gene clusters of jaw vertebrates display a tight genomic organization, which has no equivalent in any other bilateria genomes sequenced thus far. It was previously argued that such a topological consolidation toward a condensed, metagenic structure was due to the accumulation of long-range regulations flanking Hox loci, a phenomenon made possible by the successive genome duplications that occurred at the root of the vertebrate lineage, similar to gene neofunctionalization but applied to a coordinated multigenic system. Here, we propose that the emergence of such large vertebrate regulatory landscapes containing a range of global enhancers was greatly facilitated by the presence of topologically associating domains (TADs). These chromatin domains, mostly constitutive, may have been used as genomic niches where novel regulations could evolve due to both the preexistence of a structural backbone poised to integrate novel regulatory inputs, and a highly adaptive transcriptional readout. We propose a scenario for the coevolution of such TADs and the emergence of pleiotropy at ancestral vertebrate Hox loci.

Large scale genomic reorganization of topological domains at the HoxD locus

BACKGROUND

The transcriptional activation of HoxD genes during mammalian limb development involves dynamic interactions with two topologically associating domains (TADs) flanking the HoxD cluster. In particular, the activation of the most posterior HoxD genes in developing digits is controlled by regulatory elements located in the centromeric TAD (C-DOM) through long-range contacts.

RESULTS

To assess the structure-function relationships underlying such interactions, we measured compaction levels and TAD discreteness using a combination of chromosome conformation capture (4C-seq) and DNA FISH. We assessed the robustness of the TAD architecture by using a series of genomic deletions and inversions that impact the integrity of this chromatin domain and that remodel long-range contacts. We report multi-partite associations between HoxD genes and up to three enhancers. We find that the loss of native chromatin topology leads to the remodeling of TAD structure following distinct parameters.

CONCLUSIONS

Our results reveal that the recomposition of TAD architectures after large genomic re-arrangements is dependent on a boundary-selection mechanism in which CTCF mediates the gating of long-range contacts in combination with genomic distance and sequence specificity. Accordingly, the building of a recomposed TAD at this locus depends on distinct functional and constitutive parameters.

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.

Snakes: hatching of a model system for Evo-Devo?

Evo-Devo studies rely on a collection of animal model systems belonging to different phylogenetic branches to try and understand how organisms carrying a similar set of genes and pathways can develop into such a variety of shapes and sizes. The squamate clade, however, has only recently started to receive the attention it deserves in particular due to extreme morphological and metabolic aspects and, consequently, the important insights that it could bring in different fields. The recent sequencing of several squamate genomes as well as the generation of high quality trancriptomes for different snake tissues now provide the necessary tools to complement biological studies. Here, we briefly report on recent work involving developing snake embryos to illustrate their interest to assess vertebrate developmental mechanisms. We also discuss the relevance to use snake species as Evo-Devo model systems and potential ways to cross the important limitations intrinsically associated with developmental and genetic studies of these fascinating animals.

Tetrapod axial evolution and developmental constraints; Empirical underpinning by a mouse model

The tetrapod vertebral column has become increasingly complex during evolution as an adaptation to a terrestrial life. At the same time, the evolution of the vertebral formula became subject to developmental constraints acting on the size of the cervical and thoraco-lumbar regions. In the course of our studies concerning the evolution of Hox gene regulation, we produced a transgenic mouse model expressing fish Hox genes, which displayed a reduced number of thoraco-lumbar vertebrae and concurrent sacral homeotic transformations. Here, we analyze this mutant stock and conclude that the ancestral, pre-tetrapodial Hox code already possessed the capacity to induce vertebrae with sacral characteristics. This suggests that alterations in the interpretation of the Hox code may have participated to the evolution of this region in tetrapods, along with potential modifications of the HOX proteins themselves. With its reduced vertebral number, this mouse stock violates a previously described developmental constraint, which applies to the thoraco-lumbar region. The resulting offset between motor neuron morphology, vertebral patterning and the relative positioning of hind limbs illustrates that the precise orchestration of the Hox-clock in parallel with other ontogenetic pathways places constraints on the evolvability of the body plan.

•

A transgenic mouse line expressing fish Hox genes has anterior homeotic transformations.

•

Fish Hox genes are capable of inducing tetrapod specific vertebral characters.

•

A sacral Hox-code influences adult hindlimb position, yet not the position of limb budding.

Clustering of mammalian Hox genes with other H3K27me3 targets within an active nuclear domain

Embryogenesis requires the precise activation and repression of many transcriptional regulators. The Polycomb group proteins and the associated H3K27me3 histone mark are essential to maintain the inactive state of many of these genes. Mammalian Hox genes are targets of Polycomb proteins and form local 3D clusters centered on the H3K27me3 mark. More distal contacts have also been described, yet their selectivity, dynamics, and relation to other layers of chromatin organization remained elusive. We report that repressed Hox genes form mutual intra- and interchromosomal interactions with other genes located in strong domains labeled by H3K27me3. These interactions occur in a central and active nuclear environment that consists of the HiC compartment A, away from peripheral lamina-associated domains. Interactions are independent of nearby H3K27me3-marked loci and determined by chromosomal distance and cell-type-specific scaling factors, thus inducing a moderate reorganization during embryogenesis. These results provide a simplified view of nuclear organization whereby Polycomb proteins may have evolved to repress genes located in gene-dense regions whose position is restricted to central, active, nuclear environments.

Considerations when investigating lncRNA function in vivo

Although a small number of the vast array of animal long non-coding RNAs (lncRNAs) have known effects on cellular processes examined in vitro, the extent of their contributions to normal cell processes throughout development, differentiation and disease for the most part remains less clear. Phenotypes arising from deletion of an entire genomic locus cannot be unequivocally attributed either to the loss of the lncRNA per se or to the associated loss of other overlapping DNA regulatory elements. The distinction between cis- or trans-effects is also often problematic. We discuss the advantages and challenges associated with the current techniques for studying the in vivo function of lncRNAs in the light of different models of lncRNA molecular mechanism, and reflect on the design of experiments to mutate lncRNA loci. These considerations should assist in the further investigation of these transcriptional products of the genome.

DOI:http://dx.doi.org/10.7554/eLife.03058.001

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.

Evolving Hox activity profiles govern diversity in locomotor systems

The emergence of limb-driven locomotor behaviors was a key event in the evolution of vertebrates and fostered the transition from aquatic to terrestrial life. We show that the generation of limb-projecting lateral motor column (LMC) neurons in mice relies on a transcriptional autoregulatory module initiated via transient activity of multiple genes within the HoxA and HoxC clusters. Repression of this module at thoracic levels restricts expression of LMC determinants, thus dictating LMC position relative to the limbs. This suppression is mediated by a key regulatory domain that is specifically found in the Hoxc9 proteins of appendage-bearing vertebrates. The profile of Hoxc9 expression inversely correlates with LMC position in land vertebrates and likely accounts for the absence of LMC neurons in limbless species such as snakes. Thus, modulation of both Hoxc9 protein function and Hoxc9 gene expression likely contributed to evolutionary transitions between undulatory and ambulatory motor circuit connectivity programs.

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.

Topology of mammalian developmental enhancers and their regulatory landscapes

How a complex animal can arise from a fertilized egg is one of the oldest and most fascinating questions of biology, the answer to which is encoded in the genome. Body shape and organ development, and their integration into a functional organism all depend on the precise expression of genes in space and time. The orchestration of transcription relies mostly on surrounding control sequences such as enhancers, millions of which form complex regulatory landscapes in the non-coding genome. Recent research shows that high-order chromosome structures make an important contribution to enhancer functionality by triggering their physical interactions with target genes.

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.

Combined function of HoxA and HoxB clusters in neural crest cells

The evolution of chordates was accompanied by critical anatomical innovations in craniofacial development, along with the emergence of neural crest cells. The potential of these cells to implement a craniofacial program in part depends upon the (non-)expression of Hox genes. For instance, the development of jaws requires the inhibition of Hox genes function in the first pharyngeal arch. In contrast, Hox gene products induce craniofacial structures in more caudal territories. To further investigate which Hox gene clusters are involved in this latter role, we generated HoxA;HoxB cluster double mutant animals in cranial neural crest cells. We observed the appearance of a supernumerary dentary-like bone with an endochondral ossification around a neo-Meckel's cartilage matrix and an attachment of neo-muscle demonstrating that HoxB genes enhance the phenotype induced by the deletion of the HoxA cluster alone. In addition, a cervical and hypertrophic thymus was associated with the supernumerary dentary-like bone, which may reflect its ancestral position near the filtrating system. Altogether these results show that the HoxA and HoxB clusters cooperated during evolution to lead to present craniofacial diversity.

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.

Duplications of hox gene clusters and the emergence of vertebrates

The vertebrate body plan is characterized by an increased complexity relative to that of all other chordates and large-scale gene amplifications have been associated with key morphological innovations leading to their remarkable evolutionary success. Here, we use compound full Hox clusters deletions to investigate how Hox genes duplications may have contributed to the emergence of vertebrate-specific innovations. We show that the combined deletion of HoxA and HoxB leads to an atavistic heart phenotype, suggesting that the ancestral HoxA/B cluster was co-opted to help in diversifying the complex organ in vertebrates. Other phenotypic effects observed seem to illustrate the resurgence of ancestral (plesiomorphic) features. This indicates that the duplications of Hox clusters were associated with the recruitment or formation of novel cis-regulatory controls, which were key to the evolution of many vertebrate features and hence to the evolutionary radiation of this group.

Chromatin looping and organization at developmentally regulated gene loci

Developmentally regulated genes are often controlled by distant enhancers, silencers and insulators, to implement their correct transcriptional programs. In recent years, the development of 3C and derived techniques (4C, 5C, HiC, ChIA-PET, etc.) has confirmed that chromatin looping is an important mechanism for the transfer of regulatory information in mammalian cells. At many developmentally regulated gene loci, transcriptional activation is indeed accompanied by the formation of chromatin loops between genes and distant enhancers. Similarly, dynamic looping between insulator elements and changes in local 3D organization may be observed upon variation in transcriptional activity. Chromatin looping also occurs at silent gene loci, where its function remains less understood. In lineage-committed cells, partial 3D configurations are detected at loci that are activated at later stages. However, these partial configurations usually lack promoter-enhancer loops that accompany transcriptional activation, suggesting they have structural functions. Definitive evidence for a repressive role of chromatin looping is still lacking. Chromatin loops have been reported at repressed loci but, alternatively, they may act as a distraction for active loops. Together, these mechanisms allow fine-tuning of regulatory programs, thus providing further diversity in the transcriptional control of developmentally regulated gene loci.

A genetic basis for altered sexual behavior in mutant female mice

Although neural substrates of mammalian female mating behavior have been described, the association between complex courtship activity and specific underlying mechanisms remains elusive. We have isolated a mouse line that unexpectedly shows altered female social behavior with increased investigation of males and increased genital biting. We investigated adult individuals by behavioral observation and genetic and molecular neuroanatomy methods. We report exacerbated inverse pursuits and incapacitating bites directed at the genitals of stud males. This extreme deviation from wild-type female courtship segregates with a deletion of the Hoxd1 to Hoxd9 genomic region. This dominant Atypical female courtship allele (HoxD(Afc)) induces ectopic Hoxd10 gene expression in several regions in newborn forebrain transitorily and stably in a sparse subpopulation of cells in the cornu ammonis fields of adult hippocampus, which may thus lead to an abnormal modulation in the sexual behavior of mutant females. The resulting compulsive sexual solicitation behavior displayed by the most affected individuals suggests new avenues to study the genetic and molecular bases of normal and pathological mammalian affect and raises the potential involvement of the hippocampus in the control of female courtship behavior. The potential relevance to human 2q.31.1 microdeletion syndrome is discussed.

Landscapes and archipelagos: spatial organization of gene regulation in vertebrates

Vertebrate genes controlling critical developmental processes are often regulated by complex sets of global enhancer sequences, located at a distance, within neighboring gene deserts. Recent technological advances have made it possible to investigate the spatial organization of these 'regulatory landscapes'. The integration of such datasets with information on chromatin status, transcriptional activity and nuclear localization of these loci, as well as the effects of genetic modifications thereof, may bring a more comprehensive understanding of tissue- and/or stage-specific gene regulation in both normal and pathological contexts. Here, we review the impact of recent technological advances on our understanding of large-scale gene regulation in vertebrates, by focusing on paradigmatic gene loci.

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.