Dynamic evolution of precise regulatory encodings creates the clustered site signature of enhancers

Homeodomain complex formation and biomolecular condensates in Hox gene regulation

Hox genes are a family of homeodomain transcription factors that regulate specialized morphological structures along the anterior-posterior axis of metazoans. Over the past few decades, researchers have focused on defining how Hox factors with similar in vitro DNA binding activities achieve sufficient target specificity to regulate distinct cell fates in vivo. In this review, we highlight how protein interactions with other transcription factors, many of which are also homeodomain proteins, result in the formation of transcription factor complexes with enhanced DNA binding specificity. These findings suggest that Hox-regulated enhancers utilize distinct combinations of homeodomain binding sites, many of which are low-affinity, to recruit specific Hox complexes. However, low-affinity sites can only yield reproducible responses with high transcription factor concentrations. To overcome this limitation, recent studies revealed how transcription factors, including Hox factors, use intrinsically disordered domains (IDRs) to form biomolecular condensates that increase protein concentrations. Moreover, Hox factors with altered IDRs have been associated with altered transcriptional activity and human disease states, demonstrating the importance of IDRs in mediating essential Hox output. Collectively, these studies highlight how Hox factors use their DNA binding domains, protein-protein interaction domains, and IDRs to form specific transcription factor complexes that yield accurate gene expression.

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.

Nuclear microenvironments modulate transcription from low-affinity enhancers

Transcription factors bind low-affinity DNA sequences for only short durations. It is not clear how brief, low-affinity interactions can drive efficient transcription. Here, we report that the transcription factor Ultrabithorax (Ubx) utilizes low-affinity binding sites in the Drosophila melanogaster shavenbaby (svb) locus and related enhancers in nuclear microenvironments of high Ubx concentrations. Related enhancers colocalize to the same microenvironments independently of their chromosomal location, suggesting that microenvironments are highly differentiated transcription domains. Manipulating the affinity of svb enhancers revealed an inverse relationship between enhancer affinity and Ubx concentration required for transcriptional activation. The Ubx cofactor, Homothorax (Hth), was co-enriched with Ubx near enhancers that require Hth, even though Ubx and Hth did not co-localize throughout the nucleus. Thus, microenvironments of high local transcription factor and cofactor concentrations could help low-affinity sites overcome their kinetic inefficiency. Mechanisms that generate these microenvironments could be a general feature of eukaryotic transcriptional regulation.

Functional regulatory evolution outside of the minimal even-skipped stripe 2 enhancer

Transcriptional enhancers are regions of DNA that drive precise patterns of gene expression. Although many studies have elucidated how individual enhancers can evolve, most of this work has focused on what are called 'minimal' enhancers, the smallest DNA regions that drive expression that approximates an aspect of native gene expression. Here, we explore how the Drosophila erecta even-skipped (eve) locus has evolved by testing its activity in the divergent D. melanogaster genome. We found, as has been reported previously, that the D. erecta eve stripe 2 enhancer (eveS2) fails to drive appreciable expression in D. melanogaster However, we found that a large transgene carrying the entire D. erecta eve locus drives normal eve expression, including in stripe 2. We performed a functional dissection of the region upstream of the D. erecta eveS2 region and found multiple Zelda motifs that are required for normal expression. Our results illustrate how sequences outside of minimal enhancer regions can evolve functionally through mechanisms other than changes in transcription factor-binding sites that drive patterning.

Uncoupling neurogenic gene networks in the Drosophila embryo

The EGF signaling pathway specifies neuronal identities in the Drosophila embryo by regulating developmental patterning genes such as intermediate neuroblasts defective (ind). EGFR is activated in the ventral midline and neurogenic ectoderm by the Spitz ligand, which is processed by the Rhomboid protease. CRISPR/Cas9 was used to delete defined rhomboid enhancers mediating expression at each site of Spitz processing. Surprisingly, the neurogenic ectoderm, not the ventral midline, was found to be the dominant source of EGF patterning activity. We suggest that Drosophila is undergoing an evolutionary transition in central nervous system (CNS)-organizing activity from the ventral midline to the neurogenic ectoderm.

Evolving Notch polyQ tracts reveal possible solenoid interference elements

