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

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

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

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

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

Shaping gene expression and its evolution by chromatin architecture and enhancer activity

Transcriptional regulation plays a pivotal role in orchestrating the intricate genetic programs governing embryonic development. The expression of developmental genes relies on the combined activity of several cis-regulatory elements (CREs), such as enhancers and silencers, which can be located at long linear distances from the genes that they regulate and that interact with them through establishment of chromatin loops. Mutations affecting their activity or interaction with their target genes can lead to developmental disorders and are thought to have importantly contributed to the evolution of the animal body plan. The income of next-generation-sequencing approaches has allowed identifying over a million of sequences with putative regulatory potential in the human genome. Characterizing their function and establishing gene-CREs maps is essential to decode the logic governing developmental gene expression and is one of the major challenges of the post-genomic era. Chromatin 3D organization plays an essential role in determining how CREs specifically contact their target genes while avoiding deleterious off-target interactions. Our understanding of these aspects has greatly advanced with the income of chromatin conformation capture techniques and fluorescence microscopy approaches to visualize the organization of DNA elements in the nucleus. Here we will summarize relevant aspects of how the interplay between CRE activity and chromatin 3D organization regulates developmental gene expression and how it relates to pathological conditions and the evolution of animal body plan.

Regulatory changes associated with the head to trunk developmental transition

BACKGROUND

Development of vertebrate embryos is characterized by early formation of the anterior tissues followed by the sequential extension of the axis at their posterior end to build the trunk and tail structures, first by the activity of the primitive streak and then of the tail bud. Embryological, molecular and genetic data indicate that head and trunk development are significantly different, suggesting that the transition into the trunk formation stage involves major changes in regulatory gene networks.

RESULTS

We explored those regulatory changes by generating differential interaction networks and chromatin accessibility profiles from the posterior epiblast region of mouse embryos at embryonic day (E)7.5 and E8.5. We observed changes in various cell processes, including several signaling pathways, ubiquitination machinery, ion dynamics and metabolic processes involving lipids that could contribute to the functional switch in the progenitor region of the embryo. We further explored the functional impact of changes observed in Wnt signaling associated processes, revealing a switch in the functional relevance of Wnt molecule palmitoleoylation, essential during gastrulation but becoming differentially required for the control of axial extension and progenitor differentiation processes during trunk formation. We also found substantial changes in chromatin accessibility at the two developmental stages, mostly mapping to intergenic regions and presenting differential footprinting profiles to several key transcription factors, indicating a significant switch in the regulatory elements controlling head or trunk development. Those chromatin changes are largely independent of retinoic acid, despite the key role of this factor in the transition to trunk development. We also tested the functional relevance of potential enhancers identified in the accessibility assays that reproduced the expression profiles of genes involved in the transition. Deletion of these regions by genome editing had limited effect on the expression of those genes, suggesting the existence of redundant enhancers that guarantee robust expression patterns.

CONCLUSIONS

This work provides a global view of the regulatory changes controlling the switch into the axial extension phase of vertebrate embryonic development. It also revealed mechanisms by which the cellular context influences the activity of regulatory factors, channeling them to implement one of several possible biological outputs.

Three-dimensional genome structure and function

Linear DNA undergoes a series of compression and folding events, forming various three-dimensional (3D) structural units in mammalian cells, including chromosomal territory, compartment, topologically associating domain, and chromatin loop. These structures play crucial roles in regulating gene expression, cell differentiation, and disease progression. Deciphering the principles underlying 3D genome folding and the molecular mechanisms governing cell fate determination remains a challenge. With advancements in high-throughput sequencing and imaging techniques, the hierarchical organization and functional roles of higher-order chromatin structures have been gradually illuminated. This review systematically discussed the structural hierarchy of the 3D genome, the effects and mechanisms of cis-regulatory elements interaction in the 3D genome for regulating spatiotemporally specific gene expression, the roles and mechanisms of dynamic changes in 3D chromatin conformation during embryonic development, and the pathological mechanisms of diseases such as congenital developmental abnormalities and cancer, which are attributed to alterations in 3D genome organization and aberrations in key structural proteins. Finally, prospects were made for the research about 3D genome structure, function, and genetic intervention, and the roles in disease development, prevention, and treatment, which may offer some clues for precise diagnosis and treatment of related diseases.

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.

