Kinetic analysis of segmentation gene interactions in Drosophila embryos

Precisely timed regulation of enhancer activity defines the binary expression pattern of Fushi tarazu in the Drosophila embryo

The genes that drive development each typically have many different enhancers. Properly coordinating the action of these different enhancers is crucial to correctly specifying cell-fate decisions, yet it remains poorly understood how their activity is choregraphed in time. To shed light on this question, we used recently developed single-cell live imaging tools to dissect the regulation of Fushi tarazu (Ftz) in Drosophila melanogaster embryos. Ftz is a transcription factor that is expressed in asymmetric stripes by two distinct enhancers: autoregulatory and zebra. The anterior edge of each stripe needs to be sharply defined to specify essential lineage boundaries. Here, we tracked how boundary cells commit to either a high-Ftz or low-Ftz fate by measuring Ftz protein traces in real time and simultaneously quantifying transcription from the endogenous locus and individual enhancers. This revealed that the autoregulatory enhancer does not establish this fate choice. Instead, it perpetuates the decision defined by zebra. This is contrary to the prevailing view that autoregulation drives the fate decision by causing bi-stable Ftz expression. Furthermore, we showed that the autoregulatory enhancer is not activated based on a Ftz-concentration threshold but through a timing-based mechanism. We hypothesize that this is regulated by several ubiquitously expressed pioneer-like transcription factors, which have recently been shown to act as timers in the embryo. Our work provides new insight into how precisely timed enhancer activity can directly regulate the dynamics of gene regulatory networks, which may be a general mechanism for ensuring that embryogenesis runs like clockwork.

A timer gene network is spatially regulated by the terminal system in the Drosophila embryo

In insect embryos, anteroposterior patterning is coordinated by the sequential expression of the 'timer' genes caudal, Dichaete, and odd-paired, whose expression dynamics correlate with the mode of segmentation. In Drosophila, the timer genes are expressed broadly across much of the blastoderm, which segments simultaneously, but their expression is delayed in a small 'tail' region, just anterior to the hindgut, which segments during germband extension. Specification of the tail and the hindgut depends on the terminal gap gene tailless, but beyond this the regulation of the timer genes is poorly understood. We used a combination of multiplexed imaging, mutant analysis, and gene network modelling to resolve the regulation of the timer genes, identifying 11 new regulatory interactions and clarifying the mechanism of posterior terminal patterning. We propose that a dynamic Tailless expression gradient modulates the intrinsic dynamics of a timer gene cross-regulatory module, delineating the tail region and delaying its developmental maturation.

Dynamic expression of Drosophila segmental cell surface-encoding genes and their pair-rule regulators

Drosophila segmentation is regulated by a complex network of transcription factors that include products of the pair-rule genes (PRGs). PRGs are expressed in early embryos in the primorida of alternate segmental units, establishing the repeated, segmental body plan of the fly. Despite detailed analysis of the regulatory logic among segmentation genes, the relationship between these genes and the morphological formation of segments is still poorly understood, since regulation of transcription factor expression is not sufficient to explain how segments actually form and are maintained. Cell surface proteins containing Leucine rich repeats (LRR) play a variety of roles in development, and those expressed in segmental patterns likely impact segment morphogenesis. Here we explore the relationships between the PRG network and segmentally expressed LRR-encoding (sLRR) genes. We examined expression of Toll2, Toll6, Toll7, Toll8 and tartan (trn) in wild type or PRG mutant embryos. Expression of each sLRR-encoding gene is dynamic, but each has a unique register along the anterior-posterior axis. The registers for different sLRRs are off-set from one another resulting in a continually changing set of overlapping expression patterns among the sLRR-encoding genes themselves and between the sLRR-encoding genes and the PRGs. Accordingly, each sLRR-encoding gene is regulated by a unique combination of PRGs. These findings suggest that one role of the PRG network is to promote segmentation by establishing a cell surface code: each row of cells in the two-segment-wide primordia expresses a unique combination of sLRRs, thereby translating regulatory information from the PRGs to direct segment morphogenesis.

Dynamic patterning by the Drosophila pair-rule network reconciles long-germ and short-germ segmentation

