Combinatorial patterning mechanisms in the Drosophila embryo

Expression patterns of FGF and BMP pathway genes in the tardigrade Hypsibius exemplaris

A small number of conserved signaling pathways regulate development of most animals, yet we do not know where these pathways are deployed in most embryos. This includes tardigrades, a phylum with a unique body shape. We examined expression patterns of components of the BMP and FGF signaling pathways during embryonic segmentation and mesoderm development of the tardigrade Hypsibius exemplaris. Among the patterns examined, we found that an FGF ligand gene is expressed in ectodermal segment posteriors and an FGF receptor gene is expressed in underlying endomesodermal pouches, suggesting possible FGF signaling between these developing germ layers. We found that a BMP ligand gene is expressed in lateral ectoderm and dorsolateral bands along segment posteriors, while the BMP antagonist Sog gene is expressed in lateral ectoderm and also in a subset of endomesodermal cells, suggesting a possible role of BMP signaling in dorsal-ventral patterning of lateral ectoderm. In combination with known roles of these pathways during development of common model systems, we developed hypotheses for how the BMP and FGF pathways might regulate embryo segmentation and mesoderm formation of the tardigrade H. exemplaris. These results identify the expression patterns of genes from two conserved signaling pathways for the first time in the tardigrade phylum.

Temporal integration of inductive cues on the way to gastrulation

Markers for the endoderm and mesoderm germ layers are commonly expressed together in the early embryo, potentially reflecting cells' ability to explore potential fates before fully committing. It remains unclear when commitment to a single-germ layer is reached and how it is impacted by external signals. Here, we address this important question in Drosophila, a convenient model system in which mesodermal and endodermal fates are associated with distinct cellular movements during gastrulation. Systematically applying endoderm-inducing extracellular signal-regulated kinase (ERK) signals to the ventral medial embryo-which normally only receives a mesoderm-inducing cue-reveals a critical time window during which mesodermal cell movements and gene expression are suppressed by proendoderm signaling. We identify the ERK target gene huckebein (hkb) as the main cause of the ventral furrow suppression and use computational modeling to show that Hkb repression of the mesoderm-associated gene snail is sufficient to account for a broad range of transcriptional and morphogenetic effects. Our approach, pairing precise signaling perturbations with observation of transcriptional dynamics and cell movements, provides a general framework for dissecting the complexities of combinatorial tissue patterning.

Self-organized cytoskeletal alignment during Drosophila mesoderm invagination

During tissue morphogenesis, mechanical forces are propagated across tissues, resulting in tissue shape changes. These forces in turn can influence cell behaviour, leading to a feedback process that can be described as self-organizing. Here, I discuss cytoskeletal self-organization and point to evidence that suggests its role in directing force during morphogenesis. During Drosophila mesoderm invagination, the shape of the region of cells that initiates constriction creates a mechanical pattern that in turn aligns the cytoskeleton with the axis of greatest resistance to contraction. The wild-type direction of the force controls the shape and orientation of the invaginating mesoderm. Given the ability of the actomyosin cytoskeleton to self-organize, these types of feedback mechanisms are likely to play important roles in a range of different morphogenetic events. This article is part of the discussion meeting issue 'Contemporary morphogenesis'.

Regulation of spatiotemporal limits of developmental gene expression via enhancer grammar

The regulatory specificity of a gene is determined by the structure of its enhancers, which contain multiple transcription factor binding sites. A unique combination of transcription factor binding sites in an enhancer determines the boundary of target gene expression, and their disruption often leads to developmental defects. Despite extensive characterization of binding motifs in an enhancer, it is still unclear how each binding site contributes to overall transcriptional activity. Using live imaging, quantitative analysis, and mathematical modeling, we measured the contribution of individual binding sites in transcriptional regulation. We show that binding site arrangement within the Rho-GTPase component t48 enhancer mediates the expression boundary by mainly regulating the timing of transcriptional activation along the dorsoventral axis of Drosophila embryos. By tuning the binding affinity of the Dorsal (Dl) and Zelda (Zld) sites, we show that single site modulations are sufficient to induce significant changes in transcription. Yet, no one site seems to have a dominant role; rather, multiple sites synergistically drive increases in transcriptional activity. Interestingly, Dl and Zld demonstrate distinct roles in transcriptional regulation. Dl site modulations change spatial boundaries of t48, mostly by affecting the timing of activation and bursting frequency rather than transcriptional amplitude or bursting duration. However, modulating the binding site for the pioneer factor Zld affects both the timing of activation and amplitude, suggesting that Zld may potentiate higher Dl recruitment to target DNAs. We propose that such fine-tuning of dynamic gene control via enhancer structure may play an important role in ensuring normal development.

