Shadow Enhancers Are Pervasive Features of Developmental Regulatory Networks

Enhancer Domains Predict Gene Pathogenicity and Inform Gene Discovery in Complex Disease

Non-coding transcriptional regulatory elements are critical for controlling the spatiotemporal expression of genes. Here, we demonstrate that the sizes and number of enhancers linked to a gene reflect its disease pathogenicity. Moreover, genes with redundant enhancer domains are depleted of cis-acting genetic variants that disrupt gene expression, and they are buffered against the effects of disruptive non-coding mutations. Our results demonstrate that dosage-sensitive genes have evolved a robustness to the disruptive effects of genetic variation by expanding their regulatory domains. This solves a puzzle about why genes associated with human disease are depleted of cis-eQTLs (cis-expression quantitative trait loci), suggesting that this relationship might complicate gene identification in causal genome-wide association studies (GWASs) using eQTL information, and establishes a framework for identifying non-coding regulatory variation with phenotypic consequences.

An inferred fitness consequence map of the rice genome

The extent to which sequence variation impacts plant fitness is poorly understood. High-resolution maps detailing the constraint acting on the genome, especially in regulatory sites, would be beneficial as functional annotation of noncoding sequences remains sparse. Here, we present a fitness consequence (fitCons) map for rice (Oryza sativa). We inferred fitCons scores (ρ) for 246 inferred genome classes derived from nine functional genomic and epigenomic datasets, including chromatin accessibility, messenger RNA/small RNA transcription, DNA methylation, histone modifications and engaged RNA polymerase activity. These were integrated with genome-wide polymorphism and divergence data from 1,477 rice accessions and 11 reference genome sequences in the Oryzeae. We found ρ to be multimodal, with ~9% of the rice genome falling into classes where more than half of the bases would probably have a fitness consequence if mutated. Around 2% of the rice genome showed evidence of weak negative selection, frequently at candidate regulatory sites, including a novel set of 1,000 potentially active enhancer elements. This fitCons map provides perspective on the evolutionary forces associated with genome diversity, aids in genome annotation and can guide crop breeding programs.

Coacting enhancers can have complementary functions within gene regulatory networks and promote canalization

Developmental genes are often regulated by multiple enhancers exhibiting similar spatiotemporal outputs, which are generally considered redundantly acting though few have been studied functionally. Using CRISPR-Cas9, we created deletions of two enhancers, brk5' and brk3', that drive similar but not identical expression of the gene brinker (brk) in early Drosophila embryos. Utilizing both in situ hybridization and quantitative mRNA analysis, we investigated the changes in the gene network state caused by the removal of one or both of the early acting enhancers. brk5' deletion generally phenocopied the gene mutant, including expansion of the BMP ligand decapentaplegic (dpp) as well as inducing variability in amnioserosa tissue cell number suggesting a loss of canalization. In contrast, brk3' deletion presented unique phenotypes including dorsal expansion of several ventrally expressed genes and a decrease in amnioserosa cell number. Similarly, deletions were made for two enhancers associated with the gene short-gastrulation (sog), sog.int and sog.dist, demonstrating that they also exhibit distinct patterning phenotypes and affect canalization. In summary, this study shows that similar gene expression driven by coacting enhancers can support distinct, and sometimes complementary, functions within gene regulatory networks and, moreover, that phenotypes associated with individual enhancer deletion mutants can provide insight into new gene functions.

Parallel evolution of ancient, pleiotropic enhancers underlies butterfly wing pattern mimicry

Color pattern mimicry in Heliconius butterflies is a classic case study of complex trait adaptation via selection on a few large effect genes. Association studies have linked color pattern variation to a handful of noncoding regions, yet the presumptive cis-regulatory elements (CREs) that control color patterning remain unknown. Here we combine chromatin assays, DNA sequence associations, and genome editing to functionally characterize 5 cis-regulatory elements of the color pattern gene optix We were surprised to find that the cis-regulatory architecture of optix is characterized by pleiotropy and regulatory fragility, where deletion of individual cis-regulatory elements has broad effects on both color pattern and wing vein development. Remarkably, we found orthologous cis-regulatory elements associate with wing pattern convergence of distantly related comimics, suggesting that parallel coevolution of ancestral elements facilitated pattern mimicry. Our results support a model of color pattern evolution in Heliconius where changes to ancient, multifunctional cis-regulatory elements underlie adaptive radiation.

An in silico analysis of robust but fragile gene regulation links enhancer length to robustness

Organisms must ensure that expression of genes is directed to the appropriate tissues at the correct times, while simultaneously ensuring that these gene regulatory systems are robust to perturbation. This idea is captured by a mathematical concept called r-robustness, which says that a system is robust to a perturbation in up to r - 1 randomly chosen parameters. r-robustness implies that the biological system has a small number of sensitive parameters and that this number can be used as a robustness measure. In this work we use this idea to investigate the robustness of gene regulation using a sequence level model of the Drosophila melanogaster gene even-skipped. We consider robustness with respect to mutations of the enhancer sequence and with respect to changes of the transcription factor concentrations. We find that gene regulation is r-robust with respect to mutations in the enhancer sequence and identify a number of sensitive nucleotides. In both natural and in silico predicted enhancers, the number of nucleotides that are sensitive to mutation correlates negatively with the length of the sequence, meaning that longer sequences are more robust. The exact degree of robustness obtained is dependent not only on DNA sequence, but also on the local concentration of regulatory factors. We find that gene regulation can be remarkably sensitive to changes in transcription factor concentrations at the boundaries of expression features, while it is robust to perturbation elsewhere.

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.

Functional Consequences of Chromosomal Rearrangements on Gene Expression: Not So Deleterious After All?

