Enhancer sequences influence the role of the amino-terminal domain of bicoid in transcription

Bicoid (Bcd) is a Drosophila melanogaster morphogenetic gradient that controls embryonic patterning by activating target gene expression in a concentration-dependent manner. In this study we describe experiments to determine how different enhancers respond to Bcd distinctively, focusing on two natural Bcd-responsive enhancer elements, hunchback (hb) and knirps (kni). Our results show that, on the hb enhancer element, the amino-terminal domain of Bcd (residues 1 to 91) plays primarily an inhibitory role, whereas on the kni enhancer element this same Bcd domain plays a positive role at low protein concentrations. We further demonstrate that while the amino-terminal domain is largely dispensable for cooperative binding to the hb enhancer element, it is preferentially required for cooperative binding to the kni enhancer element. Alteration of the arrangement of Bcd binding sites in the kni enhancer element reduces the role of the amino-terminal domain in cooperative DNA binding but increases the effectiveness of the self-inhibitory function. In addition, elimination of symmetric pairs of Bcd binding sites in the kni enhancer element reduces both DNA binding and activation by Bcd. We propose that the amino-terminal domain of Bcd is an enhancer-specific switch that contributes to the protein's ability to activate different target genes in distinct manners.     

A Two-Dimensional Simulation Model of the Bicoid Gradient in Drosophila

Background

Bicoid (Bcd) is a Drosophila morphogenetic protein responsible for patterning the anterior structures in embryos. Recent experimental studies have revealed important insights into the behavior of this morphogen gradient, making it necessary to develop a model that can recapitulate the biological features of the system, including its dynamic and scaling properties.

Methodology/Principal Findings

We present a biologically realistic 2-D model of the dynamics of the Bcd gradient in Drosophila embryos. This model is based on equilibrium binding of Bcd molecules to non-specific, low affinity DNA sites throughout the Drosophila genome. It considers both the diffusion media within which the Bcd gradient is formed and the dynamic and other relevant properties of bcd mRNA from which Bcd protein is produced. Our model recapitulates key features of the Bcd protein gradient observed experimentally, including its scaling properties and the stability of its nuclear concentrations during development. Our simulation model also allows us to evaluate the effects of other biological activities on Bcd gradient formation, including the dynamic redistribution of bcd mRNA in early embryos. Our simulation results suggest that, in our model, Bcd protein diffusion is important for the formation of an exponential gradient in embryos.

Conclusions/Significance

The 2-D model described in this report is a simple and versatile simulation procedure, providing a quantitative evaluation of the Bcd gradient system. Our results suggest an important role of Bcd binding to non-specific, low-affinity DNA sites in proper formation of the Bcd gradient in our model. They demonstrate that highly complex biological systems can be effectively modeled with relatively few parameters.     

Mutual repression enhances the steepness and precision of gene expression boundaries

Embryonic development is driven by spatial patterns of gene expression that determine the fate of each cell in the embryo. While gene expression is often highly erratic, embryonic development is usually exceedingly precise. In particular, gene expression boundaries are robust not only against intra-embryonic fluctuations such as noise in gene expression and protein diffusion, but also against embryo-to-embryo variations in the morphogen gradients, which provide positional information to the differentiating cells. How development is robust against intra- and inter-embryonic variations is not understood. A common motif in the gene regulation networks that control embryonic development is mutual repression between pairs of genes. To assess the role of mutual repression in the robust formation of gene expression patterns, we have performed large-scale stochastic simulations of a minimal model of two mutually repressing gap genes in Drosophila, hunchback (hb) and knirps (kni). Our model includes not only mutual repression between hb and kni, but also the stochastic and cooperative activation of hb by the anterior morphogen Bicoid (Bcd) and of kni by the posterior morphogen Caudal (Cad), as well as the diffusion of Hb and Kni between neighboring nuclei. Our analysis reveals that mutual repression can markedly increase the steepness and precision of the gap gene expression boundaries. In contrast to other mechanisms such as spatial averaging and cooperative gene activation, mutual repression thus allows for gene-expression boundaries that are both steep and precise. Moreover, mutual repression dramatically enhances their robustness against embryo-to-embryo variations in the morphogen levels. Finally, our simulations reveal that diffusion of the gap proteins plays a critical role not only in reducing the width of the gap gene expression boundaries via the mechanism of spatial averaging, but also in repairing patterning errors that could arise because of the bistability induced by mutual repression.     

High bicoid levels render the terminal system dispensable for Drosophila head development

In Drosophila, the gradient of the Bicoid (Bcd) morphogen organizes the anteroposterior axis while the ends of the embryo are patterned by the maternal terminal system. At the posterior pole, expression of terminal gap genes is mediated by the local activation of the Torso receptor tyrosine kinase (Tor). At the anterior, terminal gap genes are also activated by the Tor pathway but Bcd contributes to their activation. Here we present evidence that Tor and Bcd act independently on common target genes in an additive manner. Furthermore, we show that the terminal maternal system is not required for proper head development, since high levels of Bcd activity can functionally rescue the lack of terminal system activity at the anterior pole. This observation is consistent with a recent evolution of an anterior morphogenetic center consisting of Bcd and anterior Tor function.     

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.     