The Physical Mechanisms of Drosophila Gastrulation: Mesoderm and Endoderm Invagination

A critical juncture in early development is the partitioning of cells that will adopt different fates into three germ layers: the ectoderm, the mesoderm, and the endoderm. This step is achieved through the internalization of specified cells from the outermost surface layer, through a process called gastrulation. In Drosophila, gastrulation is achieved through cell shape changes (i.e., apical constriction) that change tissue curvature and lead to the folding of a surface epithelium. Folding of embryonic tissue results in mesoderm and endoderm invagination, not as individual cells, but as collective tissue units. The tractability of Drosophila as a model system is best exemplified by how much we know about Drosophila gastrulation, from the signals that pattern the embryo to the molecular components that generate force, and how these components are organized to promote cell and tissue shape changes. For mesoderm invagination, graded signaling by the morphogen, Spätzle, sets up a gradient in transcriptional activity that leads to the expression of a secreted ligand (Folded gastrulation) and a transmembrane protein (T48). Together with the GPCR Mist, which is expressed in the mesoderm, and the GPCR Smog, which is expressed uniformly, these signals activate heterotrimeric G-protein and small Rho-family G-protein signaling to promote apical contractility and changes in cell and tissue shape. A notable feature of this signaling pathway is its intricate organization in both space and time. At the cellular level, signaling components and the cytoskeleton exhibit striking polarity, not only along the apical-basal cell axis, but also within the apical domain. Furthermore, gene expression controls a highly choreographed chain of events, the dynamics of which are critical for primordium invagination; it does not simply throw the cytoskeletal "on" switch. Finally, studies of Drosophila gastrulation have provided insight into how global tissue mechanics and movements are intertwined as multiple tissues simultaneously change shape. Overall, these studies have contributed to the view that cells respond to forces that propagate over great distances, demonstrating that cellular decisions, and, ultimately, tissue shape changes, proceed by integrating cues across an entire embryo.

Formation, interpretation, and regulation of the Drosophila Dorsal/NF-κB gradient

The morphogen gradient of the transcription factor Dorsal in the early Drosophila embryo has become one of the most widely studied tissue patterning systems. Dorsal is a Drosophila homolog of mammalian NF-κB and patterns the dorsal-ventral axis of the blastoderm embryo into several tissue types by spatially regulating upwards of 100 zygotic genes. Recent studies using fluorescence microscopy and live imaging have quantified the Dorsal gradient and its target genes, which has paved the way for mechanistic modeling of the gradient. In this review, we describe the mechanisms behind the initiation of the Dorsal gradient and its regulation of target genes. The main focus of the review is a discussion of quantitative and computational studies of the Dl gradient system, including regulation of the Dl gradient. We conclude with a discussion of potential future directions.

A facilitated diffusion mechanism establishes the Drosophila Dorsal gradient

The transcription factor NF-κB plays an important role in the immune system, apoptosis and inflammation. Dorsal, a Drosophila homolog of NF-κB, patterns the dorsal-ventral axis in the blastoderm embryo. During this stage, Dorsal is sequestered outside the nucleus by the IκB homolog Cactus. Toll signaling on the ventral side breaks the Dorsal/Cactus complex, allowing Dorsal to enter the nucleus to regulate target genes. Fluorescent data show that Dorsal accumulates on the ventral side of the syncytial blastoderm. Here, we use modeling and experimental studies to show that this accumulation is caused by facilitated diffusion, or shuttling, of the Dorsal/Cactus complex. We also show that active Toll receptors are limiting in wild-type embryos, which is a key factor in explaining global Dorsal gradient formation. Our results suggest that shuttling is necessary for viability of embryos from mothers with compromised dorsal levels. Therefore, Cactus not only has the primary role of regulating Dorsal nuclear import, but also has a secondary role in shuttling. Given that this mechanism has been found in other, independent, systems, we suggest that it might be more prevalent than previously thought.

Stepwise Progression of Embryonic Patterning

It is long established that the graded distribution of Dorsal transcription factor influences spatial domains of gene expression along the dorsoventral (DV) axis of Drosophila melanogaster embryos. However, the more recent realization that Dorsal levels also change with time raises the question of whether these dynamics are instructive. An overview of DV axis patterning is provided, focusing on new insights identified through quantitative analysis of temporal changes in Dorsal target gene expression from one nuclear cycle to the next ('steps'). Possible roles for the stepwise progression of this patterning program are discussed including (i) tight temporal regulation of signaling pathway activation, (ii) control of gene expression cohorts, and (iii) ensuring the irreversibility of the patterning and cell fate specification process.

A conserved role for Snail as a potentiator of active transcription

