A variety of evolutionary processes shape the structure of cis -regulatory elements such as enhancers. Functionally important regulatory sequences tend to be conserved as result of negative selection against deleterious mutations and positive selection for better-canalized performance. However, some forms of stabilizing selection can maintain functional conservation of cis -regulatory elements for long periods of evolutionary time despite structural transcription binding sites turnover (Ludwig et al, 1998). In addition, compensatory evolution can even accelerate the substitution process in large populations to level greater than the neutral rate of substitution (Carter and Wagner, 2002). The effects of lethal mutations in cis -regulatory regions on viability in heterozygotes and the impact of such mutations on evolutionary processes have not been substantially addressed. We ask a fundamental biological question: How does a lethal mutant in an enhancer region affect viability when carried over a dominant “normal” allele? Here cis -regulatory mutation eve ΔMSE of even-skipped gene ( eve) in D. melanogaster , that is lethal in homozygote, has been tested for relative viability in heterozygous state. In eve ΔMSE the 480-bp fragment corresponding to MSE (Minimal Stipe Element) of the eve stripe 2 enhancer was deleted and replaced with an unrelated DNA fragment (containing white + gene marker) by ends-out homologous recombination (Ludwig et al, 2011). We discovered that relative viability of flies carrying cis -regulatory mutation eve ΔMSE in heterozygote was not reduced compared to the viability of wild-type flies. In contrast, the viability of heterozygotes carrying a homozygote-lethal nonsense eve R13 mutation (Fujioka et al, 1999) was significantly impacted. Our explanation to this surprising phenomenon relays on action of eve transcriptional autoregulation. The transcription factor (TF) EVE regulates its own expression level through the autoregulatory cis -regulatory module. It could elevate the expression level from one dose to the level sufficient to restore up to 100% fitness in the heterozygote with eve ΔMSE lethal allele. This example suggests that some cis -regulatory mutants (including the deleterious ones) may exist in populations as heterozygotes with high fitness for many generations, before possibly encountering an epistatic interaction with a compensatory mutation at a different site in the genome. Our study is consistent with the compensatory mechanisms of evolution for functionally important cis -regulatory elements.

Spatiotemporally precise and robust cell fate transitions, which depend on specific signaling cues, are fundamental to the development of appropriately patterned tissues. The fidelity and precision with which photoreceptor fates are recruited in the Drosophila eye exemplifies these principles. The fly eye consists of a highly ordered array of ~750 ommatidia, each of which contains eight distinct photoreceptors, R1-R8, specified sequentially in a precise spatial pattern. Recruitment of R1-R7 fates requires reiterative receptor tyrosine kinase / mitogen activated protein kinase (MAPK) signaling mediated by the transcriptional effector Pointed (Pnt). However the overall signaling levels experienced by R2-R5 cells are distinct from those experienced by R1, R6 and R7. A relay mechanism between two Pnt isoforms initiated by MAPK activation directs the universal transcriptional response. Here we ask how the generic Pnt response is tailored to these two rounds of photoreceptor fate transitions. We find that during R2-R5 specification PntP2 is coexpressed with a closely related but previously uncharacterized isoform, PntP3. Using CRISPR/Cas9-generated isoform specific null alleles we show that under otherwise wild type conditions, R2-R5 fate specification is robust to loss of either PntP2 or PntP3, and that the two activate pntP1 redundantly; however under conditions of reduced MAPK activity, both are required. Mechanistically, our data suggest that intrinsic activity differences between PntP2 and PntP3, combined with positive and unexpected negative transcriptional auto- and cross-regulation, buffer first-round fates against conditions of compromised RTK signaling. In contrast, in a mechanism that may be adaptive to the stronger signaling environment used to specify R1, R6 and R7 fates, the Pnt network resets to a simpler topology in which PntP2 uniquely activates pntP1 and auto-activates its own transcription. We propose that differences in expression patterns, transcriptional activities and regulatory interactions between Pnt isoforms together facilitate context-appropriate cell fate specification in different signaling environments.

