Two distinct mechanisms for differential positioning of gene expression borders involving the Drosophila gap protein giant

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.

Multi-level regulation of even-skipped stripes by the ubiquitous factor Zelda

The zinc-finger protein Zelda (Zld) is a key activator of zygotic transcription in early Drosophila embryos. Here, we study Zld-dependent regulation of the seven-striped pattern of the pair-rule gene even-skipped (eve). Individual stripes are regulated by discrete enhancers that respond to broadly distributed activators; stripe boundaries are formed by localized repressors encoded by the gap genes. The strongest effects of Zld are on stripes 2, 3 and 7, which are regulated by two enhancers in a 3.8 kb genomic fragment that includes the eve basal promoter. We show that Zld facilitates binding of the activator Bicoid and the gap repressors to this fragment, consistent with its proposed role as a pioneer protein. To test whether the effects of Zld are direct, we mutated all canonical Zld sites in the 3.8 kb fragment, which reduced expression but failed to phenocopy the abolishment of stripes caused by removing Zld in trans. We show that Zld also indirectly regulates the eve stripes by establishing specific gap gene expression boundaries, which provides the embryonic spacing required for proper stripe activation.

Functional coordination between transcription factor clustering and gene activity

The prevailing view of metazoan gene regulation is that transcription is facilitated through the formation of static activator complexes at distal regulatory regions. Here, we employed quantitative single-cell live-imaging and computational analysis to provide evidence that the dynamic assembly and disassembly process of transcription factor (TF) clusters at enhancers is a major source of transcriptional bursting in developing Drosophila embryos. We further show that the regulatory connectivity between TF clustering and burst induction is highly regulated through intrinsically disordered regions (IDRs). Addition of a poly-glutamine tract to the maternal morphogen Bicoid demonstrated that extended IDR length leads to ectopic TF clustering and burst induction from its endogenous target genes, resulting in defects in body segmentation during embryogenesis. Moreover, we successfully visualized the presence of "shared" TF clusters during the co-activation of two distant genes, which provides a concrete molecular explanation for the newly proposed "topological operon" hypothesis in metazoan gene regulation.

An additive repression mechanism sets the anterior limits of anterior pair-rule stripes 1

Segments are repeated anatomical units forming the body of insects. In Drosophila, the specification of the body takes place during the blastoderm through the segmentation cascade. Pair-rule genes such as hairy (h), even-skipped (eve), runt (run), and fushi-tarazu (ftz) are of the intermediate level of the cascade and each pair-rule gene is expressed in seven transversal stripes along the antero-posterior axis of the embryo. Stripes are formed by independent cis-regulatory modules (CRMs) under the regulation of transcription factors of maternal source and of gap proteins of the first level of the cascade. The initial blastoderm of Drosophila is a syncytium and it also coincides with the mid-blastula transition when thousands of zygotic genes are transcribed and their products are able to diffuse in the cytoplasm. Thus, we anticipated a complex regulation of the CRMs of the pair-rule stripes. The CRMs of h 1, eve 1, run 1, ftz 1 are able to be activated by bicoid (bcd) throughout the anterior blastoderm and several lines of evidence indicate that they are repressed by the anterior gap genes slp1 (sloppy-paired 1), tll (tailless) and hkb (huckebein). The modest activity of these repressors led to the premise of a combinatorial mechanism regulating the expression of the CRMs of h 1, eve 1, run 1, ftz 1 in more anterior regions of the embryo. We tested this possibility by progressively removing the repression activities of slp1, tll and hkb. In doing so, we were able to expose a mechanism of additive repression limiting the anterior borders of stripes 1. Stripes 1 respond depending on their distance from the anterior end and repressors operating at different levels.

Patterning with clocks and genetic cascades: Segmentation and regionalization of vertebrate versus insect body plans

Oscillatory and sequential processes have been implicated in the spatial patterning of many embryonic tissues. For example, molecular clocks delimit segmental boundaries in vertebrates and insects and mediate lateral root formation in plants, whereas sequential gene activities are involved in the specification of regional identities of insect neuroblasts, vertebrate neural tube, vertebrate limb, and insect and vertebrate body axes. These processes take place in various tissues and organisms, and, hence, raise the question of what common themes and strategies they share. In this article, we review 2 processes that rely on the spatial regulation of periodic and sequential gene activities: segmentation and regionalization of the anterior-posterior (AP) axis of animal body plans. We study these processes in species that belong to 2 different phyla: vertebrates and insects. By contrasting 2 different processes (segmentation and regionalization) in species that belong to 2 distantly related phyla (arthropods and vertebrates), we elucidate the deep logic of patterning by oscillatory and sequential gene activities. Furthermore, in some of these organisms (e.g., the fruit fly Drosophila), a mode of AP patterning has evolved that seems not to overtly rely on oscillations or sequential gene activities, providing an opportunity to study the evolution of pattern formation mechanisms.