Polyglutamine (polyQ) tracts in regulatory proteins are extremely polymorphic. As functional elements under selection for length, triplet repeats are prone to DNA replication slippage and indel mutations. Many polyQ tracts are also embedded within intrinsically disordered domains, which are less constrained, fast evolving, and difficult to characterize. To identify structural principles underlying polyQ tracts in disordered regulatory domains, here I analyze deep evolution of metazoan Notch polyQ tracts, which can generate alleles causing developmental and neurogenic defects. I show that Notch features polyQ tract turnover that is restricted to a discrete number of conserved “polyQ insertion slots”. Notch polyQ insertion slots are: (i) identifiable by an amphipathic “slot leader” motif; (ii) conserved as an intact C-terminal array in a 1-to-1 relationship with the N-terminal solenoid-forming ankyrin repeats (ARs); and (iii) enriched in carboxamide residues (Q/N), whose sidechains feature dual hydrogen bond donor and acceptor atoms. Correspondingly, the terminal loop and β-strand of each AR feature conserved carboxamide residues, which would be susceptible to folding interference by hydrogen bonding with residues outside the ARs. I thus suggest that Notch polyQ insertion slots constitute an array of AR interference elements (ARIEs). Notch ARIEs would dynamically compete with the delicate serial folding induced by adjacent ARs. Huntingtin, which harbors solenoid-forming HEAT repeats, also possesses a similar number of polyQ insertion slots. These results suggest that intrinsically disordered interference arrays featuring carboxamide and polyQ enrichment may constitute coupled proteodynamic modulators of solenoids.

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.

Integration of Orthogonal Signaling by the Notch and Dpp Pathways in Drosophila

The transcription factor Suppressor of Hairless and its coactivator, the Notch intracellular domain, are polyglutamine (pQ)-rich factors that target enhancer elements and interact with other locally bound pQ-rich factors. To understand the functional repertoire of such enhancers, we identify conserved regulatory belts with binding sites for the pQ-rich effectors of both Notch and BMP/Dpp signaling, and the pQ-deficient tissue selectors Apterous (Ap), Scalloped (Sd), and Vestigial (Vg). We find that the densest such binding site cluster in the genome is located in the BMP-inducible nab locus, a homolog of the vertebrate transcriptional cofactors NAB1/NAB2 We report three major findings. First, we find that this nab regulatory belt is a novel enhancer driving dorsal wing margin expression in regions of peak phosphorylated Mad in wing imaginal discs. Second, we show that Ap is developmentally required to license the nab dorsal wing margin enhancer (DWME) to read out Notch and Dpp signaling in the dorsal compartment. Third, we find that the nab DWME is embedded in a complex of intronic enhancers, including a wing quadrant enhancer, a proximal wing disc enhancer, and a larval brain enhancer. This enhancer complex coordinates global nab expression via both tissue-specific activation and interenhancer silencing. We suggest that DWME integration of BMP signaling maintains nab expression in proliferating margin descendants that have divided away from Notch-Delta boundary signaling. As such, uniform expression of genes like nab and vestigial in proliferating compartments would typically require both boundary and nonboundary lineage-specific enhancers.

The Soft Touch: Low-Affinity Transcription Factor Binding Sites in Development and Evolution

Transcription factor proteins regulate gene expression by binding to specific DNA regions. Most studies of transcription factor binding sites have focused on the highest affinity sites for each factor. There is abundant evidence, however, that binding sites with a range of affinities, including very low affinities, are critical to gene regulation. Here, we present the theoretical and experimental evidence for the importance of low-affinity sites in gene regulation and development. We also discuss the implications of the widespread use of low-affinity sites in eukaryotic genomes for robustness, precision, specificity, and evolution of gene regulation.

The Nature, Extent, and Consequences of Genetic Variation in the opa Repeats of Notch in Drosophila

