Shadow enhancers: frequently asked questions about distributed cis-regulatory information and enhancer redundancy

A Distalless-responsive enhancer of the Hox gene Sex combs reduced is required for segment- and sex-specific sensory organ development in Drosophila

Hox genes are involved in the patterning of animal body parts at multiple levels of regulatory hierarchies. Early expression of Hox genes in different domains along the embryonic anterior-posterior (A/P) axis in insects, vertebrates, and other animals establishes segmental or regional identity. However, Hox gene function is also required later in development for the patterning and morphogenesis of limbs and other organs. In Drosophila, spatiotemporal modulation of Sex combs reduced (Scr) expression within the first thoracic (T1) leg underlies the generation of segment- and sex-specific sense organ patterns. High Scr expression in defined domains of the T1 leg is required for the development of T1-specific transverse bristle rows in both sexes and sex combs in males, implying that the patterning of segment-specific sense organs involves incorporation of Scr into the leg development and sex determination gene networks. We sought to gain insight into this process by identifying the cis-and trans-regulatory factors that direct Scr expression during leg development. We have identified two cis-regulatory elements that control spatially modulated Scr expression within T1 legs. One of these enhancers directs sexually dimorphic expression and is required for the formation of T1-specific bristle patterns. We show that the Distalless and Engrailed homeodomain transcription factors act through sequences in this enhancer to establish elevated Scr expression in spatially defined domains. This enhancer functions to integrate Scr into the intrasegmental gene regulatory network, such that Scr serves as a link between leg patterning, sex determination, and sensory organ development.

Reporter Analyses Reveal Redundant Enhancers that Confer Robustness on Cis-Regulatory Mechanisms

Reporter analyses of Hox1 and Brachyury (Bra) genes have revealed examples of redundant enhancers that provide regulatory robustness. Retinoic acid (RA) activates through an RA-response element the transcription of Hox1 in the nerve cord of the ascidian Ciona intestinalis. We also found a weak RA-independent neural enhancer within the second intron of Hox1. The Hox1 gene in the larvacean Oikopleura dioica is also expressed in the nerve cord. The O. dioica genome, however, does not contain the RA receptor-encoding gene, and the expression of Hox1 has become independent of RA. We have found that the upstream sequence of the O. dioica Hox1 was able to activate reporter gene expression in the nerve cord of the C. intestinalis embryo, suggesting that an RA-independent regulatory system in the nerve cord might be common in larvaceans and ascidians. This RA-independent redundant regulatory system may have facilitated the Oikopleura ancestor losing RA signaling without an apparent impact on Hox1 expression domains. On the other hand, vertebrate Bra is expressed in the ventral mesoderm and notochord, whereas its ascidian ortholog is exclusively expressed in the notochord. Fibroblast growth factor (FGF) induces Bra in the ventral mesoderm in vertebrates, whereas it induces Bra in the notochord in ascidians. Disruption of the FGF signal does not completely silence Bra expression in ascidians, suggesting that FGF-dependent and independent enhancers might comprise a redundant regulatory system in ascidians. The existence of redundant enhancers, therefore, provides regulatory robustness that may facilitate the acquisition of new expression domains.

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.

Differential tissue specific, temporal and spatial expression patterns of the Aggrecan gene is modulated by independent enhancer elements

The transcriptional mechanism through which chondrocytes control the spatial and temporal composition of the cartilage tissue has remained largely elusive. The central aim of this study was to identify whether transcriptional enhancers played a role in the organisation of the chondrocytes in cartilaginous tissue. We focused on the Aggrecan gene (Acan) as it is essential for the normal structure and function of cartilage and it is expressed developmentally in different stages of chondrocyte maturation. Using transgenic reporter studies in mice we identified four elements, two of which showed individual chondrocyte developmental stage specificity. In particular, one enhancer (-80) distinguishes itself from the others by being predominantly active in adult cartilage. Furthermore, the -62 element uniquely drove reporter activity in early chondrocytes. The remaining chondrocyte specific enhancers, +28 and -30, showed no preference to chondrocyte type. The transcription factor SOX9 interacted with all the enhancers in vitro and mutation of SOX9 binding sites in one of the enhancers (-30) resulted in a loss of its chondrocyte specificity and ectopic enhancer reporter activity. Thus, the Acan enhancers orchestrate the precise spatiotemporal expression of this gene in cartilage types at different stages of development and adulthood.

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.

Global reorganisation of cis-regulatory units upon lineage commitment of human embryonic stem cells

Long-range cis-regulatory elements such as enhancers coordinate cell-specific transcriptional programmes by engaging in DNA looping interactions with target promoters. Deciphering the interplay between the promoter connectivity and activity of cis-regulatory elements during lineage commitment is crucial for understanding developmental transcriptional control. Here, we use Promoter Capture Hi-C to generate a high-resolution atlas of chromosomal interactions involving ~22,000 gene promoters in human pluripotent and lineage-committed cells, identifying putative target genes for known and predicted enhancer elements. We reveal extensive dynamics of cis-regulatory contacts upon lineage commitment, including the acquisition and loss of promoter interactions. This spatial rewiring occurs preferentially with predicted changes in the activity of cis-regulatory elements and is associated with changes in target gene expression. Our results provide a global and integrated view of promoter interactome dynamics during lineage commitment of human pluripotent cells.

Regulatory logic driving stable levels of defective proventriculus expression during terminal photoreceptor specification in flies

