Skip to main content

Thank you for visiting You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

Dense and pleiotropic regulatory information in a developmental enhancer


Changes in gene regulation underlie much of phenotypic evolution1. However, our understanding of the potential for regulatory evolution is biased, because most evidence comes from either natural variation or limited experimental perturbations2. Using an automated robotics pipeline, we surveyed an unbiased mutation library for a developmental enhancer in Drosophila melanogaster. We found that almost all mutations altered gene expression and that parameters of gene expression—levels, location, and state—were convolved. The widespread pleiotropic effects of most mutations may constrain the evolvability of developmental enhancers. Consistent with these observations, comparisons of diverse Drosophila larvae revealed apparent biases in the phenotypes influenced by the enhancer. Developmental enhancers may encode a higher density of regulatory information than has been appreciated previously, imposing constraints on regulatory evolution.

This is a preview of subscription content

Access options

Rent or Buy article

Get time limited or full article access on ReadCube.


All prices are NET prices.

Fig. 1: Most nucleotide mutations in E3N alter gene expression.
Fig. 2: Mutational scanning identifies a Hth binding site associated with a changed evolutionary phenotype.
Fig. 3: Ubx mutations often simultaneously change levels, timing, and locations of E3N expression.
Fig. 4: E3N enhancer architecture may constrain evolution.

Data availability

The original images (cuticle preparations and embryo images, organized into zip files) and data files are available for download and are indexed at: All fly lines will be made available from the corresponding authors upon reasonable request.

Code availability

