A global control region defines a chromosomal regulatory landscape containing the HoxD cluster

The impact of inversions across 33,924 families with rare disease from a national genome sequencing project

Detection of structural variants (SVs) is currently biased toward those that alter copy number. The relative contribution of inversions toward genetic disease is unclear. In this study, we analyzed genome sequencing data for 33,924 families with rare disease from the 100,000 Genomes Project. From a database hosting >500 million SVs, we focused on 351 genes where haploinsufficiency is a confirmed disease mechanism and identified 47 ultra-rare rearrangements that included an inversion (24 bp to 36.4 Mb, 20/47 de novo). Validation utilized a number of orthogonal approaches, including retrospective exome analysis. RNA-seq data supported the respective diagnoses for six participants. Phenotypic blending was apparent in four probands. Diagnostic odysseys were a common theme (>50 years for one individual), and targeted analysis for the specific gene had already been performed for 30% of these individuals but with no findings. We provide formal confirmation of a European founder origin for an intragenic MSH2 inversion. For two individuals with complex SVs involving the MECP2 mutational hotspot, ambiguous SV structures were resolved using long-read sequencing, influencing clinical interpretation. A de novo inversion of HOXD11-13 was uncovered in a family with Kantaputra-type mesomelic dysplasia. Lastly, a complex translocation disrupting APC and involving nine rearranged segments confirmed a clinical diagnosis for three family members and resolved a conundrum for a sibling with a single polyp. Overall, inversions play a small but notable role in rare disease, likely explaining the etiology in around 1/750 families across heterogeneous clinical cohorts.

Cellular control of protein levels: A systems biology perspective

How cells regulate protein levels is a central question of biology. Over the past decades, molecular biology research has provided profound insights into the mechanisms and the molecular machinery governing each step of the gene expression process, from transcription to protein degradation. Recent advances in transcriptomics and proteomics have complemented our understanding of these fundamental cellular processes with a quantitative, systems-level perspective. Multi-omic studies revealed significant quantitative, kinetic and functional differences between the genome, transcriptome and proteome. While protein levels often correlate with mRNA levels, quantitative investigations have demonstrated a substantial impact of translation and protein degradation on protein expression control. In addition, protein-level regulation appears to play a crucial role in buffering protein abundances against undesirable mRNA expression variation. These findings have practical implications for many fields, including gene function prediction and precision medicine.

Conserved Noncoding Elements Evolve Around the Same Genes Throughout Metazoan Evolution

Conserved noncoding elements (CNEs) are DNA sequences located outside of protein-coding genes that can remain under purifying selection for up to hundreds of millions of years. Studies in vertebrate genomes have revealed that most CNEs carry out regulatory functions. Notably, many of them are enhancers that control the expression of homeodomain transcription factors and other genes that play crucial roles in embryonic development. To further our knowledge of CNEs in other parts of the animal tree, we conducted a large-scale characterization of CNEs in more than 50 genomes from three of the main branches of the metazoan tree: Cnidaria, Mollusca, and Arthropoda. We identified hundreds of thousands of CNEs and reconstructed the temporal dynamics of their appearance in each lineage, as well as determining their spatial distribution across genomes. We show that CNEs evolve repeatedly around the same genes across the Metazoa, including around homeodomain genes and other transcription factors; they also evolve repeatedly around genes involved in neural development. We also show that transposons are a major source of CNEs, confirming previous observations from vertebrates and suggesting that they have played a major role in wiring developmental gene regulatory mechanisms since the dawn of animal evolution.

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.

A distant global control region is essential for normal expression of anterior HOXA genes during mouse and human craniofacial development