How differential levels of gene expression are controlled in post-mitotic neurons is poorly understood. In the Drosophila retina, expression of the transcription factor Defective Proventriculus (Dve) at distinct cell type-specific levels is required for terminal differentiation of color- and motion-detecting photoreceptors. Here, we find that the activities of two cis-regulatory enhancers are coordinated to drive dve expression in the fly eye. Three transcription factors act on these enhancers to determine cell-type specificity. Negative autoregulation by Dve maintains expression from each enhancer at distinct homeostatic levels. One enhancer acts as an inducible backup ('dark' shadow enhancer) that is normally repressed but becomes active in the absence of the other enhancer. Thus, two enhancers integrate combinatorial transcription factor input, feedback and redundancy to generate cell type-specific levels of dve expression and stable photoreceptor fate. This regulatory logic may represent a general paradigm for how precise levels of gene expression are established and maintained in post-mitotic neurons.

Ever-Changing Landscapes: Transcriptional Enhancers in Development and Evolution

A class of cis-regulatory elements, called enhancers, play a central role in orchestrating spatiotemporally precise gene-expression programs during development. Consequently, divergence in enhancer sequence and activity is thought to be an important mediator of inter- and intra-species phenotypic variation. Here, we give an overview of emerging principles of enhancer function, current models of enhancer architecture, genomic substrates from which enhancers emerge during evolution, and the influence of three-dimensional genome organization on long-range gene regulation. We discuss intricate relationships between distinct elements within complex regulatory landscapes and consider their potential impact on specificity and robustness of transcriptional regulation.

The Functionality and Evolution of Eukaryotic Transcriptional Enhancers

Enhancers regulate precise spatial and temporal patterns of gene expression in eukaryotes and, moreover, evolutionary changes in these modular cis-regulatory elements may represent the predominant genetic basis for phenotypic evolution. Here, we review approaches to identify and functionally analyze enhancers and their transcription factor binding sites, including assay for transposable-accessible chromatin-sequencing (ATAC-Seq) and clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9, respectively. We also explore enhancer functionality, including how transcription factor binding sites combine to regulate transcription, as well as research on shadow and super enhancers, and how enhancers can act over great distances and even in trans. Finally, we discuss recent theoretical and empirical data on how transcription factor binding sites and enhancers evolve. This includes how the function of enhancers is maintained despite the turnover of transcription factor binding sites as well as reviewing studies where mutations in enhancers have been shown to underlie morphological change.

The search for cis-regulatory driver mutations in cancer genomes

With the advent of high-throughput and relatively inexpensive whole-genome sequencing technology, the focus of cancer research has begun to shift toward analyses of somatic mutations in non-coding cis-regulatory elements of the cancer genome. Cis-regulatory elements play an important role in gene regulation, with mutations in these elements potentially resulting in changes to the expression of linked genes. The recent discoveries of recurrent TERT promoter mutations in melanoma, and recurrent mutations that create a super-enhancer regulating TAL1 expression in T-cell acute lymphoblastic leukaemia (T-ALL), have sparked significant interest in the search for other somatic cis-regulatory mutations driving cancer development. In this review, we look more closely at the TERT promoter and TAL1 enhancer alterations and use these examples to ask whether other cis-regulatory mutations may play a role in cancer susceptibility. In doing so, we make observations from the data emerging from recent research in this field, and describe the experimental and analytical approaches which could be adopted in the hope of better uncovering the true functional significance of somatic cis-regulatory mutations in cancer.

Gene regulation at a distance: From remote enhancers to 3D regulatory ensembles

Large-scale identification of elements associated with gene expression revealed that many of them are located extremely far from gene transcriptional start sites. We review here the growing evidence that show that distal cis-acting elements provide key instructions to genes, as genetic variation affecting them is growingly identified as an importance source of phenotypic diversity and disease. We discuss the different mechanisms that allow these elements to exert their regulatory functions, in a robust and specific manner, despite the large genomic distances separating them from their target genes. We particularly focus on the role of the structural organization of the genome in guiding such regulatory interactions.

Stepwise Progression of Embryonic Patterning

It is long established that the graded distribution of Dorsal transcription factor influences spatial domains of gene expression along the dorsoventral (DV) axis of Drosophila melanogaster embryos. However, the more recent realization that Dorsal levels also change with time raises the question of whether these dynamics are instructive. An overview of DV axis patterning is provided, focusing on new insights identified through quantitative analysis of temporal changes in Dorsal target gene expression from one nuclear cycle to the next ('steps'). Possible roles for the stepwise progression of this patterning program are discussed including (i) tight temporal regulation of signaling pathway activation, (ii) control of gene expression cohorts, and (iii) ensuring the irreversibility of the patterning and cell fate specification process.

An ancient yet flexible cis-regulatory architecture allows localized Hedgehog tuning by patched/Ptch1

The Hedgehog signaling pathway is part of the ancient developmental-evolutionary animal toolkit. Frequently co-opted to pattern new structures, the pathway is conserved among eumetazoans yet flexible and pleiotropic in its effects. The Hedgehog receptor, Patched, is transcriptionally activated by Hedgehog, providing essential negative feedback in all tissues. Our locus-wide dissections of the cis-regulatory landscapes of fly patched and mouse Ptch1 reveal abundant, diverse enhancers with stage- and tissue-specific expression patterns. The seemingly simple, constitutive Hedgehog response of patched/Ptch1 is driven by a complex regulatory architecture, with batteries of context-specific enhancers engaged in promoter-specific interactions to tune signaling individually in each tissue, without disturbing patterning elsewhere. This structure—one of the oldest cis-regulatory features discovered in animal genomes—explains how patched/Ptch1 can drive dramatic adaptations in animal morphology while maintaining its essential core function. It may also suggest a general model for the evolutionary flexibility of conserved regulators and pathways.

DOI:http://dx.doi.org/10.7554/eLife.13550.001

Sequential construction of a model for modular gene expression control, applied to spatial patterning of theDrosophilagenehunchback