Chromosomes are folded and organized into topologically associating domains (TADs) which provide a framework for the interaction of enhancers with the promoter of their target gene(s). Structural rearrangements observed during evolution or in disease contexts suggest that changes in genome organization strongly affect gene expression and can have drastic phenotypic effects. In this review, I will discuss how recent genomic engineering experiments reveal a more contrasted picture, suggesting that TADs are important but not always essential for gene expression regulation.

Transcriptional autoregulation of zebrafish tbx6 is required for somite segmentation

The presumptive somite boundary in the presomitic mesoderm (PSM) is defined by the anterior border of the expression domain of Tbx6 protein. During somite segmentation, the expression domain of Tbx6 is regressed by Ripply-meditated degradation of Tbx6 protein. Although the expression of zebrafish tbx6 remains restricted to the PSM, the transcriptional regulation of tbx6 remains poorly understood. Here, we show that the expression of zebrafish tbx6 is maintained by transcriptional autoregulation. We find that a proximal-located cis-regulatory module, TR1, which contains two putative T-box sites, is required for somite segmentation in the intermediate body and for proper expression of segmentation genes. Embryos with deletion of TR1 exhibit significant reduction of tbx6 expression at the 12-somite stage, although its expression is initially observed. Additionally, Tbx6 is associated with TR1 and activates its own expression in the anterior PSM. Furthermore, the anterior expansion of tbx6 expression in ripply gene mutants is suppressed in a TR1-dependent manner. The results suggest that the autoregulatory loop of zebrafish tbx6 facilitates immediate removal of Tbx6 protein through termination of its own transcription at the anterior PSM.

Genomic dissection of enhancers uncovers principles of combinatorial regulation and cell type-specific wiring of enhancer-promoter contacts

Genomic binding of transcription factors, like the glucocorticoid receptor (GR), is linked to the regulation of genes. However, as we show here, GR binding is a poor predictor of GR-dependent gene regulation even when taking the 3D organization of the genome into account. To connect GR binding sites to the regulation of genes in the endogenous genomic context, we turned to genome editing. By deleting GR binding sites, individually or in combination, we uncovered how cooperative interactions between binding sites contribute to the regulation of genes. Specifically, for the GR target gene GILZ, we show that the simultaneous presence of a cluster of GR binding sites is required for the activity of an individual enhancer and that the GR-dependent regulation of GILZ depends on multiple GR-bound enhancers. Further, by deleting GR binding sites that are shared between different cell types, we show how cell type-specific genome organization and enhancer-blocking can result in cell type-specific wiring of promoter–enhancer contacts. This rewiring allows an individual GR binding site shared between different cell types to direct the expression of distinct transcripts and thereby contributes to the cell type-specific consequences of glucocorticoid signaling.

Deletion of a Csf1r enhancer selectively impacts CSF1R expression and development of tissue macrophage populations

The proliferation, differentiation and survival of mononuclear phagocytes depend on signals from the receptor for macrophage colony-stimulating factor, CSF1R. The mammalian Csf1r locus contains a highly conserved super-enhancer, the fms-intronic regulatory element (FIRE). Here we show that genomic deletion of FIRE in mice selectively impacts CSF1R expression and tissue macrophage development in specific tissues. Deletion of FIRE ablates macrophage development from murine embryonic stem cells. Csf1rΔFIRE/ΔFIRE mice lack macrophages in the embryo, brain microglia and resident macrophages in the skin, kidney, heart and peritoneum. The homeostasis of other macrophage populations and monocytes is unaffected, but monocytes and their progenitors in bone marrow lack surface CSF1R. Finally, Csf1rΔFIRE/ΔFIRE mice are healthy and fertile without the growth, neurological or developmental abnormalities reported in Csf1r-/- rodents. Csf1rΔFIRE/ΔFIRE mice thus provide a model to explore the homeostatic, physiological and immunological functions of tissue-specific macrophage populations in adult animals.

Structural variation in the 3D genome

Structural and quantitative chromosomal rearrangements, collectively referred to as structural variation (SV), contribute to a large extent to the genetic diversity of the human genome and thus are of high relevance for cancer genetics, rare diseases and evolutionary genetics. Recent studies have shown that SVs can not only affect gene dosage but also modulate basic mechanisms of gene regulation. SVs can alter the copy number of regulatory elements or modify the 3D genome by disrupting higher-order chromatin organization such as topologically associating domains. As a result of these position effects, SVs can influence the expression of genes distant from the SV breakpoints, thereby causing disease. The impact of SVs on the 3D genome and on gene expression regulation has to be considered when interpreting the pathogenic potential of these variant types.

Multi-enhancer transcriptional hubs confer phenotypic robustness

We previously showed in Drosophila melanogaster embryos that low-affinity Ultrabithorax (Ubx)-responsive shavenbaby (svb) enhancers drive expression using localized transcriptional environments and that active svb enhancers on different chromosomes tended to colocalize (Tsai et al., 2017). Here, we test the hypothesis that these multi-enhancer 'hubs' improve phenotypic resilience to stress by buffering against decreases in transcription factor concentrations and transcriptional output. Deleting a redundant enhancer from the svb locus led to reduced trichome numbers in embryos raised at elevated temperatures. Using high-resolution fluorescence microscopy, we observed lower Ubx concentration and transcriptional output in this deletion allele. Transcription sites of the full svb cis-regulatory region inserted into a different chromosome colocalized with the svb locus, increasing Ubx concentration, the transcriptional output of svb, and partially rescuing the phenotype. Thus, multiple enhancers could reinforce a local transcriptional hub to buffer against environmental stresses and genetic perturbations, providing a mechanism for phenotypical robustness.

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.

Shaping the nebulous enhancer in the era of high-throughput assays and genome editing

Since the 1st discovery of transcriptional enhancers in 1981, their textbook definition has remained largely unchanged in the past 37 years. With the emergence of high-throughput assays and genome editing, which are switching the paradigm from bottom-up discovery and testing of individual enhancers to top-down profiling of enhancer activities genome-wide, it has become increasingly evidenced that this classical definition has left substantial gray areas in different aspects. Here we survey a representative set of recent research articles and report the definitions of enhancers they have adopted. The results reveal that a wide spectrum of definitions is used usually without the definition stated explicitly, which could lead to difficulties in data interpretation and downstream analyses. Based on these findings, we discuss the practical implications and suggestions for future studies.

