Zeba Wunderlich

Depleting Gene Activities in Early Drosophila Embryos with the “Maternal-Gal4–shRNA” System

In a developing Drosophila melanogaster embryo, mRNAs have a maternal origin, a zygotic origin, or both. During the maternal–zygotic transition, maternal products are degraded and gene expression comes under the control of the zygotic genome. To interrogate the function of mRNAs that are both maternally and zygotically expressed, it is common to examine the embryonic phenotypes derived from female germline mosaics. Recently, the development of RNAi vectors based on short hairpin RNAs (shRNAs) effective during oogenesis has provided an alternative to producing germline clones. Here, we evaluate the efficacies of: (1) maternally loaded shRNAs to knockdown zygotic transcripts and (2) maternally loaded Gal4 protein to drive zygotic shRNA expression. We show that, while Gal4-driven shRNAs in the female germline very effectively generate phenotypes for genes expressed maternally, maternally loaded shRNAs are not very effective at generating phenotypes for early zygotic genes. However, maternally loaded Gal4 protein is very efficient at generating phenotypes for zygotic genes expressed during mid-embryogenesis. We apply this powerful and simple method to unravel the embryonic functions of a number of pleiotropic genes.

A Conserved Developmental Patterning Network Produces Quantitatively Different Output in Multiple Species of Drosophila

Differences in the level, timing, or location of gene expression can contribute to alternative phenotypes at the molecular and organismal level. Understanding the origins of expression differences is complicated by the fact that organismal morphology and gene regulatory networks could potentially vary even between closely related species. To assess the scope of such changes, we used high-resolution imaging methods to measure mRNA expression in blastoderm embryos of Drosophila yakuba and Drosophila pseudoobscura and assembled these data into cellular resolution atlases, where expression levels for 13 genes in the segmentation network are averaged into species-specific, cellular resolution morphological frameworks. We demonstrate that the blastoderm embryos of these species differ in their morphology in terms of size, shape, and number of nuclei. We present an approach to compare cellular gene expression patterns between species, while accounting for varying embryo morphology, and apply it to our data and an equivalent dataset for Drosophila melanogaster. Our analysis reveals that all individual genes differ quantitatively in their spatio-temporal expression patterns between these species, primarily in terms of their relative position and dynamics. Despite many small quantitative differences, cellular gene expression profiles for the whole set of genes examined are largely similar. This suggests that cell types at this stage of development are conserved, though they can differ in their relative position by up to 3–4 cell widths and in their relative proportion between species by as much as 5-fold. Quantitative differences in the dynamics and relative level of a subset of genes between corresponding cell types may reflect altered regulatory functions between species. Our results emphasize that transcriptional networks can diverge over short evolutionary timescales and that even small changes can lead to distinct output in terms of the placement and number of equivalent cells.

Different gene regulation strategies revealed by analysis of binding motifs

Coordinated regulation of gene expression relies on transcription factors (TFs) binding to specific DNA sites. Our large-scale information-theoretical analysis of > 950 TF-binding motifs demonstrates that prokaryotes and eukaryotes use strikingly different strategies to target TFs to specific genome locations. Although bacterial TFs can recognize a specific DNA site in the genomic background, eukaryotic TFs exhibit widespread, nonfunctional binding and require clustering of sites to achieve specificity. We find support for this mechanism in a range of experimental studies and in our evolutionary analysis of DNA-binding domains. Our systematic characterization of binding motifs provides a quantitative assessment of the differences in transcription regulation in prokaryotes and eukaryotes.

Shadow enhancers enable Hunchback bifunctionality in the Drosophila embryo

Significance Enhancers are regions of regulatory DNA that control gene expression and cell fate decisions during development. Enhancers compute the expression pattern of their target gene by reading the concentrations of input regulatory proteins. Many developmental genes contain multiple enhancers that control the same output pattern, but it is unclear whether these enhancers all compute the pattern in the same way. We use measurements in single cells and computational models in Drosophila embryos to demonstrate that two enhancers that encode the same gene expression pattern compute differently: the same regulatory protein represses one enhancer and activates the other. Pairs of enhancers that generate the same pattern by performing different computations may impart special properties to developmental systems. Hunchback (Hb) is a bifunctional transcription factor that activates and represses distinct enhancers. Here, we investigate the hypothesis that Hb can activate and repress the same enhancer. Computational models predicted that Hb bifunctionally regulates the even-skipped (eve) stripe 3+7 enhancer (eve3+7) in Drosophila blastoderm embryos. We measured and modeled eve expression at cellular resolution under multiple genetic perturbations and found that the eve3+7 enhancer could not explain endogenous eve stripe 7 behavior. Instead, we found that eve stripe 7 is controlled by two enhancers: the canonical eve3+7 and a sequence encompassing the minimal eve stripe 2 enhancer (eve2+7). Hb bifunctionally regulates eve stripe 7, but it executes these two activities on different pieces of regulatory DNA—it activates the eve2+7 enhancer and represses the eve3+7 enhancer. These two “shadow enhancers” use different regulatory logic to create the same pattern.

Dissecting sources of quantitative gene expression pattern divergence between Drosophila species

Gene expression patterns can diverge between species due to changes in a genes regulatory DNA or changes in the proteins, e.g., transcription factors (TFs), that regulate the gene. We developed a modeling framework to uncover the sources of expression differences in blastoderm embryos of three Drosophila species, focusing on the regulatory circuit controlling expression of the hunchback (hb) posterior stripe. Using this framework and cellular‐resolution expression measurements of hb and its regulating TFs, we found that changes in the expression patterns of hbs TFs account for much of the expression divergence. We confirmed our predictions using transgenic D. melanogaster lines, which demonstrate that this set of orthologous cis‐regulatory elements (CREs) direct similar, but not identical, expression patterns. We related expression pattern differences to sequence changes in the CRE using a calculation of the CREs TF binding site content. By applying this calculation in both the transgenic and endogenous contexts, we found that changes in binding site content affect sensitivity to regulating TFs and that compensatory evolution may occur in circuit components other than the CRE.