Gene network simulations are increasingly used to quantify mutual gene regulation in biological tissues. These are generally based on linear interactions between single-entity regulatory and target genes. Biological genes, by contrast, commonly have multiple, partially independent, cis-regulatory modules (CRMs) for regulator binding, and can produce variant transcription and translation products. We present a modeling framework to address some of the gene regulatory dynamics implied by this biological complexity. Spatial patterning of the hunchback (hb) gene in Drosophila development involves control by three CRMs producing two distinct mRNA transcripts. We use this example to develop a differential equations model for transcription which takes into account the cis-regulatory architecture of the gene. Potential regulatory interactions are screened by a genetic algorithms (GAs) approach and compared to biological expression data.

An autonomous CEBPA enhancer specific for myeloid-lineage priming and neutrophilic differentiation

Neutrophilic differentiation is dependent on CCAAT enhancer-binding protein α (C/EBPα), a transcription factor expressed in multiple organs including the bone marrow. Using functional genomic technologies in combination with clustered regularly-interspaced short palindromic repeat (CRISPR)/CRISPR-associated protein 9 genome editing and in vivo mouse modeling, we show that CEBPA is located in a 170-kb topological-associated domain that contains 14 potential enhancers. Of these, 1 enhancer located +42 kb from CEBPA is active and engages with the CEBPA promoter in myeloid cells only. Germ line deletion of the homologous enhancer in mice in vivo reduces Cebpa levels exclusively in hematopoietic stem cells (HSCs) and myeloid-primed progenitor cells leading to severe defects in the granulocytic lineage, without affecting any other Cebpa-expressing organ studied. The enhancer-deleted progenitor cells lose their myeloid transcription program and are blocked in differentiation. Deletion of the enhancer also causes loss of HSC maintenance. We conclude that a single +42-kb enhancer is essential for CEBPA expression in myeloid cells only.

Neural Progenitors Adopt Specific Identities by Directly Repressing All Alternative Progenitor Transcriptional Programs

In the vertebrate neural tube, a morphogen-induced transcriptional network produces multiple molecularly distinct progenitor domains, each generating different neuronal subtypes. Using an in vitro differentiation system, we defined gene expression signatures of distinct progenitor populations and identified direct gene-regulatory inputs corresponding to locations of specific transcription factor binding. Combined with targeted perturbations of the network, this revealed a mechanism in which a progenitor identity is installed by active repression of the entire transcriptional programs of other neural progenitor fates. In the ventral neural tube, sonic hedgehog (Shh) signaling, together with broadly expressed transcriptional activators, concurrently activates the gene expression programs of several domains. The specific outcome is selected by repressive input provided by Shh-induced transcription factors that act as the key nodes in the network, enabling progenitors to adopt a single definitive identity from several initially permitted options. Together, the data suggest design principles relevant to many developing tissues.

•

Specific vertebrate neural progenitor populations generated in vitro

•

Gene expression dynamics, transcription factor binding assessed in neural progenitors

•

Progenitor fate selected by repressors blocking entire programs of other identities

•

Repressors counteract non-selective morphogen and pan-neural activatory inputs

Enhancers as non-coding RNA transcription units: recent insights and future perspectives

Networks of regulatory enhancers dictate distinct cell identities and cellular responses to diverse signals by instructing precise spatiotemporal patterns of gene expression. However, 35 years after their discovery, enhancer functions and mechanisms remain incompletely understood. Intriguingly, recent evidence suggests that many, if not all, functional enhancers are themselves transcription units, generating non-coding enhancer RNAs. This observation provides a fundamental insight into the inter-regulation between enhancers and promoters, which can both act as transcription units; it also raises crucial questions regarding the potential biological roles of the enhancer transcription process and non-coding enhancer RNAs. Here, we review research progress in this field and discuss several important, unresolved questions regarding the roles and mechanisms of enhancers in gene regulation.

Lineage-affiliated transcription factors bind the Gata3 Tce1 enhancer to mediate lineage-specific programs

The transcription factor GATA3 is essential for the genesis and maturation of the T cell lineage, and GATA3 dysregulation has pathological consequences. Previous studies have shown that GATA3 function in T cell development is regulated by multiple signaling pathways and that the Notch nuclear effector, RBP-J, binds specifically to the Gata3 promoter. We previously identified a T cell-specific Gata3 enhancer (Tce1) lying 280 kb downstream from the structural gene and demonstrated in transgenic mice that Tce1 promoted T lymphocyte-specific transcription of reporter genes throughout T cell development; however, it was not clear if Tce1 is required for Gata3 transcription in vivo. Here, we determined that the canonical Gata3 promoter is insufficient for Gata3 transcriptional activation in T cells in vivo, precluding the possibility that promoter binding by a host of previously implicated transcription factors alone is responsible for Gata3 expression in T cells. Instead, we demonstrated that multiple lineage-affiliated transcription factors bind to Tce1 and that this enhancer confers T lymphocyte-specific Gata3 activation in vivo, as targeted deletion of Tce1 in a mouse model abrogated critical functions of this T cell-regulatory element. Together, our data show that Tce1 is both necessary and sufficient for critical aspects of Gata3 T cell-specific transcriptional activity.

Morphogen rules: design principles of gradient-mediated embryo patterning