Redundant and Cryptic Enhancer Activities of the Drosophila yellow Gene

Cis-regulatory sequences known as enhancers play a key role in regulating gene expression. Evolutionary changes in these DNA sequences contribute to phenotypic evolution. The Drosophila yellow gene, which is required for pigmentation, has emerged as a model system for understanding how cis-regulatory sequences evolve, providing some of the most detailed insights available into how activities of orthologous enhancers have diverged between species. Here, we examine the evolution of yellow cis-regulatory sequences on a broader scale, by comparing the distribution and function of yellow enhancer activities throughout the 5' intergenic and intronic sequences of Drosophila melanogaster, D. pseudoobscura, and D. willistoni We find that cis-regulatory sequences driving expression in a particular tissue are not as modular as previously described, but rather have many redundant and cryptic enhancer activities distributed throughout the regions surveyed. Interestingly, cryptic enhancer activities of sequences from one species often drove patterns of expression observed in other species, suggesting that the frequent evolutionary changes in yellow expression observed among Drosophila species may be facilitated by gaining and losing repression of preexisting cis-regulatory sequences.

Signal Integration by Shadow Enhancers and Enhancer Duplications Varies across the Drosophila Embryo

Transcription of developmental genes is controlled by multiple enhancers. Frequently, more than one enhancer can activate transcription from the same promoter in the same cells. How is regulatory information from multiple enhancers combined to determine the overall expression output? We measure nascent transcription driven by a pair of shadow enhancers, each enhancer of the pair separately, and each duplicated, using live imaging in Drosophila embryos. This set of constructs allows us to quantify the input-output function describing signal integration by two enhancers. We show that signal integration performed by these shadow enhancers and duplications varies across the expression pattern, implying that how their activities are combined depends on the transcriptional regulators bound to the enhancers in different parts of the embryo. Characterizing signal integration by multiple enhancers is a critical step in developing conceptual and computational models of gene expression at the locus level, where multiple enhancers control transcription together.

The cis-Regulatory Atlas of the Mouse Immune System

A complete chart of cis-regulatory elements and their dynamic activity is necessary to understand the transcriptional basis of differentiation and function of an organ system. We generated matched epigenome and transcriptome measurements in 86 primary cell types that span the mouse immune system and its differentiation cascades. This breadth of data enable variance components analysis that suggests that genes fall into two distinct classes, controlled by either enhancer- or promoter-driven logic, and multiple regression that connects genes to the enhancers that regulate them. Relating transcription factor (TF) expression to the genome-wide accessibility of their binding motifs classifies them as predominantly openers or closers of local chromatin accessibility, pinpointing specific cis-regulatory elements where binding of given TFs is likely functionally relevant, validated by chromatin immunoprecipitation sequencing (ChIP-seq). Overall, this cis-regulatory atlas provides a trove of information on transcriptional regulation through immune differentiation and a foundational scaffold to define key regulatory events throughout the immunological genome.

Encounters across networks: Windows into principles of genomic regulation

Gene regulatory networks account for the ability of the genome to program development in complex multi-cellular organisms. Such networks are based on principles of gene regulation by combinations of transcription factors that bind to specific cis-regulatory DNA sites to activate transcription. These cis-regulatory regions mediate logic processing at each network node, enabling progressive increases in organismal complexity with development. Gene regulatory network explanations of development have been shown to account for patterning and cell type diversification in fly and sea urchin embryonic systems, where networks are characterized by fast coupling between transcriptional inputs and changes in target gene transcription rates, and crucial cis-regulatory elements are concentrated relatively close to the protein coding sequences of the target genes, thus facilitating their identification. Stem cell-based development in post-embryonic mammalian systems also depends on gene networks, but differs from the fly and sea urchin systems. First, the number of regulatory elements per gene and the distances between regulatory elements and the genes they control are considerably larger, forcing searches via genome-wide transcription factor binding surveys rather than functional assays. Second, the intrinsic timing of network state transitions can be slowed considerably by the need to undo stem-cell chromatin configurations, which presumably add stability to stem-cell states but retard responses to transcription factor changes during differentiation. The dispersed, partially redundant cis-regulatory systems controlling gene expression and the slow state transition kinetics in these systems already reveal new insights and opportunities to extend understanding of the repertoire of gene networks and regulatory system logic.

Evolution of embryonic cis-regulatory landscapes between divergent Phallusia and Ciona ascidians

Ascidian species of the Phallusia and Ciona genera are distantly related, their last common ancestor dating several hundred million years ago. Although their genome sequences have extensively diverged since this radiation, Phallusia and Ciona species share almost identical early morphogenesis and stereotyped cell lineages. Here, we explored the evolution of transcriptional control between P. mammillata and C. robusta. We combined genome-wide mapping of open chromatin regions in both species with a comparative analysis of the regulatory sequences of a test set of 10 pairs of orthologous early regulatory genes with conserved expression patterns. We find that ascidian chromatin accessibility landscapes obey similar rules as in other metazoa. Open-chromatin regions are short, highly conserved within each genus and cluster around regulatory genes. The dynamics of chromatin accessibility and closest-gene expression are strongly correlated during early embryogenesis. Open-chromatin regions are highly enriched in cis-regulatory elements: 73% of 49 open chromatin regions around our test genes behaved as either distal enhancers or proximal enhancer/promoters following electroporation in Phallusia eggs. Analysis of this datasets suggests a pervasive use in ascidians of "shadow" enhancers with partially overlapping activities. Cross-species electroporations point to a deep conservation of both the trans-regulatory logic between these distantly-related ascidians and the cis-regulatory activities of individual enhancers. Finally, we found that the relative order and approximate distance to the transcription start site of open chromatin regions can be conserved between Ciona and Phallusia species despite extensive sequence divergence, a property that can be used to identify orthologous enhancers, whose regulatory activity can partially diverge.