CAD files and links to the software can also be found at Github:


  1. 1.

    Wittkopp, P. J. & Kalay, G. Cis-regulatory elements: Molecular mechanisms and evolutionary processes underlying divergence. Nat. Rev. Genet. 13, 59–69 (2011).

    Article  Google Scholar 

  2. 2.

    Crocker, J. & Ilsley, G. R. Using synthetic biology to study gene regulatory evolution. Curr. Opin. Genet. Dev. 47, 91–101 (2017).

    CAS  Article  Google Scholar 

  3. 3.

    Mogno, I., Kwasnieski, J. C. & Cohen, B. A. Massively parallel synthetic promoter assays reveal the in vivo effects of binding site variants. Genome Res. 23, 1908–1915 (2013).

    CAS  Article  Google Scholar 

  4. 4.

    Patwardhan, R. P. et al. Massively parallel functional dissection of mammalian enhancers in vivo. Nat. Biotechnol. 30, 265–270 (2012).

    CAS  Article  Google Scholar 

  5. 5.

    Weingarten-Gabbay, S. et al. Systematic interrogation of human promoters. Genome Res. 29, 171–183 (2019).

    CAS  Article  Google Scholar 

  6. 6.

    de Boer, C. G. et al. Deciphering eukaryotic gene-regulatory logic with 100 million random promoters. Nat. Biotechnol. 38, 56–65 (2020).

    Article  Google Scholar 

  7. 7.

    Duveau, F., Yuan, D. C., Metzger, B. P. H., Hodgins-Davis, A. & Wittkopp, P. J. Effects of mutation and selection on plasticity of a promoter activity in Saccharomyces cerevisiae. Proc. Natl Acad. Sci. USA 114, E11218–E11227 (2017).

    CAS  Article  Google Scholar 

  8. 8.

    Crocker, J. et al. Low affinity binding site clusters confer hox specificity and regulatory robustness. Cell 160, 191–203 (2015).

    CAS  Article  Google Scholar 

  9. 9.

    Payre, F. Genetic control of epidermis differentiation in Drosophila. Int. J. Dev. Biol. 48, 207–215 (2004).

    CAS  Article  Google Scholar 

  10. 10.

    Belliveau, N. M. et al. Systematic approach for dissecting the molecular mechanisms of transcriptional regulation in bacteria. Proc. Natl Acad. Sci. USA 115, E4796–E4805 (2018).

    CAS  Article  Google Scholar 

  11. 11.

    Storey, J. D., Tibshirani, R. Statistical significance for genomewide studies. Proc. Natl Acad. Sci. USA 100, 9440–9445 (2003).

    ADS  MathSciNet  CAS  Article  Google Scholar 

  12. 12.

    Pollard, K. S., Hubisz, M. J., Rosenbloom, K. R. & Siepel, A. Detection of nonneutral substitution rates on mammalian phylogenies. Genome Res. 20, 110–121 (2010).

    CAS  Article  Google Scholar 

  13. 13.

    Smith, J. M. et al. Developmental constraints and evolution: a perspective from the Mountain Lake Conference on Development and Evolution. Q. Rev. Biol. 60, 265–287 (1985).

    Article  Google Scholar 

  14. 14.

    Uller, T., Moczek, A. P., Watson, R. A., Brakefield, P. M. & Laland, K. N. Developmental bias and evolution: a regulatory network perspective. Genetics 209, 949–966 (2018).

    Article  Google Scholar 

  15. 15.

    Rastogi, C. et al. Accurate and sensitive quantification of protein-DNA binding affinity. Proc. Natl Acad. Sci. USA 115, E3692–E3701 (2018).

    CAS  Article  Google Scholar 

  16. 16.

    Chang, M. V., Chang, J. L., Gangopadhyay, A., Shearer, A. & Cadigan, K. M. Activation of wingless targets requires bipartite recognition of DNA by TCF. Curr. Biol. 18, 1877–1881 (2008).

    CAS  Article  Google Scholar 

  17. 17.

    Nagy, O. et al. Correlated Evolution of Two Copulatory Organs via a Single cis-Regulatory Nucleotide Change. Curr. Biol. 28, 3450–3457.e13 (2018).

    CAS  Article  Google Scholar 

  18. 18.

    Sabarís, G., Laiker, I., Preger-Ben Noon, E. & Frankel, N. Actors with multiple roles: pleiotropic enhancers and the paradigm of enhancer modularity. Trends Genet. 35, 423–433 (2019).

    Article  Google Scholar 

  19. 19.

    Vincent, B. J., Estrada, J. & DePace, A. H. The appeasement of Doug: a synthetic approach to enhancer biology. Integr. Biol. 8, 475–484 (2016).

    Article  Google Scholar 

  20. 20.

    Dey, S. S., Foley, J. E., Limsirichai, P., Schaffer, D. V. & Arkin, A. P. Orthogonal control of expression mean and variance by epigenetic features at different genomic loci. Mol. Syst. Biol. 11, 806 (2015).

    Article  Google Scholar 

  21. 21.

    Stern, D. L. et al. Genetic and transgenic reagents for Drosophila simulans, D. mauritiana, D. yakuba, D. santomea, and D. virilis. G3 7, 1339–1347 (2017).

    CAS  Google Scholar 

  22. 22.

    Zabidi, M. A. et al. Enhancer-core-promoter specificity separates developmental and housekeeping gene regulation. Nature 518, 556–559 (2015).

    ADS  CAS  Article  Google Scholar 

  23. 23.

    Frankel, N. et al. Phenotypic robustness conferred by apparently redundant transcriptional enhancers. Nature 466, 490–493 (2010).

    ADS  CAS  Article  Google Scholar 

  24. 24.

    Tsai, A., Alves, M. R. & Crocker, J. Multi-enhancer transcriptional hubs confer phenotypic robustness. eLife 8, e45325 (2019).

    Article  Google Scholar 

  25. 25.

    Preger-Ben Noon, E. et al. Comprehensive analysis of a cis-regulatory region reveals pleiotropy in enhancer function. Cell Rep. 22, 3021–3031 (2018).

    CAS  Article  Google Scholar 

  26. 26.

    Crocker, J., Tsai, A. & Stern, D. L. A fully synthetic transcriptional platform for a multicellular eukaryote. Cell Rep. 18, 287–296 (2017).

    CAS  Article  Google Scholar 

  27. 27.

    Jacob, F. The Possible and the Actual (Univ. Washington Press, 1994).

  28. 28.

    Stern, D. L. & Sucena, E. Preparation of cuticles from unhatched first-instar Drosophila larvae. Cold Spring Harb. Protoc. 2011, 065532 (2011).

    Google Scholar 

  29. 29.

    Tischer, C., Hilsenstein, V., Hanson, K. & Pepperkok, R. Adaptive fluorescence microscopy by online feedback image analysis. Methods Cell Biol. 123, 489–503 (2014).

    Article  Google Scholar 

  30. 30.

    Politi, A. Z. et al. Quantitative mapping of fluorescently tagged cellular proteins using FCS-calibrated four-dimensional imaging. Nat. Protoc. 13, 1445–1464 (2018).

    CAS  Article  Google Scholar 

  31. 31.

    Schindelin, J. et al. Fiji: an open-source platform for biological-image analysis. Nat. Methods 9, 676–682 (2012).

    CAS  Article  Google Scholar 

  32. 32.

    Arganda-Carreras, I. et al. in Computer Vision Approaches to Medical Image Analysis. CVAMIA 2006. Lecture Notes in Computer Science Vol. 4241 (eds Beichel, R. R. & Sonka, M.) (Springer, 2006).

  33. 33.

    Campbell, R. notBoxPlot (2020).

  34. 34.

    Jonas. Violin Plots for Plotting Multiple Distributions (distributionPlot.m) (2020).

  35. 35.

    Cock, P. J. A. et al. Biopython: freely available Python tools for computational molecular biology and bioinformatics. Bioinformatics 25, 1422–1423 (2009).

    CAS  Article  Google Scholar 

  36. 36.

    Virtanen, P. et al. SciPy 1.0: fundamental algorithms for scientific computing in Python. Nat. Methods 17, 261–272 (2020).

    CAS  Article  Google Scholar 

  37. 37.

    McKinney, W. Data structures for statistical computing in Python. Proc. 9th Python Sci. Conf. (2010).