Craniofacial abnormalities account for approximately one third of birth defects. The regulatory programs that build the face require precisely controlled spatiotemporal gene expression, achieved through tissue-specific enhancers. Clusters of coactivated enhancers and their target genes, known as superenhancers, are important in determining cell identity but have been largely unexplored in development. In this study we identified superenhancer regions unique to human embryonic craniofacial tissue. To demonstrate the importance of such regions in craniofacial development and disease, we focused on an ~600 kb noncoding region located between NPVF and NFE2L3. We identified long range interactions with this region in both human and mouse embryonic craniofacial tissue with the anterior portion of the HOXA gene cluster. Mice lacking this superenhancer exhibit perinatal lethality, and present with highly penetrant skull defects and orofacial clefts phenocopying Hoxa2-/- mice. Moreover, we identified two cases of de novo copy number changes of the superenhancer in humans both with severe craniofacial abnormalities. This evidence suggests we have identified a critical noncoding locus control region that specifically regulates anterior HOXA genes and copy number changes are pathogenic in human patients.

The chromatin signatures of enhancers and their dynamic regulation

Enhancers are cis-regulatory elements that can stimulate gene expression from distance, and drive precise spatiotemporal gene expression profiles during development. Functional enhancers display specific features including an open chromatin conformation, Histone H3 lysine 27 acetylation, Histone H3 lysine 4 mono-methylation enrichment, and enhancer RNAs production. These features are modified upon developmental cues which impacts their activity. In this review, we describe the current state of knowledge about enhancer functions and the diverse chromatin signatures found on enhancers. We also discuss the dynamic changes of enhancer chromatin signatures, and their impact on lineage specific gene expression profiles, during development or cellular differentiation.

Ancient genomic linkage couples metabolism with erythroid development

Generation of mature cells from progenitors requires tight coupling of differentiation and metabolism. During erythropoiesis, erythroblasts are required to massively upregulate globin synthesis then clear extraneous material and enucleate to produce erythrocytes1-3. Nprl3 has remained in synteny with the α-globin genes for >500 million years4, and harbours the majority of the α-globin enhancers5. Nprl3 is a highly conserved inhibitor of mTORC1, which controls cellular metabolism. However, whether Nprl3 itself serves an erythroid role is unknown. Here, we show that Nprl3 is a key regulator of erythroid metabolism. Using Nprl3-deficient fetal liver and adult competitive bone marrow - fetal liver chimeras, we show that NprI3 is required for sufficient erythropoiesis. Loss of Nprl3 elevates mTORC1 signalling, suppresses autophagy and disrupts erythroblast glycolysis and redox control. Human CD34+ progenitors lacking NPRL3 produce fewer enucleated cells and demonstrate dysregulated mTORC1 signalling in response to nutrient availability and erythropoietin. Finally, we show that the α-globin enhancers upregulate NprI3 expression, and that this activity is necessary for optimal erythropoiesis. Therefore, the anciently conserved linkage of NprI3, α-globin and their associated enhancers has enabled coupling of metabolic and developmental control in erythroid cells. This may enable erythropoiesis to adapt to fluctuating nutritional and environmental conditions.

Hand2 represses non-cardiac cell fates through chromatin remodeling at cis- regulatory elements

Developmental studies have revealed the importance of the transcription factor Hand2 in cardiac development. Hand2 promotes cardiac progenitor differentiation and epithelial maturation, while repressing other tissue types. The mechanisms underlying the promotion of cardiac fates are far better understood than those underlying the repression of alternative fates. Here, we assess Hand2-dependent changes in gene expression and chromatin remodeling in cardiac progenitors of zebrafish embryos. Cell-type specific transcriptome analysis shows a dual function for Hand2 in activation of cardiac differentiation genes and repression of pronephric pathways. We identify functional cis- regulatory elements whose chromatin accessibility are increased in hand2 mutant cells. These regulatory elements associate with non-cardiac gene expression, and drive reporter gene expression in tissues associated with Hand2-repressed genes. We find that functional Hand2 is sufficient to reduce non-cardiac reporter expression in cardiac lineages. Taken together, our data support a model of Hand2-dependent coordination of transcriptional programs, not only through transcriptional activation of cardiac and epithelial maturation genes, but also through repressive chromatin remodeling at the DNA regulatory elements of non-cardiac genes.

Hox genes in development and beyond