Context-dependent enhancer function revealed by targeted inter-TAD relocation

The expression of some genes depends on large, adjacent regions of the genome that contain multiple enhancers. These regulatory landscapes frequently align with Topologically Associating Domains (TADs), where they integrate the function of multiple similar enhancers to produce a global, TAD-specific regulation. We asked if an individual enhancer could overcome the influence of one of these landscapes, to drive gene transcription. To test this, we transferred an enhancer from its native location, into a nearby TAD with a related yet different functional specificity. We used the biphasic regulation of Hoxd genes during limb development as a paradigm. These genes are first activated in proximal limb cells by enhancers located in one TAD, which is then silenced when the neighboring TAD activates its enhancers in distal limb cells. We transferred a distal limb enhancer into the proximal limb TAD and found that its new context suppresses its normal distal specificity, even though it is bound by HOX13 transcription factors, which are responsible for the distal activity. This activity can be rescued only when a large portion of the surrounding environment is removed. These results indicate that, at least in some cases, the functioning of enhancer elements is subordinated to the host chromatin context, which can exert a dominant control over its activity.

Here the authors show that a strong enhancer sequence can be controlled by the chromatin environment provided by a topologically associated domain (TAD) located nearby. An enhancer relocated by homologous recombination takes all the hallmarks of its new neighboring enhancers located in the recipient TAD.

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.

Essay the (unusual) heuristic value of Hox gene clusters; a matter of time?

Ever since their first report in 1984, Antennapedia-type homeobox (Hox) genes have been involved in such a series of interesting observations, in particular due to their conserved clustered organization between vertebrates and arthropods, that one may legitimately wonder about the origin of this heuristic value. In this essay, I first consider different examples where Hox gene clusters have been instrumental in providing conceptual advances, taken from various fields of research and mostly involving vertebrate embryos. These examples touch upon our understanding of genomic evolution, the revisiting of 19th century views on the relationships between development and evolution and the building of a new framework to understand long-range and pleiotropic gene regulation during development. I then discuss whether the high value of the Hox gene family, when considered as an epistemic object, is related to its clustered structure (and the absence thereof in some animal species) and, if so, what is it in such particular genetic oddities that made them so generous in providing the scientific community with interesting information.

baredSC: Bayesian approach to retrieve expression distribution of single-cell data

Background

The number of studies using single-cell RNA sequencing (scRNA-seq) is constantly growing. This powerful technique provides a sampling of the whole transcriptome of a cell. However, sparsity of the data can be a major hurdle when studying the distribution of the expression of a specific gene or the correlation between the expressions of two genes.

Results

We show that the main technical noise associated with these scRNA-seq experiments is due to the sampling, i.e., Poisson noise. We present a new tool named baredSC, for Bayesian Approach to Retrieve Expression Distribution of Single-Cell data, which infers the intrinsic expression distribution in scRNA-seq data using a Gaussian mixture model. baredSC can be used to obtain the distribution in one dimension for individual genes and in two dimensions for pairs of genes, in particular to estimate the correlation in the two genes’ expressions. We apply baredSC to simulated scRNA-seq data and show that the algorithm is able to uncover the expression distribution used to simulate the data, even in multi-modal cases with very sparse data. We also apply baredSC to two real biological data sets. First, we use it to measure the anti-correlation between Hoxd13 and Hoxa11, two genes with known genetic interaction in embryonic limb. Then, we study the expression of Pitx1 in embryonic hindlimb, for which a trimodal distribution has been identified through flow cytometry. While other methods to analyze scRNA-seq are too sensitive to sampling noise, baredSC reveals this trimodal distribution.

Conclusion

baredSC is a powerful tool which aims at retrieving the expression distribution of few genes of interest from scRNA-seq data.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04507-8.

Rethinking positional information and digit identity: The role of late interdigit signaling