Download references


A.T. is supported by the German Research Foundation (Deutsche Forschungsgemeinschaft, TS 428/1-1) and EMBL. R.S.M. is supported by funding from R35GM118336. We thank J. Zaugg and J. Wirbel for statistical advice; C. J. Standley, X. Li, R. M. Galupa, G. A. Canales, M. Perkins and M. R. P. Alves for comments on the manuscript; C. Rastogi and H. Bussemaker for the NRLB algorithms; K. Richter and J. Sager for technical support; I. Jones for mounting the Drosophila species larvae; T. Laverty for Drosophila assistance; A. Milberger for help with CAD design; and J. Payne for discussions.

Author information




T.F., M.E.v.B. and J.C. performed the measurements; R.S.M., D.L.S. and J.C. were involved in planning and supervised the work; T.F. and J.C. processed the experimental data, performed the analysis, drafted the manuscript and designed the figures with input from all the authors; J.J. and P.P. designed and built the robotics with input from J.C. A.H. and C.T. wrote the automated microscopy code and analysis tools. N.A. carried out the biochemistry experiments. T.F., A.T. and J.C. analysed the data.

Corresponding authors

Correspondence to David L. Stern or Justin Crocker.

Ethics declarations

Competing interests

The authors declare no competing interests.

Additional information

Peer review information Nature thanks Robb Krumlauf and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Extended data figures and tables

Extended Data Fig. 1 Distribution of mutations in the E3N enhancer library.

a, Mutant enhancer variants of E3N were created via degenerate PCR and integrated into the placZattB plasmid, which contains a minimalized core hsp70 promoter and the lacZ reporter gene. Plasmids were integrated into the Drosophila genome at the attP2 site. b, Pie chart depicting base-pair composition of the WT E3N enhancer. c, (Left) Histogram for all 749 mutants (dark red) is approximately normal with an average of 7 mutations per mutant. Magenta bars denote lines antibody stained (117 total), and blue lines indicate lines that were also Beta-Galactosidase stained (274 total). (Right) pie chart shows probability of mutation normalized to ATCG composition (see b). d, Manhattan plot shows the summation of all mutations within the E3N library. e, Unsmoothened “footprinting scores” from Fig. 1h. Scores plotted linearly over transcription factor binding motifs (colored and shaded regions) across the E3N genomic sequence.

Extended Data Fig. 2 An automated platform for fixing, staining, and imaging Drosophila embryos.