Hox genes encode evolutionarily conserved transcription factors that are essential for the proper development of bilaterian organisms. Hox genes are unique because they are spatially and temporally regulated during development in a manner that is dictated by their tightly linked genomic organization. Although their genetic function during embryonic development has been interrogated, less is known about how these transcription factors regulate downstream genes to direct morphogenetic events. Moreover, the continued expression and function of Hox genes at postnatal and adult stages highlights crucial roles for these genes throughout the life of an organism. Here, we provide an overview of Hox genes, highlighting their evolutionary history, their unique genomic organization and how this impacts the regulation of their expression, what is known about their protein structure, and their deployment in development and beyond.

Evolution of stickleback spines through independent cis-regulatory changes at HOXDB

Understanding the mechanisms leading to new traits or additional features in organisms is a fundamental goal of evolutionary biology. We show that HOXDB regulatory changes have been used repeatedly in different fish genera to alter the length and number of the prominent dorsal spines used to classify stickleback species. In Gasterosteus aculeatus (typically 'three-spine sticklebacks'), a variant HOXDB allele is genetically linked to shortening an existing spine and adding an additional spine. In Apeltes quadracus (typically 'four-spine sticklebacks'), a variant HOXDB allele is associated with lengthening a spine and adding an additional spine in natural populations. The variant alleles alter the same non-coding enhancer region in the HOXDB locus but do so by diverse mechanisms, including single-nucleotide polymorphisms, deletions and transposable element insertions. The independent regulatory changes are linked to anterior expansion or contraction of HOXDB expression. We propose that associated changes in spine lengths and numbers are partial identity transformations in a repeating skeletal series that forms major defensive structures in fish. Our findings support the long-standing hypothesis that natural Hox gene variation underlies key patterning changes in wild populations and illustrate how different mutational mechanisms affecting the same region may produce opposite gene expression changes with similar phenotypic outcomes.

Synthetic regulatory reconstitution reveals principles of mammalian Hox cluster regulation

Precise Hox gene expression is crucial for embryonic patterning. Intra-Hox transcription factor binding and distal enhancer elements have emerged as the major regulatory modules controlling Hox gene expression. However, quantifying their relative contributions has remained elusive. Here, we introduce "synthetic regulatory reconstitution," a conceptual framework for studying gene regulation, and apply it to the HoxA cluster. We synthesized and delivered variant rat HoxA clusters (130 to 170 kilobases) to an ectopic location in the mouse genome. We found that a minimal HoxA cluster recapitulated correct patterns of chromatin remodeling and transcription in response to patterning signals, whereas the addition of distal enhancers was needed for full transcriptional output. Synthetic regulatory reconstitution could provide a generalizable strategy for deciphering the regulatory logic of gene expression in complex genomes.

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.

Inferring mammalian tissue-specific regulatory conservation by predicting tissue-specific differences in open chromatin

BACKGROUND

Evolutionary conservation is an invaluable tool for inferring functional significance in the genome, including regions that are crucial across many species and those that have undergone convergent evolution. Computational methods to test for sequence conservation are dominated by algorithms that examine the ability of one or more nucleotides to align across large evolutionary distances. While these nucleotide alignment-based approaches have proven powerful for protein-coding genes and some non-coding elements, they fail to capture conservation of many enhancers, distal regulatory elements that control spatial and temporal patterns of gene expression. The function of enhancers is governed by a complex, often tissue- and cell type-specific code that links combinations of transcription factor binding sites and other regulation-related sequence patterns to regulatory activity. Thus, function of orthologous enhancer regions can be conserved across large evolutionary distances, even when nucleotide turnover is high.

RESULTS

We present a new machine learning-based approach for evaluating enhancer conservation that leverages the combinatorial sequence code of enhancer activity rather than relying on the alignment of individual nucleotides. We first train a convolutional neural network model that can predict tissue-specific open chromatin, a proxy for enhancer activity, across mammals. Next, we apply that model to distinguish instances where the genome sequence would predict conserved function versus a loss of regulatory activity in that tissue. We present criteria for systematically evaluating model performance for this task and use them to demonstrate that our models accurately predict tissue-specific conservation and divergence in open chromatin between primate and rodent species, vastly out-performing leading nucleotide alignment-based approaches. We then apply our models to predict open chromatin at orthologs of brain and liver open chromatin regions across hundreds of mammals and find that brain enhancers associated with neuron activity have a stronger tendency than the general population to have predicted lineage-specific open chromatin.