Polyglutamine (pQ) tracts are abundant in proteins co-interacting on DNA. The lengths of these pQ tracts can modulate their interaction strengths. However, pQ tracts >40 residues are pathologically prone to amyloidogenic self-assembly. Here, we assess the extent and consequences of variation in the pQ-encoding opa repeats of Notch in Drosophila melanogaster. We use Sanger sequencing to genotype opa sequences (5′-CAX repeats), which have resisted assembly using short sequence reads. While most sampled lines carry the major allele opa31 encoding Q₁₃HQ₁₇ or the opa32 allele encoding Q₁₃HQ₁₈, many lines carry rare alleles encoding pQ tracts >32 residues: opa33a (Q₁₄HQ₁₈), opa33b (Q₁₅HQ₁₇), opa34 (Q₁₆HQ₁₇), opa35a1/opa35a2 (Q₁₃HQ₂₁), opa36 (Q₁₃HQ₂₂), and opa37 (Q₁₃HQ₂₃). Only one rare allele encodes a tract <31 residues: opa23 (Q₁₃–Q₁₀). This opa23 allele shortens the pQ tract while simultaneously eliminating the interrupting histidine. We introgressed these opa variant alleles into common backgrounds and measured the frequency of Notch-type phenotypes. Homozygotes for the short and long opa alleles have defects in embryonic survival and sensory bristle organ patterning, and sometimes show wing notching. Consistent with functional differences between Notch opa variants, we find that a scute inversion carrying the rare opa33b allele suppresses the bristle patterning defect caused by achaete/scute insufficiency, while an equivalent scute inversion carrying opa31 manifests the patterning defect. Our results demonstrate the existence of potent pQ variants of Notch and the need for long read genotyping of key repeat variables underlying gene regulatory networks.

Discovery and information-theoretic characterization of transcription factor binding sites that act cooperatively

Transcription factor binding to the surface of DNA regulatory regions is one of the primary causes of regulating gene expression levels. A probabilistic approach to model protein-DNA interactions at the sequence level is through position weight matrices (PWMs) that estimate the joint probability of a DNA binding site sequence by assuming positional independence within the DNA sequence. Here we construct conditional PWMs that depend on the motif signatures in the flanking DNA sequence, by conditioning known binding site loci on the presence or absence of additional binding sites in the flanking sequence of each site's locus. Pooling known sites with similar flanking sequence patterns allows for the estimation of the conditional distribution function over the binding site sequences. We apply our model to the Dorsal transcription factor binding sites active in patterning the Dorsal-Ventral axis of Drosophila development. We find that those binding sites that cooperate with nearby Twist sites on average contain about 0.5 bits of information about the presence of Twist transcription factor binding sites in the flanking sequence. We also find that Dorsal binding site detectors conditioned on flanking sequence information make better predictions about what is a Dorsal site relative to background DNA than detection without information about flanking sequence features.

Identifying transcriptional cis-regulatory modules in animal genomes

Gene expression is regulated through the activity of transcription factors (TFs) and chromatin-modifying proteins acting on specific DNA sequences, referred to as cis-regulatory elements. These include promoters, located at the transcription initiation sites of genes, and a variety of distal cis-regulatory modules (CRMs), the most common of which are transcriptional enhancers. Because regulated gene expression is fundamental to cell differentiation and acquisition of new cell fates, identifying, characterizing, and understanding the mechanisms of action of CRMs is critical for understanding development. CRM discovery has historically been challenging, as CRMs can be located far from the genes they regulate, have few readily identifiable sequence characteristics, and for many years were not amenable to high-throughput discovery methods. However, the recent availability of complete genome sequences and the development of next-generation sequencing methods have led to an explosion of both computational and empirical methods for CRM discovery in model and nonmodel organisms alike. Experimentally, CRMs can be identified through chromatin immunoprecipitation directed against TFs or histone post-translational modifications, identification of nucleosome-depleted 'open' chromatin regions, or sequencing-based high-throughput functional screening. Computational methods include comparative genomics, clustering of known or predicted TF-binding sites, and supervised machine-learning approaches trained on known CRMs. All of these methods have proven effective for CRM discovery, but each has its own considerations and limitations, and each is subject to a greater or lesser number of false-positive identifications. Experimental confirmation of predictions is essential, although shortcomings in current methods suggest that additional means of validation need to be developed. For further resources related to this article, please visit the WIREs website.

CONFLICT OF INTEREST

The authors have declared no conflicts of interest for this article.

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.

Regulation of broad by the Notch pathway affects timing of follicle cell development

