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

Genome organization controls transcriptional dynamics during development

Past studies offer contradictory claims for the role of genome organization in the regulation of gene activity. Here, we show through high-resolution chromosome conformation analysis that the Drosophila genome is organized by two independent classes of regulatory sequences, tethering elements and insulators. Quantitative live imaging and targeted genome editing demonstrate that this two-tiered organization is critical for the precise temporal dynamics of Hox gene transcription during development. Tethering elements mediate long-range enhancer-promoter interactions and foster fast activation kinetics. Conversely, the boundaries of topologically associating domains (TADs) prevent spurious interactions with enhancers and silencers located in neighboring TADs. These two levels of genome organization operate independently of one another to ensure precision of transcriptional dynamics and the reliability of complex patterning processes.

Epithelial Cells in 2D and 3D Cultures Exhibit Large Differences in Higher-order Genomic Interactions

Recent studies have characterized the genomic structures of many eukaryotic cells, often focusing on their relation to gene expression. However, these studies have largely investigated cells grown in 2D cultures, although the transcriptomes of 3D-cultured cells are generally closer to their in vivo phenotypes. To examine the effects of spatial constraints on chromosome conformation, we investigated the genomic architecture of mouse hepatocytes grown in 2D and 3D cultures using in situ Hi-C. Our results reveal significant differences in higher-order genomic interactions, notably in compartment identity and strength as well as in topologically associating domain (TAD)-TAD interactions, but only minor differences are found at the TAD level. Our RNA-seq analysis reveals an up-regulated expression of genes involved in physiological hepatocyte functions in the 3D-cultured cells. These genes are associated with a subset of structural changes, suggesting that differences in genomic structure are critically important for transcriptional regulation. However, there are also many structural differences that are not directly associated with changes in gene expression, whose cause remains to be determined. Overall, our results indicate that growth in 3D significantly alters higher-order genomic interactions, which may be consequential for a subset of genes that are important for the physiological functioning of the cell.

Repurposing of the enhancer-promoter communication underlies the compensation of Mesp2 by Mesp1

Organisms are inherently equipped with buffering systems against genetic perturbations. Genetic compensation, the compensatory response by upregulating another gene or genes, is one such buffering mechanism. Recently, a well-conserved compensatory mechanism was proposed: transcriptional adaptation of homologs under the nonsense-mediated mRNA decay pathways. However, this model cannot explain the onset of all compensatory events. We report a novel genetic compensation mechanism operating over the Mesp gene locus. Mesp1 and Mesp2 are paralogs located adjacently in the genome. Mesp2 loss is partially rescued by Mesp1 upregulation in the presomitic mesoderm (PSM). Using a cultured PSM induction system, we reproduced the compensatory response in vitro and found that the Mesp2-enhancer is required to promote Mesp1. We revealed that the Mesp2-enhancer directly interacts with the Mesp1 promoter, thereby upregulating Mesp1 expression upon the loss of Mesp2. Of note, this interaction is established by genomic arrangement upon PSM development independently of Mesp2 disruption. We propose that the repurposing of this established enhancer-promoter communication is the mechanism underlying this compensatory response for the upregulation of the adjacent gene.

Genetic compensation, the compensatory response by upregulating another gene or genes, is one of the inherent mechanisms against gene disruption to confer cellular fitness. However, the regulatory mechanisms are largely unknown. Nonsense-mediated mutant mRNA degradation was recently proposed as a conserved mechanism across species to upregulate homologous genes to compensate for a disrupted gene, but this cannot explain compensation events with no mutant mRNA. This study investigated the compensation mechanism operating over adjacent paralogs, Mesp1 and Mesp2, in the genome. Mesp genes encode essential transcription factors in the presomitic mesoderm for development. In general, an enhancer is considered to activate a target gene when it physically interacts with the target. The communication of the Mesp2-enhancer with the Mesp1 promoter is established upon differentiation of the presomitic mesoderm, but this communication activates Mesp1 only when Mesp2 is disrupted, leading to compensation. We revealed a novel compensation mechanism depending on the repurposing of this enhancer-promoter communication by gene disruption. Our study also provides new insight into transcriptional regulation by providing the concept that an enhancer changes its target even among its physically interacting genes in a context-dependent manner.

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.

Detection of Allele-Specific 3D Chromatin Interactions Using High-Resolution In-Nucleus 4C-seq

Chromosome conformation capture techniques are a set of methods used to determine 3D genome organization through the capture and identification of physical contacts between pairs of genomic loci. Among them, 4C-seq (circular chromosome conformation capture coupled to high-throughput sequencing) allows for the identification and quantification of the sequences interacting with a preselected locus of interest. 4C-seq has been widely used in the literature, mainly to study chromatin loops between enhancers and promoters or between CTCF binding sites and to identify chromatin domain boundaries. As 3D-contacts may be established in an allele-specific manner, we describe an up-to-date allele-specific 4C-seq protocol, starting from the selection of allele-specific viewpoints to Illumina sequencing. This protocol has mainly been optimized for cultured mammalian cells, but can be adapted for other cell types with relatively minor changes in initial steps.