Suboptimal Intermediates Underlie Evolution of the Bicoid Homeodomain

Changes in regulatory networks generate materials for evolution to create phenotypic diversity. For transcription networks, multiple studies have shown that alterations in binding sites of cis-regulatory elements correlate well with the gain or loss of specific features of the body plan. Less is known about alterations in the amino acid sequences of the transcription factors (TFs) that bind these elements. Here we study the evolution of Bicoid (Bcd), a homeodomain (HD) protein that is critical for anterior embryo patterning in Drosophila. The ancestor of Bcd (AncBcd) emerged after a duplication of a Zerknullt (Zen)-like ancestral protein (AncZB) in a suborder of flies. AncBcd diverged from AncZB, gaining novel transcriptional and translational activities. We focus on the evolution of the HD of AncBcd, which binds to DNA and RNA, and is comprised of four subdomains: an N-terminal arm (NT) and three helices; H1, H2, and Recognition Helix (RH). Using chimeras of subdomains and gene rescue assays in Drosophila, we show that robust patterning activity of the Bcd HD (high frequency rescue to adulthood) is achieved only when amino acid substitutions in three separate subdomains (NT, H1, and RH) are combined. Other combinations of subdomains also yield full rescue, but with lower penetrance, suggesting alternative suboptimal activities. Our results suggest a multistep pathway for the evolution of the Bcd HD that involved intermediate HD sequences with suboptimal activities, which constrained and enabled further evolutionary changes. They also demonstrate critical epistatic forces that contribute to the robust function of a DNA-binding domain.

Bicoid-Dependent Activation of the Target Gene hunchback Requires a Two-Motif Sequence Code in a Specific Basal Promoter

In complex genetic loci, individual enhancers interact most often with specific basal promoters. Here we investigate the activation of the Bicoid target gene hunchback (hb), which contains two basal promoters (P1 and P2). Early in embryogenesis, P1 is silent, while P2 is strongly activated. In vivo deletion of P2 does not cause activation of P1, suggesting that P2 contains intrinsic sequence motifs required for activation. We show that a two-motif code (a Zelda binding site plus TATA) is required and sufficient for P2 activation. Zelda sites are present in the promoters of many embryonically expressed genes, but the combination of Zelda plus TATA does not seem to be a general code for early activation or Bicoid-specific activation per se. Because Zelda sites are also found in Bicoid-dependent enhancers, we propose that simultaneous binding to both enhancers and promoters independently synchronizes chromatin accessibility and facilitates correct enhancer-promoter interactions.

Speeding up anterior-posterior patterning of insects by differential initialization of the gap gene cascade

Recently, it was shown that anterior-posterior patterning genes in the red flour beetle Tribolium castaneum are expressed sequentially in waves. However, in the fruit fly Drosophila melanogaster, an insect with a derived mode of embryogenesis compared to Tribolium, anterior-posterior patterning genes quickly and simultaneously arise as mature gene expression domains that, afterwards, undergo slight posterior-to-anterior shifts. This raises the question of how a fast and simultaneous mode of patterning, like that of Drosophila, could have evolved from a rather slow sequential mode of patterning, like that of Tribolium. In this paper, we propose a mechanism for this evolutionary transition based on a switch from a uniform to a gradient-mediated initialization of the gap gene cascade by maternal Hb. The model is supported by computational analyses and experiments.

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.

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.

Quantitative system drift compensates for altered maternal inputs to the gap gene network of the scuttle fly Megaselia abdita

The segmentation gene network in insects can produce equivalent phenotypic outputs despite differences in upstream regulatory inputs between species. We investigate the mechanistic basis of this phenomenon through a systems-level analysis of the gap gene network in the scuttle fly Megaselia abdita (Phoridae). It combines quantification of gene expression at high spatio-temporal resolution with systematic knock-downs by RNA interference (RNAi). Initiation and dynamics of gap gene expression differ markedly between M. abdita and Drosophila melanogaster, while the output of the system converges to equivalent patterns at the end of the blastoderm stage. Although the qualitative structure of the gap gene network is conserved, there are differences in the strength of regulatory interactions between species. We term such network rewiring 'quantitative system drift'. It provides a mechanistic explanation for the developmental hourglass model in the dipteran lineage. Quantitative system drift is likely to be a widespread mechanism for developmental evolution.

Quantitative dynamics and increased variability of segmentation gene expression in the Drosophila Krüppel and knirps mutants

