Shape and function of the Bicoid morphogen gradient in dipteran species with different sized embryos

High mobility of bicoid captured by fluorescence correlation spectroscopy: implication for the rapid establishment of its gradient

The Bicoid (Bcd) morphogen is essential for pattern formation in fruit flies. It forms an exponential concentration gradient along the embryo AP axis and turns on cascades of target genes in distinct anterior domains. The most commonly accepted model for gradient formation assumes that Bcd travels by simple diffusion and is uniformly degraded across syncytial embryos, yet several recent studies have challenged these ideas. Here, the question of Bcd mobility was investigated using fluorescence correlation spectroscopy in live Drosophila melanogaster embryos. Bcd-EGFP molecules were found to be highly mobile in the cytoplasm during cycles 12-14, with a diffusion coefficient approximately 7 microm(2)/s. This value is large enough to explain the stable establishment of the Bcd gradient simply by diffusion before cycle 8, i.e., before the onset of zygotic transcription.

Dynamic evolution of precise regulatory encodings creates the clustered site signature of enhancers

Concentration gradients of morphogenic proteins pattern the embryonic axes of Drosophila by activating different genes at different concentrations. The neurogenic ectoderm enhancers (NEEs) activate different genes at different threshold levels of the Dorsal (Dl) morphogen, which patterns the dorsal/ventral axis. NEEs share a unique arrangement of highly constrained DNA-binding sites for Dl, Twist (Twi), Snail (Sna) and Suppressor of Hairless (Su(H)), and encode the threshold variable in the precise length of DNA that separates one well-defined Dl element from a Twi element. However, NEEs also possess dense clusters of variant Dl sites. Here, we show that these increasingly variant sites are eclipsed relic elements, which were superseded by more recently evolved threshold encodings. Given the divergence in egg size during Drosophila lineage evolution, the observed characteristic clusters of divergent sites indicate a history of frequent selection for changes in threshold responses to the Dl morphogen gradient and confirm the NEE structure/function model.

The time to measure positional information: maternal hunchback is required for the synchrony of the Bicoid transcriptional response at the onset of zygotic transcription

It is widely accepted that morphogenetic gradients determine cell identity by concentration-dependent activation of target genes. How precise is each step in the gene expression process that acts downstream of morphogens, however, remains unclear. The Bicoid morphogen is a transcription factor directly activating its target genes and provides thus a simple system to address this issue in a quantitative manner. Recent studies indicate that the Bicoid gradient is precisely established in Drosophila embryos after eight nuclear divisions (cycle 9) and that target protein expression is specified five divisions later (cycle 14), with a precision that corresponds to a relative difference of Bicoid concentration of 10%. To understand how such precision was achieved, we directly analyzed nascent transcripts of the hunchback target gene at their site of synthesis. Most anterior nuclei in cycle 11 interphasic embryos exhibit efficient biallelic transcription of hunchback and this synchronous expression is specified within a 10% difference of Bicoid concentration. The fast diffusion of Bcd-EGFP (7.7 mum(2)/s) that we captured by fluorescent correlation spectroscopy in the nucleus is consistent with this robust expression at cycle 11. However, given the interruption of transcription during mitosis, it remains too slow to be consistent with precise de novo reading of Bicoid concentration at each interphase, suggesting the existence of a memorization process that recalls this information from earlier cycles. The two anterior maternal morphogens, Bicoid and Hunchback, contribute differently to this early response: whereas Bicoid provides dose-dependent positional information along the axis, maternal Hunchback is required for the synchrony of the response and is therefore likely to be involved in this memorization process.

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.

Scaling of morphogen gradients by an expansion-repression integral feedback control

Despite substantial size variations, proportions of the developing body plan are maintained with a remarkable precision. Little is known about the mechanisms that ensure this adaptation (scaling) of pattern with size. Most models of patterning by morphogen gradients do not support scaling. In contrast, we show that scaling arises naturally in a general feedback topology, in which the range of the morphogen gradient increases with the abundance of some diffusible molecule, whose production, in turn, is repressed by morphogen signaling. We term this mechanism "expansion-repression" and show that it can function within a wide range of biological scenarios. The expansion-repression scaling mechanism is analogous to an integral-feedback controller, a key concept in engineering that is likely to be instrumental also in maintaining biological homeostasis.

Precision and scaling in morphogen gradient read-out

Morphogen gradients infer cell fate as a function of cellular position. Experiments in Drosophila embryos have shown that the Bicoid (Bcd) gradient is precise and exhibits some degree of scaling. We present experimental results on the precision of Bcd target genes for embryos with a single, double or quadruple dose of bicoid demonstrating that precision is highest at mid-embryo and position dependent, rather than gene dependent. This confirms that the major contribution to precision is achieved already at the Bcd gradient formation. Modeling this dynamic process, we investigate precision for inter-embryo fluctuations in different parameters affecting gradient formation. Within our modeling framework, the observed precision can only be achieved by a transient Bcd profile. Studying different extensions of our modeling framework reveals that scaling is generally position dependent and decreases toward the posterior pole. Our measurements confirm this trend, indicating almost perfect scaling except for anterior most expression domains, which overcompensate fluctuations in embryo length.

