Transgenic analysis of Hoxd gene regulation during digit development

Hnrnpk is essential for embryonic limb bud development as a transcription activator and a collaborator of insulator protein Ctcf

Proper development of the limb bud relies on the concordance of various signals, but its molecular mechanisms have not yet been fully illustrated. Here we report that heterogeneous nuclear ribonucleoprotein K (hnRNPK) is essential for limb bud development. Its ablation in the limb bud results in limbless forelimbs and severe deformities of the hindlimbs. In terms of mechanism, hnRNPK functions as a transcription activator for the vital genes involved in the three regulatory axes of limb bud development. Simultaneously, for the first time we elucidate that hnRNPK binds to and coordinates with the insulator protein CCCTC binding factor (CTCF) to maintain a three-dimensional chromatin architecture. Ablation of hnRNPK weakens the binding strength of CTCF to topologically associating domain (TAD) boundaries, then leading to the loose TADs, and decreased interactions between promoters and enhancers, and further decreased transcription of developmental genes. Our study establishes a fundamental and novel role of hnRNPK in regulating limb bud development.

Shared retinoic acid responsive enhancers coordinately regulate nascent transcription of Hoxb coding and non-coding RNAs in the developing mouse neural tube

Signaling pathways regulate the patterns of Hox gene expression that underlie their functions in the specification of axial identity. Little is known about the properties of cis-regulatory elements and underlying transcriptional mechanisms that integrate graded signaling inputs to coordinately control Hox expression. Here, we optimized a single molecule fluorescent in situ hybridization (smFISH) technique with probes spanning introns to evaluate how three shared retinoic acid response element (RARE)-dependent enhancers in the Hoxb cluster regulate patterns of nascent transcription in vivo at the level of single cells in wild-type and mutant embryos. We predominately detect nascent transcription of only a single Hoxb gene in each cell, with no evidence for simultaneous co-transcriptional coupling of all or specific subsets of genes. Single and/or compound RARE mutations indicate that each enhancer differentially impacts global and local patterns of nascent transcription, suggesting that selectivity and competitive interactions between these enhancers is important to robustly maintain the proper levels and patterns of nascent Hoxb transcription. This implies that rapid and dynamic regulatory interactions potentiate transcription of genes through combined inputs from these enhancers in coordinating the retinoic acid response.

Shaping Hox gene activity to generate morphological diversity across vertebrate phylogeny

The importance of Hox genes for the development and evolution of the vertebrate axial skeleton and paired appendages has been recognized for already several decades. The steady growth of genomic sequence data from an increasing number of vertebrate species, together with the improvement of methods to analyze genomic structure and interactions, as well as to control gene activity in various species has refined our understanding of Hox gene activity in development and evolution. Here, I will review recent data addressing the influence of Hox regulatory processes in the evolution of the fins and the emergence of the tetrapod limb. In addition, I will discuss the involvement of posterior Hox genes in the control of vertebrate axial extension, focusing on an apparently divergent activity that Hox13 paralog group genes have on the regulation of tail bud development in mouse and zebrafish embryos.

Transcriptional Regulation and Implications for Controlling Hox Gene Expression

Hox genes play key roles in axial patterning and regulating the regional identity of cells and tissues in a wide variety of animals from invertebrates to vertebrates. Nested domains of Hox expression generate a combinatorial code that provides a molecular framework for specifying the properties of tissues along the A–P axis. Hence, it is important to understand the regulatory mechanisms that coordinately control the precise patterns of the transcription of clustered Hox genes required for their roles in development. New insights are emerging about the dynamics and molecular mechanisms governing transcriptional regulation, and there is interest in understanding how these may play a role in contributing to the regulation of the expression of the clustered Hox genes. In this review, we summarize some of the recent findings, ideas and emerging mechanisms underlying the regulation of transcription in general and consider how they may be relevant to understanding the transcriptional regulation of Hox genes.

Transcriptional enhancers and their communication with gene promoters

Transcriptional enhancers play a key role in the initiation and maintenance of gene expression programmes, particularly in metazoa. How these elements control their target genes in the right place and time is one of the most pertinent questions in functional genomics, with wide implications for most areas of biology. Here, we synthesise classic and recent evidence on the regulatory logic of enhancers, including the principles of enhancer organisation, factors that facilitate and delimit enhancer-promoter communication, and the joint effects of multiple enhancers. We show how modern approaches building on classic insights have begun to unravel the complexity of enhancer-promoter relationships, paving the way towards a quantitative understanding of gene control.

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.