Here we characterize the response of the Drosophila segmentation system to mutations in two gap genes, Kr and kni, in the form of single or double homozygotes and single heterozygotes. Segmentation gene expression in these genotypes was quantitatively monitored with cellular resolution in space and 6.5 to 13min resolution in time. As is the case with wild type, we found that gene expression domains in the posterior portion of the embryo shift to the anterior over time. In certain cases, such as the gt posterior domain in Kr mutants, the shifts are significantly larger than is seen in wild type embryos. We also investigated the effects of Kr and kni on the variability of gene expression. Mutations often produce variable phenotypes, and it is well known that the cuticular phenotype of Kr mutants is variable. We sought to understand the molecular basis of this effect. We find that throughout cycle 14A the relative levels of eve and ftz expression in stripes 2 and 3 are variable among individual embryos. Moreover, in Kr and kni mutants, unlike wild type, the variability in positioning of the posterior Hb domain and eve stripe 7 is not decreased or filtered with time. The posterior Gt domain in Kr mutants is highly variable at early times, but this variability decreases when this domain shifts in the anterior direction to the position of the neighboring Kni domain. In contrast to these findings, positional variability throughout the embryo does not decrease over time in double Kr;kni mutants. In heterozygotes the early expression patterns of segmentation genes resemble patterns seen in homozygous mutants but by the onset of gastrulation they become similar to the wild type patterns. Finally, we note that gene expression levels are reduced in Kr and kni mutant embryos and have a tendency to decrease over time. This is a surprising result in view of the role that mutual repression is thought to play in the gap gene system.

Efficient reverse-engineering of a developmental gene regulatory network

Understanding the complex regulatory networks underlying development and evolution of multi-cellular organisms is a major problem in biology. Computational models can be used as tools to extract the regulatory structure and dynamics of such networks from gene expression data. This approach is called reverse engineering. It has been successfully applied to many gene networks in various biological systems. However, to reconstitute the structure and non-linear dynamics of a developmental gene network in its spatial context remains a considerable challenge. Here, we address this challenge using a case study: the gap gene network involved in segment determination during early development of Drosophila melanogaster. A major problem for reverse-engineering pattern-forming networks is the significant amount of time and effort required to acquire and quantify spatial gene expression data. We have developed a simplified data processing pipeline that considerably increases the throughput of the method, but results in data of reduced accuracy compared to those previously used for gap gene network inference. We demonstrate that we can infer the correct network structure using our reduced data set, and investigate minimal data requirements for successful reverse engineering. Our results show that timing and position of expression domain boundaries are the crucial features for determining regulatory network structure from data, while it is less important to precisely measure expression levels. Based on this, we define minimal data requirements for gap gene network inference. Our results demonstrate the feasibility of reverse-engineering with much reduced experimental effort. This enables more widespread use of the method in different developmental contexts and organisms. Such systematic application of data-driven models to real-world networks has enormous potential. Only the quantitative investigation of a large number of developmental gene regulatory networks will allow us to discover whether there are rules or regularities governing development and evolution of complex multi-cellular organisms.

A system of repressor gradients spatially organizes the boundaries of Bicoid-dependent target genes

The homeodomain (HD) protein Bicoid (Bcd) is thought to function as a gradient morphogen that positions boundaries of target genes via threshold-dependent activation mechanisms. Here, we analyze 66 Bcd-dependent regulatory elements and show that their boundaries are positioned primarily by repressive gradients that antagonize Bcd-mediated activation. A major repressor is the pair-rule protein Runt (Run), which is expressed in an opposing gradient and is necessary and sufficient for limiting Bcd-dependent activation. Evidence is presented that Run functions with the maternal repressor Capicua and the gap protein Kruppel as the principal components of a repression system that correctly orders boundaries throughout the anterior half of the embryo. These results put conceptual limits on the Bcd morphogen hypothesis and demonstrate how the Bcd gradient functions within the gene network that patterns the embryo.

Temporal coordination of gene networks by Zelda in the early Drosophila embryo