Dichotomy in redundant enhancers points to presence of initiators of gene regulation

Background

The regulatory landscape of a gene locus often consists of several functionally redundant enhancers establishing phenotypic robustness and evolutionary stability of its regulatory program. However, it is unclear what mechanisms are employed by redundant enhancers to cooperatively orchestrate gene expression.

Results

By comparing redundant enhancers to single enhancers (enhancers present in a single copy in a gene locus), we observed that the DNA sequence encryption differs between these two classes of enhancers, suggesting a difference in their regulatory mechanisms. Initiator enhancers, which are a subset of redundant enhancers and show similar sequence encryption to single enhancers, differ from the rest of redundant enhancers in their sequence encryption, evolutionary conservation and proximity to target genes. Genes hosting initiator enhancers in their loci feature elevated levels of expression. Initiator enhancers show a high level of 3D chromatin contacts with both transcription start sites and regular enhancers, suggesting their roles as primary activators and intermediate catalysts of gene expression, through which the regulatory signals of redundant enhancers are propagated to the target genes. In addition, GWAS and eQTLs variants are significantly enriched in initiator enhancers compared to redundant enhancers, suggesting a key functional role these sequences play in gene regulation.

Conclusions

The specific characteristics and widespread abundance of initiator enhancers advocate for a possible universal hierarchical mechanism of tissue-specific gene regulation involving multiple redundant enhancers acting through initiator enhancers.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5335-0) contains supplementary material, which is available to authorized users.

Human sex reversal is caused by duplication or deletion of core enhancers upstream of SOX9

Disorders of sex development (DSDs) are conditions affecting development of the gonads or genitalia. Variants in two key genes, SRY and its target SOX9, are an established cause of 46,XY DSD, but the genetic basis of many DSDs remains unknown. SRY-mediated SOX9 upregulation in the early gonad is crucial for testis development, yet the regulatory elements underlying this have not been identified in humans. Here, we identified four DSD patients with overlapping duplications or deletions upstream of SOX9. Bioinformatic analysis identified three putative enhancers for SOX9 that responded to different combinations of testis-specific regulators. All three enhancers showed synergistic activity and together drive SOX9 in the testis. This is the first study to identify SOX9 enhancers that, when duplicated or deleted, result in 46,XX or 46,XY sex reversal, respectively. These enhancers provide a hitherto missing link by which SRY activates SOX9 in humans, and establish SOX9 enhancer mutations as a significant cause of DSD.

SRY and its target SOX9 are known key determinants in testis development. Here the authors by studying duplications and deletions upstream of SOX9 from patient samples with disorders of sex development (DSD) reveal enhancers for SOX9 critical for human sex development and DSD.

Promoter-enhancer looping and shadow enhancers of the mouse αA-crystallin locus

Gene regulation by enhancers is important for precise temporal and spatial gene expression. Enhancers can drive gene expression regardless of their location, orientation or distance from the promoter. Changes in chromatin conformation and chromatin looping occur to bring the promoter and enhancers into close proximity. αA-crystallin ranks among one of the most abundantly expressed genes and proteins in the mammalian lens. The αA-crystallin locus is characterized by a 16 kb chromatin domain marked by two distal enhancers, 5′ DCR1 and 3′ DCR3. Here we used chromatin conformation capture (3C) analysis and transgenic approaches to analyze temporal control of the mouse αA-crystallin gene. We find that DCR1 is necessary, but not sufficient alone to drive expression at E10.5 in the mouse lens pit. Chromatin looping revealed interaction between the promoter and the region 3′ to DCR1, identifying a novel enhancer region in the αA-crystallin locus. We determined that this novel enhancer region, DCR1S, recapitulates the temporal control by DCR1. Acting as shadow enhancers, DCR1 and DCR1S are able to control expression in the lens vesicle at E11.5. It remains to be elucidated however, which region of the αA-crystallin locus is responsible for expression in the lens pit at E10.5.

Summary: The αA-crystallin ranks amongst the most highly expressed tissue-specific genes. It is an advantageous model system to probe both promoter-enhancer looping and to identify distal enhancers and their temporal/spatial activities.

Multiple inputs into a posterior-specific regulatory network in the Ciona notochord

The gene regulatory networks underlying Ciona notochord fate specification and differentiation have been extensively investigated, but the regulatory basis for regionalized expression within the notochord is not understood. Here we identify three notochord-expressed genes, C11.331, C12.115 and C8.891, with strongly enriched expression in the secondary notochord cells at the posterior tip of the tail. C11.331 and C12.115 share a distinctive expression pattern that is highly enriched in the secondary notochord lineage but also graded within that lineage with the strongest expression at the posterior tip. Both genes show similar responses to pharmacological perturbations of Wnt and FGF signaling, consistent with an important role for Wnt and FGF ligands expressed at the tail tip. Reporter analysis indicates that the C11.331 cis-regulatory regions are extensively distributed, with multiple non-overlapping regions conferring posterior notochord-enriched expression. Fine-scale analysis of a minimal cis-regulatory module identifies discrete positive and negative elements including a strong silencer. Truncation of the silencer region leads to increased expression in the primary notochord, indicating that C11.331 expression is influenced by putative regulators of primary versus secondary notochord fate. The minimal CRM contains predicted ETS, GATA, LMX and Myb sites, all of which lead to reduced expression in secondary notochord when mutated. These results show that the posterior-enriched notochord expression of C11.331 depends on multiple inputs, including Wnt and FGF signals from the tip of the tail, multiple notochord-specific regulators, and yet-to-be identified regulators of regional identity within the notochord.

Modulating transcription factor activity: Interfering with protein-protein interaction networks