Seminal work from John Fallon's lab has illuminated how digit identity determination involves ongoing late regulation and occurs progressively during phalanx formation. Complementary genetic analyses in mice and several papers in this special issue have begun to flesh out how interdigit signaling accomplishes this, but major questions remain unaddressed, including how uncommitted progenitors from which phalanges arise are maintained, and what factors set limits on digit extension and phalanx number, particularly in mammals. This review summarizes what has been learned in the two decades since control of digit identity by late interdigit signals was first identified and what remains poorly understood, and will hopefully spark renewed interest in a process that is critical to evolutionary limb adaptations but nevertheless remains enigmatic.

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.

Chromatin Conformation in Development and Disease

Chromatin domains and loops are important elements of chromatin structure and dynamics, but much remains to be learned about their exact biological role and nature. Topological associated domains and functional loops are key to gene expression and hold the answer to many questions regarding developmental decisions and diseases. Here, we discuss new findings, which have linked chromatin conformation with development, differentiation and diseases and hypothesized on various models while integrating all recent findings on how chromatin architecture affects gene expression during development, evolution and disease.

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.

Hoxd13/Bmp2-mediated mechanism involved in zebrafish finfold design

The overexpression of hoxd13a during zebrafish fin development causes distal endochondral expansion and simultaneous reduction of the finfold, mimicking the major events thought to have happened during the fin-to-limb transition in Vertebrates. We investigated the effect of hoxd13a overexpression on putative downstream targets and found it to cause downregulation of proximal fin identity markers (meis1 and emx2) and upregulation of genes involved in skeletogenesis/patterning (fbn1, dacha) and AER/Finfold maintenance (bmps). We then show that bmp2b overexpression leads to finfold reduction, recapitulating the phenotype observed in hoxd13a-overexpressing fins. In addition, we show that during the development of the long finfold in leo^t1/lof^dt1 mutants, hoxd13a and bmp2b are downregulated. Our results suggest that modulation of the transcription factor Hoxd13 during evolution may have been involved in finfold reduction through regulation of the Bmp signalling that then activated apoptotic mechanisms impairing finfold elongation.

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.

Transcriptional Trajectories in Mouse Limb Buds Reveal the Transition from Anterior-Posterior to Proximal-Distal Patterning at Early Limb Bud Stage

Limb patterning relies in large part on the function of the Hox family of developmental genes. While the differential expression of Hox genes shifts from the anterior-posterior (A-P) to the proximal-distal (P-D) axis around embryonic day 11 (E11), whether this shift coincides with a more global change of A-P to P-D patterning program remains unclear. By performing and analyzing the transcriptome of the developing limb bud from E10.5 to E12.5, at single-cell resolution, we have uncovered transcriptional trajectories that revealed a general switch from A-P to P-D genetic program between E10.5 and E11.5. Interestingly, all the transcriptional trajectories at E10.5 end with cells expressing either proximal or distal markers suggesting a progressive acquisition of P-D identity. Moreover, we identified three categories of genes expressed in the distal limb mesenchyme characterized by distinct temporal expression dynamics. Among these are Hoxa13 and Hoxd13 (Hox13 hereafter), which start to be expressed around E10.5, and importantly the binding of the HOX13 factors was observed within or in the neighborhood of several of the distal limb genes. Our data are consistent with previous evidence suggesting that the transition from the early/proximal to the late/distal transcriptome of the limb mesenchyme largely relies on HOX13 function. Based on these results and the evidence that HOX13 factors restrict Hoxa11 expression to the proximal limb, in progenitor cells of the zeugopod, we propose that HOX13 act as a key determinant of P-D patterning.

Establishing the pattern of the vertebrate limb

The vertebrate limb continues to serve as an influential model of growth, morphogenesis and pattern formation. With this Review, we aim to give an up-to-date picture of how a population of undifferentiated cells develops into the complex pattern of the limb. Focussing largely on mouse and chick studies, we concentrate on the positioning of the limbs, the formation of the limb bud, the establishment of the principal limb axes, the specification of pattern, the integration of pattern formation with growth and the determination of digit number. We also discuss the important, but little understood, topic of how gene expression is interpreted into morphology.

