Transcriptional repression by the Drosophila giant protein: cis element positioning provides an alternative means of interpreting an effector gradient

Deep molecular learning of transcriptional control of a synthetic CRE enhancer and its variants

Massively parallel reporter assay measures transcriptional activities of various cis-regulatory modules (CRMs) in a single experiment. We developed a thermodynamic computational model framework that calculates quantitative levels of gene expression directly from regulatory DNA sequences. Using the framework, we investigated the molecular mechanisms of cis-regulatory mutations of a synthetic enhancer that cause abnormal gene expression. We found that, in a human cell line, competitive binding between family transcription factors (TFs) with slightly different binding preferences significantly increases the accuracy of recapitulating the transcriptional effects of thousands of single- or multi-mutations. We also discovered that even if various harmful mutations occurred in an activator binding site, CRM could stably maintain or even increase gene expression through a certain form of competitive binding between family TFs. These findings enhance understanding the effect of SNPs and indels on CRMs and would help building robust custom-designed CRMs for biologics production and gene therapy.

Global repression by tailless during segmentation

The orphan nuclear receptor Tailless (Tll) exhibits conserved roles in brain formation and maintenance that are shared, for example, with vertebrate orthologous forms (Tlx). However, the early expression of tll in two gap domains in the segmentation cascade of Drosophila is unusual even for most other insects. Here we investigate tll regulation on pair-rule stripes. With ectopic misexpression of tll we detected unexpected repression of almost all pair-rule stripes of hairy (h), even-skipped (eve), runt (run), and fushi-tarazu (ftz). Examining Tll embryonic ChIP-chip data with regions mapped as Cis-Regulatory Modules (CRMs) of pair-rule stripes we verified Tll interactions to these regions. With the ChIP-chip data we also verified Tll interactions to the CRMs of gap domains and in the misexpression assay, Tll-mediated repression on Kruppel (Kr), kni (kni) and giant (gt) according to their differential sensitivity to Tll. These results with gap genes confirmed previous data from the literature and argue against indirect repression roles of Tll in the striped pattern. Moreover, the prediction of Tll binding sites in the CRMs of eve stripes and the mathematical modeling of their removal using an experimentally validated theoretical framework shows effects on eve stripes compatible with the absence of a repressor binding to the CRMs. In addition, modeling increased tll levels in the embryo results in the differential repression of eve stripes, agreeing well with the results of the misexpression assay. In genetic assays we investigated eve 5, that is strongly repressed by the ectopic domain and representative of more central stripes not previously implied to be under direct regulation of tll. While this stripe is little affected in tll-, its posterior border is expanded in gt- but detected with even greater expansion in gt-;tll-. We end up by discussing tll with key roles in combinatorial repression mechanisms to contain the expression of medial patterns of the segmentation cascade in the extremities of the embryo.

Predictive modeling reveals that higher-order cooperativity drives transcriptional repression in a synthetic developmental enhancer

A challenge in quantitative biology is to predict output patterns of gene expression from knowledge of input transcription factor patterns and from the arrangement of binding sites for these transcription factors on regulatory DNA. We tested whether widespread thermodynamic models could be used to infer parameters describing simple regulatory architectures that inform parameter-free predictions of more complex enhancers in the context of transcriptional repression by Runt in the early fruit fly embryo. By modulating the number and placement of Runt binding sites within an enhancer, and quantifying the resulting transcriptional activity using live imaging, we discovered that thermodynamic models call for higher-order cooperativity between multiple molecular players. This higher-order cooperativity captures the combinatorial complexity underlying eukaryotic transcriptional regulation and cannot be determined from simpler regulatory architectures, highlighting the challenges in reaching a predictive understanding of transcriptional regulation in eukaryotes and calling for approaches that quantitatively dissect their molecular nature.

Transcriptional Enhancers in Drosophila

Key discoveries in Drosophila have shaped our understanding of cellular "enhancers." With a special focus on the fly, this chapter surveys properties of these adaptable cis-regulatory elements, whose actions are critical for the complex spatial/temporal transcriptional regulation of gene expression in metazoa. The powerful combination of genetics, molecular biology, and genomics available in Drosophila has provided an arena in which the developmental role of enhancers can be explored. Enhancers are characterized by diverse low- or high-throughput assays, which are challenging to interpret, as not all of these methods of identifying enhancers produce concordant results. As a model metazoan, the fly offers important advantages to comprehensive analysis of the central functions that enhancers play in gene expression, and their critical role in mediating the production of phenotypes from genotype and environmental inputs. A major challenge moving forward will be obtaining a quantitative understanding of how these cis-regulatory elements operate in development and disease.

Signal Integration by Shadow Enhancers and Enhancer Duplications Varies across the Drosophila Embryo

Transcription of developmental genes is controlled by multiple enhancers. Frequently, more than one enhancer can activate transcription from the same promoter in the same cells. How is regulatory information from multiple enhancers combined to determine the overall expression output? We measure nascent transcription driven by a pair of shadow enhancers, each enhancer of the pair separately, and each duplicated, using live imaging in Drosophila embryos. This set of constructs allows us to quantify the input-output function describing signal integration by two enhancers. We show that signal integration performed by these shadow enhancers and duplications varies across the expression pattern, implying that how their activities are combined depends on the transcriptional regulators bound to the enhancers in different parts of the embryo. Characterizing signal integration by multiple enhancers is a critical step in developing conceptual and computational models of gene expression at the locus level, where multiple enhancers control transcription together.

Inference of Transcription Factor Regulation Patterns Using Gene Expression Covariation in Natural Populations of Drosophila melanogaster