Drosophila segmentation is a well-established paradigm for developmental pattern formation. However, the later stages of segment patterning, regulated by the "pair-rule" genes, are still not well understood at the system level. Building on established genetic interactions, I construct a logical model of the Drosophila pair-rule system that takes into account the demonstrated stage-specific architecture of the pair-rule gene network. Simulation of this model can accurately recapitulate the observed spatiotemporal expression of the pair-rule genes, but only when the system is provided with dynamic "gap" inputs. This result suggests that dynamic shifts of pair-rule stripes are essential for segment patterning in the trunk and provides a functional role for observed posterior-to-anterior gap domain shifts that occur during cellularisation. The model also suggests revised patterning mechanisms for the parasegment boundaries and explains the aetiology of the even-skipped null mutant phenotype. Strikingly, a slightly modified version of the model is able to pattern segments in either simultaneous or sequential modes, depending only on initial conditions. This suggests that fundamentally similar mechanisms may underlie segmentation in short-germ and long-germ arthropods.

Odd-paired controls frequency doubling in Drosophila segmentation by altering the pair-rule gene regulatory network

The Drosophila embryo transiently exhibits a double-segment periodicity, defined by the expression of seven 'pair-rule' genes, each in a pattern of seven stripes. At gastrulation, interactions between the pair-rule genes lead to frequency doubling and the patterning of 14 parasegment boundaries. In contrast to earlier stages of Drosophila anteroposterior patterning, this transition is not well understood. By carefully analysing the spatiotemporal dynamics of pair-rule gene expression, we demonstrate that frequency-doubling is precipitated by multiple coordinated changes to the network of regulatory interactions between the pair-rule genes. We identify the broadly expressed but temporally patterned transcription factor, Odd-paired (Opa/Zic), as the cause of these changes, and show that the patterning of the even-numbered parasegment boundaries relies on Opa-dependent regulatory interactions. Our findings indicate that the pair-rule gene regulatory network has a temporally modulated topology, permitting the pair-rule genes to play stage-specific patterning roles.

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

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

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 and epigene networks: two levels in organizing the hereditary system

Intracellular governing gene networks consisting of genes and regulatory bonds among them are considered as the first level in organizing the hereditary system. We give examples of both prokaryotic gene network that controls the development of the lambda-phage and eukaryotic gene network that controls the early Drosophila ontogenesis. Using the method of generalized threshold models kinetic curves are shown for some gene products of these networks. Gene networks that govern ontogenetic processes can be envisioned as epigene networks, the networks of the subsequent level in organizing the hereditary systems. Based on the mathematical model our computer experiments show that even the simplest hypothetical two-epigene network is capable of ensuring divergent determination, conservation of determinate states and reproduction of the initial "zygotic" functional state. In addition, the experimental results are given on construction an artificial epigene.

Stripy Ftz target genes are coordinately regulated by Ftz-F1

During development, cascades of regulatory genes act in a hierarchical fashion to subdivide the embryo into increasingly specified body regions. This has been best characterized in Drosophila, where genes encoding regulatory transcription factors form a network to direct the development of the basic segmented body plan. The pair-rule genes are pivotal in this process as they are responsible for the first subdivision of the embryo into repeated metameric units. The Drosophila pair-rule gene fushi tarazu (ftz) is a derived Hox gene expressed in and required for the development of alternate parasegments. Previous studies suggested that Ftz achieves its distinct regulatory specificity as a segmentation protein by interacting with a ubiquitously expressed cofactor, the nuclear receptor Ftz-F1. However, the downstream target genes regulated by Ftz and other pair-rule genes to direct segment formation are not known. In this study, we selected candidate Ftz targets by virtue of their early expression in Ftz-like stripes. This identified two new Ftz target genes, drumstick (drm) and no ocelli (noc), and confirmed that Ftz regulates a serotonin receptor (5-HT2). These are the earliest Ftz targets identified to date and all are coordinately regulated by Ftz-F1. Engrailed (En), the best-characterized Ftz/Ftz-F1 downstream target, is not an intermediate in regulation. The drm genomic region harbors two separate seven-stripe enhancers, identified by virtue of predicted Ftz-F1 binding sites, and these sites are necessary for stripe expression in vivo. We propose that pair-rule genes, exemplified by Ftz/Ftz-F1, promote segmentation by acting at different hierarchical levels, regulating first, other segmentation genes; second, other regulatory genes that in turn control specific cellular processes such as tissue differentiation; and, third, 'segmentation realizator genes' that are directly involved in morphogenesis.