In past years, much attention has focused on the gene networks that regulate early developmental processes, but less attention has been paid to how multiple networks and processes are temporally coordinated. Recently the discovery of the transcriptional activator Zelda (Zld), which binds to CAGGTAG and related sequences present in the enhancers of many early-activated genes in Drosophila, hinted at a mechanism for how batteries of genes could be simultaneously activated. Here we use genome-wide binding and expression assays to identify Zld target genes in the early embryo with the goal of unraveling the gene circuitry regulated by Zld. We found that Zld binds to genes involved in early developmental processes such as cellularization, sex determination, neurogenesis, and pattern formation. In the absence of Zld, many target genes failed to be activated, while others, particularly the patterning genes, exhibited delayed transcriptional activation, some of which also showed weak and/or sporadic expression. These effects disrupted the normal sequence of patterning-gene interactions and resulted in highly altered spatial expression patterns, demonstrating the significance of a timing mechanism in early development. In addition, we observed prevalent overlap between Zld-bound regions and genomic "hotspot" regions, which are bound by many developmental transcription factors, especially the patterning factors. This, along with the finding that the most over-represented motif in hotspots, CAGGTA, is the Zld binding site, implicates Zld in promoting hotspot formation. We propose that Zld promotes timely and robust transcriptional activation of early-gene networks so that developmental events are coordinated and cell fates are established properly in the cellular blastoderm embryo.

Huckebein is part of a combinatorial repression code in the anterior blastoderm

The hierarchy of the segmentation cascade responsible for establishing the Drosophila body plan is composed by gap, pair-rule and segment polarity genes. However, no pair-rule stripes are formed in the anterior regions of the embryo. This lack of stripe formation, as well as other evidence from the literature that is further investigated here, led us to the hypothesis that anterior gap genes might be involved in a combinatorial mechanism responsible for repressing the cis-regulatory modules (CRMs) of hairy (h), even-skipped (eve), runt (run), and fushi-tarazu (ftz) anterior-most stripes. In this study, we investigated huckebein (hkb), which has a gap expression domain at the anterior tip of the embryo. Using genetic methods we were able to detect deviations from the wild-type patterns of the anterior-most pair-rule stripes in different genetic backgrounds, which were consistent with Hkb-mediated repression. Moreover, we developed an image processing tool that, for the most part, confirmed our assumptions. Using an hkb misexpression system, we further detected specific repression on anterior stripes. Furthermore, bioinformatics analysis predicted an increased significance of binding site clusters in the CRMs of h 1, eve 1, run 1 and ftz 1when Hkb was incorporated in the analysis, indicating that Hkb plays a direct role in these CRMs. We further discuss that Hkb and Slp1, which is the other previously identified common repressor of anterior stripes, might participate in a combinatorial repression mechanism controlling stripe CRMs in the anterior parts of the embryo and define the borders of these anterior stripes.

Combinatorial activation and concentration-dependent repression of the Drosophila even skipped stripe 3+7 enhancer

Despite years of study, the precise mechanisms that control position-specific gene expression during development are not understood. Here, we analyze an enhancer element from the even skipped (eve) gene, which activates and positions two stripes of expression (stripes 3 and 7) in blastoderm stage Drosophila embryos. Previous genetic studies showed that the JAK-STAT pathway is required for full activation of the enhancer, whereas the gap genes hunchback (hb) and knirps (kni) are required for placement of the boundaries of both stripes. We show that the maternal zinc-finger protein Zelda (Zld) is absolutely required for activation, and present evidence that Zld binds to multiple non-canonical sites. We also use a combination of in vitro binding experiments and bioinformatics analysis to redefine the Kni-binding motif, and mutational analysis and in vivo tests to show that Kni and Hb are dedicated repressors that function by direct DNA binding. These experiments significantly extend our understanding of how the eve enhancer integrates positive and negative transcriptional activities to generate sharp boundaries in the early embryo.

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.

giant is a bona fide gap gene in the intermediate germband insect, Oncopeltus fasciatus

Drosophila undergoes a form of development termed long germ segmentation, where all segments are specified nearly simultaneously so that by the blastoderm stage, the entire body plan has been determined. This mode of segmentation is evolutionarily derived. Most insects undergo short or intermediate germ segmentation, where only anterior segments are specified early, and posterior segments are sequentially specified during germband elongation. These embryological differences imply that anterior and posterior segments might rely upon different molecular mechanisms. In Drosophila, embryos mutant for giant show a gap in the anterior as well fusions of several abdominal segments. In Tribolium, a short germ beetle, giant is required for segmental identity, but not formation, in gnathal segments and also for segmentation of the entire abdomen. This raises the possibility that giant might not act as a gap gene in short and intermediate germ insects. Oncopeltus fasciatus is an intermediate germ insect that is an outgroup to the clade containing Drosophila and Tribolium. We cloned the Oncopeltus homolog of giant and determined its expression and function during segmentation. We find that Oncopeltus giant is a canonical gap gene in the maxillary and labial segments and also plays a gap-like role in the first four abdominal segments. Our results suggest that giant was a bona fide gap gene in the ancestor of these insects with this role being lost in the lineage leading towards Tribolium. This highlights the conservation of anterior patterning and evolutionary plasticity of the genetic regulation controlling posterior segmentation, even in short and intermediate germ insects.