Shaping a Morphogen Gradient for Positional Precision

Morphogen gradients, which provide positional information to cells in a developing tissue, could in principle adopt any nonuniform profile. To our knowledge, how the profile of a morphogen gradient affects positional precision has not been well studied experimentally. Here, we compare the positional precision provided by the Drosophila morphogenetic protein Bicoid (Bcd) in wild-type (wt) embryos with embryos lacking an interacting cofactor. The Bcd gradient in the latter case exhibits decreased positional precision around mid-embryo compared with its wt counterpart. The domain boundary of Hunchback (Hb), a target activated by Bcd, becomes more variable in mutant embryos. By considering embryo-to-embryo, internal, and measurement fluctuations, we dissect mathematically the relevant sources of fluctuations that contribute to the error in positional information. Using this approach, we show that the defect in Hb boundary positioning in mutant embryos is directly reflective of an altered Bcd gradient profile with increasing flatness toward mid-embryo. Furthermore, we find that noise in the Bcd input signal is dominated by internal fluctuations but, due to time and spatial averaging, the spatial precision of the Hb boundary is primarily affected by embryo-to-embryo variations. Our results demonstrate that the positional information provided by the wt Bcd gradient profile is highly precise and necessary for patterning precision.     

Bicoid - morphogen function revisited

Bicoid (Bcd) functions as a morphogen during Drosophila development. Accordingly, bcd mRNA is maternally localized to the anterior pole of the embryo, and Bcd forms an anterior/posterior gradient, which functions in a concentration dependent fashion. Thus, nuclei receiving identical amounts of Bcd should express the same target genes. However, we found that ectopic, uniform expression of Bcd causes anterior gene expression in the posterior with mirror image polarity, indicating that one or several additional factors must provide positional information. Recently, we have shown that one of these factors is Capicua (Cic), a ubiquitous maternal repressor that is down-regulated at the embryonic termini by maternal Torso, a key component of the maternal terminal system. Cic acts on Bcd dependent enhancer elements by repression and thereby controls the posterior limit of Bcd target gene expression. Based on these new findings, we propose that spatial control of gene expression in the anterior region of the embryo is not solely the result of Bcd morphogen action. Rather, it relies on a "morphogenic network" that integrates the terminal system and Bcd activities, providing both polarity and spatial information to the prospective head region of the developing embryo.     

Model for the robust establishment of precise proportions in the early Drosophila embryo

During embryonic development, a spatial pattern is formed in which proportions are established precisely. As an early pattern formation step in Drosophila embryos, an anterior-posterior gradient of Bicoid (Bcd) induces hunchback (hb) expression (Nature 337 (1989) 138; Nature 332 (1988) 281). In contrast to the Bcd gradient, the Hb profile includes information about the scale of the embryo. Furthermore, the resulting hb expression pattern shows a much lower embryo-to-embryo variability than the Bcd gradient (Nature 415 (2002) 798). An additional graded posterior repressing activity could theoretically account for the observed scaling. However, we show that such a model cannot produce the observed precision in the Hb boundary, such that a fundamentally different mechanism must be at work. We describe and simulate a model that can account for the observed precise generation of the scaled Hb profile in a highly robust manner. The proposed mechanism includes Staufen (Stau), an RNA binding protein that appears essential to precision scaling (Nature 415 (2002) 798). In the model, Stau is released from both ends of the embryo and relocalizes hb RNA by increasing its mobility. This leads to an effective transport of hb away from the respective Stau sources. The balance between these opposing effects then gives rise to scaling and precision. Considering the biological importance of robust precision scaling and the simplicity of the model, the same principle may be employed more often during development.     

Morphogen gradient formation and action

In a recent publication, we identified a novel F-box protein, encoded by fates-shifted (fsd), that plays a role in targeting Bcd for ubiquitination and degradation. Our analysis of mutant Drosophila embryos suggests that Bcd protein degradation is important for proper gradient formation and developmental fate specification. Here we describe further experiments that lead to an estimate of Bcd half-life, <15 min, in embryos during the time of gradient formation. We use our findings to evaluate different models of Bcd gradient formation. With this new estimate, we simulate the Bcd gradient formation process in our own biologically realistic 2-D model. Finally, we discuss the role of Bcd-encoded positional information in controlling the positioning and precision of developmental decisions.     

Morphogen gradient formation and action: insights from studying Bicoid protein degradation

In a recent publication, we identified a novel F-box protein, encoded by fates-shifted (fsd), that plays a role in targeting Bcd for ubiquitination and degradation. Our analysis of mutant Drosophila embryos suggests that Bcd protein degradation is important for proper gradient formation and developmental fate specification. Here we describe further experiments that lead to an estimate of Bcd half-life, < 15 min, in embryos during the time of gradient formation. We use our findings to evaluate different models of Bcd gradient formation. With this new estimate, we simulate the Bcd gradient formation process in our own biologically realistic 2-D model. Finally, we discuss the role of Bcd-encoded positional information in controlling the positioning and precision of developmental decisions.