Co-Regulated Genes and Gene Clusters

There are many co-regulated genes in eukaryotic cells. The coordinated activation or repression of such genes occurs at specific stages of differentiation, or under the influence of external stimuli. As a rule, co-regulated genes are dispersed in the genome. However, there are also gene clusters, which contain paralogous genes that encode proteins with similar functions. In this aspect, they differ significantly from bacterial operons containing functionally linked genes that are not paralogs. In this review, we discuss the reasons for the existence of gene clusters in vertebrate cells and propose that clustering is necessary to ensure the possibility of selective activation of one of several similar genes.

Gene Regulatory and Expression Differences between Mouse and Pig Limb Buds Provide Insights into the Evolutionary Emergence of Artiodactyl Traits

Digit loss/reductions are evolutionary adaptations in cursorial mammals such as pigs. To gain mechanistic insight into these processes, we performed a comparative molecular analysis of limb development in mouse and pig embryos, which revealed a loss of anterior-posterior polarity during distal progression of pig limb bud development. These alterations in pig limb buds are paralleled by changes in the mesenchymal response to Sonic hedgehog (SHH) signaling, which is altered upstream of the reduction and loss of Fgf8 expression in the ectoderm that overlaps the reduced and vestigial digit rudiments of the pig handplate, respectively. Furthermore, genome-wide open chromatin profiling using equivalent developmental stages of mouse and pig limb buds reveals the functional divergence of about one-third of the regulatory genome. This study uncovers widespread alterations in the regulatory landscapes of genes essential for limb development that likely contributed to the morphological diversion of artiodactyl limbs from the pentadactyl archetype of tetrapods.

A developmental perspective of homology and evolutionary novelty

The development and evolution of multicellular body plans is complex. Many distinct organs and body parts must be reproduced at each generation, and those that are traceable over long time scales are considered homologous. Among the most pressing and least understood phenomena in evolutionary biology is the mode by which new homologs, or "novelties" are introduced to the body plan and whether the developmental changes associated with such evolution deserve special treatment. In this chapter, we address the concepts of homology and evolutionary novelty through the lens of development. We present a series of case studies, within insects and vertebrates, from which we propose a developmental model of multicellular organ identity. With this model in hand, we make predictions regarding the developmental evolution of body plans and highlight the need for more integrative analysis of developing systems.

HOX13-dependent chromatin accessibility underlies the transition towards the digit development program

Hox genes encode transcription factors (TFs) that establish morphological diversity in the developing embryo. The similar DNA-binding motifs of the various HOX TFs contrast with the wide-range of HOX-dependent genetic programs. The influence of the chromatin context on HOX binding specificity remains elusive. Here, we used the developing limb as a model system to compare the binding specificity of HOXA13 and HOXD13 (HOX13 hereafter), which are required for digit formation, and HOXA11, involved in forearm/leg development. We find that upon ectopic expression in distal limb buds, HOXA11 binds sites normally HOX13-specific. Importantly, these sites are loci whose chromatin accessibility relies on HOX13. Moreover, we show that chromatin accessibility specific to the distal limb requires HOX13 function. Based on these results, we propose that HOX13 TFs pioneer the distal limb-specific chromatin accessibility landscape for the proper implementation of the distal limb developmental program.

Pioneer factors direct cell fate through switching inaccessible chromatin to an accessible state at specific target enhancers. Here the authors show that HOX13 transcription factors have a pioneer activity which is required for the proper implementation of the distal limb developmental program.

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.

Phenotype loss is associated with widespread divergence of the gene regulatory landscape in evolution