Spatial organization of transcribing loci during early genome activation in Drosophila

The early Drosophila embryo provides unique experimental advantages for addressing fundamental questions of gene regulation at multiple levels of organization, from individual gene loci to the entire genome. Using 1.5-h-old Drosophila embryos undergoing the first wave of genome activation,¹ we detected ∼110 discrete "speckles" of RNA polymerase II (RNA Pol II) per nucleus, two of which were larger and localized to the histone locus bodies (HLBs).^2,3 In the absence of the primary driver of Drosophila genome activation, the pioneer factor Zelda (Zld),^1,4,5 70% fewer speckles were present; however, the HLBs tended to be larger than wild-type (WT) HLBs, indicating that RNA Pol II accumulates at the HLBs in the absence of robust early-gene transcription. We observed a uniform distribution of distances between active genes in the nuclei of both WT and zld mutant embryos, indicating that early co-regulated genes do not cluster into nuclear sub-domains. However, in instances whereby transcribing genes did come into close 3D proximity (within 400 nm), they were found to have distinct RNA Pol II speckles. In contrast to the emerging model whereby active genes are clustered to facilitate co-regulation and sharing of transcriptional resources, our data support an "individualist" model of gene control at early genome activation in Drosophila. This model is in contrast to a "collectivist" model, where active genes are spatially clustered and share transcriptional resources, motivating rigorous tests of both models in other experimental systems.

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.

CTCF knockout in zebrafish induces alterations in regulatory landscapes and developmental gene expression

Coordinated chromatin interactions between enhancers and promoters are critical for gene regulation. The architectural protein CTCF mediates chromatin looping and is enriched at the boundaries of topologically associating domains (TADs), which are sub-megabase chromatin structures. In vitro CTCF depletion leads to a loss of TADs but has only limited effects over gene expression, challenging the concept that CTCF-mediated chromatin structures are a fundamental requirement for gene regulation. However, how CTCF and a perturbed chromatin structure impacts gene expression during development remains poorly understood. Here we link the loss of CTCF and gene regulation during patterning and organogenesis in a ctcf knockout zebrafish model. CTCF absence leads to loss of chromatin structure and affects the expression of thousands of genes, including many developmental regulators. Our results demonstrate the essential role of CTCF in providing the structural context for enhancer-promoter interactions, thus regulating developmental genes.

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.

Conserved Mechanisms, Novel Anatomies: The Developmental Basis of Fin Evolution and the Origin of Limbs

The transformation of paired fins into tetrapod limbs is one of the most intensively scrutinized events in animal evolution. Early anatomical and embryological datasets identified distinctive morphological regions within the appendage and posed hypotheses about how the loss, gain, and transformation of these regions could explain the observed patterns of both extant and fossil appendage diversity. These hypotheses have been put to the test by our growing understanding of patterning mechanisms that regulate formation of the appendage axes, comparisons of gene expression data from an array of phylogenetically informative taxa, and increasingly sophisticated and elegant experiments leveraging the latest molecular approaches. Together, these data demonstrate the remarkable conservation of developmental mechanisms, even across phylogenetically and morphologically disparate taxa, as well as raising new questions about the way we view homology, evolutionary novelty, and the often non-linear connection between morphology and gene expression. In this review, we present historical hypotheses regarding paired fin evolution and limb origins, summarize key aspects of central appendage patterning mechanisms in model and non-model species, address how modern comparative developmental data interface with our understanding of appendage anatomy, and highlight new approaches that promise to provide new insight into these well-traveled questions.

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

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

Understanding 3D genome organization by multidisciplinary methods

Understanding how chromatin is folded in the nucleus is fundamental to understanding its function. Although 3D genome organization has been historically difficult to study owing to a lack of relevant methodologies, major technological breakthroughs in genome-wide mapping of chromatin contacts and advances in imaging technologies in the twenty-first century considerably improved our understanding of chromosome conformation and nuclear architecture. In this Review, we discuss methods of 3D genome organization analysis, including sequencing-based techniques, such as Hi-C and its derivatives, Micro-C, DamID and others; microscopy-based techniques, such as super-resolution imaging coupled with fluorescence in situ hybridization (FISH), multiplex FISH, in situ genome sequencing and live microscopy methods; and computational and modelling approaches. We describe the most commonly used techniques and their contribution to our current knowledge of nuclear architecture and, finally, we provide a perspective on up-and-coming methods that open possibilities for future major discoveries.

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

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

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

3D genome organization during lymphocyte development and activation