The Drosophila blastoderm and the vertebrate neural tube are archetypal examples of morphogen-patterned tissues that create precise spatial patterns of different cell types. In both tissues, pattern formation is dependent on molecular gradients that emanate from opposite poles. Despite distinct evolutionary origins and differences in time scales, cell biology and molecular players, both tissues exhibit striking similarities in the regulatory systems that establish gene expression patterns that foreshadow the arrangement of cell types. First, signaling gradients establish initial conditions that polarize the tissue, but there is no strict correspondence between specific morphogen thresholds and boundary positions. Second, gradients initiate transcriptional networks that integrate broadly distributed activators and localized repressors to generate patterns of gene expression. Third, the correct positioning of boundaries depends on the temporal and spatial dynamics of the transcriptional networks. These similarities reveal design principles that are likely to be broadly applicable to morphogen-patterned tissues.

Gene regulation during development in the light of topologically associating domains

During embryonic development, complex transcriptional programs govern the precision of gene expression. Many key developmental genes are regulated via cis-regulatory elements that are located far away in the linear genome. How sequences located hundreds of kilobases away from a promoter can influence its activity has been the subject of numerous speculations, which all underline the importance of the 3D-organization of the genome. The recent advent of chromosome conformation capture techniques has put into focus the subdivision of the genome into topologically associating domains (TADs). TADs may influence regulatory activities on multiple levels. The relative invariance of TAD limits across cell types suggests that they may form fixed structural domains that could facilitate and/or confine long-range regulatory interactions. However, most recent studies suggest that interactions within TADs are more variable and dynamic than initially described. Hence, different models are emerging regarding how TADs shape the complex 3D conformations, and thereafter influence the networks of cis-interactions that govern gene expression during development. For further resources related to this article, please visit the WIREs website.

Krüppel Expression Levels Are Maintained through Compensatory Evolution of Shadow Enhancers

Many developmental genes are controlled by shadow enhancers—pairs of enhancers that drive overlapping expression patterns. We hypothesized that compensatory evolution can maintain the total expression of a gene, while individual shadow enhancers diverge between species. To test this hypothesis, we analyzed expression driven by orthologous pairs of shadow enhancers from Drosophila melanogaster, Drosophila yakuba, and Drosophila pseudoobscura that control expression of Krüppel, a transcription factor that patterns the anterior-posterior axis of blastoderm embryos. We found that the expression driven by the pair of enhancers is conserved between these three species, but expression levels driven by the individual enhancers are not. Using sequence analysis and experimental perturbation, we show that each shadow enhancer is regulated by different transcription factors. These results support the hypothesis that compensatory evolution can occur between shadow enhancers, which has implications for mechanistic and evolutionary studies of gene regulation.

Functions of noncoding sequences in mammalian genomes

Most of the mammalian genome consists of nucleotide sequences not coding for proteins. Exons of genes make up only 3% of the human genome, while the significance of most other sequences remains unknown. Recent genome studies with high-throughput methods demonstrate that the so-called noncoding part of the genome may perform important functions. This hypothesis is supported by three groups of experimental data: 1) approximately 10% of the sequences, most of which are located in noncoding parts of the genome, is evolutionarily conserved and thus can be of functional importance; 2) up to 99% of the mammalian genome is being transcribed forming short and long noncoding RNAs in addition to common mRNA; and 3) mutations in noncoding parts of the genome can be accompanied by progression of pathological states of the organism. In the light of these data, in the review we consider the functional role of numerous known sequences of noncoding parts of the genome including introns, DNA methylation regions, enhancers and locus control regions, insulators, S/MAR sequences, pseudogenes, and genes of noncoding RNAs, as well as transposons and simple repeats of centromeric and telomeric regions of chromosomes. The assumption is made that the intergenic noncoding sequences without definite/clear functions can be involved in spatial organization of genetic loci in interphase nuclei.

Genetic Changes to a Transcriptional Silencer Element Confers Phenotypic Diversity within and between Drosophila Species

The modification of transcriptional regulation has become increasingly appreciated as a major contributor to morphological evolution. However, the role of negative-acting control elements (e.g. silencers) in generating morphological diversity has been generally overlooked relative to positive-acting “enhancer” elements. The highly variable body coloration patterns among Drosophilid insects represents a powerful model system in which the molecular alterations that underlie phenotypic diversity can be defined. In a survey of pigment phenotypes among geographically disparate Japanese populations of Drosophila auraria, we discovered a remarkable degree of variation in male-specific abdominal coloration. In testing the expression patterns of the major pigment-producing enzymes, we found that phenotypes uniquely correlated with differences in the expression of ebony, a gene required for yellow-colored cuticle. Assays of ebony’s transcriptional control region indicated that a lightly pigmented strain harbored cis-regulatory mutations that caused correlated changes in its expression. Through a series of chimeric reporter constructs between light and dark strain alleles, we localized function-altering mutations to a conserved silencer that mediates a male-specific pattern of ebony repression. This suggests that the light allele was derived through the loss of this silencer’s activity. Furthermore, examination of the ebony gene of D. serrata, a close relative of D. auraria which secondarily lost male-specific pigmentation revealed the parallel loss of this silencer element. These results demonstrate how loss-of-function mutations in a silencer element resulted in increased gene expression. We propose that the mutational inactivation of silencer elements may represent a favored path to evolve gene expression, impacting morphological traits.

One of the greatest challenges in understanding the relationship between genotype and phenotype is to discern how changes in DNA affect the normal functioning of genes. Mutations may generate a new function for a gene, yet it is frequently observed that they inactivate some aspect of a gene’s normal capacity. Investigations focused on understanding the developmental basis for the evolution of anatomical structures has found a prevalent role for mutations that alter developmental gene regulation. In animals, genes are transcriptionally activated in specific tissues during development by regulatory sequences distributed across their expansive non-protein coding regions. Regulatory elements known as silencers act to prevent genes from being expressed in certain tissues, providing a mechanism for precise control. Here, we show how a silencer that prevents expression of a pigment-producing enzyme in certain Drosophila species has repeatedly been subject to inactivating mutations that increased this gene’s expression. This example illustrates how such negative-acting regulatory sequences can represent a convenient target for increasing gene expression through the loss of a genetic element.