Detecting the genomic changes underlying phenotypic changes between species is a main goal of evolutionary biology and genomics. Evolutionary theory predicts that changes in cis-regulatory elements are important for morphological changes. We combined genome sequencing, functional genomics and genome-wide comparative analyses to investigate regulatory elements in lineages that lost morphological traits. We first show that limb loss in snakes is associated with widespread divergence of limb regulatory elements. We next show that eye degeneration in subterranean mammals is associated with widespread divergence of eye regulatory elements. In both cases, sequence divergence results in an extensive loss of transcription factor binding sites. Importantly, diverged regulatory elements are associated with genes required for normal limb patterning or normal eye development and function, suggesting that regulatory divergence contributed to the loss of these phenotypes. Together, our results show that genome-wide decay of the phenotype-specific cis-regulatory landscape is a hallmark of lost morphological traits.

Cis-regulatory elements are important factors for morphological changes. Here, the authors show widespread divergence of limb and eye regulatory elements in limb loss in snakes and eye degeneration in subterranean mammals respectively.

Distal Limb Patterning Requires Modulation of cis-Regulatory Activities by HOX13

The combinatorial expression of Hox genes along the body axes is a major determinant of cell fate and plays a pivotal role in generating the animal body plan. Loss of HOXA13 and HOXD13 transcription factors (HOX13) leads to digit agenesis in mice, but how HOX13 proteins regulate transcriptional outcomes and confer identity to the distal-most limb cells has remained elusive. Here, we report on the genome-wide profiling of HOXA13 and HOXD13 in vivo binding and changes of the transcriptome and chromatin state in the transition from the early to the late-distal limb developmental program, as well as in Hoxa13^−/−; Hoxd13^−/−limbs. Our results show that proper termination of the early limb transcriptional program and activation of the late-distal limb program are coordinated by the dual action of HOX13 on cis-regulatory modules.

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.

Alterations in anterior-posterior patterning and its accompanying changes along the proximal-distal axis during the fin-to-limb transition

The evolution from fins to limbs was one of the most successful innovations for vertebrates, allowing them to vastly expand their behaviors and habitats. Fossil records suggest that morphological changes occurred not only along the proximal-distal axis included appearance of the autopod, but also occurred along the anterior-posterior axis included reductions in the size and number of basal bones and digits. This review focuses on recent progress in developmental and genetic studies aimed at elucidating the mechanisms underlying alteration of anterior-posterior patterning and its accompanying changes along the proximal-distal axis during the fin-to-limb transition.

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.

HoxD expression in the fin-fold compartment of basal gnathostomes and implications for paired appendage evolution

The role of Homeobox transcription factors during fin and limb development have been the focus of recent work investigating the evolutionary origin of limb-specific morphologies. Here we characterize the expression of HoxD genes, as well as the cluster-associated genes Evx2 and LNP, in the paddlefish Polyodon spathula, a basal ray-finned fish. Our results demonstrate a collinear pattern of nesting in early fin buds that includes HoxD14, a gene previously thought to be isolated from global Hox regulation. We also show that in both Polyodon and the catshark Scyliorhinus canicula (a representative chondrichthyan) late phase HoxD transcripts are present in cells of the fin-fold and co-localize with And1, a component of the dermal skeleton. These new data support an ancestral role for HoxD genes in patterning the fin-folds of jawed vertebrates, and fuel new hypotheses about the evolution of cluster regulation and the potential downstream differentiation outcomes of distinct HoxD-regulated compartments.

A single three-dimensional chromatin compartment in amphioxus indicates a stepwise evolution of vertebrate Hox bimodal regulation

The HoxA and HoxD gene clusters of jawed vertebrates are organized into bipartite three-dimensional chromatin structures that separate long-range regulatory inputs coming from the anterior and posterior Hox-neighboring regions. This architecture is instrumental in allowing vertebrate Hox genes to pattern disparate parts of the body, including limbs. Almost nothing is known about how these three-dimensional topologies originated. Here we perform extensive 4C-seq profiling of the Hox cluster in embryos of amphioxus, an invertebrate chordate. We find that, in contrast to the architecture in vertebrates, the amphioxus Hox cluster is organized into a single chromatin interaction domain that includes long-range contacts mostly from the anterior side, bringing distant cis-regulatory elements into contact with Hox genes. We infer that the vertebrate Hox bipartite regulatory system is an evolutionary novelty generated by combining ancient long-range regulatory contacts from DNA in the anterior Hox neighborhood with new regulatory inputs from the posterior side.

Favorable genomic environments for cis-regulatory evolution: A novel theoretical framework