The transcription factors of the Snail family are key regulators of epithelial-mesenchymal transitions, cell morphogenesis, and tumor metastasis. Since its discovery in Drosophila ∼25 years ago, Snail has been extensively studied for its role as a transcriptional repressor. Here we demonstrate that Drosophila Snail can positively modulate transcriptional activation. By combining information on in vivo occupancy with expression profiling of hand-selected, staged snail mutant embryos, we identified 106 genes that are potentially directly regulated by Snail during mesoderm development. In addition to the expected Snail-repressed genes, almost 50% of Snail targets showed an unanticipated activation. The majority of "Snail-activated" genes have enhancer elements cobound by Twist and are expressed in the mesoderm at the stages of Snail occupancy. Snail can potentiate Twist-mediated enhancer activation in vitro and is essential for enhancer activity in vivo. Using a machine learning approach, we show that differentially enriched motifs are sufficient to predict Snail's regulatory response. In silico mutagenesis revealed a likely causative motif, which we demonstrate is essential for enhancer activation. Taken together, these data indicate that Snail can potentiate enhancer activation by collaborating with different activators, providing a new mechanism by which Snail regulates development.

Size-dependent regulation of dorsal-ventral patterning in the early Drosophila embryo

How natural variation in embryo size affects patterning of the Drosophila embryo dorsal-ventral (DV) axis is not known. Here we examined quantitatively the relationship between nuclear distribution of the Dorsal transcription factor, boundary positions for several target genes, and DV axis length. Data were obtained from embryos of a wild-type background as well as from mutant lines inbred to size select embryos of smaller or larger sizes. Our data show that the width of the nuclear Dorsal gradient correlates with DV axis length. In turn, for some genes expressed along the DV axis, the boundary positions correlate closely with nuclear Dorsal levels and with DV axis length; while the expression pattern of others is relatively constant and independent of the width of the Dorsal gradient. In particular, the patterns of snail (sna) and ventral nervous system defective (vnd) correlate with nuclear Dorsal levels and exhibit scaling to DV length; while the pattern of intermediate neuroblasts defective (ind) remains relatively constant with respect to changes in Dorsal and DV length. However, in mutants that exhibit an abnormal expansion of the Dorsal gradient which fails to scale to DV length, only sna follows the Dorsal distribution and exhibits overexpansion; in contrast, vnd and ind do not overexpand suggesting some additional mechanism acts to refine the dorsal boundaries of these two genes. Thus, our results argue against the idea that the Dorsal gradient works as a global system of relative coordinates along the DV axis and suggest that individual targets respond to changes in embryo size in a gene-specific manner.

Transcription factors: from enhancer binding to developmental control

Developmental progression is driven by specific spatiotemporal domains of gene expression, which give rise to stereotypically patterned embryos even in the presence of environmental and genetic variation. Views of how transcription factors regulate gene expression are changing owing to recent genome-wide studies of transcription factor binding and RNA expression. Such studies reveal patterns that, at first glance, seem to contrast with the robustness of the developmental processes they encode. Here, we review our current knowledge of transcription factor function from genomic and genetic studies and discuss how different strategies, including extensive cooperative regulation (both direct and indirect), progressive priming of regulatory elements, and the integration of activities from multiple enhancers, confer specificity and robustness to transcriptional regulation during development.

RTK signaling modulates the Dorsal gradient

The dorsoventral (DV) axis of the Drosophila embryo is patterned by a nuclear gradient of the Rel family transcription factor, Dorsal (Dl), that activates or represses numerous target genes in a region-specific manner. Here, we demonstrate that signaling by receptor tyrosine kinases (RTK) reduces nuclear levels and transcriptional activity of Dl, both at the poles and in the mid-body of the embryo. These effects depend on wntD, which encodes a Dl antagonist belonging to the Wingless/Wnt family of secreted factors. Specifically, we show that, via relief of Groucho- and Capicua-mediated repression, the Torso and EGFR RTK pathways induce expression of WntD, which in turn limits Dl nuclear localization at the poles and along the DV axis. Furthermore, this RTK-dependent control of Dl is important for restricting expression of its targets in both contexts. Thus, our results reveal a new mechanism of crosstalk, whereby RTK signals modulate the spatial distribution and activity of a developmental morphogen in vivo.

Pattern formation by graded and uniform signals in the early Drosophila embryo