Gene regulatory networks control the complex programs that drive development. Deciphering the connections between transcription factors (TFs) and target genes is challenging, in part because TFs bind to thousands of places in the genome but control expression through a subset of these binding events. We hypothesize that we can combine natural variation of expression levels and predictions of TF binding sites to identify TF targets. We gather RNA-seq data from 71 genetically distinct F1 Drosophila melanogaster embryos and calculate the correlations between TF and potential target genes' expression levels, which we call "regulatory strength." To separate direct and indirect TF targets, we hypothesize that direct TF targets will have a preponderance of binding sites in their upstream regions. Using 14 TFs active during embryogenesis, we find that 12 TFs showed a significant correlation between their binding strength and regulatory strength on downstream targets, and 10 TFs showed a significant correlation between the number of binding sites and the regulatory effect on target genes. The general roles, e.g. bicoid's role as an activator, and the particular interactions we observed between our TFs, e.g. twist's role as a repressor of sloppy paired and odd paired, generally coincide with the literature.

Different modes of enhancer-specific regulation by Runt and Even-skipped during Drosophila segmentation

Expression of the Drosophila slp1 gene depends on nonadditive interactions between two cis-regulatory enhancers. These enhancers are repressed by preventing either Pol II recruitment or release of promoter-proximal paused Pol II in a manner that is both enhancer and transcription factor specific and can account for their nonadditive interaction.

The initial metameric expression of the Drosophila sloppy paired 1 (slp1) gene is controlled by two distinct cis-regulatory DNA elements that interact in a nonadditive manner to integrate inputs from transcription factors encoded by the pair-rule segmentation genes. We performed chromatin immunoprecipitation on reporter genes containing these elements in different embryonic genotypes to investigate the mechanism of their regulation. The distal early stripe element (DESE) mediates both activation and repression by Runt. We find that the differential response of DESE to Runt is due to an inhibitory effect of Fushi tarazu (Ftz) on P-TEFb recruitment and the regulation of RNA polymerase II (Pol II) pausing. The proximal early stripe element (PESE) is also repressed by Runt, but in this case, Runt prevents PESE-dependent Pol II recruitment and preinitiation complex (PIC) assembly. PESE is also repressed by Even-skipped (Eve), but, of interest, this repression involves regulation of P-TEFb recruitment and promoter-proximal Pol II pausing. These results demonstrate that the mode of slp1 repression by Runt is enhancer specific, whereas the mode of repression of the slp1 PESE enhancer is transcription factor specific. We propose a model based on these differential regulatory interactions that accounts for the nonadditive interactions between the PESE and DESE enhancers during Drosophila segmentation.

Quantitative perturbation-based analysis of gene expression predicts enhancer activity in early Drosophila embryo

Enhancers constitute one of the major components of regulatory machinery of metazoans. Although several genome-wide studies have focused on finding and locating enhancers in the genomes, the fundamental principles governing their internal architecture and cis-regulatory grammar remain elusive. Here, we describe an extensive, quantitative perturbation analysis targeting the dorsal-ventral patterning gene regulatory network (GRN) controlled by Drosophila NF-κB homolog Dorsal. To understand transcription factor interactions on enhancers, we employed an ensemble of mathematical models, testing effects of cooperativity, repression, and factor potency. Models trained on the dataset correctly predict activity of evolutionarily divergent regulatory regions, providing insights into spatial relationships between repressor and activator binding sites. Importantly, the collective predictions of sets of models were effective at novel enhancer identification and characterization. Our study demonstrates how experimental dataset and modeling can be effectively combined to provide quantitative insights into cis-regulatory information on a genome-wide scale.

DOI:http://dx.doi.org/10.7554/eLife.08445.001

DNA contains regions known as genes, which may be “transcribed” to produce the RNA molecules that act as templates for building proteins and regulate cell activity. Proteins called transcription factors can bind to specific sequences of DNA to influence whether nearby genes are transcribed. For example, so-called enhancer regions of DNA contain several binding sites for transcription factors, and this binding activates gene transcription. Little is known about how the transcription factor binding sites are organized in enhancer regions, which makes it difficult to use DNA sequence information alone to predict the regulation of genes.

A transcription factor called Dorsal controls the activity of a network of genes that plays a crucial role in the development of fruit fly embryos. Dorsal binds to the enhancer region of a gene called rhomboid, which has been well studied and is known to be a fairly typical example of an enhancer region.

To understand the regulatory information encoded in the DNA sequences of enhancers, Sayal, Dresch et al. have now used a technique called perturbation analysis to investigate the interactions that are likely to occur between Dorsal and other transcription factors as they bind to the rhomboid enhancer. This technique involves systematically mutating the enhancer to remove different combinations of transcription factor binding sites and quantitatively investigating the effect this has on gene activity. A large set of mathematical models were then trained using this data and shown to correctly predict the activity of a range of other gene regulatory regions. The collective predictions of the models identified new enhancer regions and revealed details about how different types of transcription factor binding sites are arranged within enhancers.

As we enter an era where the DNA sequences of entire human populations are increasingly accessible, we would like to know the functional significance of changes in gene regulatory regions. Sayal, Dresch et al. show that the regulatory properties of specific control proteins are accessible by employing quantitative experiments and mathematical models. Similar studies will be required to learn how mutations found across the genome may alter gene expression, leading to better diagnosis and treatment of disease.

DOI:http://dx.doi.org/10.7554/eLife.08445.002

Morphogen rules: design principles of gradient-mediated embryo patterning

The Drosophila blastoderm and the vertebrate neural tube are archetypal examples of morphogen-patterned tissues that create precise spatial patterns of different cell types. In both tissues, pattern formation is dependent on molecular gradients that emanate from opposite poles. Despite distinct evolutionary origins and differences in time scales, cell biology and molecular players, both tissues exhibit striking similarities in the regulatory systems that establish gene expression patterns that foreshadow the arrangement of cell types. First, signaling gradients establish initial conditions that polarize the tissue, but there is no strict correspondence between specific morphogen thresholds and boundary positions. Second, gradients initiate transcriptional networks that integrate broadly distributed activators and localized repressors to generate patterns of gene expression. Third, the correct positioning of boundaries depends on the temporal and spatial dynamics of the transcriptional networks. These similarities reveal design principles that are likely to be broadly applicable to morphogen-patterned tissues.

Sequence-based model of gap gene regulatory network

BACKGROUND