The Hippo pathway is emerging as a key evolutionarily conserved signaling mechanism that controls organ size. Three membrane-associated proteins, Kibra, Merlin, and Expanded, regulate pathway activity, but the precise molecular mechanism by which they function is still poorly understood. Here we provide evidence that Merlin and Kibra activate Hippo signaling in parallel to Expanded at a spatially distinct cellular domain, the medial apical cortex. Merlin and Kibra together recruit the adapter protein Salvador, which in turn recruits the core kinase Hippo. In addition, we show that Crumbs has a dual effect on Hippo signaling. Crumbs promotes the ability of Expanded to activate the pathway but also sequesters Kibra to downregulate Hippo signaling. Together, our findings elucidate the mechanism of Hippo pathway activation by Merlin and Kibra, identify a subcellular domain for Hippo pathway regulation, and demonstrate differential activity of upstream regulators in different subcellular domains.

The evolution of canalized traits is a central question in evolutionary biology. Natural variation in highly conserved traits can provide clues about their evolutionary potential. Here we investigate natural variation in a conserved trait—even-skipped (eve) expression at the cellular blastoderm stage of embryonic development in Drosophila melanogaster. Expression of the pair-rule gene eve was quantitatively measured in three inbred lines derived from a natural population of D. melanogaster. One line showed marked differences in the spacing, amplitude and timing of formation of the characteristic seven-striped pattern over a 50-min period prior to the onset of gastrulation. Stripe 5 amplitude and the width of the interstripe between stripes 4 and 5 were both reduced in this line, while the interstripe distance between stripes 3 and 4 was increased. Engrailed expression in stage 10 embryos revealed a statistically significant increase in the length of parasegment 6 and a decrease in the length of parasegments 8 and 9. These changes are larger than those previously reported between D. melanogaster and D. pseudoobscura, two species that are thought to have diverged from a common ancestor over 25 million years ago. This line harbors a rare 448 bp deletion in the first intron of knirps (kni). This finding suggested that reduced Kni levels caused the deviant eve expression, and indeed we observed lower levels of Kni protein at early cycle 14A in L2 compared to the other two lines. A second of the three lines displayed an approximately 20% greater level of expression for all seven eve stripes. The three lines are each viable and fertile, and none display a segmentation defect as adults, suggesting that early-acting variation in eve expression is ameliorated by developmental buffering mechanisms acting later in development. Canalization of the segmentation pathway may reduce the fitness consequences of genetic variation, thus allowing the persistence of mutations with unexpectedly strong gene expression phenotypes.

Changes in regulatory DNA contribute to phenotypic differences within and between taxa. Comparative studies show that many transcription factor binding sites (TFBS) are conserved between species whereas functional studies reveal that some mutations segregating within species alter TFBS function. Consistently, in this analysis of 13 regulatory elements in Drosophila melanogaster populations, single base and insertion/deletion polymorphism are rare in characterized regulatory elements. Experimentally defined TFBS are nearly devoid of segregating mutations and, as has been shown before, are quite conserved. For instance 8 of 11 Hunchback binding sites in the stripe 3+7 enhancer of even-skipped are conserved between D. melanogaster and Drosophila virilis. Oddly, we found a 72 bp deletion that removes one of these binding sites (Hb8), segregating within D. melanogaster. Furthermore, a 45 bp deletion polymorphism in the spacer between the stripe 3+7 and stripe 2 enhancers, removes another predicted Hunchback site. These two deletions are separated by ∼250 bp, sit on distinct haplotypes, and segregate at appreciable frequency. The Hb8Δ is at 5 to 35% frequency in the new world, but also shows cosmopolitan distribution. There is depletion of sequence variation on the Hb8Δ-carrying haplotype. Quantitative genetic tests indicate that Hb8Δ affects developmental time, but not viability of offspring. The Eve expression pattern differs between inbred lines, but the stripe 3 and 7 boundaries seem unaffected by Hb8Δ. The data reveal segregating variation in regulatory elements, which may reflect evolutionary turnover of characterized TFBS due to drift or co-evolution.