Cis-regulatory changes are arguably the primary evolutionary source of animal morphological diversity. With the recent explosion of genome-wide comparisons of the cis-regulatory content in different animal species is now possible to infer general principles underlying enhancer evolution. However, these studies have also revealed numerous discrepancies and paradoxes, suggesting that the mechanistic causes and modes of cis-regulatory evolution are still not well understood and are probably much more complex than generally appreciated. Here, we argue that the mutational mechanisms and genomic regions generating new regulatory activities must comply with the constraints imposed by the molecular properties of cis-regulatory elements (CREs) and the organizational features of long-range chromatin interactions. Accordingly, we propose a new integrative evolutionary framework for cis-regulatory evolution based on two major premises for the origin of novel enhancer activity: (i) an accessible chromatin environment and (ii) compatibility with the 3D structure and interactions of pre-existing CREs. Mechanisms and DNA sequences not fulfilling these premises, will be less likely to have a measurable impact on gene expression and as such, will have a minor contribution to the evolution of gene regulation. Finally, we discuss current comparative cis-regulatory data under the light of this new evolutionary model, and propose that the two most prominent mechanisms for the evolution of cis-regulatory changes are the overprinting of ancestral CREs and the exaptation of transposable elements.

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.

Deep conservation of wrist and digit enhancers in fish

There is no obvious morphological counterpart of the autopod (wrist/ankle and digits) in living fishes. Comparative molecular data may provide insight into understanding both the homology of elements and the evolutionary developmental mechanisms behind the fin to limb transition. In mouse limbs the autopod is built by a "late" phase of Hoxd and Hoxa gene expression, orchestrated by a set of enhancers located at the 5' end of each cluster. Despite a detailed mechanistic understanding of mouse limb development, interpretation of Hox expression patterns and their regulation in fish has spawned multiple hypotheses as to the origin and function of "autopod" enhancers throughout evolution. Using phylogenetic footprinting, epigenetic profiling, and transgenic reporters, we have identified and functionally characterized hoxD and hoxA enhancers in the genomes of zebrafish and the spotted gar, Lepisosteus oculatus, a fish lacking the whole genome duplication of teleosts. Gar and zebrafish "autopod" enhancers drive expression in the distal portion of developing zebrafish pectoral fins, and respond to the same functional cues as their murine orthologs. Moreover, gar enhancers drive reporter gene expression in both the wrist and digits of mouse embryos in patterns that are nearly indistinguishable from their murine counterparts. These functional genomic data support the hypothesis that the distal radials of bony fish are homologous to the wrist and/or digits of tetrapods.

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

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

Quantitative modeling of a gene's expression from its intergenic sequence

Modeling a gene's expression from its intergenic locus and trans-regulatory context is a fundamental goal in computational biology. Owing to the distributed nature of cis-regulatory information and the poorly understood mechanisms that integrate such information, gene locus modeling is a more challenging task than modeling individual enhancers. Here we report the first quantitative model of a gene's expression pattern as a function of its locus. We model the expression readout of a locus in two tiers: 1) combinatorial regulation by transcription factors bound to each enhancer is predicted by a thermodynamics-based model and 2) independent contributions from multiple enhancers are linearly combined to fit the gene expression pattern. The model does not require any prior knowledge about enhancers contributing toward a gene's expression. We demonstrate that the model captures the complex multi-domain expression patterns of anterior-posterior patterning genes in the early Drosophila embryo. Altogether, we model the expression patterns of 27 genes; these include several gap genes, pair-rule genes, and anterior, posterior, trunk, and terminal genes. We find that the model-selected enhancers for each gene overlap strongly with its experimentally characterized enhancers. Our findings also suggest the presence of sequence-segments in the locus that would contribute ectopic expression patterns and hence were "shut down" by the model. We applied our model to identify the transcription factors responsible for forming the stripe boundaries of the studied genes. The resulting network of regulatory interactions exhibits a high level of agreement with known regulatory influences on the target genes. Finally, we analyzed whether and why our assumption of enhancer independence was necessary for the genes we studied. We found a deterioration of expression when binding sites in one enhancer were allowed to influence the readout of another enhancer. Thus, interference between enhancer activities was a possible factor necessitating enhancer independence in our model.

Clustering of tissue-specific sub-TADs accompanies the regulation of HoxA genes in developing limbs