The detailed analysis of transcriptional regulation is crucially important for understanding biological processes. The gap gene network in Drosophila attracts large interest among researches studying mechanisms of transcriptional regulation. It implements the most upstream regulatory layer of the segmentation gene network. The knowledge of molecular mechanisms involved in gap gene regulation is far less complete than that of genetics of the system. Mathematical modeling goes beyond insights gained by genetics and molecular approaches. It allows us to reconstruct wild-type gene expression patterns in silico, infer underlying regulatory mechanism and prove its sufficiency.

RESULTS

We developed a new model that provides a dynamical description of gap gene regulatory systems, using detailed DNA-based information, as well as spatial transcription factor concentration data at varying time points. We showed that this model correctly reproduces gap gene expression patterns in wild type embryos and is able to predict gap expression patterns in Kr mutants and four reporter constructs. We used four-fold cross validation test and fitting to random dataset to validate the model and proof its sufficiency in data description. The identifiability analysis showed that most model parameters are well identifiable. We reconstructed the gap gene network topology and studied the impact of individual transcription factor binding sites on the model output. We measured this impact by calculating the site regulatory weight as a normalized difference between the residual sum of squares error for the set of all annotated sites and for the set with the site of interest excluded.

CONCLUSIONS

The reconstructed topology of the gap gene network is in agreement with previous modeling results and data from literature. We showed that 1) the regulatory weights of transcription factor binding sites show very weak correlation with their PWM score; 2) sites with low regulatory weight are important for the model output; 3) functional important sites are not exclusively located in cis-regulatory elements, but are rather dispersed through regulatory region. It is of importance that some of the sites with high functional impact in hb, Kr and kni regulatory regions coincide with strong sites annotated and verified in Dnase I footprint assays.

In silico evolution of the hunchback gene indicates redundancy in cis-regulatory organization and spatial gene expression

Biological development depends on the coordinated expression of genes in time and space. Developmental genes have extensive cis-regulatory regions which control their expression. These regions are organized in a modular manner, with different modules controlling expression at different times and locations. Both how modularity evolved and what function it serves are open questions. We present a computational model for the cis-regulation of the hunchback (hb) gene in the fruit fly (Drosophila). We simulate evolution (using an evolutionary computation approach from computer science) to find the optimal cis-regulatory arrangements for fitting experimental hb expression patterns. We find that the cis-regulatory region tends to readily evolve modularity. These cis-regulatory modules (CRMs) do not tend to control single spatial domains, but show a multi-CRM/multi-domain correspondence. We find that the CRM-domain correspondence seen in Drosophila evolves with a high probability in our model, supporting the biological relevance of the approach. The partial redundancy resulting from multi-CRM control may confer some biological robustness against corruption of regulatory sequences. The technique developed on hb could readily be applied to other multi-CRM developmental genes.

Ancestral resurrection of the Drosophila S2E enhancer reveals accessible evolutionary paths through compensatory change

Upstream regulatory sequences that control gene expression evolve rapidly, yet the expression patterns and functions of most genes are typically conserved. To address this paradox, we have reconstructed computationally and resurrected in vivo the cis-regulatory regions of the ancestral Drosophila eve stripe 2 element and evaluated its evolution using a mathematical model of promoter function. Our feed-forward transcriptional model predicts gene expression patterns directly from enhancer sequence. We used this functional model along with phylogenetics to generate a set of possible ancestral eve stripe 2 sequences for the common ancestors of 1) D. simulans and D. sechellia; 2) D. melanogaster, D. simulans, and D. sechellia; and 3) D. erecta and D. yakuba. These ancestral sequences were synthesized and resurrected in vivo. Using a combination of quantitative and computational analysis, we find clear support for functional compensation between the binding sites for Bicoid, Giant, and Krüppel over the course of 40-60 My of Drosophila evolution. We show that this compensation is driven by a coupling interaction between Bicoid activation and repression at the anterior and posterior border necessary for proper placement of the anterior stripe 2 border. A multiplicity of mechanisms for binding site turnover exemplified by Bicoid, Giant, and Krüppel sites, explains how rapid sequence change may occur while maintaining the function of the cis-regulatory element.

A synthetic biology approach to the development of transcriptional regulatory models and custom enhancer design

Synthetic biology offers novel opportunities for elucidating transcriptional regulatory mechanisms and enhancer logic. Complex cis-regulatory sequences--like the ones driving expression of the Drosophila even-skipped gene--have proven difficult to design from existing knowledge, presumably due to the large number of protein-protein interactions needed to drive the correct expression patterns of genes in multicellular organisms. This work discusses two novel computational methods for the custom design of enhancers that employ a sophisticated, empirically validated transcriptional model, optimization algorithms, and synthetic biology. These synthetic elements have both utilitarian and academic value, including improving existing regulatory models as well as evolutionary questions. The first method involves the use of simulated annealing to explore the sequence space for synthetic enhancers whose expression output fit a given search criterion. The second method uses a novel optimization algorithm to find functionally accessible pathways between two enhancer sequences. These paths describe a set of mutations wherein the predicted expression pattern does not significantly vary at any point along the path. Both methods rely on a predictive mathematical framework that maps the enhancer sequence space to functional output.

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.

Nucleosomes and the accessibility problem

Eukaryotic DNA is packaged in nucleosomes. How does this sequestration affect the ability of transcription regulators to access their sites? We cite evidence against the idea that nucleosome positioning is determined primarily by the intrinsic propensities of DNA sequences to form nucleosomes--such that, for example, regulatory sites would be 'nucleosome-free'. Instead, studies in yeast show that nucleosome positioning is primarily determined by specific DNA-binding proteins. Where nucleosomes would otherwise compete with regulatory protein binding (a modest but potentially biologically important effect), this obstacle can be relieved by at least two strategies for exposing regulatory sites. In contrast to their lack of effect on nucleosome positioning, DNA sequence differences do directly affect both the efficiencies with which nucleosomes form in regions flanking regulatory sites before induction, and the extent of their removal upon induction. These nucleosomes, evidently, inhibit basal transcription but are poised to be removed quickly upon command.

