A context-dependent bifurcation in the Pointed transcriptional effector network contributes specificity and robustness to retinal cell fate acquisition.

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.

Natural variation of the expression pattern of the segmentation gene even-skipped in melanogaster.

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.

Rearrangements of 2.5 kilobases of noncoding DNA from the Drosophila even-skipped locus define predictive rules of genomic cis-regulatory logic.

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.

Consequences of eukaryotic enhancer architecture for gene expression dynamics, development, and fitness.

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.

Canalization of segmentation and its evolution in Drosophila.

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.

Functional evolution of noncoding DNA.

Noncoding DNA in eukaryotes encodes functionally important signals for the regulation of chromosome assembly, DNA replication, and gene expression. The increasing availability of whole-genome sequences of related taxa has led to interest in the evolution of these signals, and the phylogenetic footprints they produce. Cis-regulatory sequences controlling gene expression are often conserved among related species, but are rarely conserved between distantly related taxa. Several experimentally characterized regulatory elements have failed to show sequence similarity even between closely related species.

Functional evolution of a cis-regulatory module.

Lack of knowledge about how regulatory regions evolve in relation to their structure-function may limit the utility of comparative sequence analysis in deciphering cis-regulatory sequences. To address this we applied reverse genetics to carry out a functional genetic complementation analysis of a eukaryotic cis-regulatory module-the even-skipped stripe 2 enhancer-from four Drosophila species. The evolution of this enhancer is non-clock-like, with important functional differences between closely related species and functional convergence between distantly related species. Functional divergence is attributable to differences in activation levels rather than spatiotemporal control of gene expression. Our findings have implications for understanding enhancer structure-function, mechanisms of speciation and computational identification of regulatory modules.

Nanos plays a conserved role in axial patterning outside of the Diptera.

Axial patterning is a fundamental event in early development, and molecules involved in determining the body axes provide a coordinate system for subsequent patterning. While orthologs of Drosophila bicoid and nanos play a conserved role in anteroposterior (AP) patterning within at least a subset of Diptera, conservation of this process has not yet been demonstrated outside of the flies. Indeed, it has been argued that bicoid, an instrumental "anterior" factor in Drosophila melanogaster, acquired this role during the evolution of more-derived dipterans. Interestingly, the interaction of Drosophila maternal nanos and maternal hunchback provides a system for patterning the AP axis that is partially redundant to the anterior system. Previous studies in grasshoppers suggest that hunchback may play a conserved role in axial patterning in this insect, but this function may be supplied solely by the zygotic component of hunchback expression. Here we provide evidence that the early pattern of zygotic grasshopper Hunchback expression is achieved through translational repression that may be mediated through the action of grasshopper nanos. This is consistent with the notion that an anterior gradient system is not necessary in all insects and that the posterior pole "probably conveys longitudinal polarity on the ensuing germ anlage".

Kibra and Merlin Activate the Hippo Pathway Spatially Distinct from and Independent of Expanded.

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.

Naturally occurring deletions of hunchback binding sites in the even-skipped stripe 3+7 enhancer.

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.

Genetic complexity in a Drosophila model of diabetes-associated misfolded human proinsulin.

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.

Effect of genetic variation in a Drosophila model of diabetes-associated misfolded human proinsulin.

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.

Sex-specific pattern formation during early Drosophila development.

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.

Evolution and inheritance of early embryonic patterning in Drosophila simulans and D. sechellia.

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.

Artificial selection on egg size perturbs early pattern formation in Drosophila melanogaster.

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.