CONCLUSION

The framework presented here provides a mechanism to annotate tissue-specific regulatory function across hundreds of genomes and to study enhancer evolution using predicted regulatory differences rather than nucleotide-level conservation measurements.

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.

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.

A novel microdeletion upstream of HOXD13 in a Chinese family with synpolydactyly

Synpolydactyly (SPD) is a digital malformation with the typical clinical phenotype of the webbing of 3/4 fingers and/or 4/5 toes, and combined with polydactyly. In this study, we investigated a Chinese family with SPD and genetic analysis found that all of the affected individuals in the family carry a heterozygous 11,451 bp microdeletion at chr2:176933872-176945322 (GRCh37), which is located upstream of HOXD13 gene, the known disease gene for SPD1. All the affected individuals in the family carry the heterozygous deletion variant, and the variant co-segregated with SPD in the family. Thus, we speculate that the 11,451 bp microdeletion is the disease-causing variant in the family. To date, the microdeletion associating with SPD1 which we identified is the smallest deletion upstream of the HOXD13 gene and not altering the sequence of the HOXD13 gene.

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.

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.

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.

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.

The plant endoplasmic reticulum: an organized chaos of tubules and sheets with multiple functions

The endoplasmic reticulum is a fascinating organelle at the core of the secretory pathway. It is responsible for the synthesis of one third of the cellular proteome and, in plant cells, it produces receptors and transporters of hormones as well as the proteins responsible for the biosynthesis of critical components of a cellulosic cell wall. The endoplasmic reticulum structure resembles a spider-web network of interconnected tubules and cisternae that pervades the cell. The study of the dynamics and interaction of this organelles with other cellular structures such as the plasma membrane, the Golgi apparatus and the cytoskeleton, have been permitted by the implementation of fluorescent protein and advanced confocal imaging. In this review, we report on the findings that contributed towards the understanding of the endoplasmic reticulum morphology and function with the aid of fluorescent proteins, focusing on the contributions provided by pioneering work from the lab of the late Professor Chris Hawes.

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.

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.

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.

Disparate expression specificities coded by a shared Hox-C enhancer

Can a single regulatory sequence be shared by two genes undergoing functional divergence? Here we describe a single promiscuous enhancer within the Drosophila Antennapedia Complex, EO053, that directs aspects of the expression of two adjacent genes, pb (a Hox2 ortholog) and zen2 (a divergent Hox3 paralog), with disparate spatial and temporal expression patterns. We were unable to separate the pb-like and zen2-like specificities within EO053, and we identify sequences affecting both expression patterns. Importantly, genomic deletion experiments demonstrate that EO053 cooperates with additional pb- and zen2-specific enhancers to regulate the mRNA expression of both genes. We examine sequence conservation of EO053 within the Schizophora, and show that patterns of synteny between the Hox2 and Hox3 orthologs in Arthropods are consistent with a shared regulatory relationship extending prior to the Hox3/zen divergence. Thus, EO053 represents an example of two genes having evolved disparate outputs while utilizing this shared regulatory region.

Editorial note: This article has been through an editorial process in which the authors decide how to respond to the issues raised during peer review. The Reviewing Editor's assessment is that all the issues have been addressed (see decision letter).

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

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

The regulatory landscapes of developmental genes

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

Reconstruction of the Global Neural Crest Gene Regulatory Network In Vivo

Precise control of developmental processes is encoded in the genome in the form of gene regulatory networks (GRNs). Such multi-factorial systems are difficult to decode in vertebrates owing to their complex gene hierarchies and dynamic molecular interactions. Here we present a genome-wide in vivo reconstruction of the GRN underlying development of the multipotent neural crest (NC) embryonic cell population. By coupling NC-specific epigenomic and transcriptional profiling at population and single-cell levels with genome/epigenome engineering in vivo, we identify multiple regulatory layers governing NC ontogeny, including NC-specific enhancers and super-enhancers, novel trans-factors, and cis-signatures allowing reverse engineering of the NC-GRN at unprecedented resolution. Furthermore, identification and dissection of divergent upstream combinatorial regulatory codes has afforded new insights into opposing gene circuits that define canonical and neural NC fates early during NC ontogeny. Our integrated approach, allowing dissection of cell-type-specific regulatory circuits in vivo, has broad implications for GRN discovery and investigation.