Sarcopterygian fin ontogeny elucidates the origin of hands with digits

How the hand and digits originated from fish fins during the Devonian fin-to-limb transition remains unsolved. Controversy in this conundrum stems from the scarcity of ontogenetic data from extant lobe-finned fishes. We report the patterning of an autopod-like domain by hoxa13 during fin development of the Australian lungfish, the most closely related extant fish relative of tetrapods. Differences from tetrapod limbs include the absence of digit-specific expansion of hoxd13 and hand2 and distal limitation of alx4 and pax9, which potentially evolved through an enhanced response to shh signaling in limbs. These developmental patterns indicate that the digit program originated in postaxial fin radials and later expanded anteriorly inside of a preexisting autopod-like domain during the evolution of limbs. Our findings provide a genetic framework for the transition of fins into limbs that supports the significance of classical models proposing a bending of the tetrapod metapterygial axis.

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 sterlet sturgeon genome sequence and the mechanisms of segmental rediploidization

Sturgeons seem to be frozen in time. The archaic characteristics of this ancient fish lineage place it in a key phylogenetic position at the base of the ~30,000 modern teleost fish species. Moreover, sturgeons are notoriously polyploid, providing unique opportunities to investigate the evolution of polyploid genomes. We assembled a high-quality chromosome-level reference genome for the sterlet, Acipenser ruthenus. Our analysis revealed a very low protein evolution rate that is at least as slow as in other deep branches of the vertebrate tree, such as that of the coelacanth. We uncovered a whole-genome duplication that occurred in the Jurassic, early in the evolution of the entire sturgeon lineage. Following this polyploidization, the rediploidization of the genome included the loss of whole chromosomes in a segmental deduplication process. While known adaptive processes helped conserve a high degree of structural and functional tetraploidy over more than 180 million years, the reduction of redundancy of the polyploid genome seems to have been remarkably random.

A genome assembly of the sterlet, Acipenser ruthenus, reveals a whole-genome duplication early in the evolution of the entire sturgeon lineage and provides details about the rediploidization of the genome.

The formation of the thumb requires direct modulation of Gli3 transcription by Hoxa13

In the tetrapod limb, the digits (fingers or toes) are the elements most subject to morphological diversification in response to functional adaptations. However, despite their functional importance, the mechanisms controlling digit morphology remain poorly understood. Here we have focused on understanding the special morphology of the thumb (digit 1), the acquisition of which was an important adaptation of the human hand. To this end, we have studied the limbs of the Hoxa13 mouse mutant that specifically fail to form digit 1. We show that, consistent with the role of Hoxa13 in Hoxd transcriptional regulation, the expression of Hoxd13 in Hoxa13 mutant limbs does not extend into the presumptive digit 1 territory, which is therefore devoid of distal Hox transcripts, a circumstance that can explain its agenesis. The loss of Hoxd13 expression, exclusively in digit 1 territory, correlates with increased Gli3 repressor activity, a Hoxd negative regulator, resulting from increased Gli3 transcription that, in turn, is due to the release from the negative modulation exerted by Hox13 paralogs on Gli3 regulatory sequences. Our results indicate that Hoxa13 acts hierarchically to initiate the formation of digit 1 by reducing Gli3 transcription and by enabling expansion of the 5'Hoxd second expression phase, thereby establishing anterior-posterior asymmetry in the handplate. Our work uncovers a mutual antagonism between Gli3 and Hox13 paralogs that has important implications for Hox and Gli3 gene regulation in the context of development and evolution.

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.

Pax3 cooperates with Ldb1 to direct local chromosome architecture during myogenic lineage specification

