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
Subscription info for Chinese customers
We have a dedicated website for our Chinese customers. Please go to naturechina.com to subscribe to this journal.
Rent or Buy article
Get time limited or full article access on ReadCube.
All prices are NET prices.
The original images (cuticle preparations and embryo images, organized into zip files) and data files are available for download and are indexed at: https://www.embl.de/download/crocker/Dense_and_pleiotropic_regulatory_information_in_a_developmental_enhancer/index.html. All fly lines will be made available from the corresponding authors upon reasonable request.
CAD files and links to the software can also be found at Github: https://github.com/tfuqua95/Flyspresso-CAD-files.
Wittkopp, P. J. & Kalay, G. Cis-regulatory elements: Molecular mechanisms and evolutionary processes underlying divergence. Nat. Rev. Genet. 13, 59–69 (2011).
Crocker, J. & Ilsley, G. R. Using synthetic biology to study gene regulatory evolution. Curr. Opin. Genet. Dev. 47, 91–101 (2017).
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).
Patwardhan, R. P. et al. Massively parallel functional dissection of mammalian enhancers in vivo. Nat. Biotechnol. 30, 265–270 (2012).
Weingarten-Gabbay, S. et al. Systematic interrogation of human promoters. Genome Res. 29, 171–183 (2019).
de Boer, C. G. et al. Deciphering eukaryotic gene-regulatory logic with 100 million random promoters. Nat. Biotechnol. 38, 56–65 (2020).
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).
Crocker, J. et al. Low affinity binding site clusters confer hox specificity and regulatory robustness. Cell 160, 191–203 (2015).
Payre, F. Genetic control of epidermis differentiation in Drosophila. Int. J. Dev. Biol. 48, 207–215 (2004).
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).
Storey, J. D., Tibshirani, R. Statistical significance for genomewide studies. Proc. Natl Acad. Sci. USA 100, 9440–9445 (2003).
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).
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).
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).
Rastogi, C. et al. Accurate and sensitive quantification of protein-DNA binding affinity. Proc. Natl Acad. Sci. USA 115, E3692–E3701 (2018).
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).
Nagy, O. et al. Correlated Evolution of Two Copulatory Organs via a Single cis-Regulatory Nucleotide Change. Curr. Biol. 28, 3450–3457.e13 (2018).
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).
Vincent, B. J., Estrada, J. & DePace, A. H. The appeasement of Doug: a synthetic approach to enhancer biology. Integr. Biol. 8, 475–484 (2016).
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).
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).
Zabidi, M. A. et al. Enhancer-core-promoter specificity separates developmental and housekeeping gene regulation. Nature 518, 556–559 (2015).
Frankel, N. et al. Phenotypic robustness conferred by apparently redundant transcriptional enhancers. Nature 466, 490–493 (2010).
Tsai, A., Alves, M. R. & Crocker, J. Multi-enhancer transcriptional hubs confer phenotypic robustness. eLife 8, e45325 (2019).
Preger-Ben Noon, E. et al. Comprehensive analysis of a cis-regulatory region reveals pleiotropy in enhancer function. Cell Rep. 22, 3021–3031 (2018).
Crocker, J., Tsai, A. & Stern, D. L. A fully synthetic transcriptional platform for a multicellular eukaryote. Cell Rep. 18, 287–296 (2017).
Jacob, F. The Possible and the Actual (Univ. Washington Press, 1994).
Stern, D. L. & Sucena, E. Preparation of cuticles from unhatched first-instar Drosophila larvae. Cold Spring Harb. Protoc. 2011, 065532 (2011).
Tischer, C., Hilsenstein, V., Hanson, K. & Pepperkok, R. Adaptive fluorescence microscopy by online feedback image analysis. Methods Cell Biol. 123, 489–503 (2014).
Politi, A. Z. et al. Quantitative mapping of fluorescently tagged cellular proteins using FCS-calibrated four-dimensional imaging. Nat. Protoc. 13, 1445–1464 (2018).
Schindelin, J. et al. Fiji: an open-source platform for biological-image analysis. Nat. Methods 9, 676–682 (2012).
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).
Campbell, R. notBoxPlot https://github.com/raacampbell/notBoxPlot (2020).
Jonas. Violin Plots for Plotting Multiple Distributions (distributionPlot.m) https://uk.mathworks.com/matlabcentral/fileexchange/23661-violin-plots-for-plotting-multiple-distributions-distributionplot-m (2020).
Cock, P. J. A. et al. Biopython: freely available Python tools for computational molecular biology and bioinformatics. Bioinformatics 25, 1422–1423 (2009).
Virtanen, P. et al. SciPy 1.0: fundamental algorithms for scientific computing in Python. Nat. Methods 17, 261–272 (2020).
McKinney, W. Data structures for statistical computing in Python. Proc. 9th Python Sci. Conf. (2010).
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.
The authors declare no competing interests.
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
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.
a–d, 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. e–h, 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. i–k, 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.
a, Images acquired from automated imaging are compiled into a large montage image. b–d, 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). e–g, 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).
a–s, 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 (a – s).
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). c–p, 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. s–y, 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.
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.
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). c–j, 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 (c – j).
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.
About this article
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). https://doi.org/10.1038/s41586-020-2816-5
Scientific Reports (2021)
Nature Reviews Genetics (2021)
Nature Plants (2021)