Fryns type mesomelic dysplasia of the upper limbs caused by inverted duplications of the HOXD gene cluster

The HoxD cluster is critical for vertebrate limb development. Enhancers located in both the telomeric and centromeric gene deserts flanking the cluster regulate the transcription of HoxD genes. In rare patients, duplications, balanced translocations or inversions misregulating HOXD genes are responsible for mesomelic dysplasia of the upper and lower limbs. By aCGH, whole-genome mate-pair sequencing, long-range PCR and fiber fluorescent in situ hybridization, we studied patients from two families displaying mesomelic dysplasia limited to the upper limbs. We identified microduplications including the HOXD cluster and showed that microduplications were in an inverted orientation and inserted between the HOXD cluster and the telomeric enhancers. Our results highlight the existence of an autosomal dominant condition consisting of isolated ulnar dysplasia caused by microduplications inserted between the HOXD cluster and the telomeric enhancers. The duplications likely disconnect the HOXD9 to HOXD11 genes from their regulatory sequences. This presumptive loss-of-function may have contributed to the phenotype. In both cases, however, these rearrangements brought HOXD13 closer to telomeric enhancers, suggesting that the alterations derive from the dominant-negative effect of this digit-specific protein when ectopically expressed during the early development of forearms, through the disruption of topologically associating domain structure at the HOXD locus.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

FOXO3 on the Road to Longevity: Lessons From SNPs and Chromatin Hubs

Health span is driven by a precise interplay between genes and the environment. Cell response to environmental cues is mediated by signaling cascades and genetic variants that affect gene expression by regulating chromatin plasticity. Indeed, they can promote the interaction of promoters with regulatory elements by forming active chromatin hubs.

FOXO3 encodes a transcription factor with a strong impact on aging and age-related phenotypes, as it regulates stress response, therefore affecting lifespan. A significant association has been shown between human longevity and several FOXO3 variants located in intron 2. This haplotype block forms a putative aging chromatin hub in which FOXO3 has a central role, as it modulates the physical connection and activity of neighboring genes involved in age-related processes.

Here we describe the role of FOXO3 and its single-nucleotide polymorphisms (SNPs) in healthy aging, with a focus on the enhancer region encompassing the SNP rs2802292, which upregulates FOXO3 expression and can promote the activity of the aging hub in response to different stress stimuli. FOXO3 protective effect on lifespan may be due to the accessibility of this region to transcription factors promoting its expression. This could in part explain the differences in FOXO3 association with longevity between genders, as its activity in females may be modulated by estrogens through estrogen receptor response elements located in the rs2802292-encompassing region. Altogether, the molecular mechanisms described here may help establish whether the rs2802292 SNP can be taken advantage of in predictive medicine and define the potential of targeting FOXO3 for age-related diseases.

Inference of the ancestral vertebrate phenotype through vestiges of the whole-genome duplications

Inferring the phenotype of the last common ancestor of living vertebrates is a challenging problem because of several unresolvable factors. They include the lack of reliable out-groups of living vertebrates, poor information about less fossilizable organs and specialized traits of phylogenetically important species, such as lampreys and hagfishes (e.g. secondary loss of vertebrae in adult hagfishes). These factors undermine the reliability of ancestral reconstruction by traditional character mapping approaches based on maximum parsimony. In this article, we formulate an approach to hypothesizing ancestral vertebrate phenotypes using information from the phylogenetic and functional properties of genes duplicated by genome expansions in early vertebrate evolution. We named the conjecture as ‘chronological reconstruction of ohnolog functions (CHROF)’. This CHROF conjecture raises the possibility that the last common ancestor of living vertebrates may have had more complex traits than currently thought.

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.

Two Tier Hox Collinearity Mediates Vertebrate Axial Patterning