Chromatin looping allows enhancer-bound regulatory factors to influence transcription. Large domains, referred to as topologically associated domains, participate in genome organization. However, the mechanisms underlining interactions within these domains, which control gene expression, are not fully understood. Here we report that activation of embryonic myogenesis is associated with establishment of long-range chromatin interactions centered on Pax3-bound loci. Using mass spectrometry and genomic studies, we identify the ubiquitously expressed LIM-domain binding protein 1 (Ldb1) as the mediator of looping interactions at a subset of Pax3 binding sites. Ldb1 is recruited to Pax3-bound elements independently of CTCF-Cohesin, and is necessary for efficient deposition of H3K4me1 at these sites and chromatin looping. When Ldb1 is deleted in Pax3-expressing cells in vivo, specification of migratory myogenic progenitors is severely impaired. These results highlight Ldb1 requirement for Pax3 myogenic activity and demonstrate how transcription factors can promote formation of sub-topologically associated domain interactions involved in lineage specification.

'Building a perfect body': control of vertebrate organogenesis by PBX-dependent regulatory networks

Pbx genes encode transcription factors that belong to the TALE (three-amino-acid loop extension) superclass of homeodomain proteins. We have witnessed a surge in information about the roles of this gene family as leading actors in the transcriptional control of development. PBX proteins represent a clear example of how transcription factors can regulate developmental processes by combinatorial properties, acting within multimeric complexes to implement activation or repression of transcription depending on their interaction partners. Here, we revisit long-emphasized functions of PBX transcription factors as cofactors for HOX proteins, major architects of the body plan. We further discuss new knowledge on roles of PBX proteins in different developmental contexts as upstream regulators of Hox genes-as factors that interact with non-HOX proteins and can work independently of HOX-as well as potential pioneer factors. Committed to building a perfect body, PBX proteins govern regulatory networks that direct essential morphogenetic processes and organogenesis in vertebrate development. Perturbations of PBX-dependent networks can cause human congenital disease and cancer.

Hoxa13 regulates expression of common Hox target genes involved in cartilage development to coordinate the expansion of the autopodal anlage

To elucidate the role of Hox genes in limb cartilage development, we identified the target genes of HOXA11 and HOXA13 by ChIP‐Seq. The ChIP DNA fragment contained evolutionarily conserved sequences and multiple highly conserved HOX binding sites. A substantial portion of the HOXA11 ChIP fragment overlapped with the HOXA13 ChIP fragment indicating that both factors share common targets. Deletion of the target regions neighboring Bmp2 or Tshz2 reduced their expression in the autopod suggesting that they function as the limb bud‐specific enhancers. We identified the Hox downstream genes as exhibiting expression changes in the Hoxa13 knock out (KO) and Hoxd11‐13 deletion double mutant (Hox13 dKO) autopod by Genechip analysis. The Hox downstream genes neighboring the ChIP fragment were defined as the direct targets of Hox. We analyzed the spatial expression pattern of the Hox target genes that encode two different categories of transcription factors during autopod development and Hox13dKO limb bud. (a) Bcl11a, encoding a repressor of cartilage differentiation, was expressed in the E11.5 autopod and was substantially reduced in the Hox13dKO. (b) The transcription factors Aff3, Bnc2, Nfib and Runx1t1 were expressed in the zeugopodal cartilage but not in the autopod due to the repressive or relatively weak transcriptional activity of Hox13 at E11.5. Interestingly, the expression of these genes was later observed in the autopodal cartilage at E12.5. These results indicate that Hox13 transiently suspends the cartilage differentiation in the autopodal anlage via multiple pathways until establishing the paddle‐shaped structure required to generate five digits.

Anatomic Origin of Osteochondrogenic Progenitors Impacts Sensitivity to EWS-FLI1-Induced Transformation