The early Drosophila embryo is patterned by graded distributions of maternal transcription factors. Recent studies revealed that pattern formation by these graded signals depends on uniformly expressed transcriptional activators, such as Zelda. Removal of Zelda influences both the timing and the spatial expression domains for most of the genes controlled by maternal gradients. We demonstrate that some of these patterning defects, which range from temporal delay to loss of expression, can be rationalized with the use of a mathematical model based on cooperative binding of graded and uniform factors. This model makes a number of predictions, which we confirm experimentally by analyzing the expression of short gastrulation (sog), a gene that is controlled by a combination of the Dorsal morphogen gradient and Zelda. The proposed model suggests a general mechanism for the formation of nested gene expression domains, which is a hallmark of tissue patterning by morphogen gradients. According to this mechanism, the differential effects of a morphogen on its target genes can depend on their differential sensitivity to uniform factors.

A computational statistics approach for estimating the spatial range of morphogen gradients

A crucial issue in studies of morphogen gradients relates to their range: the distance over which they can act as direct regulators of cell signaling, gene expression and cell differentiation. To address this, we present a straightforward statistical framework that can be used in multiple developmental systems. We illustrate the developed approach by providing a point estimate and confidence interval for the spatial range of the graded distribution of nuclear Dorsal, a transcription factor that controls the dorsoventral pattern of the Drosophila 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.

High resolution mapping of Twist to DNA in Drosophila embryos: Efficient functional analysis and evolutionary conservation

Cis-regulatory modules (CRMs) function by binding sequence specific transcription factors, but the relationship between in vivo physical binding and the regulatory capacity of factor-bound DNA elements remains uncertain. We investigate this relationship for the well-studied Twist factor in Drosophila melanogaster embryos by analyzing genome-wide factor occupancy and testing the functional significance of Twist occupied regions and motifs within regions. Twist ChIP-seq data efficiently identified previously studied Twist-dependent CRMs and robustly predicted new CRM activity in transgenesis, with newly identified Twist-occupied regions supporting diverse spatiotemporal patterns (>74% positive, n = 31). Some, but not all, candidate CRMs require Twist for proper expression in the embryo. The Twist motifs most favored in genome ChIP data (in vivo) differed from those most favored by Systematic Evolution of Ligands by EXponential enrichment (SELEX) (in vitro). Furthermore, the majority of ChIP-seq signals could be parsimoniously explained by a CABVTG motif located within 50 bp of the ChIP summit and, of these, CACATG was most prevalent. Mutagenesis experiments demonstrated that different Twist E-box motif types are not fully interchangeable, suggesting that the ChIP-derived consensus (CABVTG) includes sites having distinct regulatory outputs. Further analysis of position, frequency of occurrence, and sequence conservation revealed significant enrichment and conservation of CABVTG E-box motifs near Twist ChIP-seq signal summits, preferential conservation of ±150 bp surrounding Twist occupied summits, and enrichment of GA- and CA-repeat sequences near Twist occupied summits. Our results show that high resolution in vivo occupancy data can be used to drive efficient discovery and dissection of global and local cis-regulatory logic.

Deciphering the genome's regulatory code: the many languages of DNA

The generation of patterns and the diversity of cell types in a multicellular organism require differential gene regulation. At the heart of this process are enhancers or cis-regulatory modules (CRMs), genomic regions that are bound by transcription factors (TFs) that control spatio-temporal gene expression in developmental networks. To date, only a few CRMs have been studied in detail and the underlying cis-regulatory code is not well understood. Here, we review recent progress on the genome-wide identification of CRMs with chromatin immunoprecipitation of TF-DNA complexes followed by microarrays (ChIP-on-chip). We focus on two computational approaches that have succeeded in predicting the expression pattern driven by a CRM either based on TF binding site preferences and their expression levels, or quantitative analysis of CRM occupancy by key TFs. We also discuss the current limits of these methods and highlight some of the key problems that have to be solved to gain a more complete understanding of the structure and function of CRMs.

Dynamics of the Dorsal morphogen gradient

The dorsoventral (DV) patterning of the Drosophila embryo depends on the nuclear localization gradient of Dorsal (Dl), a protein related to the mammalian NF-kappaB transcription factors. Current understanding of how the Dl gradient works has been derived from studies of its transcriptional interpretation, but the gradient itself has not been quantified. In particular, it is not known whether the Dl gradient is stable or dynamic during the DV patterning of the embryo. To address this question, we developed a mathematical model of the Dl gradient and constrained its parameters by experimental data. Based on our computational analysis, we predict that the Dl gradient is dynamic and, to a first approximation, can be described as a concentration profile with increasing amplitude and constant shape. These time-dependent properties of the Dl gradient are different from those of the Bicoid and MAPK phosphorylation gradients, which pattern the anterior and terminal regions of the embryo. Specifically, the gradient of the nuclear levels of Bicoid is stable, whereas the pattern of MAPK phosphorylation changes in both shape and amplitude. We attribute these striking differences in the dynamics of maternal morphogen gradients to the differences in the initial conditions and chemistries of the anterior, DV, and terminal systems.