ad, Collecting Drosophila embryos. (a) Custom fly chambers were made, holding up to 24 different strains. (b) An explosion-view of the fly chambers. Embryo meshes (red) can attach and detach from the fly chambers and are suspended above an apple juice-agar plate. (c) Embryos are collected onto the embryo meshes and washed with saline solution and bleached. (d) Embryos are loaded into a fixation plate. eh, Components of the robot. (e) The fluid-dispensing manifold. Seven pneumo-hydraulic syringe pumps are coupled to the fluid-dispensing manifold; one pump for priming chemicals into the fluid-dispensing manifold, and six pumps for dispensing chemicals into the fixation plate. (f) The fluid-separating manifold uses 24 small syringes to aspirate fluid from the isotonic shocking attachments. (g and g’) Different components of the robot. (h) Cross-section of the fixation plate and aspiration tips and syringes. 24 small aspiration tips draw fluid from the top of each well within the isotonic shocking attachment and six main dual-purpose tips dispense and aspirate fluid into and out of the bottom of the wells. ik, The adaptive feedback imaging pipeline. (i) Samples are mounted on multi-well slides. (j) An overview tile-scan of each well is taken and x,y coordinates for embryos (green) are identified either manually or computationally. (k) For each coordinate, a fast, low-resolution confocal stack is automatically acquired. An algorithm determines the embryo’s z position and rotation, yielding a bounding box within which a high-resolution, 3D stack of the entire embryo is acquired. See also Methods. l, Control E3N WT embryos were fixed and stained on the robot. A single embryo in the same orientation and age from each well was selected and the individual nuclear fluorescence intensities were measured in AU, arbitrary units of fluorescence intensity. In plots, centre line is mean, upper and lower limits are standard deviation.

Extended Data Fig. 3 Methods of image and data analysis.

a, Images acquired from automated imaging are compiled into a large montage image. bd, Registering multiple images using fiducials. An embryo acquired during automated imaging (b) can be automatically rotated in 3D space using ELAV (teal) as a fiducial. Once properly rotated, maximum projections of the ventral half can be computed (c). Finally, the 2-D projections can be elastically registered – or deformed – to align multiple samples (d). eg, Methods of measuring expression patterns. (e) Sliding window analysis. A box is drawn between A2 and A3 and centred within A2. Multiple measurements are taken, sliding the box across the stripe. Each point on the boxplot represents one measurement within the box. In box plots, centre line is mean, upper and lower limits are standard deviation and whiskers show 95% CIs. (f) State method analysis. A row of cell-sized regions of interest are dragged down across the A2 stripe. Each point on the boxplot represents a single cell. (g) Plot profile analysis. A box is drawn from the A1 to A5 and the mean intensity is taken for each column of pixels and plotted (N = 10 embryos). Shaded areas indicate ± 1 SD, solid line is the mean expression. Scale bars, 100 μm. Embryos are matched to scale respectively in (a) and in (b-e), and (g).

Extended Data Fig. 4 Single base pair mutations and E3N conservation.

as, Example embryos carrying individual E3N::lacZ variants with single mutations. Constructs are ordered from smallest to largest effect sizes. t, u, PhyloP scores across the E3N enhancer sequence. Locations of the single mutations and their PhyloP scores are highlighted as magenta bars. v, E3N sequence alignment between 10 Drosophila species. Scale bars, 100 μm (a). Embryos are matched to scale respectively (as).

Extended Data Fig. 5 Testing additional Hth-Exd motifs in E3N.

a, b, Embryos carrying E3N::lacZ reporter constructs in a WT w1118 background (a) and hth homeodomain-less (HthHM) hth 100.1 background stained with anti-β-Galactosidase (b). cp, Embryos carrying E3N::lacZ reporter constructs stained with anti-β-Galactosidase adjacent to their respective expression plot profiles. Constructs contain mutations in Hth1 (CTGGCA → CCCCCC), Hth2 (TGACAA → CCCCCC), Hth3 (TTGTCG → CCCCCC), and Hth4 (TGAGAG → CCCCCC). (c and d) E3N WT. (e and f) E3N with Hth1 site changed. (g and h) E3N with Hth2 site changed. (i and j) E3N with Hth1 and Hth2 sites changed. (k and l) E3N with Ubx3 site changed (CATAATTTGT → CAGGGTTTGT). (m and n) E3N with Hth3 and Hth4 sites changed. (o and p) E3N with Hth1, Hth2, Hth3, and Hth4 sites changed. In all plots, the black and magenta lines denote the average expression driven by the wild-type and modified enhancers, respectively (n = 10 for each genotype). Shaded areas indicate ± 1 s.d. AU, arbitrary units of fluorescence intensity. q, top, Schematic for the E3N enhancer, denoting binding sites and possible protein-to-protein interactions. q, bottom, Schematic for different E3N fragments tested. r, Multiple-species alignment of Hth1, 2 and the UBX-Exd site. sy, Electromobility shift assays (EMSA) for different fragments of E3N denoted in (q). All EMSAs were run on native (non-denatured) gels. HthHM/Exd is the homeodomain-less (HthHM) isoform of Hth incubated with Exd. HthFL/Exd is the Hth isoform with a homeodomain, incubated with Exd. Fragments tested with the WT Hth binding site and a mutated form. (s) EMSA for fragment-f with Hth2 mutated (t) additionally with increasing Ubx concentrations. (u) EMSA for fragment-a with Hth1 and 2 mutated. (v) EMSA for fragment-a and fragment-b with Hth3 and Hth4 mutated. (w) EMSA for fragment-a, fragment-c, and fragment-d. (x) EMSA for fragment-e with Hth1 mutated. (y) EMSA for fragment-g with Hth2 mutated. Scale bars, 100 μm. Embryos are matched to scale.