Stochastic Spatio-Temporal Dynamic Model for Gene/Protein Interaction Network in Early Drosophila Development

In order to investigate the possible mechanisms for eve stripe formation of Drosophila embryo, a spatio-temporal gene/protein interaction network model is proposed to mimic dynamic behaviors of protein synthesis, protein decay, mRNA decay, protein diffusion, transcription regulations and autoregulation to analyze the interplay of genes and proteins at different compartments in early embryogenesis. In this study, we use the maximum likelihood (ML) method to identify the stochastic 3-D Embryo Space-Time (3-DEST) dynamic model for gene/protein interaction network via 3-D mRNA and protein expression data and then use the Akaike Information Criterion (AIC) to prune the gene/protein interaction network. The identified gene/protein interaction network allows us not only to analyze the dynamic interplay of genes and proteins on the border of eve stripes but also to infer that eve stripes are established and maintained by network motifs built by the cooperation between transcription regulations and diffusion mechanisms in early embryogenesis. Literature reference with the wet experiments of gene mutations provides a clue for validating the identified network. The proposed spatio-temporal dynamic model can be extended to gene/protein network construction of different biological phenotypes, which depend on compartments, e.g. postnatal stem/progenitor cell differentiation.

Cerebral cortex development: From progenitors patterning to neocortical size during evolution

The central nervous system is composed of thousands of distinct neurons that are assembled in a highly organized structure. In order to form functional neuronal networks, distinct classes of cells have to be generated in a precise number, in a spatial and temporal hierarchy and to be positioned at specific coordinates. An exquisite coordination of appropriate growth of competent territories and their patterning is required for regionalization and neurogenesis along both the anterior-posterior and dorso-ventral axis of the developing nervous system. The neocortex represents the brain territory that has undergone a major increase in its relative size during the course of mammalian evolution. In this review we will discuss how the fine tuning of growth and cell fate patterning plays a crucial role in the achievement of the final size of central nervous system structures and how divergence might have contributed to the surface increase of the cerebral cortex in mammals. In particular, we will describe how lack of precision might have been instrumental to neocortical evolution.

Anterior-posterior positional information in the absence of a strong Bicoid gradient

The Bicoid (Bcd) transcription factor is distributed as a long-range concentration gradient along the anterior posterior (AP) axis of the Drosophila embryo. Bcd is required for the activation of a series of target genes, which are expressed at specific positions within the gradient. Here we directly tested whether different concentration thresholds within the Bcd gradient establish the relative positions of its target genes by flattening the gradient and systematically varying expression levels. Genome-wide expression profiles were used to estimate the total number of Bcd target genes, and a general correlation was found between the Bcd concentration required for activation and the positions where target genes are expressed in wild-type embryos. However, concentrations required for target gene activation in embryos with flattened Bcd were consistently lower than those present at each target gene's position in the wild-type gradient, suggesting that Bcd is in excess at every position along the AP axis. Also, several Bcd target genes were positioned in correctly ordered stripes in embryos with flattened Bcd, and we suggest that these stripes are normally regulated by interactions between Bcd and the terminal patterning system. Our findings argue strongly against the strict interpretation of the Bcd morphogen hypothesis, and support the idea that target gene positioning involves combinatorial interactions that are mediated by the binding site architecture of each gene's cis-regulatory elements.

Precise registration of gene expression boundaries by a repressive morphogen in Drosophila

BACKGROUND

Morphogen gradients are thought to create concentration thresholds that differentially position the expression boundaries of multiple target genes. Despite intensive study, it is still unclear how the concentration profiles within gradients are spatially related to the critical patterning thresholds they generate.

RESULTS

Here we use a combination of quantitative measurements and ectopic-misexpression experiments to examine the transcriptional-repression activities of the Hunchback (Hb) protein gradient in Drosophila embryos. Our results define five expression boundaries that are set primarily by differences in Hb concentration and two boundaries that are set by combinatorial mechanisms involving Hb and at least one other repressor.

CONCLUSIONS

Hb functions as a repressive morphogen, but only within a specific range of concentrations ( approximately 40% to approximately 4.4% of maximum Hb concentration), within which there are at least four distinct concentration thresholds. The lower limit of the range reflects a position where the slope of the gradient becomes too shallow for resolution by specific target genes. Concentrations above the upper limit do not contribute directly to differential-repression mechanisms, but they provide a robust source that permits proper functioning of the gradient in heterozygous embryos that contain only one functional hb gene.

Deriving structure from evolution: metazoan segmentation