A gene expression atlas of a bicoid -depleted Drosophila embryo reveals early canalization of cell fate

In developing embryos, gene regulatory networks drive cells towards discrete terminal fates, a process called canalization. We studied the behavior of the anterior-posterior segmentation network in Drosophila melanogaster embryos by depleting a key maternal input, bicoid (bcd), and measuring gene expression patterns of the network at cellular resolution. This method results in a gene expression atlas containing the levels of mRNA or protein expression of 13 core patterning genes over six time points for every cell of the blastoderm embryo. This is the first cellular resolution dataset of a genetically perturbed Drosophila embryo that captures all cells in 3D. We describe the technical developments required to build this atlas and how the method can be employed and extended by others. We also analyze this novel dataset to characterize the degree and timing of cell fate canalization in the segmentation network. We find that in two layers of this gene regulatory network, following depletion of bcd, individual cells rapidly canalize towards normal cell fates. This result supports the hypothesis that the segmentation network directly canalizes cell fate, rather than an alternative hypothesis whereby cells are initially mis-specified and later eliminated by apoptosis. Our gene expression atlas provides a high resolution picture of a classic perturbation and will enable further computational modeling of canalization and gene regulation in this transcriptional network. Summary: Drosophila bicoid mutant embryos show severe patterning phenotypes, but individual cells retain wild-type fates - as judged by expression of key patterning genes at single-cell resolution.

Comparing mRNA levels using in situ hybridization of a target gene and co-stain

In situ hybridization is an important technique for measuring the spatial expression patterns of mRNA in cells, tissues, and whole animals. However, mRNA levels cannot be compared across experiments using typical protocols. Here we present a semi-quantitative method to compare mRNA levels of a gene across multiple samples. This method yields an estimate of the error in the measurement to allow statistical comparison. Our method uses a typical in situ hybridization protocol to stain for a target gene and an internal standard, which we refer to as a co-stain. As a proof of concept, we apply this method to multiple lines of transgenic Drosophila embryos, harboring constructs that express reporter genes to different levels. We generated this test set by mutating enhancer sequences to contain different numbers of binding sites for Zelda, a transcriptional activator. We demonstrate that using a co-stain with in situ hybridization is an effective method to compare mRNA levels across samples. This method requires only minor modifications to existing in situ hybridization protocols and uses straightforward analysis techniques. This strategy can be broadly applied to detect quantitative, spatially resolved changes in mRNA levels.

How gene order is influenced by the biophysics of transcription regulation

What are the forces that shape the structure of prokaryotic genomes: the order of genes, their proximity, and their orientation? Coregulation and coordinated horizontal gene transfer are believed to promote the proximity of functionally related genes and the formation of operons. However, forces that influence the structure of the genome beyond the level of a single operon remain unknown. Here, we show that the biophysical mechanism by which regulatory proteins search for their sites on DNA can impose constraints on genome structure. Using simulations, we demonstrate that rapid and reliable gene regulation requires that the transcription factor (TF) gene be close to the site on DNA the TF has to bind, thus promoting the colocalization of TF genes and their targets on the genome. We use parameters that have been measured in recent experiments to estimate the relevant length and times scales of this process and demonstrate that the search for a cognate site may be prohibitively slow if a TF has a low copy number and is not colocalized. We also analyze TFs and their sites in a number of bacterial genomes, confirm that they are colocalized significantly more often than expected, and show that this observation cannot be attributed to the pressure for coregulation or formation of selfish gene clusters, thus supporting the role of the biophysical constraint in shaping the structure of prokaryotic genomes. Our results demonstrate how spatial organization can influence timing and noise in gene expression.

Modeling transcriptional networks in Drosophila development at multiple scales

Quantitative models of developmental processes can provide insights at multiple scales. Ultimately, models may be particularly informative for key questions about network level behavior during development such as how does the system respond to environmental perturbation, or operate reliably in different genetic backgrounds? The transcriptional networks that pattern the Drosophila embryo have been the subject of numerous quantitative experimental studies coupled to modeling frameworks in recent years. In this review, we describe three studies that consider these networks at different levels of molecular detail and therefore result in different types of insights. We also discuss other developmental transcriptional networks operating in Drosophila, with the goal of highlighting what additional insights they may provide.

Krüppel Expression Levels Are Maintained through Compensatory Evolution of Shadow Enhancers

Many developmental genes are controlled by shadow enhancers—pairs of enhancers that drive overlapping expression patterns. We hypothesized that compensatory evolution can maintain the total expression of a gene, while individual shadow enhancers diverge between species. To test this hypothesis, we analyzed expression driven by orthologous pairs of shadow enhancers from Drosophila melanogaster, Drosophila yakuba, and Drosophila pseudoobscura that control expression of Krüppel, a transcription factor that patterns the anterior-posterior axis of blastoderm embryos. We found that the expression driven by the pair of enhancers is conserved between these three species, but expression levels driven by the individual enhancers are not. Using sequence analysis and experimental perturbation, we show that each shadow enhancer is regulated by different transcription factors. These results support the hypothesis that compensatory evolution can occur between shadow enhancers, which has implications for mechanistic and evolutionary studies of gene regulation.