How to make stripes: deciphering the transition from non-periodic to periodic patterns in Drosophila segmentation

The generation of metameric body plans is a key process in development. In Drosophila segmentation, periodicity is established rapidly through the complex transcriptional regulation of the pair-rule genes. The 'primary' pair-rule genes generate their 7-stripe expression through stripe-specific cis-regulatory elements controlled by the preceding non-periodic maternal and gap gene patterns, whereas 'secondary' pair-rule genes are thought to rely on 7-stripe elements that read off the already periodic primary pair-rule patterns. Using a combination of computational and experimental approaches, we have conducted a comprehensive systems-level examination of the regulatory architecture underlying pair-rule stripe formation. We find that runt (run), fushi tarazu (ftz) and odd skipped (odd) establish most of their pattern through stripe-specific elements, arguing for a reclassification of ftz and odd as primary pair-rule genes. In the case of run, we observe long-range cis-regulation across multiple intervening genes. The 7-stripe elements of run, ftz and odd are active concurrently with the stripe-specific elements, indicating that maternal/gap-mediated control and pair-rule gene cross-regulation are closely integrated. Stripe-specific elements fall into three distinct classes based on their principal repressive gap factor input; stripe positions along the gap gradients correlate with the strength of predicted input. The prevalence of cis-elements that generate two stripes and their genomic organization suggest that single-stripe elements arose by splitting and subfunctionalization of ancestral dual-stripe elements. Overall, our study provides a greatly improved understanding of how periodic patterns are established in the Drosophila embryo.

Gene length may contribute to graded transcriptional responses in the Drosophila embryo

An important question in developmental biology is how relatively shallow gradients of morphogens can reliably establish a series of distinct transcriptional readouts. Current models emphasize interactions between transcription factors binding in distinct modes to cis-acting sequences of target genes. Another recent idea is that the cis-acting interactions may amplify preexisting biases or prepatterns to establish robust transcriptional responses. In this study, we examine the possible contribution of one such source of prepattern, namely gene length. We developed quantitative imaging tools to measure gene expression levels for several loci at a time on a single-cell basis and applied these quantitative imaging tools to dissect the establishment of a gene expression border separating the mesoderm and neuroectoderm in the early Drosophila embryo. We first characterized the formation of a transient ventral-to-dorsal gradient of the Snail (Sna) repressor and then examined the relationship between this gradient and repression of neural target genes in the mesoderm. We found that neural genes are repressed in a nested pattern within a zone of the mesoderm abutting the neuroectoderm, where Sna levels are graded. While several factors may contribute to the transient graded response to the Sna gradient, our analysis suggests that gene length may play an important, albeit transient, role in establishing these distinct transcriptional responses. One prediction of the gene-length-dependent transcriptional patterning model is that the co-regulated genes knirps (a short gene) and knirps-related (a long gene) should be transiently expressed in domains of differing widths, which we confirmed experimentally. These findings suggest that gene length may contribute to establishing graded responses to morphogen gradients by providing transient prepatterns that are subsequently amplified and stabilized by traditional cis-regulatory interactions.

The gap gene network

Gap genes are involved in segment determination during the early development of the fruit fly Drosophila melanogaster as well as in other insects. This review attempts to synthesize the current knowledge of the gap gene network through a comprehensive survey of the experimental literature. I focus on genetic and molecular evidence, which provides us with an almost-complete picture of the regulatory interactions responsible for trunk gap gene expression. I discuss the regulatory mechanisms involved, and highlight the remaining ambiguities and gaps in the evidence. This is followed by a brief discussion of molecular regulatory mechanisms for transcriptional regulation, as well as precision and size-regulation provided by the system. Finally, I discuss evidence on the evolution of gap gene expression from species other than Drosophila. My survey concludes that studies of the gap gene system continue to reveal interesting and important new insights into the role of gene regulatory networks in development and evolution.

Deciphering a transcriptional regulatory code: modeling short-range repression in the Drosophila embryo

A well-defined set of transcriptional regulatory modules was created and analyzed in the Drosophila embryo.

Fractional occupancy-based models were developed to explain the interaction of short range transcriptional repressors with endogenous activators by using quantitative data from these modules.

Our fractional occupancy-based modeling uncovered specific quantitative features of short-range repressors; a complex nonlinear quenching relationship, similar quenching efficiencies for different activators, and modest levels of cooperativity

The extension of the study to endogenous enhancers highlighted several features of enhancer architecture design in Drosophila embryos.

Transcriptional regulatory information, represented by patterns of protein-binding sites on DNA, comprises an important portion of genetic coding. Despite the abundance of genomic sequences now available, identifying and characterizing this information remain a major challenge. Minor changes in protein-binding sites can have profound effects on gene expression, and such changes have been shown to underlie important aspects of disease and evolution. Thus, an important aim in contemporary systems biology is to develop a global understanding of the transcriptional regulatory code, allowing prediction of gene output based on DNA sequence information. Recent studies have focused on endogenous transcriptional regulatory sequences (Janssens et al, 2006; Zinzen et al, 2006; Segal et al, 2008); however, distinct enhancers differ in many features, including transcription factor activity, spacing, and cooperativity, making it difficult to learn the effects of individual features and generalize them to other cis-regulatory elements. We have pursued a bottom up approach to understand the mechanistic processing of regulatory elements by the transcriptional machinery, using a well-defined and characterized set of repressors and activators in Drosophila blastoderm embryos. The study focuses on the Giant, Krüppel, Knirps, and Snail proteins, which have been characterized as short-range repressors, able to act locally to interfere with activator function (quenching) (Gray et al, 1994; Arnosti et al, 1996a). Such repressors have central functions in development.