The identification and validation of gene-gene interactions is a major challenge in human studies. Here, we explore an approach for studying epistasis in humans using a Drosophila melanogaster model of neonatal diabetes mellitus. Expression of the mutant preproinsulin (hINS(C96Y)) in the eye imaginal disc mimics the human disease: it activates conserved stress-response pathways and leads to cell death (reduction in eye area). Dominant-acting variants in wild-derived inbred lines from the Drosophila Genetics Reference Panel produce a continuous, highly heritable distribution of eye-degeneration phenotypes in a hINS(C96Y) background. A genome-wide association study (GWAS) in 154 sequenced lines identified a sharp peak on chromosome 3L, which mapped to a 400-bp linkage block within an intron of the gene sulfateless (sfl). RNAi knockdown of sfl enhanced the eye-degeneration phenotype in a mutant-hINS-dependent manner. RNAi against two additional genes in the heparan sulfate (HS) biosynthetic pathway (ttv and botv), in which sfl acts, also modified the eye phenotype in a hINS(C96Y)-dependent manner, strongly suggesting a novel link between HS-modified proteins and cellular responses to misfolded proteins. Finally, we evaluated allele-specific expression difference between the two major sfl-intronic haplotypes in heterozygtes. The results showed significant heterogeneity in marker-associated gene expression, thereby leaving the causal mutation(s) and its mechanism unidentified. In conclusion, the ability to create a model of human genetic disease, map a QTL by GWAS to a specific gene, and validate its contribution to disease with available genetic resources and the potential to experimentally link the variant to a molecular mechanism demonstrate the many advantages Drosophila holds in determining the genetic underpinnings of human disease.

Rearrangements of about 2.5 kilobases of regulatory DNA located 5' of the transcription start site of the Drosophila even-skipped locus generate large-scale changes in the expression of even-skipped stripes 2, 3, and 7. The most radical effects are generated by juxtaposing the minimal stripe enhancers MSE2 and MSE3 for stripes 2 and 3 with and without small "spacer" segments less than 360 bp in length. We placed these fusion constructs in a targeted transformation site and obtained quantitative expression data for these transformants together with their controlling transcription factors at cellular resolution. These data demonstrated that the rearrangements can alter expression levels in stripe 2 and the 2-3 interstripe by a factor of more than 10. We reasoned that this behavior would place tight constraints on possible rules of genomic cis-regulatory logic. To find these constraints, we confronted our new expression data together with previously obtained data on other constructs with a computational model. The model contained representations of thermodynamic protein-DNA interactions including steric interference and cooperative binding, short-range repression, direct repression, activation, and coactivation. The model was highly constrained by the training data, which it described within the limits of experimental error. The model, so constrained, was able to correctly predict expression patterns driven by enhancers for other Drosophila genes; even-skipped enhancers not included in the training set; stripe 2, 3, and 7 enhancers from various Drosophilid and Sepsid species; and long segments of even-skipped regulatory DNA that contain multiple enhancers. The model further demonstrated that elevated expression driven by a fusion of MSE2 and MSE3 was a consequence of the recruitment of a portion of MSE3 to become a functional component of MSE2, demonstrating that cis-regulatory "elements" are not elementary objects.

Drosophila melanogaster has been widely used as a model of human Mendelian disease, but its value in modeling complex disease has received little attention. Fly models of complex disease would enable high-resolution mapping of disease-modifying loci and the identification of novel targets for therapeutic intervention. Here, we describe a fly model of permanent neonatal diabetes mellitus and explore the complexity of this model. The approach involves the transgenic expression of a misfolded mutant of human preproinsulin, hINS(C96Y), which is a cause of permanent neonatal diabetes. When expressed in fly imaginal discs, hINS(C96Y) causes a reduction of adult structures, including the eye, wing, and notum. Eye imaginal discs exhibit defects in both the structure and the arrangement of ommatidia. In the wing, expression of hINS(C96Y) leads to ectopic expression of veins and mechano-sensory organs, indicating disruption of wild-type signaling processes regulating cell fates. These readily measurable "disease" phenotypes are sensitive to temperature, gene dose, and sex. Mutant (but not wild-type) proinsulin expression in the eye imaginal disc induces IRE1-mediated XBP1 alternative splicing, a signal for endoplasmic reticulum stress response activation, and produces global change in gene expression. Mutant hINS transgene tester strains, when crossed to stocks from the Drosophila Genetic Reference Panel, produce F1 adults with a continuous range of disease phenotypes and large broad-sense heritability. Surprisingly, the severity of mutant hINS-induced disease in the eye is not correlated with that in the notum in these crosses, nor with eye reduction phenotypes caused by the expression of two dominant eye mutants acting in two different eye development pathways, Drop (Dr) or Lobe (L), when crossed into the same genetic backgrounds. The tissue specificity of genetic variability for mutant hINS-induced disease has, therefore, its own distinct signature. The genetic dominance of disease-specific phenotypic variability in our model of misfolded human proinsulin makes this approach amenable to genome-wide association study in a simple F1 screen of natural variation.