During Drosophila oogenesis, activation of Notch signaling in the follicular epithelium (FE) around stage 6 of oogenesis is essential for entry into the endocycle and a series of other changes such as cell differentiation and migration of subsets of the follicle cells. Notch induces the expression of zinc finger protein Hindsight and suppresses homeodomain protein Cut to regulate the mitotic/endocycle (ME) switch. Here we report that broad (br), encoding a small group of zinc-finger transcription factors resulting from alternative splicing, is a transcriptional target of Notch nuclear effector Suppressor of Hairless (Su(H)). The early pattern of Br in the FE, uniformly expressed except in the polar cells, is established by Notch signaling around stage 6, through the binding of Su(H) to the br early enhancer (brE) region. Mutation of the Su(H) binding site leads to a significant reduction of brE reporter expression in follicle cells undergoing the endocycle. Chromatin immunoprecipitation results further confirm Su(H) binding to the br early enhancer. Consistent with its expression in follicle cells during midoogenesis, loss of br function results in a delayed entry into the endocycle. Our findings suggest an important role of br in the timing of follicle cell development, and its transcriptional regulation by the Notch pathway.

Drosophila functional elements are embedded in structurally constrained sequences

Modern functional genomics uncovered numerous functional elements in metazoan genomes. Nevertheless, only a small fraction of the typical non-exonic genome contains elements that code for function directly. On the other hand, a much larger fraction of the genome is associated with significant evolutionary constraints, suggesting that much of the non-exonic genome is weakly functional. Here we show that in flies, local (30–70 bp) conserved sequence elements that are associated with multiple regulatory functions serve as focal points to a pattern of punctuated regional increase in G/C nucleotide frequencies. We show that this pattern, which covers a region tenfold larger than the conserved elements themselves, is an evolutionary consequence of a shift in the balance between gain and loss of G/C nucleotides and that it is correlated with nucleosome occupancy across multiple classes of epigenetic state. Evidence for compensatory evolution and analysis of SNP allele frequencies show that the evolutionary regime underlying this balance shift is likely to be non-neutral. These data suggest that current gaps in our understanding of genome function and evolutionary dynamics are explicable by a model of sparse sequence elements directly encoding for function, embedded into structural sequences that help to define the local and global epigenomic context of such functional elements.

A key challenge in functional genomics is to predict evolutionary dynamics from functional annotation of the genome and vice versa. Modern epigenomic studies helped assign function to numerous new sequence elements, but left most of the genome essentially uncharacterized. Evolutionary genomics, on the other hand, consistently suggests that a much larger fraction of the un-annotated genome evolves under selective pressure. We hypothesize that this function-selection gap can be attributed to sequences that facilitate the physical organization of functional elements, such as transcription factor binding sites, within chromosomes. We exemplify this by studying in detail the sequences embedding small conserved elements (CEs) in Drosophila. We show that, while CEs have typically high AT content, high GC content levels around them are maintained by a non-neutral evolutionary balance between gain and loss of GC nucleotides. This non-uniform pattern is highly correlated with nucleosome organization around CEs, potentially imposing an evolutionary constraint on as much as one quarter of the genome. We suggest this can at least partly explain the above function-selection gap. Weak evolutionary constraints on “structural” sequences (at scales ranging from one nucleosome to recently described multi-megabase topological domains) may affect genome evolution just like structural motifs shape protein evolution.

Variation in the dorsal gradient distribution is a source for modified scaling of germ layers in Drosophila

Specification of germ layers along the dorsoventral axis by morphogenetic gradients is an ideal model to study scaling properties of gradients and cell fate changes during evolution. Classical anatomical studies in divergent insects (e.g., flies and grasshoppers) revealed that the neuroectodermal size is conserved and originates similar numbers of neuroblasts of homologous identity. In contrast, mesodermal domains vary significantly in closely related Drosophila species. To further investigate the underlying mechanisms of scaling of germ layers across Drosophila species, we quantified the Dorsal (Dl)/NF-κB gradient, the main morphogenetic gradient that initiates separation of the mesoderm, neuroectoderm, and ectoderm. We discovered a variable range of Toll activation across species and found that Dl activates mesodermal genes at the same threshold levels in melanogaster sibling species. We also show that the Dl gradient distribution can be modulated by nuclear size and packing densities. We propose that variation in mesodermal size occurs at a fast evolutionary rate and is an important mechanism to define the ventral boundary of the neuroectoderm.