Modeling the dynamics of transcriptional gene regulatory networks for animal development

The dynamic process of cell fate specification is regulated by networks of regulatory genes. The architecture of the network defines the temporal order of specification events. To understand the dynamic control of the developmental process, the kinetics of mRNA and protein synthesis and the response of the cis-regulatory modules to transcription factor concentration must be considered. Here we review mathematical models for mRNA and protein synthesis kinetics which are based on experimental measurements of the rates of the relevant processes. The model comprises the response functions of cis-regulatory modules to their transcription factor inputs, by incorporating binding site occupancy and its dependence on biologically measurable quantities. We use this model to simulate gene expression, to distinguish between cis-regulatory execution of "AND" and "OR" logic functions, rationalize the oscillatory behavior of certain transcriptional auto-repressors and to show how linked subcircuits can be dealt with. Model simulations display the effects of mutation of binding sites, or perturbation of upstream gene expression. The model is a generally useful tool for understanding gene regulation and the dynamics of cell fate specification.

Global gene expression analysis identifies PDEF transcriptional networks regulating cell migration during cancer progression

Prostate derived ETS factor (PDEF) is an ETS (epithelial-specific E26 transforming sequence) family member that has been identified as a potential tumor suppressor. In multiple invasive breast cancer cells, PDEF expression inhibits cell migration by preventing the acquisition of directional morphological polarity conferred by changes in cytoskeleton organization. In this study, microarray analysis was used to identify >200 human genes that displayed a common differential expression pattern in three invasive breast cancer cell lines after expression of exogenous PDEF protein. Gene ontology associations and data mining analysis identified focal adhesion, adherens junctions, cell adhesion, and actin cytoskeleton regulation as cell migration-associated interaction pathways significantly impacted by PDEF expression. Validation experiments confirmed the differential expression of four cytoskeleton-associated genes with known functional associations with these pathways: uPA, uPAR, LASP1, and VASP. Significantly, chromatin immunoprecipitation studies identified PDEF as a direct negative regulator of the metastasis-associated gene uPA and phenotypic rescue experiments demonstrate that exogenous urokinase plasminogen activator (uPA) expression can restore the migratory ability of invasive breast cancer cells expressing PDEF. Furthermore, immunofluorescence studies identify the subcellular relocalization of urokinase plasminogen activator receptor (uPAR), LIM and SH3 protein (LASP1), and vasodilator-stimulated protein (VASP) as a possible mechanism accounting for the loss of morphological polarity observed upon PDEF expression.

Specification and positioning of parasegment grooves in Drosophila

Developmental boundaries ensure that cells fated to participate in a particular structure are brought together or maintained at the appropriate locale within developing embryos. Parasegment grooves mark the position of boundaries that separate every segment of the Drosophila embryo into anterior and posterior compartments. Here, we dissect the genetic hierarchy that controls the formation of this morphological landmark. We report that primary segment polarity genes (engrailed, hedgehog and wingless) are not involved in specifying the position of parasegment grooves. Wingless signalling plays only a permissive role by triggering the formation of grooves at cellular interfaces defined by the ON/OFF state of expression of the earlier acting pair-rule genes eve and ftz. We suggest that the transcription factors encoded by these genes activate two programmes in parallel: a cell fate programme mediated by segment polarity genes and a boundary/epithelial integrity programme mediated by unknown target genes.

Controlling gene expression in Drosophila using engineered zinc finger protein transcription factors

Zinc finger protein transcription factors (ZFP TFs) have been designed to control the expression of endogenous genes in a variety of cells. However, thus far the use of engineered ZFP TFs in germline transgenic settings has been restricted to plants. Here we report that ZFP TFs can regulate gene expression in transgenic Drosophila. To demonstrate this, we targeted the promoter of the well-characterized fushi tarazu (ftz) gene with a ZFP TF activator using the VP16 activation domain from Herpes simplex virus, and ZFP TF repressors using the Drosophila methyl-CpG binding domain (MBD)-like Delta protein. Heat-shock-inducible expression of the ZFP TF activator and repressors resulted in reciprocal effects on ftz regulation, as deduced from changes in the staining pattern and intensity of ftz and en gene expression, and from the cuticular analysis of first instar larvae. These data demonstrate the utility of ZFP TFs as tools for controlling gene expression in the context of a metazoan organism.