Ewing sarcomas predominantly arise in pelvic and stylopod bones (i.e., femur and humerus), likely as a consequence of EWS-FLI1 oncogene-induced transformation of mesenchymal stem/progenitor cells (MSCs). MSCs located in the embryonic superficial zone cells (eSZ) of limbs express anatomically distinct posterior Hox genes. Significantly, high expression of posterior HOXD genes, especially HOXD13, is a hallmark of Ewing sarcoma. These data drove our hypothesis that Hox genes in posterior skeleton MSCs contribute to Ewing sarcoma tumorigenesis. We isolated eSZ cells from stylopod and zeugopod (i.e., tibia/fibula, radius/ulna) bones, from wild-type and Hoxd13 mutant embryos, and tested the impact of EWS-FLI1 transduction on cell proliferation, gene expression, and tumorigenicity. Our data demonstrate that both stylopod and zeugopod eSZ cells tolerate EWS-FLI1 but that stylopod eSZ cells are relatively more susceptible, demonstrating changes in proliferation and gene expression consistent with initiation of malignant transformation. Significantly, loss of Hoxd13 had no impact, showing that it is dispensable for the initiation of EWS-FLI1-induced transformation in mouse MSCs. These findings show that MSCs from anatomically distinct sites are differentially susceptible to EWS-FLI1-induced transformation, supporting the premise that the dominant presentation of Ewing sarcoma in pelvic and stylopod bones is attributable to anatomically-defined differences in MSCs.

Reassessing the Role of Hox Genes during Vertebrate Development and Evolution

Since their discovery Hox genes have been at the core of the established models explaining the development and evolution of the vertebrate body plan as well as its paired appendages. Recent work brought new light to their role in the patterning processes along the main body axis. These studies show that Hox genes do not control the basic layout of the vertebrate body plan but carry out region-specific patterning instructions loaded on the derivatives of axial progenitors by Hox-independent processes. Furthermore, the finding that Hox clusters are embedded in functional chromatin domains, which critically impacts their expression, has significantly altered our understanding of the mechanisms of Hox gene regulation. This new conceptual framework has broadened our understanding of both limb development and the evolution of vertebrate paired appendages.

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.

Dynamic 3D chromatin architecture contributes to enhancer specificity and limb morphogenesis

The regulatory specificity of enhancers and their interaction with gene promoters is thought to be controlled by their sequence and the binding of transcription factors. By studying Pitx1, a regulator of hindlimb development, we show that dynamic changes in chromatin conformation can restrict the activity of enhancers. Inconsistent with its hindlimb-restricted expression, Pitx1 is controlled by an enhancer (Pen) that shows activity in forelimbs and hindlimbs. By Capture Hi-C and three-dimensional modeling of the locus, we demonstrate that forelimbs and hindlimbs have fundamentally different chromatin configurations, whereby Pen and Pitx1 interact in hindlimbs and are physically separated in forelimbs. Structural variants can convert the inactive into the active conformation, thereby inducing Pitx1 misexpression in forelimbs, causing partial arm-to-leg transformation in mice and humans. Thus, tissue-specific three-dimensional chromatin conformation can contribute to enhancer activity and specificity in vivo and its disturbance can result in gene misexpression and disease.

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.

Control of directionality of chromatin folding for the inter- and intra-domain contacts at the Tfap2c-Bmp7 locus

BACKGROUND

Contact domains of chromatin serve as a fundamental unit to regulate action of enhancers for target genes. Looping between a pair of CCCTC-binding factor (CTCF)-binding sites in convergent orientations underlies the formation of contact domains, while those in divergent orientations establish domain boundaries. However, every CTCF site is not necessarily engaged in loop or boundary structures, leaving functions of CTCF in varied genomic contexts still elusive. The locus containing Tfap2c and Bmp7 encompasses two contact domains separated by a region between the two genes, termed transition zone (TZ), characterized by two arrays of CTCF sites in divergent configuration. In this study, we created deletion and inversion alleles of these and other regions across the locus and investigated how they impinge on the conformation.

RESULTS

Deletion of the whole two CTCF arrays with the CRISPR/Cas9 system resulted in impairment of blocking of chromatin contacts by the TZ, as assessed by the circular chromatin conformation capture assay (4C-seq). Deletion and inversion of either of the two arrays similarly, but less pronouncedly, led to reduction in the blocking activity. Thus, the divergent configuration provides the TZ with the strong boundary activity. Uniquely, we show the TZ harbors a 50-kb region within one of the two arrays that contacts broadly with the both flanking intervals, regardless of the presence or orientation of the other CTCF array. Further, we show the boundary CTCF array has little impact on intra-domain folding; instead, locally associating CTCF sites greatly affect it.