What does it take to evolve an enhancer? A simulation-based study of factors influencing the emergence of combinatorial regulation

There is widespread interest today in understanding enhancers, which are regulatory elements typically harboring several transcription factor binding sites and mediating the combinatorial effect of transcription factors on gene expression. The evolution of enhancers poses interesting unanswered questions, for example, the evolutionary time taken for a typical enhancer to emerge or the factors shaping its evolution. Existing approaches to cis-regulatory evolution have often ignored the combinatorial nature and varied biochemical mechanisms of gene regulation encoded in enhancers. We report on our investigation of enhancer evolution through the use of PEBCRES, a framework for evolutionary simulation of enhancers that employs a mechanistic and well-supported sequence-to-expression model to assign fitness to the evolving enhancer genotype. We estimated the time necessary to evolve, from genomic background, enhancers capable of driving complex gene expression patterns similar to those involved in early development in Drosophila. We found the time-to-evolve to range between 0.5 and 10 Myr, and to vary greatly with the target expression pattern, complexity of the real enhancer known to encode that pattern, and the strength of input from specific transcription factors. To our knowledge, this is the first estimate of waiting times for realistic enhancers to evolve. The in silico evolved enhancers had, with a few interesting exceptions, site compositions similar to those seen in real enhancers for the same patterns. Our simulations also revealed that certain features of an enhancer might evolve not due to their biological function but as aids to the evolutionary process itself.

Investigating the transcriptional control of cardiovascular development

Transcriptional regulation of thousands of genes instructs complex morphogenetic and molecular events for heart development. Cardiac transcription factors choreograph gene expression at each stage of differentiation by interacting with cofactors, including chromatin-modifying enzymes, and by binding to a constellation of regulatory DNA elements. Here, we present salient examples relevant to cardiovascular development and heart disease, and review techniques that can sharpen our understanding of cardiovascular biology. We discuss the interplay between cardiac transcription factors, cis-regulatory elements, and chromatin as dynamic regulatory networks, to orchestrate sequential deployment of the cardiac gene expression program.

The selection and function of cell type-specific enhancers

The human body contains several hundred cell types, all of which share the same genome. In metazoans, much of the regulatory code that drives cell type-specific gene expression is located in distal elements called enhancers. Although mammalian genomes contain millions of potential enhancers, only a small subset of them is active in a given cell type. Cell type-specific enhancer selection involves the binding of lineage-determining transcription factors that prime enhancers. Signal-dependent transcription factors bind to primed enhancers, which enables these broadly expressed factors to regulate gene expression in a cell type-specific manner. The expression of genes that specify cell type identity and function is associated with densely spaced clusters of active enhancers known as super-enhancers. The functions of enhancers and super-enhancers are influenced by, and affect, higher-order genomic organization.

Partially redundant enhancers cooperatively maintain Mammalian pomc expression above a critical functional threshold

Cell-specific expression of many genes is conveyed by multiple enhancers, with each individual enhancer controlling a particular expression domain. In contrast, multiple enhancers drive similar expression patterns of some genes involved in embryonic development, suggesting regulatory redundancy. Work in Drosophila has indicated that functionally overlapping enhancers canalize development by buffering gene expression against environmental and genetic disturbances. However, little is known about regulatory redundancy in vertebrates and in genes mainly expressed during adulthood. Here we study nPE1 and nPE2, two phylogenetically conserved mammalian enhancers that drive expression of the proopiomelanocortin gene (Pomc) to the same set of hypothalamic neurons. The simultaneous deletion of both enhancers abolished Pomc expression at all ages and induced a profound metabolic dysfunction including early-onset extreme obesity. Targeted inactivation of either nPE1 or nPE2 led to very low levels of Pomc expression during early embryonic development indicating that both enhancers function synergistically. In adult mice, however, Pomc expression is controlled additively by both enhancers, with nPE1 being responsible for ∼80% and nPE2 for ∼20% of Pomc transcription. Consequently, nPE1 knockout mice exhibit mild obesity whereas nPE2-deficient mice maintain a normal body weight. These results suggest that nPE2-driven Pomc expression is compensated by nPE1 at later stages of development, essentially rescuing the earlier phenotype of nPE2 deficiency. Together, these results reveal that cooperative interactions between the enhancers confer robustness of Pomc expression against gene regulatory disturbances and preclude deleterious metabolic phenotypes caused by Pomc deficiency in adulthood. Thus, our study demonstrates that enhancer redundancy can be used by genes that control adult physiology in mammals and underlines the potential significance of regulatory sequence mutations in common diseases.

The stability of animal form and function in the face of genetic and environmental variation relies on consistent gene expression. Multiple enhancers, each specifying a unique regulatory domain, control the precise spatiotemporal expression of many genes. However, in some genes apparently redundant enhancers regulate expression in overlapping cell-specific patterns. Although this arrangement has been shown to be important for developmental robustness in invertebrates, the role of apparently redundant enhancers in vertebrate species and in genes functioning in adulthood is poorly understood. Here, we show that expression of the mammalian Pomc gene is controlled in a tissue-specific manner by two such apparently redundant enhancers. We used targeted deletion of the individual enhancers to delineate their respective contributions to Pomc expression in the brain. Since Pomc expression from its intact locus exceeds the sum of the individual enhancer contributions to Pomc mRNA levels in embryonic mice, we infer a synergistic action between the enhancers during development. In contrast, the interaction between the enhancers is additive in adult mice. Deletion of both enhancers simultaneously almost completely abolished Pomc expression and the mutant mice displayed extreme obesity and metabolic dysfunction, while deletion of the individual enhancers had a modest or no phenotypic effect. Together, our results demonstrate that the two enhancers cooperatively maintain Pomc expression above a critical functional threshold.