A Schnurri/Mad/Medea complex attenuates the dorsal-twist gradient readout at vnd

Morphogen gradients are used in developing embryos, where they subdivide a field of cells into territories characterized by distinct cell fate potentials. Such systems require both a spatially-graded distribution of the morphogen, and an ability to encode different responses at different target genes. However, the potential for different temporal responses is also present because morphogen gradients typically provide temporal cues, which may be a potential source of conflict. Thus, a low threshold response adapted for an early temporal onset may be inappropriate when the desired spatial response is a spatially-limited, high-threshold expression pattern. Here, we identify such a case with the Drosophila vnd locus, which is a target of the dorsal (dl) nuclear concentration gradient that patterns the dorsal/ventral (D/V) axis of the embryo. The vnd gene plays a critical role in the "ventral dominance" hierarchy of vnd, ind, and msh, which individually specify distinct D/V neural columnar fates in increasingly dorsal ectodermal compartments. The role of vnd in this regulatory hierarchy requires early temporal expression, which is characteristic of low-threshold responses, but its specification of ventral neurogenic ectoderm demands a relatively high-threshold response to dl. We show that the Neurogenic Ectoderm Enhancer (NEE) at vnd takes additional input from the complementary Dpp gradient via a conserved Schnurri/Mad/Medea silencer element (SSE) unlike NEEs at brk, sog, rho, and vn. These results show how requirements for conflicting temporal and spatial responses to the same gradient can be solved by additional inputs from complementary gradients.

Lost in the map

Organismal development and evolution are complex, multifaceted processes that depend intimately on context. They are subject to environmental influences, chance appearance and fixation of mutations, and numerous other idiosyncrasies. Genomics is detailing the molecular signature of effects of these mechanisms on phenotypes, but because numerous distinct evolutionary explanations can produce a given genomic pattern, the molecular details, rather than elucidating process, typically distract from explanatory insight and contribute little to predictive capability. While genomic research has burgeoned, direct study of evolutionary and developmental processes has lagged. We advocate for reinvigoration of direct study of process, along with refocusing of attention on questions of broad biological import, as more productive of urgently needed insights, which genomic approaches are not providing.

Analysis of variation at transcription factor binding sites in Drosophila and humans

BACKGROUND

Advances in sequencing technology have boosted population genomics and made it possible to map the positions of transcription factor binding sites (TFBSs) with high precision. Here we investigate TFBS variability by combining transcription factor binding maps generated by ENCODE, modENCODE, our previously published data and other sources with genomic variation data for human individuals and Drosophila isogenic lines.

RESULTS

We introduce a metric of TFBS variability that takes into account changes in motif match associated with mutation and makes it possible to investigate TFBS functional constraints instance-by-instance as well as in sets that share common biological properties. We also take advantage of the emerging per-individual transcription factor binding data to show evidence that TFBS mutations, particularly at evolutionarily conserved sites, can be efficiently buffered to ensure coherent levels of transcription factor binding.

CONCLUSIONS

Our analyses provide insights into the relationship between individual and interspecies variation and show evidence for the functional buffering of TFBS mutations in both humans and flies. In a broad perspective, these results demonstrate the potential of combining functional genomics and population genetics approaches for understanding gene regulation.

A machine learning approach for identifying novel cell type-specific transcriptional regulators of myogenesis