Modelling the Drosophila embryo

I provide a historical overview on the use of mathematical models to gain insight into pattern formation during early development of the fruit fly Drosophila melanogaster. It is my intention to illustrate how the aims and methodology of modelling have changed from the early beginnings of a theoretical developmental biology in the 1960s to modern-day systems biology. I show that even early modelling attempts addressed interesting and relevant questions, which were not tractable by experimental approaches. Unfortunately, their validation was severely hampered by a lack of specificity and appropriate experimental evidence. There is a simple lesson to be learned from this: we cannot deduce general rules for pattern formation from first principles or spurious reproduction of developmental phenomena. Instead, we must infer such rules (if any) from detailed and accurate studies of specific developmental systems. To achieve this, mathematical modelling must be closely integrated with experimental approaches. I report on progress that has been made in this direction in the past few years and illustrate the kind of novel insights that can be gained from such combined approaches. These insights demonstrate the great potential (and some pitfalls) of an integrative, systems-level investigation of pattern formation.

Morphogen gradient formation

How morphogen gradients are formed in target tissues is a key question for understanding the mechanisms of morphological patterning. Here, we review different mechanisms of morphogen gradient formation from theoretical and experimental points of view. First, a simple, comprehensive overview of the underlying biophysical principles of several mechanisms of gradient formation is provided. We then discuss the advantages and limitations of different experimental approaches to gradient formation analysis.

Stable, precise, and reproducible patterning of bicoid and hunchback molecules in the early Drosophila embryo

Precise patterning of morphogen molecules and their accurate reading out are of key importance in embryonic development. Recent experiments have visualized distributions of proteins in developing embryos and shown that the gradient of concentration of Bicoid morphogen in Drosophila embryos is established rapidly after fertilization and remains stable through syncytial mitoses. This stable Bicoid gradient is read out in a precise way to distribute Hunchback with small fluctuations in each embryo and in a reproducible way, with small embryo-to-embryo fluctuation. The mechanisms of such stable, precise, and reproducible patterning through noisy cellular processes, however, still remain mysterious. To address these issues, here we develop the one- and three-dimensional stochastic models of the early Drosophila embryo. The simulated results show that the fluctuation in expression of the hunchback gene is dominated by the random arrival of Bicoid at the hunchback enhancer. Slow diffusion of Hunchback protein, however, averages out this intense fluctuation, leading to the precise patterning of distribution of Hunchback without loss of sharpness of the boundary of its distribution. The coordinated rates of diffusion and transport of input Bicoid and output Hunchback play decisive roles in suppressing fluctuations arising from the dynamical structure change in embryos and those arising from the random diffusion of molecules, and give rise to the stable, precise, and reproducible patterning of Bicoid and Hunchback distributions.

Formation of the bicoid morphogen gradient: an mRNA gradient dictates the protein gradient

The Bicoid (Bcd) protein gradient is generally believed to be established in pre-blastoderm Drosophila embryos by the diffusion of Bcd protein after translation of maternal mRNA, which serves as a strictly localized source of Bcd at the anterior pole. However, we previously published evidence that the Bcd gradient is preceded by a bcd mRNA gradient. Here, we have revisited and extended this observation by showing that the bcd mRNA and Bcd protein gradient profiles are virtually identical at all times. This confirms our previous conclusion that the Bcd gradient is produced by a bcd mRNA gradient rather than by diffusion. Based on our observation that bcd mRNA colocalizes with Staufen (Stau), we propose that the bcd mRNA gradient forms by a novel mechanism involving quasi-random active transport of a Stau-bcd mRNA complex through a nonpolar microtubular network, which confines the bcd mRNA to the cortex of the embryo.

Regulative feedback in pattern formation: towards a general relativistic theory of positional information

Positional specification by morphogen gradients is traditionally viewed as a two-step process. A gradient is formed and then interpreted, providing a spatial metric independent of the target tissue, similar to the concept of space in classical mechanics. However, the formation and interpretation of gradients are coupled, dynamic processes. We introduce a conceptual framework for positional specification in which cellular activity feeds back on positional information encoded by gradients, analogous to the feedback between mass-energy distribution and the geometry of space-time in Einstein's general theory of relativity. We discuss how such general relativistic positional information (GRPI) can guide systems-level approaches to pattern formation.

Egg size, embryonic development time and ovoviviparity in Drosophila species

Lengths, widths and volumes of eggs from 11 species of Drosophila whose genomes have been fully sequenced exhibit significant variation that is not explained by their phylogenetic relationships. Furthermore, egg size differences are unrelated to embryonic development time in these species. In addition, two of the species, Drosophila sechellia and, to a lesser degree, D. yakuba, both ecological specialists, exhibit ovoviviparity, suggesting that female control over oviposition in these species differs from what is observed in D. melanogaster. The interspecific differences in these reproductive characters, coupled with the availability of whole genome sequences for each, provide an unprecedented opportunity to examine their evolution.