The aim our study was to enable ab initio predictions of enhancer function, given defined quantities of regulatory proteins and the sequence of the enhancer (Figure 1). We have generated a large quantitative data set using fluorescent confocal laser scanning microscopy to determine the inputs (Giant, Krüppel, and Knirps protein levels) and outputs (lacZ mRNA levels) of the regulatory elements introduced into Drosophila by transgenesis. We analyzed the effect of altering specific features of a set of related gene modules, designed to uncover critical aspects of repression, including quenching distance, cooperativity, and overall factor potency.

We generated specific descriptions for each regulatory element using fractional occupancy-based modeling and identified quantitative values for parameters affecting transcriptional regulation in vivo, and these parameters were used to build and test the model. Through this process, we uncovered earlier unknown features that allow correct predictions of regulation by short-range repressors, including a non-monotonic distance function for quenching, which implicates possible phasing effects, a modest contribution for repressor–repressor cooperativity, and similarity in repression of disparate activators.

By applying these parameters to a model of the endogenous rhomboid enhancer, we uncovered novel insights into the architecture of this enhancer (Figure 8). Our study provides essential quantitative elements of a transcriptional regulatory code that will allow extensive analysis of genomic information in Drosophila melanogaster and related organisms. Extension of these predictive models should facilitate the development of more sophisticated computational algorithms for the identification and functional characterization of novel regulatory elements. The development of such quantitative modeling tools will change our understanding of the genome from essentially a parts list to a dynamically regulated system, and will greatly facilitate studies in disease, population genetics, and evolutionary biology.

Systems biology seeks a genomic-level interpretation of transcriptional regulatory information represented by patterns of protein-binding sites. Obtaining this information without direct experimentation is challenging; minor alterations in binding sites can have profound effects on gene expression, and underlie important aspects of disease and evolution. Quantitative modeling offers an alternative path to develop a global understanding of the transcriptional regulatory code. Recent studies have focused on endogenous regulatory sequences; however, distinct enhancers differ in many features, making it difficult to generalize to other cis-regulatory elements. We applied a systematic approach to simpler elements and present here the first quantitative analysis of short-range transcriptional repressors, which have central functions in metazoan development. Our fractional occupancy-based modeling uncovered unexpected features of these proteins' activity that allow accurate predictions of regulation by the Giant, Knirps, Krüppel, and Snail repressors, including modeling of an endogenous enhancer. This study provides essential elements of a transcriptional regulatory code that will allow extensive analysis of genomic information in Drosophila melanogaster and related organisms.

Groucho corepressor functions as a cofactor for the Knirps short-range transcriptional repressor

Despite the pervasive roles for repressors in transcriptional control, the range of action of these proteins on cis regulatory elements remains poorly understood. Knirps has essential roles in patterning the Drosophila embryo by means of short-range repression, an activity that is essential for proper regulation of complex transcriptional control elements. Short-range repressors function in a local fashion to interfere with the activity of activators or basal promoters within approximately 100 bp. In contrast, long-range repressors such as Hairy act over distances >1 kb. The functional distinction between these two classes of repressors has been suggested to stem from the differential recruitment of the CtBP corepressor to short-range repressors and Groucho to long-range repressors. Contrary to this differential recruitment model, we report that Groucho is a functional part of the Knirps short-range repression complex. The corepressor interaction is mediated via an eh-1 like motif present in the N terminus and a conserved region present in the central portion of Knirps. We also show that this interaction is important for the CtBP-independent repression activity of Knirps and is required for regulation of even-skipped. Our study uncovers a previously uncharacterized interaction between proteins previously thought to function in distinct repression pathways, and indicates that the Groucho corepressor can be differentially harnessed to execute short- and long-range repression.

Image processing and analysis for quantifying gene expression from early Drosophila embryos

Correlation of quantities of transcriptional activators and repressors with the mRNA output of target genes is a central issue for modeling gene regulation. In multicellular organisms, both spatial and temporal differences in gene expression must be taken into account; this can be achieved by use of in situ hybridization followed by confocal laser scanning microscopy (CLSM). Here we present a method to correlate the protein levels of the short-range repressor Giant with lacZ mRNA produced by reporter genes using images of Drosophila blastoderm embryos taken by CLSM. The image stacks from CLSM are processed using a semiautomatic algorithm to produce correlations between the repressor levels and lacZ mRNA reporter genes. We show that signals derived from CLSM are proportional to actual mRNA levels. Our analysis reveals that a suggested parabolic form of the background fluorescence in confocal images of early Drosophila embryos is evident most prominently in flattened specimens, with intact embryos exhibiting a more linear background. The data extraction described in this paper is primarily conceived for analysis of synthetic reporter genes that are designed to decipher cis-regulatory grammar, but the techniques are generalizable for quantitative analysis of other engineered or endogenous genes in embryos.

Gene regulation: boundaries within limits

Quantitative measurements of the Hunchback transcription factor in Drosophila embryos show that its target genes can respond with a high degree of precision to the exact level of the protein, simulating a continuous, analog readout, while other target genes show a combinatorial effect, resembling a Boolean logic element.

The Drosophila gap gene giant has an anterior segment identity function mediated through disconnected and teashirt

The C2H2 zinc-finger-containing transcription factors encoded by the disconnected (disco) and teashirt (tsh) genes contribute to the regionalization of the Drosophila embryo by establishing fields in which specific Homeotic complex (Hom-C) proteins can function. In Drosophila embryos, disco and the paralogous disco-related (disco-r) are expressed throughout most of the epidermis of the head segments, but only in small patches in the trunk segments. Conversely, tsh is expressed extensively in the trunk segments, with little or no accumulation in the head segments. Little is known about the regulation of these genes; for example, what limits their expression to these domains? Here, we report the regulatory effects of gap genes on the spatial expression of disco, disco-r, and tsh during Drosophila embryogenesis. The data shed new light on how mutations in giant (gt) affect patterning within the anterior gt domain, demonstrating homeotic function in this domain. However, the homeosis does not occur through altered expression of the Hom-C genes but through changes in the regulation of disco and tsh.

Spreading of a corepressor linked to action of long-range repressor hairy