Segmentation is a common feature of disparate clades of metazoans, and its evolution is a central problem of evolutionary developmental biology. We evolved in silico regulatory networks by a mutation/selection process that just rewards the number of segment boundaries. For segmentation controlled by a static gradient, as in long-germ band insects, a cascade of adjacent repressors reminiscent of gap genes evolves. For sequential segmentation controlled by a moving gradient, similar to vertebrate somitogenesis, we invariably observe a very constrained evolutionary path or funnel. The evolved state is a cell autonomous 'clock and wavefront' model, with the new attribute of a separate bistable system driven by an autonomous clock. Early stages in the evolution of both modes of segmentation are functionally similar, and simulations suggest a possible path for their interconversion. Our computation illustrates how complex traits can evolve by the incremental addition of new functions on top of pre-existing traits.

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.

Enhancer responses to similarly distributed antagonistic gradients in development

Formation of spatial gene expression patterns in development depends on transcriptional responses mediated by gene control regions, enhancers. Here, we explore possible responses of enhancers to overlapping gradients of antagonistic transcriptional regulators in the Drosophila embryo. Using quantitative models based on enhancer structure, we demonstrate how a pair of antagonistic transcription factor gradients with similar or even identical spatial distributions can lead to the formation of distinct gene expression domains along the embryo axes. The described mechanisms are sufficient to explain the formation of the anterior and the posterior knirps expression, the posterior hunchback expression domain, and the lateral stripes of rhomboid expression and of other ventral neurogenic ectodermal genes. The considered principles of interaction between antagonistic gradients at the enhancer level can also be applied to diverse developmental processes, such as domain specification in imaginal discs, or even eyespot pattern formation in the butterfly wing.

The early development of the fruit fly embryo depends on an intricate but well-studied gene regulatory network. In fly eggs, maternally deposited gene products—morphogenes—form spatial concentration gradients. The graded distribution of the maternal morphogenes initiates a cascade of gene interactions leading to embryo development. Gradients of activators and repressors regulating common target genes may produce different outcomes depending on molecular mechanisms, mediating their function. Here, we describe quantitative mathematical models for the interplay between gradients of positive and negative transcriptional regulators—proteins, activating or repressing their target genes through binding the gene's regulatory DNA sequences. We predict possible spatial outcomes of the transcriptional antagonistic interactions in fly development and consider examples where the predicted cases may take place.

Multiple Promoter Targeting Sequences exist in Abdominal-B to regulate long-range gene activation

In complex genomes, insulators set up chromatin domain boundaries and protect promoters from inappropriate activation by enhancers from neighboring genes. The Drosophila Abdominal-B locus uses insulator elements to organize its large regulatory region into several body segment-specific chromatin domains. This organization leads to a problem in enhancer-promoter communication, that is, how do distal enhancers activate the Abd-B promoter when there are several insulators in between? This issue is partially resolved by the Promoter Targeting Sequence, which can overcome the enhancer blocking effect of an insulator. In this study, we describe a new Promoter Targeting Sequence, PTS-6, from the Abd-B 3' regulatory region. PTS-6, comprised of approximately 200 bp, was found to bypass both homologous Abdominal-B insulators, such as Fab-7 and Fab-8, and a heterologous insulator, suHw. Most importantly, it also overcomes a combination of two insulators such as Fab-7/Fab-8. Thus, PTS-6 could, in principle, target remote enhancers that are separated from the Abd-B promoter by multiple insulators. In addition, PTS-6 selectively targets the distal enhancer to only one transgenic promoter, and it strongly facilitates Abd-B enhancers. These results suggest that promoter targeting is necessary for long-range enhancer-promoter communication in Abd-B, and PTS elements could be a common occurrence in large, complex genetic loci.

The Promoter Targeting Sequence mediates epigenetically heritable transcription memory

Large gene complexes frequently use "specialized" DNA elements to ensure proper regulation of gene activities. The Promoter Targeting Sequence (PTS) from the Abdominal-B locus of the Drosophila Bithorax complex overcomes an insulator, and facilitates, yet restricts, distant enhancers to a single promoter. We found that this promoter-targeting activity is independent of an enhancer's tissue or temporal specificity, and can be remembered in all somatic cells in the absence of promoter activation. It requires an insulator for its establishment, but can be maintained by the PTS in the absence of an insulator. More importantly, the promoter-targeting activity can be remembered after the transgene is translocated to new chromosomal locations. These results suggest that promoter targeting is established independent of enhancer activity, and is maintained epigenetically throughout development and subsequent generations.

Groucho-dependent repression by Sloppy-paired 1 differentially positions anterior pair-rule stripes in the Drosophila embryo