Biophysical parameters that govern transcription factors activity are binding locations across the genome, dwelling time at these regulatory elements and specific protein-protein interactions. Most molecular strategies used to develop small compounds that block transcription factors activity have been based on biochemistry and cell biology methods that that do not take into consideration these key biophysical features. Here, we review the advance in the field of transcription factor biology and describe how their interactome and transcriptional regulation on a genome wide scale have been deciphered. We suggest that this new knowledge has the potential to be used to implement innovative research drug discovery program.

Cooperation, cis-interactions, versatility and evolutionary plasticity of multiple cis-acting elements underlie krox20 hindbrain regulation

Cis-regulation plays an essential role in the control of gene expression, and is particularly complex and poorly understood for developmental genes, which are subject to multiple levels of modulation. In this study, we performed a global analysis of the cis-acting elements involved in the control of the zebrafish developmental gene krox20. krox20 encodes a transcription factor required for hindbrain segmentation and patterning, a morphogenetic process highly conserved during vertebrate evolution. Chromatin accessibility analysis reveals a cis-regulatory landscape that includes 6 elements participating in the control of initiation and autoregulatory aspects of krox20 hindbrain expression. Combining transgenic reporter analyses and CRISPR/Cas9-mediated mutagenesis, we assign precise functions to each of these 6 elements and provide a comprehensive view of krox20 cis-regulation. Three important features emerged. First, cooperation between multiple cis-elements plays a major role in the regulation. Cooperation can surprisingly combine synergy and redundancy, and is not restricted to transcriptional enhancer activity (for example, 4 distinct elements cooperate through different modes to maintain autoregulation). Second, several elements are unexpectedly versatile, which allows them to be involved in different aspects of control of gene expression. Third, comparative analysis of the elements and their activities in several vertebrate species reveals that this versatility is underlain by major plasticity across evolution, despite the high conservation of the gene expression pattern. These characteristics are likely to be of broad significance for developmental genes.

Diverse Cis-Regulatory Mechanisms Contribute to Expression Evolution of Tandem Gene Duplicates

Pairs of duplicated genes generally display a combination of conserved expression patterns inherited from their unduplicated ancestor and newly acquired domains. However, how the cis-regulatory architecture of duplicated loci evolves to produce these expression patterns is poorly understood. We have directly examined the gene-regulatory evolution of two tandem duplicates, the Drosophila Ly6 genes CG9336 and CG9338, which arose at the base of the drosophilids between 40 and 60 Ma. Comparing the expression patterns of the two paralogs in four Drosophila species with that of the unduplicated ortholog in the tephritid Ceratitis capitata, we show that they diverged from each other as well as from the unduplicated ortholog. Moreover, the expression divergence appears to have occurred close to the duplication event and also more recently in a lineage-specific manner. The comparison of the tissue-specific cis-regulatory modules (CRMs) controlling the paralog expression in the four Drosophila species indicates that diverse cis-regulatory mechanisms, including the novel tissue-specific enhancers, differential inactivation, and enhancer sharing, contributed to the expression evolution. Our analysis also reveals a surprisingly variable cis-regulatory architecture, in which the CRMs driving conserved expression domains change in number, location, and specificity. Altogether, this study provides a detailed historical account that uncovers a highly dynamic picture of how the paralog expression patterns and their underlying cis-regulatory landscape evolve. We argue that our findings will encourage studying cis-regulatory evolution at the whole-locus level to understand how interactions between enhancers and other regulatory levels shape the evolution of gene expression.

The Conserved Sonic Hedgehog Limb Enhancer Consists of Discrete Functional Elements that Regulate Precise Spatial Expression

Expression of sonic hedgehog (Shh) in the limb bud is regulated by an enhancer called the zone of polarizing activity regulatory sequence (ZRS), which, in evolution, belongs to an ancient group of highly conserved cis regulators found in all classes of vertebrates. Here, we examined the endogenous ZRS in mice, using genome editing to establish the relationship between enhancer composition and embryonic phenotype. We show that enhancer activity is a consolidation of distinct activity domains. Spatial restriction of Shh expression is mediated by a discrete repressor module, whereas levels of gene expression are controlled by large overlapping domains containing varying numbers of HOXD binding sites. The number of HOXD binding sites regulates expression levels incrementally. Substantial portions of conserved sequence are dispensable, indicating the presence of sequence redundancy. We propose a collective model for enhancer activity in which function is an integration of discrete expression activities and redundant components that drive robust expression.

•

The ancient vertebrate enhancer, the ZRS, shows sequence plasticity

•

Discrete regulatory activities are assigned to specific sites in the enhancer

•

The number of HOXD binding sites determines the level of Shh expression

•

Robust expression is a collective of regulatory and redundant information

Disrupting the three-dimensional regulatory topology of the Pitx1 locus results in overtly normal development

Developmental gene expression patterns are orchestrated by thousands of distant-acting transcriptional enhancers. However, identifying enhancers essential for the expression of their target genes has proven challenging. Maps of long-range regulatory interactions may provide the means to identify enhancers crucial for developmental gene expression. To investigate this hypothesis, we used circular chromosome conformation capture coupled with interaction maps in the mouse limb to characterize the regulatory topology of Pitx1, which is essential for hindlimb development. We identified a robust hindlimb-specific interaction between Pitx1 and a putative hindlimb-specific enhancer. To interrogate the role of this interaction in Pitx1 regulation, we used genome editing to delete this enhancer in mouse. Although deletion of the enhancer completely disrupts the interaction, Pitx1 expression in the hindlimb is only mildly affected, without any detectable compensatory interactions between the Pitx1 promoter and potentially redundant enhancers. Pitx1 enhancer null mice did not exhibit any of the characteristic morphological defects of the Pitx1-/- mutant. Our results suggest that robust, tissue-specific physical interactions at essential developmental genes have limited predictive value for identifying enhancer mutations with strong loss-of-function phenotypes.