Extended Data Fig. 6 The effects of Ubx affinity on morphology.

a, b, Schematic output from NRLB15 shows predicted binding affinity for Exd::Ubx heterodimers across the E3N sequence, where black peaks are on the 5′ strand and red peaks respectively on the 3′ strand. Affinity plots are shown for Drosophila melanogaster (a) and Drosophila virilis (b). c, d, Drosophila cuticle preps for flies with WT E3N driving svb cDNA (c), or E3N with increased Ubx binding affinity driving svb cDNA (d). Trichomes were counted within a region of interest (teal box) defined by anatomic epithelial sensory cells (*). Arrows and brackets demarcate ectopic trichomes. e, Boxplots comparing trichome numbers in the A1 segment in the region of interest from panels (c) and (d) (n = 13, P < 0.02), see also Tsai et al., 201824. In box plots, centre line is mean, upper and lower limits are standard deviation and whiskers show 95% CIs. Scale bars, 25 μm each.

Extended Data Fig. 7 Extensive pleiotropic effects across the E3N enhancer.

a, b, Plot comparing the percent of lines with pleiotropic or ectopic expression versus the number of mutations based on antibody staining (a) and Beta-Galactosidase staining (b). cj, A subset of mutants with pleiotropic effects. (c) Line 145-2 drives ectopic expression in the developing wing and haltere discs (7/7 embryos). (d) Line 139-6 drives wider stripes and increased expression, as well as ectopic expression between the stripes and (e) on the dorsal side, (5/5 embryos). (f) Line 40-8 drives a split stripe pattern, where the middle row of nuclei within the stripes is not active (6/6 embryos). (g) Line 93-4 expression varies along the anterior-posterior axis (5/5 embryos). (h) Line 77-9 drives ectopic expression in the salivary glands (5/5 embryos). (i) Line 81-7 drives expression in the developing mouth hooks (5/5 embryos.) (j) Line 15-2v activates expression at stage 10 and drives ectopic expression throughout the embryo in multiple developmental stages (14/14 embryos). k, Plot of footprinting scores versus E3N sequence. Magenta is the footprinting score (σi, see methods). The higher the peak, the higher probability that a mutation will change expression. Gray plots are the mutation coverage for the number of lines screened per base (Mi, see methods). l, EWAC scores represent p-values from a log of odds ratio test for the association of a mutation changing expression. Dashed lines denote p- and q-values11, respectively. See Materials and Methods. m, Plot comparing the percent of lines with changed expression for mutations in the overlapping Pan/Hth site. n, Quantification of the staining intensities in the stripe and naked domains with the indicated reporter construct using the “sliding window” technique (Extended Data Fig. 3e). N = the number of embryos per line measured, from left to right: 10, 10, 7, 9, 10, 7, 10, 6, 10, 8, 4, 8, 10, 10, 10. In box plots, centre line is mean, upper and lower limits are standard deviation and whiskers show 95% CIs. Scale bars, 100 μm. Embryos are matched to scale (cj).

Extended Data Fig. 8 Cuticle preps from 60 Drosophila species across approximately 100 million years of evolution.

a, Phylogenetic tree of Drosophila species studied here, spanning approximately 150 million years of evolution. Red indicates a loss of trichomes. b, Representative cuticle preps for Drosophila species. See also Fig. 4. Scale bars, 25 μm each.

Supplementary information

Rights and permissions

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Fuqua, T., Jordan, J., van Breugel, M.E. et al. Dense and pleiotropic regulatory information in a developmental enhancer. Nature 587, 235–239 (2020).

Download citation

Further reading


By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.


Quick links

Nature Briefing

Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.

Get the most important science stories of the day, free in your inbox. Sign up for Nature Briefing