A screen for genes that interact with the Drosophila pair-rule segmentation gene fushi tarazu

The pair-rule gene fushi tarazu (ftz) of Drosophila is expressed at the blastoderm stage in seven stripes that serve to define the even-numbered parasegments. ftz encodes a DNA-binding homeodomain protein and is known to regulate genes of the segment polarity, homeotic, and pair-rule classes. Despite intensive analysis in a number of laboratories, how ftz is regulated and how it controls its targets are still poorly understood. To help understand these processes, we conducted a screen to identify dominant mutations that enhance the lethality of a ftz temperature-sensitive mutant. Twenty-six enhancers were isolated, which define 21 genes. All but one of the mutations recovered show a maternal effect in their interaction with ftz. Three of the enhancers proved to be alleles of the known ftz protein cofactor gene ftz-f1, demonstrating the efficacy of the screen. Four enhancers are alleles of Atrophin (Atro), the Drosophila homolog of the human gene responsible for the neurodegenerative disease dentatorubral-pallidoluysian atrophy. Embryos from Atro mutant germ-line mothers lack the even-numbered (ftz-dependent) engrailed stripes and show strong ftz-like segmentation defects. These defects likely result from a reduction in Even-skipped (Eve) repression ability, as Atro has been shown to function as a corepressor for Eve. In this study, we present evidence that Atro is also a member of the trithorax group (trxG) of Hox gene regulators. Atro appears to be particularly closely related in function to the trxG gene osa, which encodes a component of the brahma chromatin remodeling complex. One additional gene was identified that causes pair-rule segmentation defects in embryos from homozygous mutant germ-line mothers. The single allele of this gene, called bek, also causes nuclear abnormalities similar to those caused by alleles of the Trithorax-like gene, which encodes the GAGA factor.

Drawing lines in the sand: even skipped et al. and parasegment boundaries

The pair-rule segmentation gene even skipped (eve) is required to activate engrailed stripes and to organize odd-numbered parasegments (PSs). The protein product Eve has been shown to be an active repressor of transcription, and recent models for Eve function suggest that activation of engrailed is indirect, but these models have not been fully tested. Here we identify the forkhead domain transcription factor Sloppy-paired as the key intermediate in the initial activation of engrailed by Eve in odd-numbered parasegments. We also analyze the roles of the transcription factors Runt and Odd-skipped in this process. Detailed analysis of engrailed and pair-rule gene expression in various mutant combinations shows how eve activates engrailed by repressing these engrailed repressors, and further indicates that mutual repression among pair-rule genes plays an important role in establishing parasegment boundaries. We present a new model of pair-rule gene function that explains the response of these boundaries to the relative levels of Eve and Fushi Tarazu.

Segmenting the fly embryo: a logical analysis of the pair-rule cross-regulatory module

This manuscript reports a dynamical analysis of the pair-rule cross-regulatory module controlling segmentation in Drosophila melanogaster. We propose a logical model accounting for the ability of the pair-rule module to determine the formation of alternate juxtaposed Engrailed- and Wingless-expressing cells that form the (para)segmental boundaries. This module has the intrinsic capacity to generate four distinct expression states, each characterized by the expression of a particular combination of pair-rule genes or expression mode. The selection of one of these expression modes depends on the maternal and gap inputs, but also crucially on cross-regulations among pair-rule genes. The latter are instrumental in the interpretation of the maternal-gap pre-pattern. Our logical model allows the qualitative reproduction of the patterns of pair-rule gene expressions corresponding to the wild type situation, to loss-of-function and cis-regulatory mutations, and to ectopic pair-rule expressions. Furthermore, this model provides a formal explanation for the morphogenetic role of the initial bell-shaped expression of the gene even-skipped, i.e. for the distinct effects of different levels of the Even-skipped protein on its target pair-rule genes. It also accounts for the requirement of Even-skipped for the formation of all Engrailed-stripes. Finally, it provides new insights into the roles and evolutionary origins of the apparent redundancies in the regulatory architecture of the pair-rule module.