CONCLUSIONS

Our results show that the TZ not only separates the two domains, but also bears a wide interval that shows isotropic behavior of chromatin folding, indicating a potentially complex nature of actual boundaries in the genome. We also show that CTCF-binding sites inside a domain greatly contribute to the intra-domain folding of chromatin. Thus, the study reveals diverse and context-dependent roles of CTCF in organizing chromatin conformation at different levels.

Problems in Fish-to-Tetrapod Transition: Genetic Expeditions Into Old Specimens

The fish-to-tetrapod transition is one of the fundamental problems in evolutionary biology. A significant amount of paleontological data has revealed the morphological trajectories of skeletons, such as those of the skull, vertebrae, and appendages in vertebrate history. Shifts in bone differentiation, from dermal to endochondral bones, are key to explaining skeletal transformations during the transition from water to land. However, the genetic underpinnings underlying the evolution of dermal and endochondral bones are largely missing. Recent genetic approaches utilizing model organisms—zebrafish, frogs, chickens, and mice—reveal the molecular mechanisms underlying vertebrate skeletal development and provide new insights for how the skeletal system has evolved. Currently, our experimental horizons to test evolutionary hypotheses are being expanded to non-model organisms with state-of-the-art techniques in molecular biology and imaging. An integration of functional genomics, developmental genetics, and high-resolution CT scanning into evolutionary inquiries allows us to reevaluate our understanding of old specimens. Here, we summarize the current perspectives in genetic programs underlying the development and evolution of the dermal skull roof, shoulder girdle, and appendages. The ratio shifts of dermal and endochondral bones, and its underlying mechanisms, during the fish-to-tetrapod transition are particularly emphasized. Recent studies have suggested the novel cell origins of dermal bones, and the interchangeability between dermal and endochondral bones, obscuring the ontogenetic distinction of these two types of bones. Assimilation of ontogenetic knowledge of dermal and endochondral bones from different structures demands revisions of the prevalent consensus in the evolutionary mechanisms of vertebrate skeletal shifts.

Regulatory integration of Hox factor activity with T-box factors in limb development

In tetrapods, Tbx4, Tbx5 and Hox cluster genes are crucial for forelimb and hindlimb development and mutations in these genes are responsible for congenital limb defects. The molecular basis of their integrated mechanisms of action in the context of limb development remains poorly understood. We studied Tbx4 and Hoxc10 owing to their overlapping loss-of-function phenotypes and colocalized expression in mouse hindlimb buds. We report an extensive overlap between Tbx4 and Hoxc10 genome occupancy and their putative target genes. Tbx4 and Hoxc10 interact directly with each other, have the ability to bind to a previously unrecognized T-box-Hox composite DNA motif and show synergistic activity when acting on reporter genes. Pitx1, the master regulator for hindlimb specification, also shows extensive genomic colocalization with Tbx4 and Hoxc10. Genome occupancy by Tbx4 in hindlimb buds is similar to Tbx5 occupancy in forelimbs. By contrast, another Hox factor, Hoxd13, also interacts with Tbx4/Tbx5 but antagonizes Tbx4/Tbx5-dependent transcriptional activity. Collectively, the modulation of Tbx-dependent activity by Hox factors acting on common DNA targets may integrate different developmental processes for the balanced formation of proportionate limbs.

Fin-fold development in paddlefish and catshark and implications for the evolution of the autopod