A two tier mechanism mediates Hox collinearity. Besides the familiar collinear chromatin modification within each Hox cluster (nanocollinearity), there is also a macrocollinearity tier. Individual Hox clusters and individual cells are coordinated and synchronized to generate multiscale (macro and nano) collinearity in the early vertebrate embryo. Macro-collinearity is mediated by three non-cell autonomous Hox–Hox interactions. These mediate temporal collinearity in early NOM (non-organizer mesoderm), time space translation where temporal collinearity is translated to spatial collinearity along the early embryo’s main body axis and neural transformation, where Hox expression is copied monospecifically from NOM mesoderm to overlying neurectoderm in the late gastrula. Unlike nanocollinearity, which is Hox cluster restricted, axial macrocollinearity extends into the head and EAD domains, thus covering the whole embryonic anterior-posterior (A-P) axis. EAD: extreme anterior domain, the only A-P axial domain anterior to the head. The whole time space translation mechanism interacts with A-P signaling pathways via “decision points,” separating different domains on the axis.

Disruption of Chromosomal Architecture of cox2 Locus Sensitizes Lung Cancer Cells to Radiotherapy

Despite treatment of lung cancer with radiotherapy and chemotherapy, the survival rate of lung cancer patients remains poor. Previous studies demonstrated the importance of upregulation of inflammatory factors, such as cyclooxygenase 2 (cox2), in tumor tolerance. In the present study, we investigated the role of cox2 in radiosensitivity of lung cancer. Our results showed that the combination treatment of radiation with aspirin, an anti-inflammatory drug, induced a synergistic reduction of cell survival in A549 and H1299 lung cancer cells. In comparison with normal human lung fibroblasts (NHLFs), the cell viability was significantly decreased and the level of apoptosis was remarkably enhanced in A549 cells. Mechanistic studies revealed that the reduction of cox2 by aspirin in A549 and H1299 was caused by disruption of the chromosomal architecture of the cox2 locus. Moreover, the disruption of chromatin looping was mediated by the inhibition of nuclear translocation of p65 and decreased enrichment of p65 at cox2-regulatory elements. Importantly, disorganization of the chromosomal architecture of cox2 triggered A549 cells sensitive to γ-radiation by the induction of apoptosis. In conclusion, we present evidence of an effective therapeutic treatment targeting the epigenetic regulation of lung cancer and a potential strategy to overcome radiation resistance in cancer cells.

Mutations in LNPK, Encoding the Endoplasmic Reticulum Junction Stabilizer Lunapark, Cause a Recessive Neurodevelopmental Syndrome

The dynamic shape of the endoplasmic reticulum (ER) is a reflection of its wide variety of critical cell biological functions. Consequently, perturbation of ER-shaping proteins can cause a range of human phenotypes. Here, we describe three affected children (from two consanguineous families) who carry homozygous loss-of-function mutations in LNPK (previously known as KIAA1715); this gene encodes lunapark, which is proposed to serve as a curvature-stabilizing protein within tubular three-way junctions of the ER. All individuals presented with severe psychomotor delay, intellectual disability, hypotonia, epilepsy, and corpus callosum hypoplasia, and two of three showed mild cerebellar hypoplasia and atrophy. Consistent with a proposed role in neurodevelopmental disease, LNPK was expressed during brain development in humans and mice and was present in neurite-like processes in differentiating human neural progenitor cells. Affected cells showed the absence of full-length lunapark, aberrant ER structures, and increased luminal mass density. Together, our results implicate the ER junction stabilizer lunapark in establishing the corpus callosum.

ER-phagy requires Lnp1, a protein that stabilizes rearrangements of the ER network

The endoplasmic reticulum (ER) forms a contiguous network of tubules and sheets that is predominantly associated with the cell cortex in yeast. Upon treatment with rapamycin, the ER undergoes degradation by selective autophagy. This process, termed ER-phagy, requires Atg40, a selective autophagy receptor that localizes to the cortical ER. Here we report that ER-phagy also requires Lnp1, an ER membrane protein that normally resides at the three-way junctions of the ER network, where it serves to stabilize the network as it is continually remodeled. Rapamycin treatment increases the expression of Atg40, driving ER domains marked by Atg40 puncta to associate with Atg11, a scaffold protein needed to form autophagosomes. Although Atg40 largely localizes to the cortical ER, the autophagy machinery resides in the cell interior. The localization of Atg40 to sites of autophagosome formation is blocked in an lnp1Δ mutant or upon treatment of wild-type cells with the actin-depolymerizing drug Latrunculin A. This prevents the association of Atg40 with Atg11 and the packaging of the ER into autophagosomes. We propose that Lnp1 is needed to stabilize the actin-dependent remodeling of the ER that is essential for ER-phagy.