The deleterious effects of different X-chromosome dosage in males and females are buffered by a process called dosage compensation, which in Drosophila is achieved through a doubling of X-linked transcription in males. The male-specific lethal complex mediates this process, but is known to act only after gastrulation. Recent work has shown that the transcription of X-linked genes is also upregulated in males prior to gastrulation; whether it results in functional dosage compensation is not known. Absent or partial early dosage compensation raises the possibility of sex-biased expression of key developmental genes, such as the segmentation genes controlling anteroposterior patterning. We assess the functional output of early dosage compensation by measuring the expression of even-skipped (eve) with high spatiotemporal resolution in male and female embryos. We show that eve has a sexually dimorphic pattern, suggesting an interaction with either X-chromosome dose or the sex determination system. By manipulating the gene copy number of an X-linked transcription factor, giant (gt), we traced sex-biased eve patterning to gt dose, indicating that early dosage compensation is functionally incomplete. Despite sex-biased eve expression, the gene networks downstream of eve are able to produce sex-independent segmentation, a point that we establish by measuring the proportions of segments in elongated germ-band embryos. Finally, we use a whole-locus eve transgene with modified cis regulation to demonstrate that segment proportions have a sex-dependent sensitivity to subtle changes in Eve expression. The sex independence of downstream segmentation despite this sensitivity to Eve expression implies that additional autosomal gene- or pathway-specific mechanisms are required to ameliorate the effects of partial early dosage compensation.

Pattern formation in Drosophila embryogenesis has been widely investigated as a developmental and evolutionary model of robustness. To ask whether genetic variation for pattern formation is suppressed in this system, artificial selection for divergent egg size was used to challenge the scaling of even-skipped (eve) pattern formation in mitotic cycle 14 (stage 5) embryos of Drosophila melanogaster. Three-dimensional confocal imaging revealed shifts in the allometry of eve pair-rule stripes along both anterior-posterior (A-P) and dorsoventral (D-V) axes as a correlated response to egg size selection, indicating the availability of genetic variation for this buffered trait. Environmental perturbation was not required for the manifestation of this variation. The number of nuclei at the cellular blastoderm stage also changed in response to selection, with large-egg selected lines having more than 1000 additional nuclei relative to small-egg lines. This increase in nuclear number in larger eggs does not scale with egg size, however, as nuclear density is inversely correlated with egg length. Nuclear density varies along the A-P axis but does not correlate with the shift in eve stripe allometry between the selection treatments. Despite its macroevolutionary conservation, both eve stripe patterning and blastoderm cell number vary genetically both within and between closely related species.

The regulatory logic of time- and tissue-specific gene expression has mostly been dissected in the context of the smallest DNA fragments that, when isolated, recapitulate native expression in reporter assays. It is not known if the genomic sequences surrounding such fragments, often evolutionarily conserved, have any biological function or not. Using an enhancer of the even-skipped gene of Drosophila as a model, we investigate the functional significance of the genomic sequences surrounding empirically identified enhancers. A 480 bp long “minimal stripe element” is able to drive even-skipped expression in the second of seven stripes but is embedded in a larger region of 800 bp containing evolutionarily conserved binding sites for required transcription factors. To assess the overall fitness contribution made by these binding sites in the native genomic context, we employed a gene-replacement strategy in which whole-locus transgenes, capable of rescuing even-skipped- lethality to adulthood, were substituted for the native gene. The molecular phenotypes were characterized by tagging Even-skipped with a fluorescent protein and monitoring gene expression dynamics in living embryos. We used recombineering to excise the sequences surrounding the minimal enhancer and site-specific transgenesis to create co-isogenic strains differing only in their stripe 2 sequences. Remarkably, the flanking sequences were dispensable for viability, proving the sufficiency of the minimal element for biological function under normal conditions. These sequences are required for robustness to genetic and environmental perturbation instead. The mutant enhancers had measurable sex- and dose-dependent effects on viability. At the molecular level, the mutants showed a destabilization of stripe placement and improper activation of downstream genes. Finally, we demonstrate through live measurements that the peripheral sequences are required for temperature compensation. These results imply that seemingly redundant regulatory sequences beyond the minimal enhancer are necessary for robust gene expression and that “robustness” itself must be an evolved characteristic of the wild-type enhancer.