The evolutionary origin of the autopod involved a loss of the fin-fold and associated dermal skeleton with a concomitant elaboration of the distal endoskeleton to form a wrist and digits. Developmental studies, primarily from teleosts and amniotes, suggest a model for appendage evolution in which a delay in the AER-to-fin-fold conversion fuelled endoskeletal expansion by prolonging the function of AER-mediated regulatory networks. Here, we characterize aspects of paired fin development in the paddlefish Polyodon spathula (a non-teleost actinopterygian) and catshark Scyliorhinus canicula (chondrichthyan) to explore aspects of this model in a broader phylogenetic context. Our data demonstrate that in basal gnathostomes, the autopod marker HoxA13 co-localizes with the dermoskeleton component And1 to mark the position of the fin-fold, supporting recent work demonstrating a role for HoxA13 in zebrafish fin ray development. Additionally, we show that in paddlefish, the proximal fin and fin-fold mesenchyme share a common mesodermal origin, and that components of the Shh/LIM/Gremlin/Fgf transcriptional network critical to limb bud outgrowth and patterning are expressed in the fin-fold with a profile similar to that of tetrapods. Together these data draw contrast with hypotheses of AER heterochrony and suggest that limb-specific morphologies arose through evolutionary changes in the differentiation outcome of conserved early distal patterning compartments.

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.

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.

The three-dimensional genome: regulating gene expression during pluripotency and development

The precise expression of genes in time and space during embryogenesis is largely influenced by communication between enhancers and promoters, which is propagated and governed by the physical proximity of these elements in the nucleus. Here, we review how chromatin domains organize the genome by guiding enhancers to their target genes thereby preventing non-specific interactions with other neighboring regions. We also discuss the dynamics of chromatin interactions between enhancers and promoters, as well as the consequent changes in gene expression, that occur in pluripotent cells and during development. Finally, we evaluate how genomic changes such as deletions, inversions and duplications affect 3D chromatin configuration overall and lead to ectopic enhancer-promoter contacts, and thus gene misexpression, which can contribute to abnormal development and disease.

Fins into limbs: Recent insights from sarcopterygian fish

Limbs with digits evolved as sarcopterygian fish transitioned to a terrestrial life, giving rise to modern tetrapods. Since the Devonian, most of the sarcopterygian fish diversity became extinct, with the only surviving representatives being two coelacanth and six lungfish species. As the sister group of tetrapods, sarcopterygian fish constitute the ideal models to address questions regarding the transition of vertebrates from water to land. However, distantly related yet experimentally amenable teleost fish species have instead become the organisms of choice for geneticists and developmental biologists. Comparative studies using teleosts, such as zebrafish, have greatly enriched our knowledge on the genetic and regulatory mechanisms underlying fin development. However, their highly derived fin anatomy provides limited insights on the origin of limbs and the teleost-specific whole-genome duplication represents a further complication to comparisons of gene function and regulation. In recent years, refined morphological and behavioral studies and access to lungfish embryos and availability of genetic resources have underscored the value of nontetrapod sarcopterygians as model organisms. Here we review recent studies using lungfishes and coelacanths that shed light on various aspects of the evolution of the tetrapod limb, including the origin of the tetrapod pelvis, limb musculature, digits, locomotion, and regenerative capacity.

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.

On the Formation of Digits and Joints during Limb Development

Critical steps in forming the vertebrate limb include the positioning of digits and the positioning of joints within each digit. Recent studies have proposed that the iterative series of digits is established by a Turing-like mechanism generating stripes of chondrogenic domains. However, re-examination of available data suggest that digits are actually patterned as evenly spaced spots, not stripes, which then elongate into rod-shaped digit rays by incorporating new cells at their tips. Moreover, extension of the digit rays and the patterning of the joints occur simultaneously at the distal tip, implying that an integrated model is required to fully understand these processes.

HOX13 proteins: the molecular switcher in Hoxd bimodal regulation

Ros discusses the report by Beccari et al. on the mechanism of HOX13 regulation of HoxD gene expression in the developing tetrapod limb to create a region of low Hox expression that ultimately gives rise to the wrist as well as the evolutionary implications of this finding.

The striking correlation between the genomic arrangement of Hox genes and their temporal and spatial pattern of expression during embryonic development has been a source of fascination since its discovery. This correspondence has been used as a privileged example in the investigation of the connection between genomic architecture and function. In this issue of Genes & Development, Beccari and colleagues (pp. 1172–1186) make a big step forward in understanding Hox gene regulation during limb development by showing the pivotal role of HOXA13 and HOXD13 proteins in the transition from a proximal to a distal type of Hoxd transcriptional regulation.