Cell-Specific PEAR1 Methylation Studies Reveal a Locus that Coordinates Expression of Multiple Genes

Chromosomal interactions connect distant enhancers and promoters on the same chromosome, activating or repressing gene expression. PEAR1 encodes the Platelet-Endothelial Aggregation Receptor 1, a contact receptor involved in platelet function and megakaryocyte and endothelial cell proliferation. PEAR1 expression during megakaryocyte differentiation is controlled by DNA methylation at its first CpG island. We identified a PEAR1 cell-specific methylation sensitive region in endothelial cells and megakaryocytes that showed strong chromosomal interactions with ISGL20L2, RRNAD1, MRLP24, HDGF and PRCC, using available promoter capture Hi-C datasets. These genes are involved in ribosome processing, protein synthesis, cell cycle and cell proliferation. We next studied the methylation and expression profile of these five genes in Human Umbilical Vein Endothelial Cells (HUVECs) and megakaryocyte precursors. While cell-specific PEAR1 methylation corresponded to variability in expression for four out of five genes, no methylation change was observed in their promoter regions across cell types. Our data suggest that PEAR1 cell-type specific methylation changes may control long distance interactions with other genes. Further studies are needed to show whether such interaction data might be relevant for the genome-wide association data that showed a role for non-coding PEAR1 variants in the same region and platelet function, platelet count and cardiovascular risk.

Assessing sufficiency and necessity of enhancer activities for gene expression and the mechanisms of transcription activation

Enhancers are important genomic regulatory elements directing cell type-specific transcription. They assume a key role during development and disease, and their identification and functional characterization have long been the focus of scientific interest. The advent of next-generation sequencing and clustered regularly interspaced short palindromic repeat (CRISPR)/Cas9-based genome editing has revolutionized the means by which we study enhancer biology. In this review, we cover recent developments in the prediction of enhancers based on chromatin characteristics and their identification by functional reporter assays and endogenous DNA perturbations. We discuss that the two latter approaches provide different and complementary insights, especially in assessing enhancer sufficiency and necessity for transcription activation. Furthermore, we discuss recent insights into mechanistic aspects of enhancer function, including findings about cofactor requirements and the role of post-translational histone modifications such as monomethylation of histone H3 Lys4 (H3K4me1). Finally, we survey how these approaches advance our understanding of transcription regulation with respect to promoter specificity and transcriptional bursting and provide an outlook covering open questions and promising developments.

Distinct epigenetic signatures elucidate enhancer-gene relationships that delineate CIMP and non-CIMP colorectal cancers

Epigenetic changes, like DNA methylation, affect gene expression and in colorectal cancer (CRC), a distinct phenotype called the CpG island methylator phenotype (“CIMP”) has significantly higher levels of DNA methylation at so-called “Type C loci” within the genome. We postulate that enhancer-gene pairs are coordinately controlled through DNA methylation in order to regulate the expression of key genes/biomarkers for a particular phenotype.

Firstly, we found 24 experimentally-validated enhancers (VISTA enhancer browser) that contained statistically significant (FDR-adjusted q-value of <0.01) differentially methylated regions (DMRs) (1000bp) in a study of CIMP versus non-CIMP CRCs. Of these, the methylation of 2 enhancers, 1702 and 1944, were found to be very well correlated with the methylation of the genes Wnt3A and IGDCC3, respectively, in two separate and independent datasets.