Chromosomes have a complex three-dimensional (3D) architecture comprising A/B compartments, topologically associating domains and promoter-enhancer interactions. At all these levels, the 3D genome has functional consequences for gene transcription and therefore for cellular identity. The development and activation of lymphocytes involves strict control of gene expression by transcription factors (TFs) operating in a three-dimensionally organized chromatin landscape. As lymphocytes are indispensable for tissue homeostasis and pathogen defense, and aberrant lymphocyte activity is involved in a wide range of human morbidities, acquiring an in-depth understanding of the molecular mechanisms that control lymphocyte identity is highly relevant. Here we review current knowledge of the interplay between 3D genome organization and transcriptional control during B and T lymphocyte development and antigen-dependent activation, placing special emphasis on the role of TFs.

Topologically Associating Domains and Regulatory Landscapes in Development, Evolution and Disease

Animal genomes are folded in topologically associating domains (TADs) that have been linked to the regulation of the genes they contain by constraining regulatory interactions between cis-regulatory elements and promoters. Therefore, TADs are proposed as structural scaffolds for the establishment of regulatory landscapes (RLs). In this review, we discuss recent advances in the connection between TADs and gene regulation, their relationship with gene RLs and their dynamics during development and differentiation. Moreover, we describe how restructuring TADs may lead to pathological conditions, which explains their high evolutionary conservation, but at the same time it provides a substrate for the emergence of evolutionary innovations that lay at the origin of vertebrates and other phylogenetic clades.

Of Necks, Trunks and Tails: Axial Skeletal Diversity among Vertebrates

The axial skeleton of all vertebrates is composed of individual units known as vertebrae. Each vertebra has individual anatomical attributes, yet they can be classified in five different groups, namely cervical, thoracic, lumbar, sacral and caudal, according to shared characteristics and their association with specific body areas. Variations in vertebral number, size, morphological features and their distribution amongst the different regions of the vertebral column are a major source of the anatomical diversity observed among vertebrates. In this review I will discuss the impact of those variations on the anatomy of different vertebrate species and provide insights into the genetic origin of some remarkable morphological traits that often serve to classify phylogenetic branches or individual species, like the long trunks of snakes or the long necks of giraffes.

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.

Two Enhancers Regulate HoxB Genes Expression During Retinoic Acid-Induced Early Embryonic Stem Cells Differentiation Through Long-Range Chromatin Interactions

Homeobox B cluster (HoxB) genes play important roles in retinoic acid (RA)-induced early embryonic stem cells (ESCs) differentiation. Knowledge of regulation network of HoxB is important to further unveil the mechanism of ESCs differentiation. In this study, we identified two enhancers that were activated by RA treatment and 4C data showed long-range interactions between HoxB genes and the two enhancers. CRISPR/Cas9-mediated individual or compound deletion of the two enhancers significantly inhibits HoxB gene expression, and transcriptome analysis revealed that RA-induced early ESCs differentiation was blocked in the enhancer KO cells. We propose new mechanism by which two enhancers regulate HoxB gene expression by different regulation modes during RA-induced early ESCs differentiation through long-range chromatin interactions.

Comparison of the Hi-C, GAM and SPRITE methods using polymer models of chromatin

Hi-C, split-pool recognition of interactions by tag extension (SPRITE) and genome architecture mapping (GAM) are powerful technologies utilized to probe chromatin interactions genome wide, but how faithfully they capture three-dimensional (3D) contacts and how they perform relative to each other is unclear, as no benchmark exists. Here, we compare these methods in silico in a simplified, yet controlled, framework against known 3D structures of polymer models of murine and human loci, which can recapitulate Hi-C, GAM and SPRITE experiments and multiplexed fluorescence in situ hybridization (FISH) single-molecule conformations. We find that in silico Hi-C, GAM and SPRITE bulk data are faithful to the reference 3D structures whereas single-cell data reflect strong variability among single molecules. The minimal number of cells required in replicate experiments to return statistically similar contacts is different across the technologies, being lowest in SPRITE and highest in GAM under the same conditions. Noise-to-signal levels follow an inverse power law with detection efficiency and grow with genomic distance differently among the three methods, being lowest in GAM for genomic separations >1 Mb.

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.

Mechanisms of enhancer action: the known and the unknown

Differential gene expression mechanisms ensure cellular differentiation and plasticity to shape ontogenetic and phylogenetic diversity of cell types. A key regulator of differential gene expression programs are the enhancers, the gene-distal cis-regulatory sequences that govern spatiotemporal and quantitative expression dynamics of target genes. Enhancers are widely believed to physically contact the target promoters to effect transcriptional activation. However, our understanding of the full complement of regulatory proteins and the definitive mechanics of enhancer action is incomplete. Here, we review recent findings to present some emerging concepts on enhancer action and also outline a set of outstanding questions.

Insight Into Chromatin-Enriched RNA: A Key Chromatin Regulator in Tumors