The Drosophila body plan is composed of a linear array of cephalic, thoracic, and abdominal segments along the anterior posterior axis. The number and positions of individual segments are established by a transcriptional network comprised of maternal effect, gap, pair-rule, and segment polarity genes. The sloppy-paired (slp) locus contains two genes (slp1 and slp2) that are expressed in overlapping striped patterns in the presumptive thorax and abdomen. Previous studies suggest that these genes function at the pair-rule and segment polarity levels to establish the spacing and polarity of thoracic and abdominal segments. One of these genes (slp1) is also expressed in a broad anterior domain that appears before the striped patterns. There are severe cephalic defects in slp1 mutants, including the complete loss of the mandibular segment, but the molecular roles played by Slp1 in anterior patterning are not clear. Here, we present evidence that the anterior Slp1 domain acts as a gradient to differentially repress the anteriormost stripes of several different pair-rule genes. This repressive gradient contributes to the precise spatial arrangement of anterior pair-rule stripe borders required for expression of the first engrailed stripe and the formation of the mandibular segment. These results suggest that Slp1 functions as a gap gene-like repressor, in addition to its roles at the pair-rule and segment polarity levels of the hierarchy. The Slp1 protein contains a protein motif (EH1) which mediates binding to the transcriptional corepressor Groucho (Gro). We show that this domain is required for Slp1-mediated repression in vivo.

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.

Ventral dominance governs sequential patterns of gene expression across the dorsal–ventral axis of the neuroectoderm in the Drosophila embryo

A nuclear concentration gradient of the maternal transcription factor Dorsal establishes three tissues across the dorsal-ventral axis of precellular Drosophila embryos: mesoderm, neuroectoderm, and dorsal ectoderm. Subsequent interactions among Dorsal target genes subdivide the mesoderm and dorsal ectoderm. Here we investigate the subdivision of the neuroectoderm by three conserved homeobox genes, ventral nervous system defective (vnd), intermediate neuroblasts defective (ind), and muscle segment homeobox (msh). These genes divide the ventral nerve cord into three columns along the dorsal-ventral axis. Sequential patterns of vnd, ind, and msh expression are established prior to gastrulation and evidence is presented that these genes respond to distinct thresholds of the Dorsal gradient. Maintenance of these patterns depends on cross-regulatory interactions, whereby genes expressed in ventral regions repress those expressed in more dorsal regions. This "ventral dominance" includes regulatory genes that are expressed in the mesectoderm and mesoderm. At least some of these regulatory interactions are direct. For example, the misexpression of vnd in transgenic embryos represses ind and msh, and the addition of Vnd binding sites to a heterologous enhancer is sufficient to mediate repression. The N-terminal domain of Vnd contains a putative eh1 repression domain that binds Groucho in vitro. Mutations in this domain diminish Groucho binding and also attenuate repression in vivo. We discuss the significance of ventral dominance with respect to the patterning of the vertebrate neural tube, and compare it with the previously observed phenomenon of posterior prevalence, which governs sequential patterns of Hox gene expression across the anterior-posterior axis of metazoan embryos.

BEST: a novel computational approach for comparing gene expression patterns from early stages of Drosophila melanogaster development

Embryonic gene expression patterns are an indispensable part of modern developmental biology. Currently, investigators must visually inspect numerous images containing embryonic expression patterns to identify spatially similar patterns for inferring potential genetic interactions. The lack of a computational approach to identify pattern similarities is an impediment to advancement in developmental biology research because of the rapidly increasing amount of available embryonic gene expression data. Therefore, we have developed computational approaches to automate the comparison of gene expression patterns contained in images of early stage Drosophila melanogaster embryos (prior to the beginning of germ-band elongation); similarities and differences in gene expression patterns in these early stages have extensive developmental effects. Here we describe a basic expression search tool (BEST) to retrieve best matching expression patterns for a given query expression pattern and a computational device for gene interaction inference using gene expression pattern images and information on the associated genotypes and probes. Analysis of a prototype collection of Drosophila gene expression pattern images is presented to demonstrate the utility of these methods in identifying biologically meaningful matches and inferring gene interactions by direct image content analysis. In particular, the use of BEST searches for gene expression patterns is akin to that of BLAST searches for finding similar sequences. These computational developmental biology methodologies are likely to make the great wealth of embryonic gene expression pattern data easily accessible and to accelerate the discovery of developmental networks.

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.

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.

A homeodomain protein code specifies progenitor cell identity and neuronal fate in the ventral neural tube