Transcriptional repressor proteins play key roles in the control of gene expression in development. For the Drosophila embryo, the following two functional classes of repressors have been described: short-range repressors such as Knirps that locally inhibit the activity of enhancers and long-range repressors such as Hairy that can dominantly inhibit distal elements. Several long-range repressors interact with Groucho, a conserved corepressor that is homologous to mammalian TLE proteins. Groucho interacts with histone deacetylases and histone proteins, suggesting that it may effect repression by means of chromatin modification; however, it is not known how long-range effects are mediated. Using embryo chromatin immunoprecipitation, we have analyzed a Hairy-repressible gene in the embryo during activation and repression. When inactivated, repressors, activators, and coactivators cooccupy the promoter, suggesting that repression is not accomplished by the displacement of activators or coactivators. Strikingly, the Groucho corepressor is found to be recruited to the transcribed region of the gene, contacting a region of several kilobases, concomitant with a loss of histone H3 and H4 acetylation. Groucho has been shown to form higher-order complexes in vitro; thus, our observations suggest that long-range effects may be mediated by a "spreading" mechanism, modifying chromatin over extensive regions to inhibit transcription.

Permissive and instructive anterior patterning rely on mRNA localization in the wasp embryo

The long-germ mode of embryogenesis, in which segments arise simultaneously along the anteriorposterior axis, has evolved several times in different lineages of the holometabolous, or fully metamorphosing, insects. Drosophila's long-germ fate map is established largely by the activity of the dipteran-specific Bicoid (Bcd) morphogen gradient, which operates both instructively and permissively to accomplish anterior patterning. By contrast, all nondipteran long-germ insects must achieve anterior patterning independently of bcd. We show that bcd's permissive function is mimicked in the wasp by a maternal repression system in which anterior localization of the wasp ortholog of giant represses anterior expression of the trunk gap genes so that head and thorax can properly form.

Functional interaction between the Drosophila knirps short range transcriptional repressor and RPD3 histone deacetylase

Knirps and other short range transcriptional repressors play critical roles in patterning the early Drosophila embryo. These repressors are known to bind the C-terminal binding protein corepressor, but their mechanism of action is poorly understood. We purified functional recombinant Knirps protein from transgenic embryos to identify possible cofactors that contribute to the activity of this protein. The protein migrates in a complex of approximately 450 kDa and was found to copurify with the Rpd3 histone deacetylase protein during a double affinity purification procedure. Association of Rpd3 with Knirps was dependent on the presence of the C-terminal binding protein-dependent repression domain of Knirps. Previous studies of an rpd3 mutant had not shown defects in the pattern of expression of even-skipped, a target of the Knirps repressor. However, in embryos doubly heterozygous for knirps and rpd3, a marked increase in the frequency of defects in the Knirps-regulated posterior domain of even-skipped expression was found, indicating that Rpd3 contributes to Knirps repression activity in vivo. This finding implicates deacetylation in the mechanism of short range repression in Drosophila.

cis-regulatory logic of short-range transcriptional repression in Drosophila melanogaster

Bioinformatics analysis of transcriptional control is guided by knowledge of the characteristics of cis-regulatory regions or enhancers. Features such as clustering of binding sites and co-occurrence of binding sites have aided enhancer identification, but quantitative predictions of enhancer function are not yet generally feasible. To facilitate the analysis of regulatory sequences in Drosophila melanogaster, we identified quantitative parameters that affect the activity of short-range transcriptional repressors, proteins that play key roles in development. In addition to the previously noted distance dependence, repression is strongly influenced by the stoichiometry, affinity, spacing, and arrangement of activator binding sites. Repression is insensitive to the type of activation domain, suggesting that short-range repression may primarily affect activators at the level of DNA binding. The activity of several short-range, but not long-range, repressors is circumscribed by the same quantitative parameters. This cis-regulatory "grammar" may aid the identification of enhancers regulated by short-range repressors and facilitate bioinformatic prediction of the functional output of transcriptional regulatory sequences.

A high-throughput method for quantifying gene expression data from early Drosophila embryos

We describe an automated high-throughput method to measure protein levels in single nuclei in blastoderm embryos of Drosophila melanogaster by means of immunofluorescence. The method consists of a chain of specific algorithms assembled into an image processing pipeline. This pipeline transforms a confocal scan of an embryo stained with fluorescently tagged antibodies into a text file. This text file contains a numerical identifier for each nucleus, the coordinates of its centroid, and the average concentrations of three proteins in that nucleus. The central algorithmic component of the method is the automatic identification of nuclei by edge detection with the use of watersheds as an error-correction step. This method provides high-throughput quantification at cellular resolution.

Transcriptional enhancers: Intelligent enhanceosomes or flexible billboards?

In higher eukaryotes, transcriptional enhancers play critical roles in the integration of cellular signaling information, but apart from a few well-studied model enhancers, we lack a general picture of transcriptional information processing by most enhancers. Here we discuss recent studies that have provided fresh insights on information processing that occurs on enhancers, and propose that in addition to the highly cooperative and coordinate action of "enhanceosomes", a less integrative, but more flexible form of information processing is mediated by information display or "billboard" enhancers. Application of these models has important ramifications not only for the biochemical analysis of transcription, but also for the wider fields of bioinformatics and evolutionary biology.

CtBP contributes quantitatively to Knirps repression activity in an NAD binding-dependent manner

Transcriptional repressors often employ multiple activities, but the molecular mechanisms and physiological relevance of this functional diversity remain obscure. The Drosophila melanogaster Knirps repressor uses CtBP corepressor-dependent and -independent pathways. To separately analyze the components of Knirps repression activity, we elucidated the specific repression properties of CtBP and of Knirps subdomains. Like Knirps, CtBP represses adjacent transcriptional activators; but unlike Knirps, CtBP is unable to repress basal promoter elements. We determined that the ability of CtBP to recapitulate only a subset of Knirps activities is due to a quantitative, rather than qualitative, deficiency in repression activity. The CtBP-dependent portion of Knirps synergizes with the CtBP-independent repression activity to potently repress promoter elements from enhancer- or promoter-proximal positions. This result indicates that multiple repression activities are combined to exceed critical thresholds on target genes. CtBP mutant proteins unable to bind NAD fail to interact with DNA-bound factors. We show that DNA-binding Gal4-CtBP fusion proteins also require NAD binding for activity, indicating that NAD plays a role in repression at a step subsequent to CtBP recruitment to the promoter.