Chromatin-enriched RNAs (cheRNAs) constitute a special class of long noncoding RNAs (lncRNAs) that are enriched around chromatin and function to activate neighboring or distal gene transcription. Recent studies have shown that cheRNAs affect chromatin structure and gene expression by recruiting chromatin modifiers or acting as bridges between distal enhancers and promoters. The abnormal transcription of cheRNAs plays an important role in the occurrence of many diseases, particularly tumors. The critical effect of cancer stem cells (CSCs) on the formation and development of tumors is well known, but the function of cheRNAs in tumorigenesis, especially in CSC proliferation and stemness maintenance, is not yet fully understood. This review focuses on the mechanisms of cheRNAs in epigenetic regulation and chromatin conformation and discusses the way cheRNAs function in CSCs to deepen the understanding of tumorigenesis and provide novel insight to advance tumor-targeting therapy.

The giant axolotl genome uncovers the evolution, scaling, and transcriptional control of complex gene loci

The axolotl is an important model organism because it is a tetrapod with a similar body plan to humans. Unlike humans, the axolotl regenerates limbs and other complex tissues. Therefore, the axolotl contributes to understanding evolution, development, and regeneration. With sophisticated tools for gene modification and tissue labeling, a fully assembled genome sequence was a sorely missing resource. Assembly was difficult because the genome size is 10× that of humans. Here, we use a cross-linking strategy called Hi-C to link together fragmented genome sequences to chromosome scale. We show that gene regulation occurs over very large genomic distances and that mitotic chromosomes are packaged efficiently.

Vertebrates harbor recognizably orthologous gene complements but vary 100-fold in genome size. How chromosomal organization scales with genome expansion is unclear, and how acute changes in gene regulation, as during axolotl limb regeneration, occur in the context of a vast genome has remained a riddle. Here, we describe the chromosome-scale assembly of the giant, 32 Gb axolotl genome. Hi-C contact data revealed the scaling properties of interphase and mitotic chromosome organization. Analysis of the assembly yielded understanding of the evolution of large, syntenic multigene clusters, including the Major Histocompatibility Complex (MHC) and the functional regulatory landscape of the Fibroblast Growth Factor 8 (Axfgf8) region. The axolotl serves as a primary model for studying successful regeneration.

Dmrt2 promotes transition of endochondral bone formation by linking Sox9 and Runx2

Endochondral bone formation is fundamental for skeletal development. During this process, chondrocytes undergo multiple steps of differentiation and coordinated transition from a proliferating to a hypertrophic stage, which is critical to advance skeletal development. Here, we identified the transcription factor Dmrt2 (double-sex and mab-3 related transcription factor 2) as a Sox9-inducible gene that promotes chondrocyte hypertrophy in pre-hypertrophic chondrocytes. Epigenetic analysis further demonstrated that Sox9 regulates Dmrt2 expression through an active enhancer located 18 kb upstream of the Dmrt2 gene and that this enhancer's chromatin status is progressively activated through chondrocyte differentiation. Dmrt2-knockout mice exhibited a dwarf phenotype with delayed initiation of chondrocyte hypertrophy. Dmrt2 augmented hypertrophic chondrocyte gene expression including Ihh through physical and functional interaction with Runx2. Furthermore, Dmrt2 deficiency reduced Runx2-dependent Ihh expression. Our findings suggest that Dmrt2 is critical for sequential chondrocyte differentiation during endochondral bone formation and coordinates the transcriptional network between Sox9 and Runx2.

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.

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.

CTCF loss has limited effects on global genome architecture in Drosophila despite critical regulatory functions

Vertebrate genomes are partitioned into contact domains defined by enhanced internal contact frequency and formed by two principal mechanisms: compartmentalization of transcriptionally active and inactive domains, and stalling of chromosomal loop-extruding cohesin by CTCF bound at domain boundaries. While Drosophila has widespread contact domains and CTCF, it is currently unclear whether CTCF-dependent domains exist in flies. We genetically ablate CTCF in Drosophila and examine impacts on genome folding and transcriptional regulation in the central nervous system. We find that CTCF is required to form a small fraction of all domain boundaries, while critically controlling expression patterns of certain genes and supporting nervous system function. We also find that CTCF recruits the pervasive boundary-associated factor Cp190 to CTCF-occupied boundaries and co-regulates a subset of genes near boundaries together with Cp190. These results highlight a profound difference in CTCF-requirement for genome folding in flies and vertebrates, in which a large fraction of boundaries are CTCF-dependent and suggest that CTCF has played mutable roles in genome architecture and direct gene expression control during metazoan evolution.

Developmental hourglass and heterochronic shifts in fin and limb development

