Cooperative DNA-binding by Bicoid provides a mechanism for threshold-dependent gene activation in the Drosophila embryo

The Bicoid morphogen directs pattern formation along the anterior-posterior (A-P) axis of the Drosophila embryo. Bicoid is distributed in a concentration gradient that decreases exponentially from the anterior pole, however, it transcribes target genes such as hunchback in a step-function-like pattern; the expression domain is uniform and has a sharply defined posterior boundary. A 'gradient-affinity' model proposed to explain Bicoid action states that (i) cooperative gene activation by Bicoid generates the sharp on/off switch for target gene transcription and (ii) target genes with different affinities for Bicoid are expressed at different positions along the A-P axis. Using an in vivo yeast assay and in vitro methods, we show that Bicoid binds DNA with pairwise cooperativity; Bicoid bound to a strong site helps Bicoid bind to a weak site. These results support the first aspect of the model, providing a mechanism by which Bicoid generates sharp boundaries of gene expression. However, contrary to the second aspect of the model, we find no significant difference between the affinity of Bicoid for the anterior gene hunchback and the posterior gene knirps, We propose, instead, that the arrangement of Bicoids bound to the target gene presents a unique signature to the transcription machinery that, in combination with overall affinity, regulates the extent of gene transcription along the A-P axis

Expression and regulation of caudal in the lower cyclorrhaphan fly Megaselia

The homeobox gene caudal (cad) regulates posterior development in Drosophila. In early embryos, the cad protein (CAD) is expressed in a posterior-to-anterior concentration gradient, which contributes polarity to the developing embryo. The CAD gradient is complementary to and dependent on the anterior pattern organizer Bicoid (BCD), which represses the translation of ubiquitous maternal cad transcripts in the anterior embryo through a direct interaction with the cad 3′ untranslated region (UTR). Here, we show that early embryos of the lower cyclorrhaphan fly Megaselia express the putative cad orthologue Mab-cad throughout the posterior three quarters of the blastoderm but lack maternal transcripts. In transgenic blastoderm embryos of Drosophila, Mab-cad cis-regulatory DNA drives the expression of a reporter gene in a similar pattern, while Mab-cad 3′ UTR fails to mediate translational repression of a ubiquitously transcribed reporter. For another lower cyclorrhaphan fly (Lonchoptera) and two related outgroup taxa of Cyclorrhapha (Empis, Haematopota), we report maternal cad expression in ovarian follicles. Together, our results suggest that BCD is not required for the translational repression of Mab-cad, and that maternal cad expression was lost in the Megaselia lineage

The Drosophila Gap Gene Network Is Composed of Two Parallel Toggle Switches

Drosophila “gap” genes provide the first response to maternal gradients in the early fly embryo. Gap genes are expressed in a series of broad bands across the embryo during first hours of development. The gene network controlling the gap gene expression patterns includes inputs from maternal gradients and mutual repression between the gap genes themselves. In this study we propose a modular design for the gap gene network, involving two relatively independent network domains. The core of each network domain includes a toggle switch corresponding to a pair of mutually repressive gap genes, operated in space by maternal inputs. The toggle switches present in the gap network are evocative of the phage lambda switch, but they are operated positionally (in space) by the maternal gradients, so the synthesis rates for the competing components change along the embryo anterior-posterior axis. Dynamic model, constructed based on the proposed principle, with elements of fractional site occupancy, required 5–7 parameters to fit quantitative spatial expression data for gap gradients. The identified model solutions (parameter combinations) reproduced major dynamic features of the gap gradient system and explained gap expression in a variety of segmentation mutants

Dual-functioning transcription factors in the developmental gene network of Drosophila melanogaster

Quantitative models for transcriptional regulation have shown great promise for advancing our understanding of the biological mechanisms underlying gene regulation. However, all of the models to date assume a transcription factor (TF) to have either activating or repressing function towards all the genes it is regulating.In this paper we demonstrate, on the example of the developmental gene network in D. melanogaster, that the data-fit can be improved by up to 40% if the model is allowing certain TFs to have dual function, that is, acting as activator for some genes and as repressor for others. We demonstrate that the improvement is not due to additional flexibility in the model but rather derived from the data itself. We also found no evidence for the involvement of other known site-specific TFs in regulating this network. Finally, we propose SUMOylation as a candidate biological mechanism allowing TFs to switch their role when a small ubiquitin-like modifier (SUMO) is covalently attached to the TF. We strengthen this hypothesis by demonstrating that the TFs predicted to have dual function also contain the known SUMO consensus motif, while TFs predicted to have only one role lack this motif.We argue that a SUMOylation-dependent mechanism allowing TFs to have dual function represents a promising area for further research and might be another step towards uncovering the biological mechanisms underlying transcriptional regulation

Evolutionary Mirages: Selection on Binding Site Composition Creates the Illusion of Conserved Grammars in Drosophila Enhancers