Shadow enhancers enable Hunchback bifunctionality in the Drosophila embryo

Hunchback (Hb) is a bifunctional transcription factor that activates and represses distinct enhancers. Here, we investigate the hypothesis that Hb can activate and repress the same enhancer. Computational models predicted that Hb bifunctionally regulates the even-skipped (eve) stripe 3+7 enhancer (eve3+7) in Drosophila blastoderm embryos. We measured and modeled eve expression at cellular resolution under multiple genetic perturbations and found that the eve3+7 enhancer could not explain endogenous eve stripe 7 behavior. Instead, we found that eve stripe 7 is controlled by two enhancers: the canonical eve3+7 and a sequence encompassing the minimal eve stripe 2 enhancer (eve2+7). Hb bifunctionally regulates eve stripe 7, but it executes these two activities on different pieces of regulatory DNA--it activates the eve2+7 enhancer and represses the eve3+7 enhancer. These two "shadow enhancers" use different regulatory logic to create the same pattern.

Neural precursor-specific expression of multiple Drosophila genes is driven by dual enhancer modules with overlapping function

Transcriptional cis-regulatory modules (CRMs), or enhancers, are responsible for directing gene expression in specific territories and cell types during development. In some instances, the same gene may be served by two or more enhancers with similar specificities. Here we show that the utilization of dual, or "shadow", enhancers is a common feature of genes that are active specifically in neural precursor (NP) cells in Drosophila. By genome-wide computational discovery of statistically significant clusters of binding motifs for both proneural activator (P) proteins and basic helix-loop-helix (bHLH) repressor (R) factors (a "P+R" regulatory code), we have identified NP-specific enhancer modules associated with multiple genes expressed in this cell type. These CRMs are distinct from those previously identified for the corresponding gene, establishing the existence of a dual-enhancer arrangement in which both modules reside close to the gene they serve. Using wild-type and mutant reporter gene constructs in vivo, we show that P sites in these modules mediate activation by proneural factors in "proneural cluster" territories, whereas R sites mediate repression by bHLH repressors, which serves to restrict expression specifically to NP cells. To our knowledge, our results identify the first direct targets of these bHLH repressors. Finally, using genomic rescue constructs for neuralized (neur), we demonstrate that each of the gene's two NP-specific enhancers is sufficient to rescue neur function in the lateral inhibition process by which adult sensory organ precursor (SOP) cells are specified, but that deletion of both enhancers results in failure of this event.

Co-regulation of invected and engrailed by a complex array of regulatory sequences in Drosophila

invected (inv) and engrailed (en) form a gene complex that extends about 115 kb. These two genes encode highly related homeodomain proteins that are co-regulated in a complex manner throughout development. Our dissection of inv/en regulatory DNA shows that most enhancers are spread throughout a 62 kb region. We used two types of constructs to analyze the function of this DNA: P-element based reporter constructs with small pieces of DNA fused to the en promoter driving lacZ expression and large constructs with HA-tagged en and inv inserted in the genome with the phiC31 system. In addition, we generated deletions of inv and en DNA in situ and assayed their effects on inv/en expression. Our results support and extend our knowledge of inv/en regulation. First, inv and en share regulatory DNA, most of which is flanking the en transcription unit. In support of this, a 79-kb HA-en transgene can rescue inv en double mutants to viable, fertile adults. In contrast, an 84-kb HA-inv transgene lacks most of the enhancers for inv/en expression. Second, there are multiple enhancers for inv/en stripes in embryos; some of these may be redundant but others play discrete roles at different stages of embryonic development. Finally, no small reporter construct gave expression in the posterior compartment of imaginal discs, a hallmark of inv/en expression. Robust expression of HA-en in the posterior compartment of imaginal discs is evident from the 79-kb HA-en transgene, while a 45-kb HA-en transgene gives weaker, variable imaginal disc expression. We suggest that the activity of the imaginal disc enhancer(s) is dependent on the chromatin structure of the inv/en domain.

Analytic approaches to stochastic gene expression in multicellular systems

Deterministic thermodynamic models of the complex systems, which control gene expression in metazoa, are helping researchers identify fundamental themes in the regulation of transcription. However, quantitative single cell studies are increasingly identifying regulatory mechanisms that control variability in expression. Such behaviors cannot be captured by deterministic models and are poorly suited to contemporary stochastic approaches that rely on continuum approximations, such as Langevin methods. Fortunately, theoretical advances in the modeling of transcription have assembled some general results that can be readily applied to systems being explored only through a deterministic approach. Here, I review some of the recent experimental evidence for the importance of genetically regulating stochastic effects during embryonic development and discuss key results from Markov theory that can be used to model this regulation. I then discuss several pairs of regulatory mechanisms recently investigated through a Markov approach. In each case, a deterministic treatment predicts no difference between the mechanisms, but the statistical treatment reveals the potential for substantially different distributions of transcriptional activity. In this light, features of gene regulation that seemed needlessly complex evolutionary baggage may be appreciated for their key contributions to reliability and precision of gene expression.

The embryo as a laboratory: quantifying transcription in Drosophila