How genetic changes are linked to morphological novelties and developmental constraints remains elusive. Here, we investigate genetic apparatuses that distinguish fish fins from tetrapod limbs by analyzing transcriptomes and open-chromatin regions (OCRs). Specifically, we compared mouse forelimb buds with the pectoral fin buds of an elasmobranch, the brown-banded bamboo shark (Chiloscyllium punctatum). A transcriptomic comparison with an accurate orthology map revealed both a mass heterochrony and hourglass-shaped conservation of gene expression between fins and limbs. Furthermore, open-chromatin analysis suggested that access to conserved regulatory sequences is transiently increased during mid-stage limb development. During this stage, stage-specific and tissue-specific OCRs were also enriched. Together, early and late stages of fin/limb development are more permissive to mutations than middle stages, which may have contributed to major morphological changes during the fin-to-limb evolution. We hypothesize that the middle stages are constrained by regulatory complexity that results from dynamic and tissue-specific transcriptional controls.

Animals come in all shapes and sizes. This diversity arose through genetic mutations during evolution, but it is unclear exactly how these variations led to the formation of new shapes. There is increasing evidence to suggest that not all shapes are possible and that variability between animals is limited by a phenomenon known as “developmental constraint”. These limitations direct parts of the body towards a specific shape as they develop in the embryo. Therefore, understanding the mechanisms underlying these developmental constraints could help explain how different body shapes evolved.

The limbs of humans and other mammals evolved from the fins of fish, and this transition is often used to study the role developmental constraints play in evolution. This is an ideal model as there is already a detailed fossil record mapping this evolutionary event, and data pinpointing some of the genes involved in the development of limbs and fins. But this data is incomplete, and a full comparison between the genes activated in the fin and the limb during embryonic development had not been achieved. This is because most fish used for research have undergone recent genetic changes, making it hard to spot which genetic differences are linked to the evolution of the limb.

To overcome this barrier, Onimaru et al. compared genetic data from the developing limbs of mice to the developing fins of the brown-banded bamboo shark, which evolves much slower than other fish. This revealed that although many genes commonly played a role in the development of the fin and the limb in the embryo, the activity of these shared genes was not the same. For example, genes that switched on in the late stages of limb development, switched off in the late stages of fin development. But in the middle of development, those differences were relatively small and both species activated very similar sets of genes. Many of these genes were pleiotropic, which means they have important roles in other tissues and therefore mutate less often. This suggests that the mid-stage of limb development is under the strongest level of constraint.

Darwin’s theory of natural selection explains that mutations drive evolution. But the theory cannot predict what kinds of new body shapes new mutations will produce. Understanding how the activity levels of different genes affect development could help to fill this knowledge gap. This has potential medical applications, for example, understanding why some genetic changes cause more serious problems than others. This work suggests that mutations in genes that are active during the mid-stage of limb development may have the most serious impact.

Giant lungfish genome elucidates the conquest of land by vertebrates

Lungfishes belong to lobe-fined fish (Sarcopterygii) that, in the Devonian period, 'conquered' the land and ultimately gave rise to all land vertebrates, including humans1-3. Here we determine the chromosome-quality genome of the Australian lungfish (Neoceratodus forsteri), which is known to have the largest genome of any animal. The vast size of this genome, which is about 14× larger than that of humans, is attributable mostly to huge intergenic regions and introns with high repeat content (around 90%), the components of which resemble those of tetrapods (comprising mainly long interspersed nuclear elements) more than they do those of ray-finned fish. The lungfish genome continues to expand independently (its transposable elements are still active), through mechanisms different to those of the enormous genomes of salamanders. The 17 fully assembled lungfish macrochromosomes maintain synteny to other vertebrate chromosomes, and all microchromosomes maintain conserved ancient homology with the ancestral vertebrate karyotype. Our phylogenomic analyses confirm previous reports that lungfish occupy a key evolutionary position as the closest living relatives to tetrapods4,5, underscoring the importance of lungfish for understanding innovations associated with terrestrialization. Lungfish preadaptations to living on land include the gain of limb-like expression in developmental genes such as hoxc13 and sall1 in their lobed fins. Increased rates of evolution and the duplication of genes associated with obligate air-breathing, such as lung surfactants and the expansion of odorant receptor gene families (which encode proteins involved in detecting airborne odours), contribute to the tetrapod-like biology of lungfishes. These findings advance our understanding of this major transition during vertebrate evolution.

Folding Keratin Gene Clusters during Skin Regional Specification

Regional specification is critical for skin development, regeneration, and evolution. The contribution of epigenetics in this process remains unknown. Here, using avian epidermis, we find two major strategies regulate β-keratin gene clusters. (1) Over the body, macro-regional specificities (scales, feathers, claws, etc.) established by typical enhancers control five subclusters located within the epidermal differentiation complex on chromosome 25; (2) within a feather, micro-regional specificities are orchestrated by temporospatial chromatin looping of the feather β-keratin gene cluster on chromosome 27. Analyses suggest a three-factor model for regional specification: competence factors (e.g., AP1) make chromatin accessible, regional specifiers (e.g., Zic1) target specific genome regions, and chromatin regulators (e.g., CTCF and SATBs) establish looping configurations. Gene perturbations disrupt morphogenesis and histo-differentiation. This chicken skin paradigm advances our understanding of how regulation of big gene clusters can set up a two-dimensional body surface map.