The clustering of transcription factor binding sites in developmental enhancers and the apparent preferential conservation of clustered sites have been widely interpreted as proof that spatially constrained physical interactions between transcription factors are required for regulatory function. However, we show here that selection on the composition of enhancers alone, and not their internal structure, leads to the accumulation of clustered sites with evolutionary dynamics that suggest they are preferentially conserved. We simulated the evolution of idealized enhancers from Drosophila melanogaster constrained to contain only a minimum number of binding sites for one or more factors. Under this constraint, mutations that destroy an existing binding site are tolerated only if a compensating site has emerged elsewhere in the enhancer. Overlapping sites, such as those frequently observed for the activator Bicoid and repressor Krüppel, had significantly longer evolutionary half-lives than isolated sites for the same factors. This leads to a substantially higher density of overlapping sites than expected by chance and the appearance that such sites are preferentially conserved. Because D. melanogaster (like many other species) has a bias for deletions over insertions, sites tended to become closer together over time, leading to an overall clustering of sites in the absence of any selection for clustered sites. Since this effect is strongest for the oldest sites, clustered sites also incorrectly appear to be preferentially conserved. Following speciation, sites tend to be closer together in all descendent species than in their common ancestors, violating the common assumption that shared features of species' genomes reflect their ancestral state. Finally, we show that selection on binding site composition alone recapitulates the observed number of overlapping and closely neighboring sites in real D. melanogaster enhancers. Thus, this study calls into question the common practice of inferring “cis-regulatory grammars” from the organization and evolutionary dynamics of developmental enhancers

Finding regulatory elements and regulatory motifs: a general probabilistic framework

Over the last two decades a large number of algorithms has been developed for regulatory motif finding. Here we show how many of these algorithms, especially those that model binding specificities of regulatory factors with position specific weight matrices (WMs), naturally arise within a general Bayesian probabilistic framework. We discuss how WMs are constructed from sets of regulatory sites, how sites for a given WM can be discovered by scanning of large sequences, how to cluster WMs, and more generally how to cluster large sets of sites from different WMs into clusters. We discuss how 'regulatory modules', clusters of sites for subsets of WMs, can be found in large intergenic sequences, and we discuss different methods for ab initio motif finding, including expectation maximization (EM) algorithms, and motif sampling algorithms. Finally, we extensively discuss how module finding methods and ab initio motif finding methods can be extended to take phylogenetic relations between the input sequences into account, i.e. we show how motif finding and phylogenetic footprinting can be integrated in a rigorous probabilistic framework. The article is intended for readers with a solid background in applied mathematics, and preferably with some knowledge of general Bayesian probabilistic methods. The main purpose of the article is to elucidate that all these methods are not a disconnected set of individual algorithmic recipes, but that they are just different facets of a single integrated probabilistic theory

Drosophila TIEG Is a Modulator of Different Signalling Pathways Involved in Wing Patterning and Cell Proliferation

Acquisition of a final shape and size during organ development requires a regulated program of growth and patterning controlled by a complex genetic network of signalling molecules that must be coordinated to provide positional information to each cell within the corresponding organ or tissue. The mechanism by which all these signals are coordinated to yield a final response is not well understood. Here, I have characterized the Drosophila ortholog of the human TGF-β Inducible Early Gene 1 (dTIEG). TIEG are zinc-finger proteins that belong to the Krüppel-like factor (KLF) family and were initially identified in human osteoblasts and pancreatic tumor cells for the ability to enhance TGF-β response. Using the developing wing of Drosophila as “in vivo” model, the dTIEG function has been studied in the control of cell proliferation and patterning. These results show that dTIEG can modulate Dpp signalling. Furthermore, dTIEG also regulates the activity of JAK/STAT pathway suggesting a conserved role of TIEG proteins as positive regulators of TGF-β signalling and as mediators of the crosstalk between signalling pathways acting in a same cellular context

The Formation of the Bicoid Morphogen Gradient Requires Protein Movement from Anteriorly Localized mRNA

New quantitative data show that the Bicoid morphogen gradient is generated from a dynamic localized source and that protein gradient formation requires protein movement along the anterior-posterior axis

A Cis-Regulatory Map of the Drosophila Genome

Systematic annotation of gene regulatory elements is a major challenge in genome science. Direct mapping of chromatin modification marks and transcriptional factor binding sites genome-wide1, 2 has successfully identified specific subtypes of regulatory elements3. In Drosophila several pioneering studies have provided genome-wide identification of Polycomb response elements4, chromatin states5, transcription factor binding sites6, 7, 8, 9, RNA polymerase II regulation8 and insulator elements10; however, comprehensive annotation of the regulatory genome remains a significant challenge. Here we describe results from the modENCODE cis-regulatory annotation project. We produced a map of the Drosophila melanogaster regulatory genome on the basis of more than 300 chromatin immunoprecipitation data sets for eight chromatin features, five histone deacetylases and thirty-eight site-specific transcription factors at different stages of development. Using these data we inferred more than 20,000 candidate regulatory elements and validated a subset of predictions for promoters, enhancers and insulators in vivo. We identified also nearly 2,000 genomic regions of dense transcription factor binding associated with chromatin activity and accessibility. We discovered hundreds of new transcription factor co-binding relationships and defined a transcription factor network with over 800 potential regulatory relationships