Insect basic leucine zipper proteins and their role in cyclic AMP-dependent regulation of gene expression

The cAMP-protein kinase A (PKA) pathway is an important intracellular signal transduction cascade that can be activated by a large variety of stimuli. Activation or inhibition of this pathway will ultimately affect the transcriptional regulation of various genes through distinct responsive sites. In vertebrates, the best- characterized nuclear targets of PKA are the cyclic AMP response element-binding (CREB) proteins. It is now well established that CREB is not only regulated by PKA, but many other kinases can exert an effect as well. Since CREB-like proteins were also discovered in invertebrates, several studies unraveling their physiological functions in this category of metazoans have been performed. This review will mainly focus on the presence and regulation of CREB proteins in insects. Differences in transcriptional responses to the PKA pathway and other CREB-regulating stimuli between cells, tissues, and even organisms can be partially attributed to the presence of different CREB isoforms. In addition, the regulation of CREB appears to show some important differences between insects and vertebrates. Since CREB is a basic leucine zipper (bZip) protein, other insect members of this important family of transcriptional regulators will be briefly discussed as well.

Design and function of transcriptional switches in Drosophila

Extensive genetic and biochemical analysis of Drosophila melanogaster has made this system an important model for characterization of transcriptional regulatory elements and factors. Given the striking conservation of transcriptional controls in metazoans, general principles derived from studies of Drosophila are expected to continue to illuminate transcriptional regulation in other systems, including vertebrates. With improvement in technologies for genetic manipulation of insects, research in Drosophila will also aid the design of systems for controlled expression of genes in other hosts. This review focuses on recent advances from Drosophila in analysis of the functional components of transcriptional switches, including basal promoters, enhancers, boundary elements, and maintenance elements.

Role of CtBP in transcriptional repression by the Drosophila giant protein

The giant protein is a short-range transcriptional repressor that refines the expression pattern of gap and pair-rule genes in the Drosophila blastoderm embryo. Short-range repressors including knirps, Krüppel, and snail utilize the CtBP cofactor for repression, but it is not known whether a functional interaction with CtBP is a general property of all short-range repressors. We studied giant repression activity in a CtBP mutant and find that this cofactor is required for giant repression of some, but not all, genes. While targets of giant such as the even-skipped stripe 2 enhancer and a synthetic lacZ reporter show clear derepression in the CtBP mutant, another giant target, the hunchback gene, is expressed normally. A more complex situation is seen with regulation of the Krüppel gene, in which one enhancer is repressed by giant in a CtBP-dependent manner, while another is repressed in a CtBP-independent manner. These results demonstrate that giant can repress both via CtBP-dependent and CtBP-independent pathways, and that promoter context is critical for determining giant-CtBP functional interaction. To initiate mechanistic studies of the giant repression activity, we have identified a minimal repression domain within giant that encompasses residues 89-205, including an evolutionarily conserved region bearing a putative CtBP binding motif.

Cell-type specificity of short-range transcriptional repressors

Transcriptional repressors can be classified as short- or long-range, according to their range of activity. Functional analysis of identified short-range repressors has been carried out largely in transgenic Drosophila, but it is not known whether general properties of short-range repressors are evident in other types of assays. To study short-range transcriptional repressors in cultured cells, we created chimeric tetracycline repressors based on Drosophila transcriptional repressors Giant, Drosophila C-terminal-binding protein (dCtBP), and Knirps. We find that Giant and dCtBP are efficient repressors in Drosophila and mammalian cells, whereas Knirps is active only in insect cells. The restricted activity of Knirps, in contrast to that of Giant, suggests that not all short-range repressors possess identical activities, consistent with recent findings showing that short-range repressors act through multiple pathways. The mammalian repressor Kid is more effective than either Giant or dCtBP in mammalian cells but is inactive in Drosophila cells. These results indicate that species-specific factors are important for the function of the Knirps and Kid repressors. Giant and dCtBP repress reporter genes in a variety of contexts, including genes that were introduced by transient transfection, carried on episomal elements, or stably integrated. This broad activity indicates that the context of the target gene is not critical for the ability of short-range repressors to block transcription, in contrast to other repressors that act only on stably integrated genes.

Analysis and function of transcriptional regulatory elements: insights from Drosophila

Analysis of gene expression is assuming an increasingly important role in elucidating the molecular basis of insect biology. Transcriptional regulation of gene expression is directed by a variety of cis-acting DNA elements that control spatial and temporal patterns of expression. This review summarizes current knowledge about properties of transcriptional regulatory elements, based largely on research in Drosophila melanogaster, and outlines ways that new technologies are providing tools to facilitate the study of transcriptional regulatory elements in other insects.

Thoracic patterning by the Drosophila gap gene hunchback

Localized gene expression patterns are critical for establishing body plans in all multicellular animals. In Drosophila, the gap gene hunchback (hb) is expressed in a dynamic pattern in anterior regions of the embryo. Hb protein is first detected as a shallow maternal gradient that prevents expression of posterior gap genes in anterior regions. hb mRNA is also expressed zygotically, first as a broad anterior domain controlled by the Bicoid (Bcd) morphogen, and then in a stripe at the position of parasegment 4 (PS4). Here, we show that the PS4-hb stripe changes the profile of the anterior Hb gradient by generating a localized peak of protein that persists until after the broad domain has started to decline. This peak is required specifically for the formation of the mesothoracic (T2) segment. At the molecular level, the PS4-hb stripe is critical for activation of the homeotic gene Antennapedia, but does not affect a gradient of Hb repressive activity formed by the combination of maternal and Bcd-dependent Hb. The repressive gradient is critical for establishing the positions of several target genes, including the gap genes Kruppel (Kr), knirps (kni), and giant (gt), and the homeotic gene Ultrabithorax (Ubx). Different Hb concentrations are sufficient for repression of gt, kni, and Ubx, but a very high level of Hb, or a combinatorial mechanism, is required for repression of Kr. These results suggest that the individual phases of hb transcription, which overlap temporally and spatially, contribute specific patterning functions in early embryogenesis.