Of numbers and movement - understanding transcription factor pathogenesis by advanced microscopy

Transcription factors (TFs) are life-sustaining and, therefore, the subject of intensive research. By regulating gene expression, TFs control a plethora of developmental and physiological processes, and their abnormal function commonly leads to various developmental defects and diseases in humans. Normal TF function often depends on gene dosage, which can be altered by copy-number variation or loss-of-function mutations. This explains why TF haploinsufficiency (HI) can lead to disease. Since aberrant TF numbers frequently result in pathogenic abnormalities of gene expression, quantitative analyses of TFs are a priority in the field. In vitro single-molecule methodologies have significantly aided the identification of links between TF gene dosage and transcriptional outcomes. Additionally, advances in quantitative microscopy have contributed mechanistic insights into normal and aberrant TF function. However, to understand TF biology, TF-chromatin interactions must be characterised in vivo, in a tissue-specific manner and in the context of both normal and altered TF numbers. Here, we summarise the advanced microscopy methodologies most frequently used to link TF abundance to function and dissect the molecular mechanisms underlying TF HIs. Increased application of advanced single-molecule and super-resolution microscopy modalities will improve our understanding of how TF HIs drive disease.

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.

The relationship between genome structure and function

Precise patterns of gene expression in metazoans are controlled by three classes of regulatory elements: promoters, enhancers and boundary elements. During differentiation and development, these elements form specific interactions in dynamic higher-order chromatin structures. However, the relationship between genome structure and its function in gene regulation is not completely understood. Here we review recent progress in this field and discuss whether genome structure plays an instructive role in regulating gene expression or is a reflection of the activity of the regulatory elements of the genome.

Senescence-activated enhancer landscape orchestrates the senescence-associated secretory phenotype in murine fibroblasts

The three-dimensional configuration of the chromatin architecture is known to be crucial for alterations in the transcriptional network; however, the underlying mechanisms of epigenetic control of senescence-related gene expression by modulating the chromatin architecture remain unknown. Here, we demonstrate frequent chromosomal compartment switching during mouse embryonic fibroblasts (MEFs) replicative senescence as characterized by senescence-inactivated (SIAEs) and -activated enhancers (SAEs) in topologically associated domains (TADs). Mechanistically, SAEs are closely correlated with senescence-associated secretory phenotype (SASP) genes, which are a key transcriptional feature of an aging microenvironment that contributes to tumor progression, aging acceleration, and immunoinflammatory responses. Moreover, SAEs can positively regulate robust changes in SASP expression. The transcription factor CCAAT/enhancer binding protein α (C/EBPα) is capable of enhancing SAE activity, which accelerates the emergence of SAEs flanking SASPs and the secretion of downstream factors, contributing to the progression of senescence. Our results provide novel insight into the TAD-related control of SASP gene expression by revealing hierarchical roles of the chromatin architecture, transcription factors, and enhancer activity in the regulation of cellular senescence.

Mammalian-specific ectodermal enhancers control the expression of Hoxc genes in developing nails and hair follicles

Vertebrate Hox genes are critical for the establishment of structures during the development of the main body axis. Subsequently, they play important roles either in organizing secondary axial structures such as the appendages, or during homeostasis in postnatal stages and adulthood. Here, we set up to analyze their elusive function in the ectodermal compartment, using the mouse limb bud as a model. We report that the HoxC gene cluster was co-opted to be transcribed in the distal limb ectoderm, where it is activated following the rule of temporal colinearity. These ectodermal cells subsequently produce various keratinized organs such as nails or claws. Accordingly, deletion of the HoxC cluster led to mice lacking nails (anonychia), a condition stronger than the previously reported loss of function of Hoxc13, which is the causative gene of the ectodermal dysplasia 9 (ECTD9) in human patients. We further identified two mammalian-specific ectodermal enhancers located upstream of the HoxC gene cluster, which together regulate Hoxc gene expression in the hair and nail ectodermal organs. Deletion of these regulatory elements alone or in combination revealed a strong quantitative component in the regulation of Hoxc genes in the ectoderm, suggesting that these two enhancers may have evolved along with the mammalian taxon to provide the level of HOXC proteins necessary for the full development of hair and nail.

HOX gene cluster (de)regulation in brain: from neurodevelopment to malignant glial tumours

HOX genes encode a family of evolutionarily conserved homeodomain transcription factors that are crucial both during development and adult life. In humans, 39 HOX genes are arranged in four clusters (HOXA, B, C, and D) in chromosomes 7, 17, 12, and 2, respectively. During embryonic development, particular epigenetic states accompany their expression along the anterior-posterior body axis. This tightly regulated temporal-spatial expression pattern reflects their relative chromosomal localization, and is critical for normal embryonic brain development when HOX genes are mainly expressed in the hindbrain and mostly absent in the forebrain region. Epigenetic marks, mostly polycomb-associated, are dynamically regulated at HOX loci and regulatory regions to ensure the finely tuned HOX activation and repression, highlighting a crucial epigenetic plasticity necessary for homeostatic development. HOX genes are essentially absent in healthy adult brain, whereas they are detected in malignant brain tumours, namely gliomas, where HOX genes display critical roles by regulating several hallmarks of cancer. Here, we review the major mechanisms involved in HOX genes (de)regulation in the brain, from embryonic to adult stages, in physiological and oncologic conditions. We focus particularly on the emerging causes of HOX gene deregulation in glioma, as well as on their functional and clinical implications.