Pattern formation in Drosophila is a widely studied example of a robust developmental system. Such robust systems pose a challenge to adaptive evolution, as they mask variation that selection may otherwise act upon. Yet we find variation in the localization of expression domains (henceforth "stripe allometry") in the pattern formation pathway. Specifically, we characterize differences in the gap genes giant and Kruppel, and the pair-rule gene even-skipped, which differ between the sibling species Drosophila simulans and D. sechellia. In a double-backcross experiment, stripe allometry is consistent with maternal inheritance of stripe positioning and multiple genetic factors, with a distinct genetic basis from embryo length. Embryos produced by F1 and F2 backcross mothers exhibit novel spatial patterns of gene expression relative to the parental species, with no measurable increase in positional variance among individuals. Buffering of novel spatial patterns in the backcross genotypes suggests that robustness need not be disrupted in order for the trait to evolve, and perhaps the system is incapable of evolving to prevent the expression of all genetic variation. This limitation, and the ability of natural selection to act on minute genetic differences that are within the "margin of error" for the buffering mechanism, indicates that developmentally buffered traits can evolve without disruption of robustness.

Segmentation in Drosophila embryogenesis occurs through a hierarchical cascade of regulatory gene expression driven by the establishment of a diffusion-mediated morphogen gradient. Here, we investigate the response of this pattern formation process to genetic variation and evolution in egg size. Specifically, we ask whether spatial localization of gap genes Kruppel (Kr) and giant (gt) and the pair-rule gene even-skipped (eve) during cellularization is robust to genetic variation in embryo length in three Drosophila melanogaster isolines and two closely related species. We identified two wild-derived strains of D. melanogaster whose eggs differ by approximately 25% in length when reared under identical conditions. These two lines, a D. melanogaster laboratory stock (w1118), and offspring from crosses between the lines all exhibit precise scaling in the placement of gap and pair-rule gene expression along the anterior-posterior axis in relation to embryo length. Genetic analysis indicates that this scaling is maternally controlled. Maternal regulation of scaling must be required for consistent localization of segmentation gene expression because embryo size, a genetically variable and adaptive trait, is maternally inherited. We also investigated spatial scaling between these D. melanogaster lines and single lines of Drosophila sechellia and Drosophila simulans, the latter two differing by approximately 25% in egg length. In contrast to the robust scaling we observed within species, localization of gene expression relative to embryo length differs significantly between the three species. Thus, the developmental mechanism that assures robust scaling within a species does not prevent rapid evolution between species.

Eukaryotic gene expression is mediated by compact cis-regulatory modules, or enhancers, which are bound by specific sets of transcription factors1. The combinatorial interaction of these bound transcription factors determines time- and tissue-specific gene activation or repression. The even-skipped stripe 2 element controls the expression of the second transverse stripe of even-skipped messenger RNA in Drosophila melanogaster embryos, and is one of the best characterized eukaryotic enhancers2,3,4. Although even-skipped stripe 2 expression is strongly conserved in Drosophila, the stripe 2 element itself has undergone considerable evolutionary change in its binding-site sequences and the spacing between them. We have investigated this apparent contradiction, and here we show that two chimaeric enhancers, constructed by swapping the 5′ and 3′ halves of the native stripe 2 elements of two species, no longer drive expression of a reporter gene in the wild-type pattern. Sequence differences between species have functional consequences, therefore, but they are masked by other co-evolved differences. On the basis of these results, we present a model for the evolution of eukaryotic regulatory sequences.

PEB-me is a predominant protein of mature Drosophila melanogaster ejaculatory bulbs. It is resolved into four or five closely spaced subfractions (apparent molecular weight 35-39 kD) by sodium dodecyl sulfate-polyacrylamide gel electrophoresis. Four electrophoretic variants of PEB-me differing in apparent molecular weight by 200-800 daltons were found. These appear to be controlled by four alleles of a gene (peb) located by recombination and deletion mapping to the 60F1-2 region of chromosome 2. A minor ejaculatory bulb protein of ca. 80 kD (hPEB) was found to be immunochemically related to PEB and possibly encoded by peb. PEB is not detected by immunoblotting techniques in virgin females, in male tissues other than the ejaculatory bulb, or during developmental stages preceding the formation of this organ. The results of transplantations of genital imaginal discs and of immature ejaculatory bulbs between two strains having different PEB alleles suggest that the ejaculatory bulb is the site of PEB synthesis. In flies mutant for tra, tra-2, dsx, or ix, tissue specificity of PEB localization is retained and the protein is found whenever the ejaculatory bulb is formed, regardless of the chromosomal sex of the fly. The protein is transferred into the female genital duct during mating, where it can be detected for up to 12 hr. Possible functions of PEB in Drosophila reproduction are discussed.