CtBP-dependent activities of the short-range Giant repressor in the Drosophila embryo

There are at least three short-range gap repressors in the precellular Drosophila embryo: Krüppel, Knirps, and Giant. Krüppel and Knirps contain related repression motifs, PxDLSxH and PxDLSxK, respectively, which mediate interactions with the dCtBP corepressor protein. Here, we present evidence that Giant might also interact with dCtBP. The misexpression of Giant in ventral regions of transgenic embryos results in the selective repression of eve stripe 5. A stripe5-lacZ transgene exhibits an abnormal staining pattern in dCtBP mutants that is consistent with attenuated repression by Giant. The analysis of Gal4-Giant fusion proteins identified a minimal repression domain that contains a sequence motif, VLDLS, which is conserved in at least two other sequence-specific repressors. Removal of this sequence from the native Giant protein does not impair its repression activity in transgenic embryos. We propose that Giant-dCtBP interactions might be indirect and mediated by an unknown bZIP subunit that forms a heteromeric complex with Giant. We also suggest that the VLDLS motif recruits an as yet unidentified corepressor protein.

Complex organization of promoter and enhancer elements regulate the tissue- and developmental stage-specific expression of the Drosophila melanogaster Gld gene

The Drosophila melanogaster Gld gene has multiple and diverse developmental and physiological functions. We report herein that interactions among proximal promoter elements and a cluster of intronically located enhancers and silencers specify the complex regulation of Gld that underlies its diverse functions. Gld expression in nonreproductive tissues is largely determined by proximal promoter elements with the exception of the embryonic labium where Gld is activated by an enhancer within the first intron. A nuclear protein, GPAL, has been identified that binds the Gpal elements in the proximal promoter region. Regulation of Gld in the reproductive organs is particularly complex, involving interactions among the Gpal proximal promoter elements, a unique TATA box, three distinct enhancer types, and one or more silencer elements. The three somatic reproductive organ enhancers each activate expression in male and female pairs of reproductive organs. One of these pairs, the male ejaculatory duct and female oviduct, are known to be developmentally homologous. We report evidence that the other two pairs of organs are developmentally homologous as well. A comprehensive model to explain the full developmental regulation of Gld and its evolution is presented.

dCtBP-dependent and -independent repression activities of the Drosophila Knirps protein

Transcriptional repressor proteins play essential roles in controlling the correct temporal and spatial patterns of gene expression in Drosophila melanogaster embryogenesis. Repressors such as Knirps, Krüppel, and Snail mediate short-range repression and interact with the dCtBP corepressor. The mechanism by which short-range repressors block transcription is not well understood; therefore, we have undertaken a detailed structure-function analysis of the Knirps protein. To provide a physiological setting for measurement of repression, the activities of endogenous or chimeric Knirps repressor proteins were assayed on integrated reporter genes in transgenic embryos. Two distinct repression functions were identified in Knirps. One repression activity depends on dCtBP binding, and this function maps to a C-terminal region of Knirps that contains a dCtBP binding motif. In addition, an N-terminal region was identified that represses in a CtBP mutant background and does not bind to the dCtBP protein in vitro. Although the dCtBP protein is important for Knirps activity on some genes, one endogenous target of the Knirps protein, the even-skipped stripe 3 enhancer, is not derepressed in a CtBP mutant. These results indicate that Knirps can utilize two different pathways to mediate transcriptional repression and suggest that the phenomenon of short-range repression may be a combination of independent activities.

Evidence for stabilizing selection in a eukaryotic enhancer element

Eukaryotic gene expression is mediated by compact cis-regulatory modules, or enhancers, which are bound by specific sets of transcription factors. 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 enhancers. 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 wildtype 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.

Discrete enhancer elements mediate selective responsiveness of enhancer of split complex genes to common transcriptional activators

In Drosophila, genes of the Enhancer of split Complex [E(spl)-C] are important components of the Notch (N) cell-cell signaling pathway, which is utilized in imaginal discs to effect a series of cell fate decisions during adult peripheral nervous system development. Seven genes in the complex encode basic helix-loop-helix (bHLH) transcriptional repressors, while 4 others encode members of the Bearded family of small proteins. A striking diversity is observed in the imaginal disc expression patterns of the various E(spl)-C genes, suggestive of a diversity of function, but the mechanistic basis of this variety has not been elucidated. Here we present strong evidence from promoter-reporter transgene experiments that regulation at the transcriptional level is primarily responsible. Certain E(spl)-C genes were known previously to be direct targets of transcriptional activation both by the N-signal-dependent activator Suppressor of Hairless [Su(H)] and by the proneural bHLH proteins achaete and scute. Our extensive sequence analysis of the promoter-proximal upstream regions of 12 transcription units in the E(spl)-C reveals that such dual transcriptional activation is likely to be the rule for at least 10 of the 12 genes. We next show that the very different wing imaginal disc expression patterns of E(spl)m4 and E(spl)mgamma are a property of small (200-300 bp), evolutionarily conserved transcriptional enhancer elements, which can confer these distinct patterns on a heterologous promoter despite their considerable structural similarity [each having three Su(H) and two proneural protein binding sites]. We also demonstrate that the characteristic inactivity of the E(spl)mgamma enhancer in the notum and margin territories of the wing disc can be overcome by elevated activity of the N receptor. We conclude that the distinctive expression patterns of E(spl)-C genes in imaginal tissues depend to a significant degree on the capacity of their transcriptional cis-regulatory apparatus to respond selectively to direct proneural- and Su(H)-mediated activation, often in only a subset of the territories and cells in which these modes of regulation are operative.