Marchantia TCP transcription factor activity correlates with three-dimensional chromatin structure

Information in the genome is not only encoded within sequence or epigenetic modifications, but is also found in how it folds in three-dimensional space. The formation of self-interacting genomic regions, named topologically associated domains (TADs), is known as a key feature of genome organization beyond the nucleosomal level. However, our understanding of the formation and function of TADs in plants is extremely limited. Here we show that the genome of Marchantia polymorpha, a member of a basal land plant lineage, exhibits TADs with epigenetic features similar to those of higher plants. By analysing various epigenetic marks across Marchantia TADs, we find that these regions generally represent interstitial heterochromatin and their borders are enriched with Marchantia transcription factor TCP1. We also identify a type of TAD that we name 'TCP1-rich TAD', in which genomic regions are highly accessible and are densely bound by TCP1 proteins. Transcription of TCP1 target genes differs on the basis gene location, and those in TCP1-rich TADs clearly show a lower expression level. In tcp1 mutant lines, neither TCP1-bound TAD borders nor TCP1-rich TADs display drastically altered chromatin organization patterns, suggesting that, in Marchantia, TCP1 is dispensable for TAD formation. However, we find that in tcp1 mutants, genes residing in TCP1-rich TADs have a greater extent of expression fold change as opposed to genes that do not belong to these TADs. Our results suggest that, besides standing as spatial chromatin-packing modules, plant TADs function as nuclear microcompartments associated with transcription factor activities.

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.

Chronic loss of STAG2 leads to altered chromatin structure contributing to de-regulated transcription in AML

Background

The cohesin complex plays a major role in folding the human genome into 3D structural domains. Mutations in members of the cohesin complex are known early drivers of myelodysplastic syndromes (MDS) and acute myeloid leukaemia (AML), with STAG2 the most frequently mutated complex member.

Methods

Here we use functional genomics (RNA-seq, ChIP-seq and HiChIP) to investigate the impact of chronic STAG2 loss on three-dimensional genome structure and transcriptional programming in a clinically relevant model of chronic STAG2 loss.

Results

The chronic loss of STAG2 led to loss of smaller loop domains and the maintenance/formation of large domains that, in turn, led to altered genome compartmentalisation. These changes in genome structure resulted in altered gene expression, including deregulation of the HOXA locus and the MAPK signalling pathway, resulting in increased sensitivity to MEK inhibition.

Conclusions

The altered genomic architecture driven by the chronic loss of STAG2 results in altered gene expression that may contribute to leukaemogenesis and may be therapeutically targeted.

Tidying-up the plant nuclear space: domains, functions, and dynamics

Understanding how the packaging of chromatin in the nucleus is regulated and organized to guide complex cellular and developmental programmes, as well as responses to environmental cues is a major question in biology. Technological advances have allowed remarkable progress within this field over the last years. However, we still know very little about how the 3D genome organization within the cell nucleus contributes to the regulation of gene expression. The nuclear space is compartmentalized in several domains such as the nucleolus, chromocentres, telomeres, protein bodies, and the nuclear periphery without the presence of a membrane around these domains. The role of these domains and their possible impact on nuclear activities is currently under intense investigation. In this review, we discuss new data from research in plants that clarify functional links between the organization of different nuclear domains and plant genome function with an emphasis on the potential of this organization for gene regulation.

Polycomb Repressive Complexes in Hox Gene Regulation: Silencing and Beyond: The Functional Dynamics of Polycomb Repressive Complexes in Hox Gene Regulation

The coordinated expression of the Hox gene family encoding transcription factors is critical for proper embryonic development and patterning. Major efforts have thus been dedicated to understanding mechanisms controlling Hox expression. In addition to the temporal and spatial sequential activation of Hox genes, proper embryonic development requires that Hox genes get differentially silenced in a cell-type specific manner as development proceeds. Factors contributing to Hox silencing include the polycomb repressive complexes (PRCs), which control gene expression through epigenetic modifications. This review focuses on PRC-dependent regulation of the Hox genes and is aimed at integrating the growing complexity of PRC functional properties in the context of Hox regulation. In particular, mechanisms underlying PRC binding dynamics as well as a series of studies that have revealed the impact of PRC on the 3D organization of the genome is discussed, which has a significant role on Hox regulation during development.

Many facades of CTCF unified by its coding for three-dimensional genome architecture