HoxA genes exhibit central roles during development and causal mutations have been found in several human syndromes including limb malformation. Despite their importance, information on how these genes are regulated is lacking. Here, we report on the first identification of bona fide transcriptional enhancers controlling HoxA genes in developing limbs and show that these enhancers are grouped into distinct topological domains at the sub-megabase scale (sub-TADs). We provide evidence that target genes and regulatory elements physically interact with each other through contacts between sub-TADs rather than by the formation of discreet “DNA loops”. Interestingly, there is no obvious relationship between the functional domains of the enhancers within the limb and how they are partitioned among the topological domains, suggesting that sub-TAD formation does not rely on enhancer activity. Moreover, we show that suppressing the transcriptional activity of enhancers does not abrogate their contacts with HoxA genes. Based on these data, we propose a model whereby chromatin architecture defines the functional landscapes of enhancers. From an evolutionary standpoint, our data points to the convergent evolution of HoxA and HoxD regulation in the fin-to-limb transition, one of the major morphological innovations in vertebrates.

Hox genes encode transcription factors with crucial roles during development. These genes are grouped in four different clusters names HoxA, B, C, and D. Mutations in genes of the HoxA and D clusters have been found in several human syndromes, affecting in some cases limb development. Despite their essential role and contrary to the genes of the HoxD cluster, little is known about how the HoxA genes are regulated. Here, we identified a large set of regulatory elements controlling HoxA genes during limb development. By studying spatial chromatin organization at the HoxA region, we found that the regulatory elements are spatially clustered regardless of their activity. Clustering of enhancers define tissue-specific chromatin domains that interact specifically with each other and with active genes in the limb. Our findings give support to the emerging concept that chromatin architecture defines the functional properties of genomes. Additionally, our study suggests a common constraint of the chromatin topology in the evolution of HoxA and HoxD regulation in the emergence of the hand/foot, which is one of the major morphological innovations in vertebrates.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Hox genes regulation in vertebrates

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

The genomic landscape of cohesin-associated chromatin interactions

Cohesin is implicated in establishing tissue-specific DNA loops that target enhancers to promoters, and also localizes to sites bound by the insulator protein CTCF, which blocks enhancer-promoter communication. However, cohesin-associated interactions have not been characterized on a genome-wide scale. Here we performed chromatin interaction analysis with paired-end tag sequencing (ChIA-PET) of the cohesin subunit SMC1A in developing mouse limb. We identified 2264 SMC1A interactions, of which 1491 (65%) involved sites co-occupied by CTCF. SMC1A participates in tissue-specific enhancer-promoter interactions and interactions that demarcate regions of correlated regulatory output. In contrast to previous studies, we also identified interactions between promoters and distal sites that are maintained in multiple tissues but are poised in embryonic stem cells and resolve to tissue-specific activated or repressed chromatin states in the mouse embryo. Our results reveal the diversity of cohesin-associated interactions in the genome and highlight their role in establishing the regulatory architecture of development.

Chromatin organization and global regulation of Hox gene clusters

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

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

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

An integrated holo-enhancer unit defines tissue and gene specificity of the Fgf8 regulatory landscape

Fgf8 encodes a key signaling factor, and its precise regulation is essential for embryo patterning. Here, we identified the regulatory modules that control Fgf8 expression during mammalian embryogenesis. These enhancers are interspersed with unrelated genes along a large region of 220 kb; yet they act on Fgf8 only. Intriguingly, this region also contains additional genuine enhancer activities that are not transformed into gene expression. Using genomic engineering strategies, we showed that these multiple and distinct regulatory modules act as a coherent unit and influence genes depending on their position rather than on their promoter sequence. These findings highlight how the structure of a locus regulates the autonomous intrinsic activities of the regulatory elements it contains and contributes to their tissue and target specificities. We discuss the implications of such regulatory systems regarding the evolution of gene expression and the impact of human genomic structural variations.

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.

Making limbs from fins

Our limbs evolved from fish fins as vertebrates colonized aquatic shallows and land. If we understood the genetics underlying this transition, could we build a limb on a fish? In this issue of Developmental Cell, Freitas et al. (2012) show that boosting Hoxd13 expression can bring zebrafish one step closer to terra firma.

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.

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.