Homotypic regulatory clusters in Drosophila

Cis-regulatory modules (CRMs) are transcription regulatory DNA segments (approximately 1 Kb range) that control the expression of developmental genes in higher eukaryotes. We analyzed clustering of known binding motifs for transcription factors (TFs) in over 60 known CRMs from 20 Drosophila developmental genes, and we present evidence that each type of recognition motif forms significant clusters within the regulatory regions regulated by the corresponding TF. We demonstrate how a search with a single binding motif can be applied to explore gene regulatory networks and to discover coregulated genes in the genome. We also discuss the potential of the clustering method in interpreting the differential response of genes to various levels of transcriptional regulators.

Evolution of Ftz protein function in insects

The Drosophila gene fushi tarazu (ftz) encodes a homeodomain-containing transcriptional regulator (Ftz) required at several stages during development. Drosophila melanogaster ftz (Dm-ftz) is first expressed in seven stripes defining alternate parasegments of the embryo--a "pair-rule" segmentation function [1, 2]. It is then expressed in specific neural precursor cells in the central nervous system and finally in the developing hindgut [3]. An Orthopteran ortholog of ftz (Sg-ftz, formally Dax) has been isolated from the grasshopper Schistocerca gregaria [4]. The pattern of Sg-ftz expression in Schistocerca embryos suggests that some developmental roles of the ftz gene are likely to be conserved between these two species (e.g., CNS functions) while others may have diverged (e.g., segmentation functions). To test whether the function of the Ftz protein itself differs between these two species, here we compare the functions of Sg-Ftz and Dm-Ftz proteins by expressing both in Drosophila embryos. Sg-ftz mimics only poorly several segmentation roles of Dm-ftz (engrailed activation, wingless repression, and embryonic cuticle transformation). However, the two proteins are similarly active in the rescue of a CNS-specific ftz mutant. These findings argue that this ftz CNS function is mediated by conserved parts of the protein, while efficient pair-rule function requires sequences present specifically in the Drosophila protein.

Establishment and maintenance of parasegmental compartments

Embryos of higher metazoans are divided into repeating units early in development. In Drosophila, the earliest segmental units to form are the parasegments. Parasegments are initially defined by alternating stripes of expression of the fushi-tarazu and even-skipped genes. How fushi-tarazu and even-skipped define the parasegment boundaries, and how parasegments are lost when fushi-tarazu or even-skipped fail to function correctly, have never been fully or properly explained. Here we show that parasegment widths are defined early by the relative levels of fushi-tarazu and even-skipped at stripe junctions. Changing these levels results in alternating wide and narrow parasegments. When shifted by 30% or more, the enlarged parasegments remain enlarged and the reduced parasegments are lost. Loss of the reduced parasegments occurs in three steps; delamination of cells from the epithelial layer, apoptosis of the delaminated cells and finally apoptosis of inappropriate cells remaining at the surface. The establishment and maintenance of vertebrate metameres may be governed by similar processes and properties.

Mechanisms regulating target gene selection by the homeodomain-containing protein Fushi tarazu

Homeodomain proteins are DNA-binding transcription factors that control major developmental patterning events. Although DNA binding is mediated by the homeodomain, interactions with other transcription factors play an unusually important role in the selection and regulation of target genes. A major question in the field is whether these cofactor interactions select target genes by modulating DNA binding site specificity (selective binding model), transcriptional activity (activity regulation model) or both. A related issue is whether the number of target genes bound and regulated is a small or large percentage of genes in the genome. In this study, we have addressed these issues using a chimeric protein that contains the strong activation domain of the viral VP16 protein fused to the Drosophila homeodomain-containing protein Fushi tarazu (Ftz). We find that genes previously thought not to be direct targets of Ftz remain unaffected by FtzVP16. Addition of the VP16 activation domain to Ftz does, however, allow it to regulate previously identified target genes at times and in regions that Ftz alone cannot. It also changes Ftz into an activator of two genes that it normally represses. Taken together, the results suggest that Ftz binds and regulates a relatively limited number of target genes, and that cofactors affect target gene specificity primarily by controlling binding site selection. Activity regulation then fine-tunes the temporal and spatial domains of promoter responses, the magnitude of these responses, and whether they are positive or negative.