CCCTC-binding factor (CTCF) is a multifunctional zinc finger protein that is conserved in metazoan species. CTCF is consistently found to play an important role in many diverse biological processes. CTCF/cohesin-mediated active chromatin 'loop extrusion' architects three-dimensional (3D) genome folding. The 3D architectural role of CTCF underlies its multifarious functions, including developmental regulation of gene expression, protocadherin (Pcdh) promoter choice in the nervous system, immunoglobulin (Ig) and T-cell receptor (Tcr) V(D)J recombination in the immune system, homeobox (Hox) gene control during limb development, as well as many other aspects of biology. Here, we review the pleiotropic functions of CTCF from the perspective of its essential role in 3D genome architecture and topological promoter/enhancer selection. We envision the 3D genome as an enormous complex architecture, with tens of thousands of CTCF sites as connecting nodes and CTCF proteins as mysterious bonds that glue together genomic building parts with distinct articulation joints. In particular, we focus on the internal mechanisms by which CTCF controls higher order chromatin structures that manifest its many façades of physiological and pathological functions. We also discuss the dichotomic role of CTCF sites as intriguing 3D genome nodes for seemingly contradictory 'looping bridges' and 'topological insulators' to frame a beautiful magnificent house for a cell's nuclear home.

Exonuclease combinations reduce noises in 3D genomics technologies

Chromosome conformation-capture technologies are widely used in 3D genomics; however, experimentally, such methods have high-noise limitations and, therefore, require significant bioinformatics efforts to extract reliable distal interactions. Miscellaneous undesired linear DNAs, present during proximity-ligation, represent a main noise source, which needs to be minimized or eliminated. In this study, different exonuclease combinations were tested to remove linear DNA fragments from a circularized DNA preparation. This method efficiently removed linear DNAs, raised the proportion of annulation and increased the valid-pairs ratio from ∼40% to ∼80% for enhanced interaction detection in standard Hi-C. This strategy is applicable for development of various 3D genomics technologies, or optimization of Hi-C sequencing efficiency.

Molecular Basis of the Function of Transcriptional Enhancers

Transcriptional enhancers are major genomic elements that control gene activity in eukaryotes. Recent studies provided deeper insight into the temporal and spatial organization of transcription in the nucleus, the role of non-coding RNAs in the process, and the epigenetic control of gene expression. Thus, multiple molecular details of enhancer functioning were revealed. Here, we describe the recent data and models of molecular organization of enhancer-driven transcription.

Active and repressed biosynthetic gene clusters have spatially distinct chromosome states

While colocalization within a bacterial operon enables coexpression of the constituent genes, the mechanistic logic of clustering of nonhomologous monocistronic genes in eukaryotes is not immediately obvious. Biosynthetic gene clusters that encode pathways for specialized metabolites are an exception to the classical eukaryote rule of random gene location and provide paradigmatic exemplars with which to understand eukaryotic cluster dynamics and regulation. Here, using 3C, Hi-C, and Capture Hi-C (CHi-C) organ-specific chromosome conformation capture techniques along with high-resolution microscopy, we investigate how chromosome topology relates to transcriptional activity of clustered biosynthetic pathway genes in Arabidopsis thaliana Our analyses reveal that biosynthetic gene clusters are embedded in local hot spots of 3D contacts that segregate cluster regions from the surrounding chromosome environment. The spatial conformation of these cluster-associated domains differs between transcriptionally active and silenced clusters. We further show that silenced clusters associate with heterochromatic chromosomal domains toward the periphery of the nucleus, while transcriptionally active clusters relocate away from the nuclear periphery. Examination of chromosome structure at unrelated clusters in maize, rice, and tomato indicates that integration of clustered pathway genes into distinct topological domains is a common feature in plant genomes. Our results shed light on the potential mechanisms that constrain coexpression within clusters of nonhomologous eukaryotic genes and suggest that gene clustering in the one-dimensional chromosome is accompanied by compartmentalization of the 3D chromosome.

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

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

Transcriptional activation during cell reprogramming correlates with the formation of 3D open chromatin hubs

Chromosome structure is a crucial regulatory factor for a wide range of nuclear processes. Chromosome conformation capture (3C)-based experiments combined with computational modelling are pivotal for unveiling 3D chromosome structure. Here, we introduce TADdyn, a tool that integrates time-course 3C data, restraint-based modelling, and molecular dynamics to simulate the structural rearrangements of genomic loci in a completely data-driven way. We apply TADdyn on in situ Hi-C time-course experiments studying the reprogramming of murine B cells to pluripotent cells, and characterize the structural rearrangements that take place upon changes in the transcriptional state of 21 genomic loci of diverse expression dynamics. By measuring various structural and dynamical properties, we find that during gene activation, the transcription starting site contacts with open and active regions in 3D chromatin domains. We propose that these 3D hubs of open and active chromatin may constitute a general feature to trigger and maintain gene transcription.

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.