Distinct classes of neurons are generated at defined positions in the ventral neural tube in response to a gradient of Sonic Hedgehog (Shh) activity. A set of homeodomain transcription factors expressed by neural progenitors act as intermediaries in Shh-dependent neural patterning. These homeodomain factors fall into two classes: class I proteins are repressed by Shh and class II proteins require Shh signaling for their expression. The profile of class I and class II protein expression defines five progenitor domains, each of which generates a distinct class of postmitotic neurons. Cross-repressive interactions between class I and class II proteins appear to refine and maintain these progenitor domains. The combinatorial expression of three of these proteins--Nkx6.1, Nkx2.2, and Irx3--specifies the identity of three classes of neurons generated in the ventral third of the neural tube.

A novel cis-regulatory element, the PTS, mediates an anti-insulator activity in the Drosophila embryo

The Abd-B Hox gene contains an extended 3' cis-regulatory region that is subdivided into a series of separate lab domains. The lab-7 domain activates Abd-B in parasegment 12 (ps12), whereas lab-8 controls expression in ps13. iab-7 is flanked by two insulators, Fab-7 and Fab-8, which are thought to prevent regulatory factors, such as Polycomb silencers, from influencing neighboring iab domains. This organization poses a potential paradox, since insulator DNAs can work in a dominant fashion to block enhancer-promoter interactions over long distances. Here, we present evidence for a novel cis-regulatory sequence located within lab-7, the promoter targeting sequence (PTS), which permits distal enhancers to overcome the blocking effects of Fab-8 and the heterologous su(Hw) insulator. We propose that the PTS converts dominant, long-range insulators into local regulatory elements that separate neighboring lab domains.

Analysis of an even-skipped rescue transgene reveals both composite and discrete neuronal and early blastoderm enhancers, and multi-stripe positioning by gap gene repressor gradients

The entire functional even-skipped locus of Drosophila melanogaster is contained within a 16 kilobase region. As a transgene, this region is capable of rescuing even-skipped mutant flies to fertile adulthood. Detailed analysis of the 7.7 kb of regulatory DNA 3′ of the transcription unit revealed ten novel, independently regulated patterns. Most of these patterns are driven by non-overlapping regulatory elements, including ones for syncytial blastoderm stage stripes 1 and 5, while a single element specifies both stripes 4 and 6. Expression analysis in gap gene mutants showed that stripe 5 is restricted anteriorly by Kruppel and posteriorly by giant, the same repressors that regulate stripe 2. Consistent with the coregulation of stripes 4 and 6 by a single cis-element, both the anterior border of stripe 4 and the posterior border of stripe 6 are set by zygotic hunchback, and the region between the two stripes is ‘carved out’ by knirps. Thus the boundaries of stripes 4 and 6 are set through negative regulation by the same gap gene domains that regulate stripes 3 and 7 (Small, S., Blair, A. and Levine, M. (1996) Dev. Biol. 175, 314–24), but at different concentrations. The 3′ region also contains a single element for neurogenic expression in ganglion mother cells 4–2a and 1–1a, and neurons derived from them (RP2, a/pCC), suggesting common regulators in these lineages. In contrast, separable elements were found for expression in EL neurons, U/CQ neurons and the mesoderm. The even-skipped 3′ untranslated region is required to maintain late stage protein expression in RP2 and a/pCC neurons, and appears to affect protein levels rather than mRNA levels. Additionally, a strong pairing-sensitive repression element was localized to the 3′ end of the locus, but was not found to contribute to efficient functional rescue.

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

Early developmental patterning of the Drosophila embryo is driven by the activities of a diverse set of maternally and zygotically derived transcription factors, including repressors encoded by gap genes such as Kruppel, knirps, giant and the mesoderm-specific snail. The mechanism of repression by gap transcription factors is not well understood at a molecular level. Initial characterization of these transcription factors suggests that they act as short-range repressors, interfering with the activity of enhancer or promoter elements 50 to 100 bp away. To better understand the molecular mechanism of short-range repression, we have investigated the properties of the Giant gap protein. We tested the ability of endogenous Giant to repress when bound close to the transcriptional initiation site and found that Giant effectively represses a heterologous promoter when binding sites are located at −55 bp with respect to the start of transcription. Consistent with its role as a short-range repressor, as the binding sites are moved to more distal locations, repression is diminished. Rather than exhibiting a sharp ‘step-function’ drop-off in activity, however, repression is progressively restricted to areas of highest Giant concentration. Less than a two-fold difference in Giant protein concentration is sufficient to determine a change in transcriptional status of a target gene. This effect demonstrates that Giant protein gradients can be differentially interpreted by target promoters, depending on the exact location of the Giant binding sites within the gene. Thus, in addition to binding site affinity and number, cis element positioning within a promoter can affect the response of a gene to a repressor gradient. We also demonstrate that a chimeric Gal4-Giant protein lacking the basic/zipper domain can specifically repress reporter genes, suggesting that the Giant effector domain is an autonomous repression domain.