Evolutionary emergence of the rac3b/rfng/sgca regulatory cluster refined mechanisms for hindbrain boundaries formation

Developmental programs often rely on parallel morphogenetic mechanisms that guarantee precise tissue architecture. While redundancy constitutes an obvious selective advantage, little is known on how novel morphogenetic mechanisms emerge during evolution. In zebrafish, rhombomeric boundaries behave as an elastic barrier, preventing cell intermingling between adjacent compartments. Here, we identify the fundamental role of the small-GTPase Rac3b in actomyosin cable assembly at hindbrain boundaries. We show that the novel rac3b/rfng/sgca regulatory cluster, which is specifically expressed at the boundaries, emerged in the Ostariophysi superorder by chromosomal rearrangement that generated new cis-regulatory interactions. By combining 4C-seq, ATAC-seq, transgenesis, and CRISPR-induced deletions, we characterized this regulatory domain, identifying hindbrain boundary-specific cis-regulatory elements. Our results suggest that the capacity of boundaries to act as an elastic mesh for segregating rhombomeric cells evolved by cooption of critical genes to a novel regulatory block, refining the mechanisms for hindbrain segmentation.

Dynamic evolution of regulatory element ensembles in primate CD4+ T cells

How evolutionary changes at enhancers affect the transcription of target genes remains an important open question. Previous comparative studies of gene expression have largely measured the abundance of messenger RNA, which is affected by post-transcriptional regulatory processes, hence limiting inferences about the mechanisms underlying expression differences. Here, we directly measured nascent transcription in primate species, allowing us to separate transcription from post-transcriptional regulation. We used precision run-on and sequencing to map RNA polymerases in resting and activated CD4+ T cells in multiple human, chimpanzee and rhesus macaque individuals, with rodents as outgroups. We observed general conservation in coding and non-coding transcription, punctuated by numerous differences between species, particularly at distal enhancers and non-coding RNAs. Genes regulated by larger numbers of enhancers are more frequently transcribed at evolutionarily stable levels, despite reduced conservation at individual enhancers. Adaptive nucleotide substitutions are associated with lineage-specific transcription and at one locus, SGPP2, we predict and experimentally validate that multiple substitutions contribute to human-specific transcription. Collectively, our findings suggest a pervasive role for evolutionary compensation across ensembles of enhancers that jointly regulate target genes.

A damped oscillator imposes temporal order on posterior gap gene expression in Drosophila

Insects determine their body segments in two different ways. Short-germband insects, such as the flour beetle Tribolium castaneum, use a molecular clock to establish segments sequentially. In contrast, long-germband insects, such as the vinegar fly Drosophila melanogaster, determine all segments simultaneously through a hierarchical cascade of gene regulation. Gap genes constitute the first layer of the Drosophila segmentation gene hierarchy, downstream of maternal gradients such as that of Caudal (Cad). We use data-driven mathematical modelling and phase space analysis to show that shifting gap domains in the posterior half of the Drosophila embryo are an emergent property of a robust damped oscillator mechanism, suggesting that the regulatory dynamics underlying long- and short-germband segmentation are much more similar than previously thought. In Tribolium, Cad has been proposed to modulate the frequency of the segmentation oscillator. Surprisingly, our simulations and experiments show that the shift rate of posterior gap domains is independent of maternal Cad levels in Drosophila. Our results suggest a novel evolutionary scenario for the short- to long-germband transition and help explain why this transition occurred convergently multiple times during the radiation of the holometabolan insects.

Different insect species exhibit one of two distinct modes of determining their body segments (known as segmentation) during development: they either use a molecular oscillator to position segments sequentially, or they generate segments simultaneously through a hierarchical gene-regulatory cascade. The sequential mode is ancestral, while the simultaneous mode has been derived from it independently several times during evolution. In this paper, we present evidence suggesting that simultaneous segmentation also involves an oscillator in the posterior end of the embryo of the vinegar fly, Drosophila melanogaster. This surprising result indicates that both modes of segment determination are much more similar than previously thought. Such similarity provides an important step towards our understanding of the frequent evolutionary transitions observed between sequential and simultaneous segmentation.

Comparative Analysis of Immune Cells Reveals a Conserved Regulatory Lexicon

Most well-characterized enhancers are deeply conserved. In contrast, genome-wide comparative studies of steady-state systems showed that only a small fraction of active enhancers are conserved. To better understand conservation of enhancer activity, we used a comparative genomics approach that integrates temporal expression and epigenetic profiles in an innate immune system. We found that gene expression programs diverge among mildly induced genes, while being highly conserved for strongly induced genes. The fraction of conserved enhancers varies greatly across gene expression programs, with induced genes and early-response genes, in particular, being regulated by a higher fraction of conserved enhancers. Clustering of conserved accessible DNA sequences within enhancers resulted in over 60 sequence motifs including motifs for known factors, as well as many with unknown function. We further show that the number of instances of these motifs is a strong predictor of the responsiveness of a gene to pathogen detection.

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.

Enhancer redundancy provides phenotypic robustness in mammalian development

Distant-acting tissue-specific enhancers vastly outnumber protein-coding genes in mammalian genomes, but the functional significance of this regulatory complexity remains insufficiently understood^1,2. Here we show that the pervasive presence of multiple enhancers with similar activities near the same gene confers phenotypic robustness to loss-of-function mutations in individual enhancers. We used genome editing to create 23 mouse deletion lines and inter-crosses, including both single and combinatorial enhancer deletions at seven distinct loci required for limb development. Surprisingly, none of ten deletions of individual enhancers caused noticeable changes in limb morphology. In contrast, removal of pairs of limb enhancers near the same gene resulted in discernible phenotypes, indicating that enhancers function redundantly in establishing normal morphology. In a genetic background sensitized by reduced baseline expression of the target gene, even single enhancer deletions caused limb abnormalities, suggesting that functional redundancy is conferred by additive effects of enhancers on gene expression levels. A genome-wide analysis integrating epigenomic and transcriptomic data from 29 developmental mouse tissues revealed that mammalian genes are very commonly associated with multiple enhancers that have similar spatiotemporal activity. Systematic exploration of three representative developmental structures (limb, brain, heart) uncovered more than a thousand cases in which five or more enhancers with redundant activity patterns were found near the same gene. Taken together, our data indicate that enhancer redundancy is a remarkably widespread feature of mammalian genomes and provides an effective regulatory buffer preventing deleterious phenotypic consequences upon loss of individual enhancers.