Transcriptional regulation of gene expression is fundamental to most cellular processes, including determination of cellular fates. Quantitative studies of transcription in cultured cells have led to significant advances in identifying mechanisms underlying transcriptional control. Recent progress allowed implementation of these same quantitative methods in multicellular organisms to ask how transcriptional regulation unfolds both in vivo and at the single molecule level in the context of embryonic development. Here we review some of these advances in early Drosophila development, which bring the embryo on par with its single celled counterparts. In particular, we discuss progress in methods to measure mRNA and protein distributions in fixed and living embryos, and we highlight some initial applications that lead to fundamental new insights about molecular transcription processes. We end with an outlook on how to further exploit the unique advantages that come with investigating transcriptional control in the multicellular context of development.

Microsatellite repeat instability fuels evolution of embryonic enhancers in Hawaiian Drosophila

For ∼30 million years, the eggs of Hawaiian Drosophila were laid in ever-changing environments caused by high rates of island formation. The associated diversification of the size and developmental rate of the syncytial fly embryo would have altered morphogenic gradients, thus necessitating frequent evolutionary compensation of transcriptional responses. We investigate the consequences these radiations had on transcriptional enhancers patterning the embryo to see whether their pattern of molecular evolution is different from non-Hawaiian species. We identify and functionally assay in transgenic D. melanogaster the Neurogenic Ectoderm Enhancers from two different Hawaiian Drosophila groups: (i) the picture wing group, and (ii) the modified mouthparts group. We find that the binding sites in this set of well-characterized enhancers are footprinted by diverse microsatellite repeat (MSR) sequences. We further show that Hawaiian embryonic enhancers in general are enriched in MSR relative to both Hawaiian non-embryonic enhancers and non-Hawaiian embryonic enhancers. We propose embryonic enhancers are sensitive to Activator spacing because they often serve as assembly scaffolds for the aggregation of transcription factor activator complexes. Furthermore, as most indels are produced by microsatellite repeat slippage, enhancers from Hawaiian Drosophila lineages, which experience dynamic evolutionary pressures, would become grossly enriched in MSR content.

In pursuit of design principles of regulatory sequences

Instructions for when, where and to what level each gene should be expressed are encoded within regulatory sequences. The importance of motifs recognized by DNA-binding regulators has long been known, but their extensive characterization afforded by recent technologies only partly accounts for how regulatory instructions are encoded in the genome. Here, we review recent advances in our understanding of regulatory sequences that influence transcription and go beyond the description of motifs. We discuss how understanding different aspects of the sequence-encoded regulation can help to unravel the genotype-phenotype relationship, which would lead to a more accurate and mechanistic interpretation of personal genome sequences.

The architecture of gene expression: integrating dispersed cis-regulatory modules into coherent regulatory domains

Specificity and precision of expression are essential for the genes that regulate developmental processes. The specialized cis-acting modules, such as enhancers, that define gene expression patterns can be distributed across large regions, raising questions about the nature of the mechanisms that underline their action. Recent data has exposed the structural 3D context in which these long-range enhancers are operating. Here, we present how these studies shed new light on principles driving long-distance regulatory relationships. We discuss the molecular mechanisms that enable and accompany the action of long-range acting elements and the integration of multiple distributed regulatory inputs into the coherent and specific regulatory programs that are key to embryonic development.

Spurious transcription factor binding: non-functional or genetically redundant?

Transcription factor binding sites (TFBSs) on the DNA are generally accepted as the key nodes of gene control. However, the multitudes of TFBSs identified in genome-wide studies, some of them seemingly unconstrained in evolution, have prompted the view that in many cases TF binding may serve no biological function. Yet, insights from transcriptional biochemistry, population genetics and functional genomics suggest that rather than segregating into 'functional' or 'non-functional', TFBS inputs to their target genes may be generally cumulative, with varying degrees of potency and redundancy. As TFBS redundancy can be diminished by mutations and environmental stress, some of the apparently 'spurious' sites may turn out to be important for maintaining adequate transcriptional regulation under these conditions. This has significant implications for interpreting the phenotypic effects of TFBS mutations, particularly in the context of genome-wide association studies for complex traits.

The chromatin landscape and transcription factors in T cell programming

T cell development from multipotent progenitors to specialized effector subsets of mature T cells is guided by the iterative action of transcription factors. At each stage, transcription factors interact not only with an existing landscape of histone modifications and nucleosome packing, but also with other bound factors, while they modify the landscape for later-arriving factors in ways that fundamentally affect the control of gene expression. This review covers insights from genome-wide analyses of transcription factor binding and resulting chromatin conformation changes that reveal roles of cytokine signaling in effector T cell programming, the ways in which one factor can completely transform the impacts of previously bound factors, and the ways in which the baseline chromatin landscape is established during early T cell lineage commitment.

In silico evolution of the hunchback gene indicates redundancy in cis-regulatory organization and spatial gene expression

Biological development depends on the coordinated expression of genes in time and space. Developmental genes have extensive cis-regulatory regions which control their expression. These regions are organized in a modular manner, with different modules controlling expression at different times and locations. Both how modularity evolved and what function it serves are open questions. We present a computational model for the cis-regulation of the hunchback (hb) gene in the fruit fly (Drosophila). We simulate evolution (using an evolutionary computation approach from computer science) to find the optimal cis-regulatory arrangements for fitting experimental hb expression patterns. We find that the cis-regulatory region tends to readily evolve modularity. These cis-regulatory modules (CRMs) do not tend to control single spatial domains, but show a multi-CRM/multi-domain correspondence. We find that the CRM-domain correspondence seen in Drosophila evolves with a high probability in our model, supporting the biological relevance of the approach. The partial redundancy resulting from multi-CRM control may confer some biological robustness against corruption of regulatory sequences. The technique developed on hb could readily be applied to other multi-CRM developmental genes.

A Gli silencer is required for robust repression of gremlin in the vertebrate limb bud