Transcriptional enhancers integrate the contributions of multiple classes of transcription factors (TFs) to orchestrate the myriad spatio-temporal gene expression programs that occur during development. A molecular understanding of enhancers with similar activities requires the identification of both their unique and their shared sequence features. To address this problem, we combined phylogenetic profiling with a DNA-based enhancer sequence classifier that analyzes the TF binding sites (TFBSs) governing the transcription of a co-expressed gene set. We first assembled a small number of enhancers that are active in Drosophila melanogaster muscle founder cells (FCs) and other mesodermal cell types. Using phylogenetic profiling, we increased the number of enhancers by incorporating orthologous but divergent sequences from other Drosophila species. Functional assays revealed that the diverged enhancer orthologs were active in largely similar patterns as their D. melanogaster counterparts, although there was extensive evolutionary shuffling of known TFBSs. We then built and trained a classifier using this enhancer set and identified additional related enhancers based on the presence or absence of known and putative TFBSs. Predicted FC enhancers were over-represented in proximity to known FC genes; and many of the TFBSs learned by the classifier were found to be critical for enhancer activity, including POU homeodomain, Myb, Ets, Forkhead, and T-box motifs. Empirical testing also revealed that the T-box TF encoded by org-1 is a previously uncharacterized regulator of muscle cell identity. Finally, we found extensive diversity in the composition of TFBSs within known FC enhancers, suggesting that motif combinatorics plays an essential role in the cellular specificity exhibited by such enhancers. In summary, machine learning combined with evolutionary sequence analysis is useful for recognizing novel TFBSs and for facilitating the identification of cognate TFs that coordinate cell type-specific developmental gene expression patterns.

Evolutionary origins of transcription factor binding site clusters

Empirical studies have revealed that regulatory DNA sequences such as enhancers or promoters often harbor multiple binding sites for the same transcription factor. Such "homotypic site clustering" has been hypothesized as arising out of functional requirements of the sequences. Here, we propose an alternative explanation of this phenomenon that multisite enhancers are common because they are favored by evolutionary sampling of the genotype-phenotype landscape. To test this hypothesis, we developed a new computational framework specialized for population genetic simulations of enhancer evolution. It uses a thermodynamics-based model of enhancer function, integrating information from strong as well as weak binding sites, to determine the strength of selection. Using this framework, we found that even when simpler genotypes exist for a desired strength of regulation, relatively complex genotypes (enhancers with more sites) are more readily reached by the simulated evolutionary process. We show that there are more ways to "build" a fit genotype with many weak sites than with a few strong sites, and this is why evolution finds complex genotypes more often. Our claims are consistent with an empirical analysis of binding site content in enhancers characterized in Drosophila melanogaster and their orthologs in other Drosophila species. We also characterized a subtle but significant difference between genotypes likely to be sampled by evolution and equally fit genotypes one would obtain by uniform sampling of the fitness landscape, that is, an "evolutionary signature" in enhancer sequences. Finally, we investigated potential effects of other factors, such as rugged fitness landscapes, short local duplications, and noise characteristics of enhancers, on the emergence of homotypic site clustering. Homotypic site clustering is an important contributor to the complexity and function of cis-regulatory sequences. This work provides a simple null hypothesis for its origin, against which alternative adaptationist explanations may be evaluated, and cautions against "evolutionary mirages" present in common features of genomic sequence. The quantitative framework we develop here can be used more generally to understand how mechanisms of enhancer action influence their composition and evolution.

Rapid evolutionary rewiring of a structurally constrained eye enhancer

BACKGROUND

Enhancers are genomic cis-regulatory sequences that integrate spatiotemporal signals to control gene expression. Enhancer activity depends on the combination of bound transcription factors as well as-in some cases-the arrangement and spacing of binding sites for these factors. Here, we examine evolutionary changes to the sequence and structure of sparkling, a Notch/EGFR/Runx-regulated enhancer that activates the dPax2 gene in cone cells of the developing Drosophila eye.

RESULTS

Despite functional and structural constraints on its sequence, sparkling has undergone major reorganization in its recent evolutionary history. Our data suggest that the relative strengths of the various regulatory inputs into sparkling change rapidly over evolutionary time, such that reduced input from some factors is compensated by increased input from different regulators. These gains and losses are at least partly responsible for the changes in enhancer structure that we observe. Furthermore, stereotypical spatial relationships between certain binding sites ("grammar elements") can be identified in all sparkling orthologs-although the sites themselves are often recently derived. We also find that low binding affinity for the Notch-regulated transcription factor Su(H), a conserved property of sparkling, is required to prevent ectopic responses to Notch in noncone cells.

CONCLUSIONS

Rapid DNA sequence turnover does not imply either the absence of critical cis-regulatory information or the absence of structural rules. Our findings demonstrate that even a severely constrained cis-regulatory sequence can be significantly rewired over a short evolutionary timescale.