Regulatory Principles Governing Tissue Specificity of Developmental Enhancers

Transcriptional enhancers are short segments of genomic DNA (50 bp to 1 kb in length) that can work over long distances (≥1 Mb) to regulate gene expression in specific cells and tissues. Genomic assays have identified on the order of 400,000 to one million putative enhancers in the human genome (e.g., ENCODE Consortium). This suggests that a typical gene is regulated by tens of enhancers, ensuring stringent regulation of gene expression in response to a variety of intrinsic and external signals. Despite the discovery of the first transcriptional enhancer more than 30 years ago, we know surprisingly little about how enhancers regulate gene expression. In particular, the relationship between primary DNA sequence and enhancer specificity remains obscure. Here we summarize recent high-throughput studies in whole embryos aimed at the systematic identification of the sequence and organizational constraints underlying enhancer function and specificity.

Transcription factor-dependent 'anti-repressive' mammalian enhancers exclude H3K27me3 from extended genomic domains

Compacted chromatin and nucleosomes are known barriers to gene expression; the nature and relative importance of other transcriptional constraints remain unclear, especially at distant enhancers. Polycomb repressor complex 2 (PRC2) places the histone mark H3K27me3 predominantly at promoters, where its silencing activity is well documented. In adult tissues, enhancers lack H3K27me3, and it is unknown whether intergenic H3K27me3 deposits affect nearby genes. In primary intestinal villus cells, we identified hundreds of tissue-restricted enhancers that require the transcription factor (TF) CDX2 to prevent the incursion of H3K27me3 from adjoining areas of elevated basal marking into large well-demarcated genome domains. Similarly, GATA1-dependent enhancers exclude H3K27me3 from extended regions in erythroid blood cells. Excess intergenic H3K27me3 in both TF-deficient tissues is associated with extreme mRNA deficits, which are significantly rescued in intestinal cells lacking PRC2. Explaining these observations, enhancers show TF-dependent binding of the H3K27 demethylase KDM6A. Thus, in diverse cell types, certain genome regions far from promoters accumulate H3K27me3, and optimal gene expression depends on enhancers clearing this repressive mark. These findings reveal new "anti-repressive" function for hundreds of tissue-specific enhancers.

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.

Multiple enhancer regions govern the transcription of CCN2 during embryonic development

CCN2 is a critical matricellular protein that is expressed in several cells with major implications in physiology and different pathologies. However, the transcriptional regulation of this gene remains obscure. We used the Encyclopaedia of DNA Elements browser (ENCODE) to visualise the region spanning from 300 kb upstream to the CCN2 start site in silico in order to identify enhancer regions that regulate transcription of this gene. Selection was based on three criteria associated with enhancer regions: 1) H3K4me1 and H3K27ac histone modifications, 2) DNase I hypersensitivity of chromatin and 3) inter-species conservation. Reporter constructs were created with sequences spanning each of the regions of interest placed upstream of an Hsp68 silent proximal promoter sequence in order to drive the expression of β-galactosidase transgene. Each of these constructs was subsequently used to create transgenic mice in which reporter gene production was assessed at the E15.5 developmental stage. Four functional enhancers were identified, with each driving distinct, tissue-specific patterns of transgene expression. An enhancer located -100 kb from the CCN2 transcription start site facilitated expression within vascular tissue. An enhancer -135 kb upstream of CCN2 drove expression within the articular chondrocytes of synovial joints. The other two enhancers, located at -198 kb and -229 kb, mediated transgene expression within dermal fibroblasts, however the most prevalent activity was found within hypertrophic chondrocytes and periosteal tissue, respectively. These findings suggest that the global expression of CCN2 during development results from the activity of several tissue-specific enhancer regions in addition to proximal regulatory elements that have previously been demonstrated to drive transcription of the gene during development.

Transcriptional precision and accuracy in development: from measurements to models and mechanisms

During development, genes are transcribed at specific times, locations and levels. In recent years, the emergence of quantitative tools has significantly advanced our ability to measure transcription with high spatiotemporal resolution in vivo. Here, we highlight recent studies that have used these tools to characterize transcription during development, and discuss the mechanisms that contribute to the precision and accuracy of the timing, location and level of transcription. We attempt to disentangle the discrepancies in how physicists and biologists use the term ‘precision' to facilitate interactions using a common language. We also highlight selected examples in which the coupling of mathematical modeling with experimental approaches has provided important mechanistic insights, and call for a more expansive use of mathematical modeling to exploit the wealth of quantitative data and advance our understanding of animal transcription.

Summary: This Review highlights how high-resolution quantitative tools and theoretical models have formed our current view of the mechanisms determining precision and accuracy in the timing, location and level of transcription in the Drosophila embryo.

Multiplex Enhancer Interference Reveals Collaborative Control of Gene Regulation by Estrogen Receptor α-Bound Enhancers