The transcriptional response to the Hedgehog (Hh) pathway is mediated by Gli proteins, which function as context-dependent transcriptional activators or repressors. However, the mechanism by which Gli proteins regulate their target genes is poorly understood. Here, we have performed the first genetic characterization of a Gli-dependent cis-regulatory module (CRM), focusing on its regulation of Grem1 in the mouse limb bud. The CRM, termed GRE1 (Gli responsive element 1), can act as both an enhancer and a silencer. The enhancer activity requires sustained Hh signaling. As a Gli-dependent silencer, GRE1 prevents ectopic transcription of Grem1 driven through additional CRMs. In doing so, GRE1 works with additional GREs to robustly regulate Grem1. We suggest that multiple Gli CRMs may be a general mechanism for mediating a robust transcriptional response to the Hh pathway.

A promoter-level mammalian expression atlas

Regulated transcription controls the diversity, developmental pathways and spatial organization of the hundreds of cell types that make up a mammal. Using single-molecule cDNA sequencing, we mapped transcription start sites (TSSs) and their usage in human and mouse primary cells, cell lines and tissues to produce a comprehensive overview of mammalian gene expression across the human body. We find that few genes are truly 'housekeeping', whereas many mammalian promoters are composite entities composed of several closely separated TSSs, with independent cell-type-specific expression profiles. TSSs specific to different cell types evolve at different rates, whereas promoters of broadly expressed genes are the most conserved. Promoter-based expression analysis reveals key transcription factors defining cell states and links them to binding-site motifs. The functions of identified novel transcripts can be predicted by coexpression and sample ontology enrichment analyses. The functional annotation of the mammalian genome 5 (FANTOM5) project provides comprehensive expression profiles and functional annotation of mammalian cell-type-specific transcriptomes with wide applications in biomedical research.

An atlas of active enhancers across human cell types and tissues

Enhancers control the correct temporal and cell-type-specific activation of gene expression in multicellular eukaryotes. Knowing their properties, regulatory activity and targets is crucial to understand the regulation of differentiation and homeostasis. Here we use the FANTOM5 panel of samples, covering the majority of human tissues and cell types, to produce an atlas of active, in vivo-transcribed enhancers. We show that enhancers share properties with CpG-poor messenger RNA promoters but produce bidirectional, exosome-sensitive, relatively short unspliced RNAs, the generation of which is strongly related to enhancer activity. The atlas is used to compare regulatory programs between different cells at unprecedented depth, to identify disease-associated regulatory single nucleotide polymorphisms, and to classify cell-type-specific and ubiquitous enhancers. We further explore the utility of enhancer redundancy, which explains gene expression strength rather than expression patterns. The online FANTOM5 enhancer atlas represents a unique resource for studies on cell-type-specific enhancers and gene regulation.

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

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

Evolution of transcriptional enhancers and animal diversity

Deciphering the genetic bases that drive animal diversity is one of the major challenges of modern biology. Although four decades ago it was proposed that animal evolution was mainly driven by changes in cis-regulatory DNA elements controlling gene expression rather than in protein-coding sequences, only now are powerful bioinformatics and experimental approaches available to accelerate studies into how the evolution of transcriptional enhancers contributes to novel forms and functions. In the introduction to this Theme Issue, we start by defining the general properties of transcriptional enhancers, such as modularity and the coexistence of tight sequence conservation with transcription factor-binding site shuffling as different mechanisms that maintain the enhancer grammar over evolutionary time. We discuss past and current methods used to identify cell-type-specific enhancers and provide examples of how enhancers originate de novo, change and are lost in particular lineages. We then focus in the central part of this Theme Issue on analysing examples of how the molecular evolution of enhancers may change form and function. Throughout this introduction, we present the main findings of the articles, reviews and perspectives contributed to this Theme Issue that together illustrate some of the great advances and current frontiers in the field.

Simulations of enhancer evolution provide mechanistic insights into gene regulation

There is growing interest in models of regulatory sequence evolution. However, existing models specifically designed for regulatory sequences consider the independent evolution of individual transcription factor (TF)-binding sites, ignoring that the function and evolution of a binding site depends on its context, typically the cis-regulatory module (CRM) in which the site is located. Moreover, existing models do not account for the gene-specific roles of TF-binding sites, primarily because their roles often are not well understood. We introduce two models of regulatory sequence evolution that address some of the shortcomings of existing models and implement simulation frameworks based on them. One model simulates the evolution of an individual binding site in the context of a CRM, while the other evolves an entire CRM. Both models use a state-of-the art sequence-to-expression model to predict the effects of mutations on the regulatory output of the CRM and determine the strength of selection. We use the new framework to simulate the evolution of TF-binding sites in 37 well-studied CRMs belonging to the anterior-posterior patterning system in Drosophila embryos. We show that these simulations provide accurate fits to evolutionary data from 12 Drosophila genomes, which includes statistics of binding site conservation on relatively short evolutionary scales and site loss across larger divergence times. The new framework allows us, for the first time, to test hypotheses regarding the underlying cis-regulatory code by directly comparing the evolutionary implications of the hypothesis with the observed evolutionary dynamics of binding sites. Using this capability, we find that explicitly modeling self-cooperative DNA binding by the TF Caudal (CAD) provides significantly better fits than an otherwise identical evolutionary simulation that lacks this mechanistic aspect. This hypothesis is further supported by a statistical analysis of the distribution of intersite spacing between adjacent CAD sites. Experimental tests confirm direct homodimeric interaction between CAD molecules as well as self-cooperative DNA binding by CAD. We note that computational modeling of the D. melanogaster CRMs alone did not yield significant evidence to support CAD self-cooperativity. We thus demonstrate how specific mechanistic details encoded in CRMs can be revealed by modeling their evolution and fitting such models to multispecies data.