We show for the first time that there are indeed distinct and dynamic changes in the methylation pattern of specific enhancer-gene pairs in CRCs. Such a coordinated epigenetic event could be indicative of an interaction between (1) enhancer 1702 and Wnt3A and (2) enhancer 1944 and IGDCC3. Moreover, our study shows that the methylation patterns of these 2 enhancer-gene pairs can potentially be used as biomarkers to delineate CIMP from non-CIMP CRCs.

Deletion of a Long-Range Dlx5 Enhancer Disrupts Inner Ear Development in Mice

Distal enhancers are thought to play important roles in the spatiotemporal regulation of gene expression during embryonic development, but few predicted enhancer elements have been shown to affect transcription of their endogenous genes or to alter phenotypes when disrupted. Here, we demonstrate that a 123.6-kb deletion within the mouse Slc25a13 gene is associated with reduced transcription of Dlx5, a gene located 660 kb away. Mice homozygous for the Slc25a13 deletion mutation [named hyperspin (hspn)] have malformed inner ears and are deaf with balance defects, whereas previously reported Slc25a13 knockout mice showed no phenotypic abnormalities. Inner ears of Slc25a13^hspn/hspn mice have malformations similar to those of Dlx5^-/- embryos, and Dlx5 expression is severely reduced in the otocyst but not the branchial arches of Slc25a13^hspn/hspn embryos, indicating that the Slc25a13^hspn deletion affects otic-specific enhancers of Dlx5 In addition, transheterozygous Slc25a13^+/hspn Dlx5^+/- mice exhibit noncomplementation with inner ear dysmorphologies similar to those of Slc25a13^hspn/hspn and Dlx5^-/-embryos, verifying a cis-acting effect of the Slc25a13^hspn deletion on Dlx5 expression. CRISPR/Cas9-mediated deletions of putative enhancer elements located within the Slc25a13^hspn deleted region failed to phenocopy the defects of Slc25a13^hspn/hspn mice, suggesting the possibility of multiple enhancers with redundant functions. Our findings in mice suggest that analogous enhancer elements in the human SLC25A13 gene may regulate DLX5 expression and underlie the hearing loss that is associated with split-hand/-foot malformation 1 syndrome. Slc25a13^hspn/hspn mice provide a new animal model for studying long-range enhancer effects on Dlx5 expression in the developing inner ear.

Developmentally regulated higher-order chromatin interactions orchestrate B cell fate commitment

Genome organization in 3D nuclear-space is important for regulation of gene expression. However, the alterations of chromatin architecture that impinge on the B cell-fate choice of multi-potent progenitors are still unclear. By integrating in situ Hi-C analyses with epigenetic landscapes and genome-wide expression profiles, we tracked the changes in genome architecture as the cells transit from a progenitor to a committed state. We identified the genomic loci that undergo developmental switch between A and B compartments during B-cell fate determination. Furthermore, although, topologically associating domains (TADs) are stable, a significant number of TADs display structural alterations that are associated with changes in cis-regulatory interaction landscape. Finally, we demonstrate the potential roles for Ebf1 and its downstream factor, Pax5, in chromatin reorganization and transcription regulation. Collectively, our studies provide a general paradigm of the dynamic relationship between chromatin reorganization and lineage-specific gene expression pattern that dictates cell-fate determination.

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.

SBE6: a novel long-range enhancer involved in driving sonic hedgehog expression in neural progenitor cells

The expression of genes with key roles in development is under very tight spatial and temporal control, mediated by enhancers. A classic example of this is the sonic hedgehog gene (Shh), which plays a pivotal role in the proliferation, differentiation and survival of neural progenitor cells both in vivo and in vitro. Shh expression in the brain is tightly controlled by several known enhancers that have been identified through genetic, genomic and functional assays. Using chromatin profiling during the differentiation of embryonic stem cells to neural progenitor cells, here we report the identification of a novel long-range enhancer for Shh—Shh-brain-enhancer-6 (SBE6)—that is located 100 kb upstream of Shh and that is required for the proper induction of Shh expression during this differentiation programme. This element is capable of driving expression in the vertebrate brain. Our study illustrates how a chromatin-focused approach, coupled to in vivo testing, can be used to identify new cell-type specific cis-regulatory elements, and points to yet further complexity in the control of Shh expression during embryonic brain development.

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.