Multiple regulatory regions have the potential to regulate a single gene, yet how these elements combine to affect gene expression remains unclear. To uncover the combinatorial relationships between enhancers, we developed Enhancer-interference (Enhancer-i), a CRISPR interference-based approach that uses 2 different repressive domains, KRAB and SID, to prevent enhancer activation simultaneously at multiple regulatory regions. We applied Enhancer-i to promoter-distal estrogen receptor α binding sites (ERBS), which cluster around estradiol-responsive genes and therefore may collaborate to regulate gene expression. Targeting individual sites revealed predominant ERBS that are completely required for the transcriptional response, indicating a lack of redundancy. Simultaneous interference of different ERBS combinations identified supportive ERBS that contribute only when predominant sites are active. Using mathematical modeling, we find strong evidence for collaboration between predominant and supportive ERBS. Overall, our findings expose a complex functional hierarchy of enhancers, where multiple loci bound by the same transcription factor combine to fine-tune the expression of target genes.

Speed regulation of genetic cascades allows for evolvability in the body plan specification of insects

During the anterior-posterior fate specification of insects, anterior fates arise in a nonelongating tissue (called the "blastoderm"), and posterior fates arise in an elongating tissue (called the "germband"). However, insects differ widely in the extent to which anterior-posterior fates are specified in the blastoderm versus the germband. Here we present a model in which patterning in both the blastoderm and germband of the beetle Tribolium castaneum is based on the same flexible mechanism: a gradient that modulates the speed of a genetic cascade of gap genes, resulting in the induction of sequential kinematic waves of gap gene expression. The mechanism is flexible and capable of patterning both elongating and nonelongating tissues, and hence converting blastodermal to germband fates and vice versa. Using RNAi perturbations, we found that blastodermal fates could be shifted to the germband, and germband fates could be generated in a blastoderm-like morphology. We also suggest a molecular mechanism underlying our model, in which gradient levels regulate the switch between two enhancers: One enhancer is responsible for sequential gene activation, and the other is responsible for freezing temporal rhythms into spatial patterns. This model is consistent with findings in Drosophila melanogaster, where gap genes were found to be regulated by two nonredundant "shadow" enhancers.

Molecular and functional genetics of the proopiomelanocortin gene, food intake regulation and obesity

A specter is haunting the world, the specter of obesity. During the last decade, this pandemia has skyrocketed threatening children, adolescents and lower income families worldwide. Although driven by an increase in the consumption of ultraprocessed edibles of poor nutritional value, the obesogenic changes in contemporary human lifestyle affect people differently, revealing that some individuals are more prone to develop increased adiposity. During the last years, we performed a variety of genetic, evolutionary, biochemical and behavioral experiments that allowed us to understand how a group of neurons present in the arcuate nucleus of the hypothalamus regulate the expression of the proopiomelanocortin (Pomc) gene and induce satiety. We disentangled the neuronal transcriptional code of Pomc by identifying the cis-acting regulatory elements and primary transcription factors controlling hypothalamic Pomc expression and determined their functional importance in the regulation of food intake and adiposity. Altogether, our studies reviewed here shed light on the power and limitations of the mammalian central satiety pathways and may contribute to the development of individual and collective strategies to reduce the debilitating effects of the self-induced obesity pandemia.

Composition and dosage of a multipartite enhancer cluster control developmental expression of Ihh (Indian hedgehog)

Copy number variations (CNVs) often include non-coding sequence and putative enhancers but how these rearrangements induce disease is poorly understood. Here we investigate CNVs involving the regulatory landscape of Indian hedgehog (IHH), causing multiple, highly localised phenotypes including craniosynostosis and synpolydactyly^1,2. We show through transgenic reporter and genome editing studies in mice that Ihh is regulated by a constellation of at least 9 enhancers with individual tissue specificities in the digit anlagen, growth plates, skull sutures and fingertips. Consecutive deletions show that they function in an additive manner resulting in growth defects of the skull and long bones. Duplications, in contrast, cause not only dose-dependent upregulation but also misexpression of Ihh, leading to abnormal phalanges, fusion of sutures and syndactyly. Thus, precise spatio-temporal control of developmental gene expression is achieved by complex multipartite enhancer ensembles. Alterations in the composition of such clusters can result in gene misexpression and disease.

Uncoupling evolutionary changes in DNA sequence, transcription factor occupancy and enhancer activity

Sequence variation within enhancers plays a major role in both evolution and disease, yet its functional impact on transcription factor (TF) occupancy and enhancer activity remains poorly understood. Here, we assayed the binding of five essential TFs over multiple stages of embryogenesis in two distant Drosophila species (with 1.4 substitutions per neutral site), identifying thousands of orthologous enhancers with conserved or diverged combinatorial occupancy. We used these binding signatures to dissect two properties of developmental enhancers: (1) potential TF cooperativity, using signatures of co-associations and co-divergence in TF occupancy. This revealed conserved combinatorial binding despite sequence divergence, suggesting protein-protein interactions sustain conserved collective occupancy. (2) Enhancer in-vivo activity, revealing orthologous enhancers with conserved activity despite divergence in TF occupancy. Taken together, we identify enhancers with diverged motifs yet conserved occupancy and others with diverged occupancy yet conserved activity, emphasising the need to functionally measure the effect of divergence on enhancer activity.

Large scale genomic reorganization of topological domains at the HoxD locus

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

A sequence level model of an intact locus predicts the location and function of nonadditive enhancers

Metazoan gene expression is controlled through the action of long stretches of noncoding DNA that contain enhancers-shorter sequences responsible for controlling a single aspect of a gene's expression pattern. Models built on thermodynamics have shown how enhancers interpret protein concentration in order to determine specific levels of gene expression, but the emergent regulatory logic of a complete regulatory locus shows qualitative and quantitative differences from isolated enhancers. Such differences may arise from steric competition limiting the quantity of DNA that can simultaneously influence the transcription machinery. We incorporated this competition into a mechanistic model of gene regulation, generated efficient algorithms for this computation, and applied it to the regulation of Drosophila even-skipped (eve). This model finds the location of enhancers and identifies which factors control the boundaries of eve expression. This model predicts a new enhancer that, when assayed in vivo, drives expression in a non-eve pattern. Incorporation of chromatin accessibility eliminates this inconsistency.