- Open Access
Seeing is believing: the Bicoid protein reveals its path
Hereditas volume 155, Article number: 28 (2018)
In this commentary, I will review the latest findings on the Bicoid (Bcd) morphogen in Drosophila, a paradigm for gradient formation taught to biology students for more than two decades. “Seeing is believing” also summarizes the erroneous steps that were needed to elucidate the mechanisms of gradient formation and the path of movement of Bcd. Initially proclaimed as a dogma in 1988 and later incorporated into the SDD model where the broad diffusion of Bcd throughout the embryo was the predominant step leading to gradient formation, the SDD model was irrefutable for more than two decades until first doubts were raised in 2007 regarding the diffusion properties of Bcd associated with the SDD model. This led to re-thinking of the issue and the definition of a new model, termed the ARTS model which could explain most of the physical constraints that were inherently associated with the SDD model. In the ARTS model, gradient formation is mediated by the mRNA which is redistributed along cortical microtubules to form a mRNA gradient which is translated to form the protein gradient. Contrary to the SDD model, there is no Bcd diffusion from the tip. The ARTS model is also compatible with the observed cortical movement of Bcd. I will critically compare the SDD and the ARTS models as well as other models, analyze the major differences, and highlight the path where Bcd is localized during early nuclear cycles.
History of mechanisms to explain the occurrence of the Bicoid morphogen gradient
Proclamation of the SDD model as a dogma in the 80’s
The Bicoid (Bcd) protein and its gradient (Fig. 1f) is one of the fascinating observations in nature. Discovered in the fruit fly, Drosophila during the late 80’s [1, 2], it was illustrated in textbooks as a paradigm for morphogen gradient formation. Not only was the morphogen gradient beautiful by appearance, it was also remarkable how the embryo managed to generate an anterior-posterior (A-P) gradient and a coordinate system based on information which initially is stored as a point source as maternal mRNA (Fig. 1a). Within 3 h of development, a precise gradient is generated revealing an exponential decay (Fig. 1f) which strictly follows mathematical rules . bcd resides at the top of the hierarchy of segmentation genes  to pattern the A-P axis and instructs the gap genes  via its concentration gradient from the anterior side. The gap genes, in turn, are expressed in broader domains along the A-P axis and control the so-called pair-rule genes which are usually expressed in 7 stripes . These, in turn, control the segment polarity genes (also referred to segmentation genes)  that further subdivide the embryo into smaller units. Finally, homeotic genes are thought the maintain the status of the established segments .
The SDD model (S standing for “synthesis”, D for “diffusion” and D for uniform “degradation”) was initially used to explain the establishment of the gradient [2, 9]. The model proposed that the bcd mRNA stays strictly at the tip during all developmental stages (Fig. 1d) serving as template for the production of Bcd which itself diffuses throughout the embryo in all directions (Fig. 1e), followed by uniform degradation. The nuclei do not contribute to the shape of the gradient, but rather function as a tool to interpret the gradient .
Attempts to validate the model were inconclusive for a long time and it took almost two decades to realize that the diffusion constant of Bcd was far too low to move to the posterior  and that SDD model might contain conceptual errors. One of the major drawbacks of the SSD model was its presumption that the Bcd protein diffused throughout the embryo (Fig. 1e). This was fueled by studies where fluorescently labeled dextrane particles injected at the anterior pole were used to simulate the diffusion of Bcd . In retrospect, it was a daring proposal, however, the constraints of this approach were clear from the very beginning. Nevertheless, the approach was too simplistic to assume that a protein would behave like a dextrane particle. Subsequently, other reports measured higher diffusion rates [13,14,15], calculated to be high enough to explain the SDD diffusion model, corroborated by a recent biophysical model analysis .
The ARTS model as a solution for most physical constraints
In 2009, a new model for Bcd gradient formation was proposed solving most of the physical constraints associated with the apparent slow diffusion constant of Bcd: the ARTS model [17,18,19,20,21], standing for “active RNA transport and synthesis”. The ARTS model incorporates the formation of a bcd mRNA gradient and involves active transport of the mRNA by means of microtubules forming a mRNA gradient first (Fig. 1b), followed by synthesis of the protein based on the mRNA gradient (Fig. 1c). One report showed that the mRNA briefly enters the yolk at nuclear cycle 4 . In summary, while the start and end points of the two proposed mechanisms are identical (Fig. 1a, f), it is the mechanism in between where the ARTS model differs completely from that of the SDD model .
In 2007, a model was proposed to involve Bcd diffusion combined with nucleocytoplasmic shuttling, but no Bcd degradation . The nuclei would serve as reversible traps that affect and slow down Bcd diffusion, while their increase in number with time would counteract the diffusive spread of Bcd. Notably, the Bcd gradient was predicted to be established before the nuclei migrate to the periphery, i. e. before nc 10. Moreover, the gradient was proposed to remain stationary during the remaining 4 nuclear cycles up to nc 14 when cellularization was reached.
Analysis of Bicoid movement in the early Drosophila embryo
To precisely monitor the path of Bcd movement during early development, a sensitive approach was developed that allowed for the study of the spatial Bcd movement during the early nuclear cycles using single confocal sections . This study revealed, for the first time, that Bcd moved at the cortex of the egg but never entered the inner portion filled with yolk (Fig. 2a, b), a feature that could not be observed in the past because live imaging of Bcd was largely impossible due to the elevated background from the yolk. Notably, the inner part of the egg, i. e. the yolk acted as a non-permissive territory. The results from this study immediately refuted the SDD model which predicted that Bcd would move throughout the yolk. If the egg was exposed to hypoxia, i. e. by water submersion, it fell into “sleep” and development was halted. However, development resumed when oxygen was restored. The hypoxia technique permitted the monitoring of the movement of Bcd under prolonged sleeping conditions, i. e. for 3 h (Fig. 2c) or 7 h (Fig. 2d). During this sleeping phase, Bcd still moved slowly at the cortex to the posterior. Again, the speed and location of Bcd movement was not compatible with the diffusion properties claimed by the SDD model , nor with that from the nucleocytoplasmic shuttle model . Notably, the inner part of the egg, i. e. the yolk still acted as a non-permissive territory (Fig. 2c, d; . This finding was unexpected, based on the previous knowledge of the diffusion of fluorescent dextrane particles which easily entered the yolk . Data published concurrently by  reiterated the notion that the mechanism for Bcd movement was presumably far more complex than previously anticipated. Experiments demonstrated that embryos exposed to smaller drugs directed against major cytoarchitectural proteins such as microtubules (MTs) or actin, also affected Bcd movement (Fig. 2e-f; . Intact MTs appeared to be indispensable for maintaining a non-permissive territory (Fig. 2e). If actin was compromised, Bcd movement as well as its stability was affected leading to a substantially-altered Bcd pattern (Fig. 2f). This data suggested that actin has a dual function for Bcd.
State of the art
The SDD model is possible under artificial conditions
The drug treatment data of  revealed that the inner part of the egg, i. e. the yolk became permissive for Bcd if embryos were bathed in substances affecting the cytoarchitecture. In vinblastine-treated embryos, the Bcd protein movement behaved exactly as the SDD model would have predicted, i. e. it moved to the posterior in a broad front (Fig. 2e) and the yolk became permissive for Bcd. However, these are artificial conditions since embryos are not exposed to vinblastine in nature. This observation also revealed several apparent weaknesses of the SDD model that were never before discussed. If the model was correct, why should the embryo translate a protein at the tip while the majority remains in the interior of the embryo, never reaching the blastoderm nuclei? Secondly, how would an insect egg three times the size of Drosophila, e. g. that of the blow fly Lucilia sericata [18, 23] solve the problem of Bcd movement through the yolk where the expansive distance in this large embryo creates further physical constraints? Thirdly, why would nature choose such a difficult path from the tip through the yolk, and then back to the cortex to enter the blastoderm nuclei?
Why does Bcd not enter the yolk?
Data on the content and structure of the inner yolk is limited, primarily due to microscopic constraints, because the laser of a confocal microscope cannot penetrate deeply into the optically dense yolk layer. Actin microfilaments in the inner yolk were described [24, 25], but require a more detailed description. For MTs, only the spindle apparatus during the nuclear cycles were bright enough to become visible . Further, attempts to stain an internal MT network for visualizing axial expansion and cortical migration had apparently failed [26, 27]. It is possible that the appearance of a MT-network involved in axial expansion and cortical expansion is nuclear cycle-dependent and consequently may be visible only for a fraction of the cycle, as was the case for the cortical MT network transporting the bcd mRNA along the cortex . It is plausible to assume that the yolk contains a plethora of cytoskeletal elements that so far have escaped detection. For example,  described the existence of long fibrous materials, presumably of MT-origin, that were observed when embryos were squashed under certain salt conditions and the yolk content examined.
Existence of different Bcd isoforms
Currently, little data is known regarding the alternative splicing of the bcd gene resulting in 5 different isoforms, some of which have been characterized [29, 30]. Of these, isoform A  represents a small homeobox-less isoform, 149 amino acids (aa) in size which is expressed at vanishingly low levels  and therefore was likely undetected by the assay of . The four other isoforms differ in differential use of a splice-acceptor site at exon 3 which is 15 nucleotides apart giving rise to a protein with 5 extra amino acids (+ 5 aa), combined with the alternative use of a translation start point further downstream. The combination of these events allows for the creation of 4 larger Bcd isoforms of 413, 418, 489 and 494 aa, respectively . Interestingly, the + 5 aa form is expressed as strongly as the commonly-used Bcd isoform lacking these 5 amino acids (unpublished data). However, given the prevalence of alternative splicing immediately upstream of the homeodomain, it was striking to learn that this splicing event was not found to occur in a close relative of D. melanogaster, D. pseudoobscura , that gave rise to the − 5 aa isoform only. All studies using a bcd cDNA in the past were done based on the original c53.46.6 bcd cDNA clone  encoding the − 5 aa form. This form allowed rescue of a bcd-null phenotype , however, it was unknown to what extent any changes of the anlagen under the rescue conditions were corrected, in analogy to the active fate map repair system in bcd+ 5 + 8 embryos that allowed the rescue of embryos with severely altered Even-skipped stripe patterns and thus substantially changed segmental anlagen [20, 34]. A likely scenario could be that the Bcd protein movement as well as its function could be largely dependent on isoforms, a question which has not been addressed. Furthermore, the smallest isoform with 149 aa likely has the capacity to travel faster to the posterior and may be the factor responsible for suppressing pole cell formation as demonstrated in . Consistent with this observation is data that over-expression of the smallest Bcd isoform during oogenesis alters the segmental anlagen in the posterior and concomitantly suppresses pole cell-formation (unpublished data). The smallest Bcd isoform may also be responsible for the surprisingly robust signal resulting from hub activity of Bcd observed at the posterior end , despite the fact that only low concentrations of Bcd were measured at the posterior pole . Unfortunately, the experimental setup of  did not allow to discriminate between the different Bcd isoforms to give a hint which isoform was actually detected. To resolve the issue of Bcd isoform movement and function, sensitive isoform-specific antibodies would be required, experiments that would be technically challenging, but not impossible.
Is active transport of Bcd possible?
This question was discussed by  in an effort to explain the low diffusion constant of Bcd. An interaction database of FlyBase  lists a number of Bcd-interacting proteins such as Bin 1 and Bin 3 [37, 38], Chip , Dampened , elF4EHP  and Fate-shifted . Of these proteins, none were directly associated with a motor protein, thus, any interaction partner may require further proteins linking Bcd to a motor machinery which so far has not been detected. In my view, the strongest argument against active long-range transport of Bcd is the existence of the bcd mRNA gradient as the template for the protein gradient [1, 17]. Moreover, the establishment of the mRNA gradient is fast enough to allow for the formation of the Bcd gradient which makes the involvement of active transport of Bcd largely superfluous.
Why does it make sense to observe Bcd residing and moving along the cortex and not moving towards the interior? Arguably, the most important finding is the cortical location of the mRNA (Fig. 1b) and the definition of the yolk and the cortex as territories which have an impact on the movement of Bcd. From an energy point of view, it would be too costly to establish a path associated with the yolk first which would pose further problems on controlling the movement of Bcd. These are the main reasons why today I believe in what I see.
Active mRNA transport, synthesis
Synthesis, diffusion, degradation
Frigerio G, Burri M, Bopp D, Baumgartner S, Noll M. Structure of the segmentation gene paired and the Drosophila PRD gene set as part of a gene network. Cell. 1986;47(5):735–46.
Driever W, Nüsslein-Volhard C. A gradient of bicoid protein in Drosophila embryos. Cell. 1988;54(1):83–93.
Houchmandzadeh B, Wieschaus E, Leibler S. Establishment of developmental precision and proportions in the early Drosophila embryo. Nature. 2002;415(6873):798–802.
Nüsslein-Volhard C, Wieschaus E. Mutations affecting segment number and polarity in Drosophila. Nature. 1980;287(5785):795–801.
Crombach A, Wotton KR, Jimenez-Guri E, Jaeger J. Gap gene regulatory dynamics evolve along a genotype network. Mol Biol Evol. 2016;33(5):1293–307.
Clark E. Dynamic patterning by the Drosophila pair-rule network reconciles long-germ and short-germ segmentation. PLoS Biol. 2017;15(9):e2002439.
Ingham PW. Drosophila segment polarity mutants and the rediscovery of the hedgehog pathway genes. Curr Top Dev Biol. 2016;116:477–88.
Pick L. Hox genes, evo-devo, and the case of the ftz gene. Chromosoma. 2016;125(3):535–51.
Driever W, Nüsslein-Volhard C. The bicoid protein determines position in the Drosophila embryo in a concentration-dependent manner. Cell. 1988;54(1):95–104.
Grimm O, Wieschaus E. The Bicoid gradient is shaped independently of nuclei. Development. 2010;137(17):2857–62.
Gregor T, Wieschaus EF, McGregor AP, Bialek W, Tank DW. Stability and nuclear dynamics of the bicoid morphogen gradient. Cell. 2007;130(1):141–52.
Gregor T, Bialek W, de Ruyter van Steveninck RR, Tank DW, Wieschaus EF. Diffusion and scaling during early embryonic pattern formation. Proc Natl Acad Sci U S A. 2005;102(51):18403–7.
Castle BT, Howard SA, Odde DJ. Assessment of transport mechanisms underlying the Bicoid Morphogen gradient. Cell Mol Bioeng. 2011;4(1):116–21.
Abu-Arish A, Porcher A, Czerwonka A, Dostatni N, Fradin C. High mobility of bicoid captured by fluorescence correlation spectroscopy: implication for the rapid establishment of its gradient. Biophys J. 2010;99(4):L33–5.
Durrieu L, Kirrmaier D, Schneidt T, Kats I, Raghavan S, Hufnagel L, Saunders TE, Knop M. Bicoid gradient formation mechanism and dynamics revealed by protein lifetime analysis. Mol Syst Biol. 2018;14:e8355.
Sigaut L, Pearson JE, Colman-Lerner A, Ponce Dawson S. Messages do diffuse faster than messengers: reconciling disparate estimates of the morphogen bicoid diffusion coefficient. PLoS Comput Biol. 2014;10(6):e1003629.
Spirov A, Fahmy K, Schneider M, Frei E, Noll M, Baumgartner S. Formation of the bicoid morphogen gradient: an mRNA gradient dictates the protein gradient. Development. 2009;136(4):605–14.
Fahmy K, Akber M, Cai X, Koul A, Hayder A, Baumgartner S. alphaTubulin 67C and Ncd are essential for establishing a cortical microtubular network and formation of the Bicoid mRNA gradient in Drosophila. PLoS One. 2014;9(11):e112053.
Ali-Murthy Z, Kornberg TB. Bicoid gradient formation and function in the Drosophila pre-syncytial blastoderm. Elife. 2016;5. https://doi.org/10.7554/eLife.13222.
Cai X, Akber M, Spirov A, Baumgartner S. Cortical movement of Bicoid in early Drosophila embryos is actin- and microtubule-dependent and disagrees with the SDD diffusion model. PLoS One. 2017;12(10):e0185443.
Lipshitz HD. Follow the mRNA: a new model for Bicoid gradient formation. Nat Rev Mol Cell Biol. 2009;10(8):509–12.
Coppey M, Berezhkovskii AM, Kim Y, Boettiger AN, Shvartsman SY. Modeling the bicoid gradient: diffusion and reversible nuclear trapping of a stable protein. Dev Biol. 2007;312(2):623–30.
Gregor T, McGregor AP, Wieschaus EF. Shape and function of the Bicoid morphogen gradient in dipteran species with different sized embryos. Dev Biol. 2008;316(2):350–8.
Warn RM, Magrath R, Webb S. Distribution of F-actin during cleavage of the Drosophila syncytial blastoderm. J Cell Biol. 1984;98(1):156–62.
Hatanaka K, Okada M. Retarded nuclear migration in Drosophila embryos with aberrant F-actin reorganization caused by maternal mutations and by cytochalasin treatment. Development. 1991;111(4):909–20.
Zalokar M, Erk I. Division and migration of nuclei during early embryogenesis of Drosophila melanogaster. J Microscopie Biol Cell. 1976;25:97–106.
Baker J, Theurkauf WE, Schubiger G. Dynamic changes in microtubule configuration correlate with nuclear migration in the preblastoderm Drosophila embryo. J Cell Biol. 1993;122(1):113–21.
Foe VE, Alberts BM. Studies of nuclear and cytoplasmic behaviour during the five mitotic cycles that precede gastrulation in Drosophila embryogenesis. J Cell Sci. 1983;61:31–70.
Rodel CJ, Gilles AF, Averof M. MicroRNAs act as cofactors in bicoid-mediated translational repression. Curr Biol. 2013;23(16):1579–84.
Berleth T, Burri M, Thoma G, Bopp D, Richstein S, Frigerio G, Noll M, Nüsslein-Volhard C. The role of localization of bicoid RNA in organizing the anterior pattern of the Drosophila embryo. EMBO J. 1988;7(6):1749–56.
Attrill H, Falls K, Goodman JL, Millburn GH, Antonazzo G, Rey AJ, Marygold SJ, FlyBase C. FlyBase: establishing a Gene Group resource for Drosophila melanogaster. Nucleic Acids Res. 2016;44(D1):D786–92.
Seeger MA, Kaufman TC. Molecular analysis of the bicoid gene from Drosophila pseudoobscura: identification of conserved domains within coding and noncoding regions of the bicoid mRNA. EMBO J. 1990;9(9):2977–87.
Driever W, Ma J, Nusslein-Volhard C, Ptashne M. Rescue of bicoid mutant Drosophila embryos by bicoid fusion proteins containing heterologous activating sequences. Nature. 1989;342(6246):149–54.
Namba R, Pazdera TM, Cerrone RL, Minden JS. Drosophila embryonic pattern repair: how embryos respond to bicoid dosage alteration. Development. 1997;124(7):1393–403.
Mir M, Reimer A, Haines JE, Li XY, Stadler M, Garcia H, Eisen MB, Darzacq X. Dense Bicoid hubs accentuate binding along the morphogen gradient. Genes Dev. 2017;31(17):1784–94.
Morrison AH, Scheeler M, Dubuis J, Gregor T. Quantifying the Bicoid morphogen gradient in living fly embryos. Cold Spring Harb Protoc. 2012;2012(4):398–406.
Zhu W, Hanes SD. Identification of drosophila bicoid-interacting proteins using a custom two-hybrid selection. Gene. 2000;245(2):329–39.
Singh N, Zhu W, Hanes SD. Sap18 is required for the maternal gene bicoid to direct anterior patterning in Drosophila melanogaster. Dev Biol. 2005;278(1):242–54.
Torigoi E, Bennani-Baiti IM, Rosen C, Gonzalez K, Morcillo P, Ptashne M, Dorsett D. Chip interacts with diverse homeodomain proteins and potentiates bicoid activity in vivo. Proc Natl Acad Sci U S A. 2000;97(6):2686–91.
Liu J, Ma J. Dampened regulates the activating potency of Bicoid and the embryonic patterning outcome in Drosophila. Nat Commun. 2013;4:2968.
Cho PF, Poulin F, Cho-Park YA, Cho-Park IB, Chicoine JD, Lasko P, Sonenberg N. A new paradigm for translational control: inhibition via 5′-3' mRNA tethering by Bicoid and the eIF4E cognate 4EHP. Cell. 2005;121(3):411–23.
Liu J, Ma J. Fates-shifted is an F-box protein that targets Bicoid for degradation and regulates developmental fate determination in Drosophila embryos. Nat Cell Biol. 2011;13(1):22–9.
S. B. thanks the Swedish Research Council, the Ekhaga, Nilsson-Ehle and Erik Philip-Sörensen Stiftelse for support, and Linda Wei and Thomas Blom for comments on the manuscript.
Swedish Research Council, award number 2010–4358.
Erik Philip-Sörensen Foundation.
Nilsson Ehle Foundation.
Availability of data and materials
The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.
Ethics approval and consent to participate
Consent for publication
The author declares that he has no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.