Distance measurements via the morphogen gradient of Bicoid in Drosophila embryos
© He et al; licensee BioMed Central Ltd. 2010
Received: 13 April 2010
Accepted: 2 August 2010
Published: 2 August 2010
Patterning along the anterior-posterior (A-P) axis in Drosophila embryos is instructed by the morphogen gradient of Bicoid (Bcd). Despite extensive studies of this morphogen, how embryo geometry may affect gradient formation and target responses has not been investigated experimentally.
In this report, we systematically compare the Bcd gradient profiles and its target expression patterns on the dorsal and ventral sides of the embryo. Our results support a hypothesis that proper distance measurement and the encoded positional information of the Bcd gradient are along the perimeter of the embryo. Our results also reveal that the dorsal and ventral sides of the embryo have a fundamentally similar relationship between Bcd and its target Hunchback (Hb), suggesting that Hb expression properties on the two sides of the embryo can be directly traced to Bcd gradient properties. Our 3-D simulation studies show that a curvature difference between the two sides of an embryo is sufficient to generate Bcd gradient properties that are consistent with experimental observations.
The findings described in this report provide a first quantitative, experimental evaluation of embryo geometry on Bcd gradient formation and target responses. They demonstrate that the physical features of an embryo, such as its shape, are integral to how pattern is formed.
A cornerstone in developmental biology is the concept of morphogens. These are chemical substances that form concentration gradients to provide positional information to a developing embryo or tissue [1–5]. Despite significant advances in the understanding of morphogens in recent years, how the geometry of an embryo or tissue may affect morphogen gradient formation has not been well investigated experimentally. This question, like other quantitative questions about morphogens, requires an accurate measurement of the gradient concentrations, a task that poses significant technical and analytical challenges [6–8]. Recent efforts, including live-imaging studies  and our own embryo staining studies [8, 9] of the gradient of the Drosophila morphogen Bicoid (Bcd), have overcome such challenges and revealed important new insights into the behavior of this morphogen gradient (see also ). Here we extend our studies to investigate the effects of embryo geometry on Bcd gradient formation and the interpretation of its encoded positional information.
Bcd, a morphogen which forms a gradient along the anterior-posterior (A-P) axis in Drosophila embryos, is responsible for patterning the anterior structures, including the head [11, 12]. It is a transcription factor that directly activates the expression of its downstream target genes in a concentration-dependent manner [13–15]. Bcd gradient formation takes place before the physical boundaries of cells (i.e., cell membranes) are formed in the embryo. It is a process  thought to involve localized synthesis of Bcd protein molecules in the anterior coupled with diffusion and degradation of these molecules throughout the embryo (see also ). At the syncytial blastoderm stage when Bcd begins to establish and refine its concentration gradient and activate its target genes, the embryo already possesses two distinct diffusion media for Bcd molecules: the inner part (yolk) of the embryo and the nuclear layer on the cortex [18, 19]. Since the cortex of the embryo is a curved surface, embryo geometry may affect how the Bcd gradient is formed and, as a consequence, how its target genes are activated during development.
In this study, we systematically analyze the behaviors of both the Bcd gradient and its target responses on the dorsal (D) and ventral (V) sides of the embryo to investigate the effects of embryo geometry on morphogen gradient formation and action. We sought to address the question, through quantitative studies of Bcd and its direct gap gene targets, of how segmentation gene expression patterns exhibit differences between the dorsal and ventral sides. Previous studies show that the expression boundaries of gap genes and expression stripes of pair-rule genes are not parallel to the D-V axis, but rather present a splay on the ventral side of the embryo [20–23]. How such a splay is established remains controversial because several models, ranging from an interaction between the A-P and D-V patterning systems  to embryo geometry  and nuclear density/movements , have been proposed. Our experimental results reveal important insights into how distance is properly measured and positional information is encoded by the Bcd gradient in embryos. We also present 3-D simulation results demonstrating that a curvature difference between the two sides of an embryo is sufficient to generate Bcd gradient properties that are consistent with experimental observations. Our systematic analyses of a direct Bcd target, Hunchback (Hb), show that both dorsal and ventral sides of the embryo are fundamentally similar to each other in terms of Hb response to the Bcd gradient. These results suggest an important role of embryo geometry, operating through its effects on the Bcd gradient, in establishing differences of segmentation gene expression patterns between the dorsal and ventral sides of the embryo. Together, our findings demonstrate that the physical features of an embryo, such as its shape, can directly affect the outcome of the patterning process.
Embryo shape and design of study
In our current study, we take advantage of the geometric differences between the dorsal and ventral sides of the embryo to facilitate our investigation of how the Bcd gradient is formed and how its direct gap gene targets are activated during development. We analyze Bcd gradient profiles and its target responses using A-P positions that are defined in two distinct ways: projected distance x from the anterior pole or contour distance c from the identical anterior pole (Fig. 1B). Technically, while x is measured by projecting the coordinate of the center of an intensity scanning window onto the A-P axis, contour distance c is measured as the route along the perimeter of either the dorsal or ventral side of the embryo (cD and cV, respectively; see Methods for details on c determination). In Fig. 1C, we plot the average difference between the two contour distances, Δcembryo/CD (Δcembryo = cD - cV), against fractional embryo length (i.e., the normalized projected distance) x/L in wt embryos. This plot reveals the following information regarding the geometric differences between the dorsal and ventral sides of the embryo: 1) cV is longer than cD at same x/L positions in most parts of the embryo, except near the anterior pole (x/L < ~0.12); 2) there are two broad regions of the embryo, ~0.12 < x/L < ~0.36 and ~0.66 < x/L < ~0.96, where cV expands/accumulates its lead over cD; 3) in the mid-embryo section along the A-P axis (~0.36 < x/L < ~0.66), Δcembryo is stable, reflective of a small curvature difference between the two sides of the embryo (see additional file 1 for a theoretical 2-D consideration of geometry and positional information). These observed differences between the dorsal and ventral sides represent important tools to facilitate our study of Bcd gradient formation and target gene activation during development. Here we compare the behaviors of the Bcd gradient profiles and its target responses on both dorsal and ventral sides of the embryo. An important premise of our work is that the observed Bcd gradient properties can directly shed light on how distance is measured and positional information encoded during development. Thus, if the positional information encoded by the Bcd gradient is based on the projected distance x in embryos, the profiles from the dorsal and ventral sides are expected to be more similar to each other as a function of x than as a function of c. The converse is expected to be true if the encoded positional information is based on the contour distance.
Evaluating two different distance measurements for the Bcd gradient
In our analysis, embryos at early nuclear cycle 14 were triple-stained with anti-Bcd and anti-Hb antibodies as well as DAPI (see Methods for details). We captured high-resolution digital images on the midsagittal plane of the embryos and measured fluorescence intensities along the cortical, nuclear layer (Fig. 1B). As further discussed in our recent work , we used raw Bcd intensities B in our analysis, with background subtracted when necessary (see Methods for details); our experimental conditions yielded Bcd intensity profiles with a reproducibility that is comparable to live-imaging studies [6, 8]. The high precision in Bcd intensity measurements and the use of raw Bcd intensities make the analyses described in this report possible. As discussed in detail previously , normalization of Bcd intensities can mask or distort the system properties. Using a pixel-by-pixel scanning approach (see Methods for details), we recorded both the intensity values and the x-y coordinates in microns in order to quantify the intensity profiles using two distinct distance measurement methods: x and c as defined above.
Measurements of Hb expression boundary positions using two methods
To determine whether the positional information encoded by the Bcd profiles on the dorsal and ventral sides of the wt embryos are properly interpreted by its downstream target responses, we measured Hb expression profiles in wild-type (wt) embryos. To ensure a direct evaluation of the relationship between Bcd and Hb, we extracted Bcd and Hb intensities in wt embryos that were simultaneously stained with both Bcd and Hb antibodies (see Methods for details). Again, we analyze Hb expression profiles in two distinct ways: as a function of either projected distance x from the anterior or contour distance c. This allows us to measure the Hb expression boundary position (the location where Hb intensity level is half maximal) as either projected distance xHb from the anterior or contour distance cHb. Fig. 2C and 2D show the average Hb intensity profiles from the dorsal and ventral sides as a function of either x or c, respectively. While these two Hb profiles are visually different from each other when measured as a function of x (Fig. 2C), they converge near mid-embryo when measured as a function of c (Fig. 2D). Quantitatively, while the xHb values on the dorsal and ventral sides of wt embryos (xHbD and xHbV, respectively) are significantly different from each other (ΔxHb = xHbD - xHbV = 248.7 - 232.2 = 16.5 μm, p = 10-8; Fig. 2C), the difference between the cHb values on these two sides of the embryos becomes much smaller (ΔcHb = cHbD - cHbV = 284.6 - 279.1 = 5.5 μm, p = 0.07; Fig. 2D). The observed value of ΔcHb, but not ΔxHb, is comparable to the diameter of a single nucleus (~6.5 μm ) or the embryo-to-embryo variation of the Hb expression boundary position (0.014 L or 7.8 μm ). These results show that Hb, as a biological readout of the Bcd gradient, exhibits expression properties that are consistent with those of the Bcd gradient in terms of how distance is measured (i.e., positional information encoded) in embryos. Together, these results suggest that both Bcd profiles and its encoded positional information (as monitored by Hb expression boundary positions) follow a function of contour distance c from the anterior, as opposed to projected distance x.
Input-output relationship on dorsal and ventral sides of the embryo
where n is the Hill coefficient, which defines the steepness of the Hb on/off switch in response to the Bcd gradient (in Fig. 3B the green profile is, for comparative purposes, the Hill curve of n = 5, a value comparable to what is observed experimentally in DNA binding studies ). Since these two parameters provide quantitative measures of the Bcd-Hb relationship, we compared their values calculated from both sides of individual embryos. Our results show that neither Kd values (expressed as arbitrary intensity units) on the dorsal and ventral sides of the embryo (5.0 ± 1.2 and 4.8 ± 1.2, respectively; p = 0.53) nor their n values (5.1 ± 2.7 and 4.9 ± 2.7; p = 0.73) exhibit a significant difference. Finally, to evaluate the precision of the Hb readout machinery in response to the Bcd gradient, we calculated the normalized standard deviations of Hb expression (δH/Hmax) on both sides of the embryo as a function of Bcd concentration (B/Kd). Our results (Fig. 3C) reveal a globally similar profile of Hb expression noise in response to Bcd on both sides of the embryo. Together, these results demonstrate a fundamentally similar Bcd-Hb relationship on the dorsal and ventral sides of the embryo.
Analyzing precision of the Bcd gradient on dorsal and ventral sides
To further evaluate the Bcd and Hb behaviors on the dorsal and ventral sides of the embryo, we compared their Bcd intensity noise δB/B (standard deviation normalized by the mean) profiles. Previous studies including our own focused primarily on the reproducibility of the Bcd gradient on the dorsal or lateral side of the embryos [6, 8, 30–32]. In all those studies, Bcd intensity profiles were measured as a function of projected distance x. Here we measure them as a function of c and directly compare the measured Bcd intensity noise on dorsal and ventral sides of the embryo. Our results show that the ventral side has a higher δB/B around cHb than the dorsal side (Fig. 3D). Although the precise mechanism(s) for this intensity noise difference remains to be elucidated, here we simply take advantage of this difference to further help us evaluate the Bcd-Hb relationship. To directly compare the precision of positional information provided by the Bcd gradient, we converted Bcd intensity errors δB (standard deviation) to its positional errors σc . Our results reveal a higher σc on the ventral side (39 ± 1.2 μm at cHb; standard deviation calculated by bootstrapping) than the dorsal side (27 ± 0.6 μm at cHb). If Hb expression properties are dictated by Bcd gradient properties as suggested above, Hb expression is expected to be more variable on the ventral side than on the dorsal side. To specifically test this prediction, we compared the precision of cHb on the two sides of the embryos. Our results show that the variability of Hb boundary on the ventral side (σcHb = 11.3 ± 1.3 μm) is higher than that on the dorsal side (σcHb = 9.1 ± 0.9 μm). These results provide further support to our conclusion that the dorsal and ventral sides of the embryo are fundamentally the same in regard to the Bcd-Hb relationship.
Contributions of Bcd gradient profiles to dorsal-ventral differences
Our results shown in Fig. 4A make four specific, experimentally testable predictions if the Bcd profiles are responsible for the expression slanting of its direct gap gene targets in embryos. First, as seen with the Hb expression boundary positions in wt embryos (Fig. 2C), targets at other boundary positions along the A-P axis (except the most anterior part of the embryo) should also exhibit a slanting effect, i.e., ΔxTarget = xTargetD - xTargetV should be a positive value. Second, while targets with a boundary position closer to mid-embryo should have a smaller ΔxTarget thus with a subdued slanting, those with a more anterior boundary position (expect near the anterior pole) should have a larger ΔxTarget with an enhanced slanting. Third, the distribution of the lines connecting the boundary positions on the two sides of the embryo should exhibit a global similarity to that of Bcd profiles shown in Fig. 4B. Finally, as seen in Fig. 4A where ΔxBcd is larger than ΔcBcd except at the highest Bcd concentrations (which correspond to locations near the anterior pole of the embryo), ΔxTarget should be larger than ΔcTarget.
To directly evaluate these predictions, we measured xHb in embryos that have either reduced (1×-bcd) or increased (3×-bcd) Bcd concentrations, thus shifting the xHb position toward the anterior and posterior, respectively. We also analyzed the expression boundary of another direct Bcd target, Orthodenticle (Otd), in wt embryos. Similar to our analysis of Hb expression in wt embryos, we measured the expression boundary positions on both dorsal and ventral sides. Our results (Fig. 4C) show that, consistent with our predictions outlined above, ΔxTarget is a positive value in all the cases tested. Furthermore, this value is highest for Otd, which has the most anterior expression boundary position among those tested, and lowest for xHb in 3×-bcd embryos, which has the expression boundary position closest to the center of the embryo. The line distributions in Fig. 4D and Fig. 4B also exhibit a global similarity where data are available. Finally, as predicted and shown in Fig. 4C, ΔcTarget (red) is smaller than ΔxTarget (blue) in all cases tested (see additional file 1 for a theoretical discussion relevant to this issue). Together, these results show that the differences in Bcd gradient profiles on the dorsal and ventral sides of the embryos are responsible for the observed slanting of its direct gap gene targets.
3-D simulation of an embryo with asymmetric curvatures
We now return to the hypothesis that the Bcd gradient is formed in a process where contour distance from the anterior is measured, as suggested by our observed Bcd gradient properties in wt embryos at early nuclear cycle 14 (see also additional files 2, 3 and 4 for results in 1×-bcd and 3×-bcd embryos). Bcd gradient formation is a complex process, which takes place during a developmental period when the embryo itself is undergoing an exponential growth in terms of the nuclear number. There are also migrations of the nuclei to the cortex of the embryo at nuclear cycle 10, leading to the formation of the cortical (nuclear) layer where Bcd activates transcription of its downstream target genes [6, 19]. How distance is measured and positional information encoded by the Bcd gradient in the embryo is governed by physical properties of both Bcd molecules and the embryo acting as the Bcd diffusion medium. We recently developed a biologically realistic 2-D model to simulate the dynamics of Bcd gradient formation in embryos . This model was based on the equilibrium binding of Bcd molecules to low-affinity, non-specific DNA bindings sites throughout the Drosophila genome. Here we expand this 2-D model to a new 3-D model (see Methods for further details).
In this study we systematically analyze the Bcd and Hb behaviors on both the dorsal and ventral sides of the embryo. We reasoned that, since these two sides of the embryos have distinct geometric features (Fig. 1, see also additional file 1 for a theoretical discussion), a direct comparison between the Bcd and Hb profiles on these two sides should shed new light on how embryo geometry may affect both Bcd gradient formation and target gene activation during development. In our analyses we use two distinct methods to measure positions along the A-P axis: the projected distance x from the anterior pole and the contour distance c along the perimeters of the dorsal or ventral sides of the embryo. Our findings support a hypothesis that the Bcd gradient profile establishes and exerts the positional information by following the contour distance c. Our simulation studies (Fig. 5) show that a curvature difference as seen on the dorsal and ventral sides of the embryo is sufficient to generate Bcd gradient properties consistent with experimental data. Mechanistically, our findings underscore the importance of the nuclear layer on the cortex of a Drosophila syncytial blastoderm as a critical medium for Bcd gradient formation. However, as further discussed below, Bcd gradient formation cannot be viewed as a simple 1-D process in the cortical layer as an isolated diffusion medium, suggesting that an interaction (exchange of Bicoid molecules) between the cortical layer and the inner part of the embryo is important in proper Bcd gradient formation.
Our quantitative studies presented in this report show that, at early nuclear cycle 14, the dorsal and ventral sides of the embryo have a fundamentally similar input-output relationship between Bcd and Hb (see additional file 9 for experimentally measured parameter values). These results provide important insights into the well-documented splay of the segmentation gene expression stripes toward the poles on the ventral side [20–23]. Previous studies did not have a quantitative consideration of both the embryo geometry and morphogen profiles and, thus, were not able to explain how such a splay is achieved mechanistically. Our current work probes, in a quantitative way, the earliest decisions an embryo has to make along the A-P axis, namely, expression of gap genes in response to the Bcd gradient in the anterior. Our results show that the slanting of the Hb and Otd expression domains can be directly traced to the differences in Bcd gradient properties between the two sides of the embryo (Fig. 4). Our input-output analyses further support the notion that Bcd gradient properties are primarily responsible for target expression differences on the dorsal and ventral sides of the embryo. Additional support to this conclusion also comes from previous findings that a Bcd-responsive reporter containing only Bcd binding sites (thus responding to Bcd as the sole activator input), or even reporters of artificially-generated A-P activator gradients, also exhibit a slant of the expression boundary . We suggest that the documented, extensive interactions between the A-P and D-V systems [26–28] may play important roles in maintaining or refining the expression of A-P genes that are either further downstream of the gene activation cascade or at times subsequent to early nuclear cycle 14, the stage of embryos analyzed in our work. We note that, unlike in the anterior (Fig. 4B), the iso-concentration contour lines of the Bcd gradient in the posterior part of the embryo do not exhibit any obvious ventral splay toward the posterior pole (data not shown). It is possible that the D-V system may play a particularly important role in A-P gene expression, in terms of ventral splay, in the posterior. This suggestion is consistent with recent quantitative findings that mutations in the D-V system preferentially reduce the splay of segmentation gene expression stripes in the posterior . An alternative explanation is that posterior gradients such as Nanos may have D-V differences analogous to the Bcd gradient. We also note that recent studies suggest that the terminal system and other gap gene products can interact with Bcd to regulate target gene expression [35–40]. However, these interactions, which are along the A-P axis, do not affect our interpretations of the sources of D-V differences regarding the observed slanting of Bcd target expression patterns.
Our experimental and simulation results suggest that proper Bcd gradient formation may require an interaction between the cortical layer and the inner part of the embryo. In particular, we found that the Bcd gradient profiles on the two sides of the same embryo are not identical even when measured as a function of contour distance c (Fig. 2B, see also simulated results in Fig. 5E). Ideally, if Bcd diffusion within the cortical layer were a simple 1-D process, we would have expected Bcd profiles on both sides of the embryo to be identical when measured as a function of contour distance c. Our experimental results show that, compared with the dorsal side, the Bcd gradient on the ventral side has a lower Bmax (see Fig. 2A and 2B) and a larger length constant λ (calculated by fitting to B x = B0e-x/λ; Fig. 2A values measured in x: λD = 99.2 ± 7.4 μm, λV = 112.6 ± 11.5 μm; Fig. 2B values measured in c: λD = 104.8 ± 6.9 μm, λV = 115.5 ± 9.0 μm). These results indicate that simple 1-D diffusion models, even when incorporating the appropriate consideration that the cortical layer is a critical medium for Bcd diffusion, are inadequate to address the geometric issues in such a complex biological process. In contrast, our 3-D simulation with an asymmetric embryo revealed results that are strikingly similar to these experimental findings including a lower Bmax and longer λ on the ventral side than the dorsal side (Fig. 5D values measured in x: λD = 92.3 μm, λV = 98.2 μm; Fig. 5E values measured in c: λD = 95.5 μm, λV = 105.8 μm). These differences are also obtained in a 2-D simulation of an asymmetric embryo in which, similar to the simulated 3-D embryo, Bcd molecules can diffuse within the inner part of the embryo and are allowed to exchange between the inner part and the cortical layer (not shown). However, in a 3-D simulation where the inner part of the embryo is impermeable to Bcd molecules (thus no exchange between the inner part and the cortical layer), the simulated λV is not larger than λD (not shown). We note that, although our 3-D simulation can replicate our experimentally observed Bcd gradient properties, it is possible that other physical features of the embryo may exert additional effects on Bcd gradient formation. For example, it has been shown that the ventral side of wt embryos has a lower nuclear density than the dorsal side [25, 41]. It is possible that the nuclear density difference may further affect Bcd gradient formation in ways (e.g., Bcd molecule diffusion) that are beyond those caused by the geometric differences, an issue that will require future investigations. Furthermore, nuclei on the cortex of the embryo at nuclear cycle 14 make systematic movements referred to as the nuclear flow  but such a flow was shown to have a relatively minor effect on the anterior Hb expression pattern. We emphasize that, since our current 3-D model can fully recapitulate our experimentally observed properties regarding the differences between the dorsal and ventral sides of the embryo, it likely has incorporated the key physical quantities most relevant to the system behavior.
The experimental and theoretical studies described in this report represent a first, systematic investigation of the effects of embryo geometry on the positional information encoded by the Bcd gradient. Our results suggest that the ventral splay of segmentation genes along the A-P axis may be directly traced to the distinct properties of the Bcd gradient (and its gap target gene responses) on the dorsal and ventral sides of the embryo. Our findings also suggest that, while the cortical layer is critical for Bcd gradient formation, Bcd diffusion likely requires an interaction between the cortical layer and the inner part of the embryo. These results demonstrate that modeling a morphogen gradient requires a realistic consideration of the physical properties, including the geometry, of the relevant biological system.
Immunostaining and microscopy
In our experiments, embryos were collected, fixed and stained with antibodies as described previously . 1×-bcd embryos were from bcd E1 /+ females and 3×-bcd embryos were from CyO, bcd/+ females (gift of Gary Struhl). For Bcd-Hb-DAPI triple staining, we used rabbit anti-Bcd  and guinea pig anti-Hb  as primary antibodies, and Alexa-594-conjugated goat anti-rabbit and Alexa-488-conjugated goat anti-guinea pig antibodies (Molecular Probes) as secondary antibodies. DAPI (Sigma) was used to stain the nuclei to facilitate imaging and data processing. For Otd-Hb-DAPI triple staining, we used guinea pig anti-Otd  and rat anti-Hb  as primary antibodies, and Alexa-488-conjugated goat anti-guinea pig and Alexa-555-conjugated goat anti-rat antibodies (Molecular Probes) as secondary antibodies. High-resolution images (1388 × 1040 8 bits/pixel) were captured on Zeiss Imager Z1 ApoTome microscope with a Zeiss Plan Apochromat 20× objective in linear settings (i.e, without any normalization during the imaging stage). Throughout this work, we use raw Bcd intensities that were captured within a linear range (i.e, without any saturation), and there were no further adjustments, except background subtraction when stated (see additional file 5 for raw Bcd intensity data). To ensure that background was measured under identical experimental conditions, we mixed experimental embryos with those from bcd E1 females (thus lacking Bcd). In our study, the mean Bcd intensities of bcd E1 embryos were simply subtracted as background from the intensity data of the experimental embryos. All stained embryos were mounted under cover slips with bridges to avoid geometric distortions, and images were taken on the midsagittal plane. The estimated noise in L measurement caused by focal plane selection is <0.1% , which is negligible compared with the measured differences between L and CD (~5%) or between CD and CV (~13%) in our current study.
Measurements of distance and fluorescence intensities
where j is the index of the scanning windows. Raw Bcd intensities were obtained by averaging the pixel intensities within the scanning windows, without further adjustments except background subtraction when necessary (Fig. 3B and 3C). The dorsal and ventral sides have similar background Bcd intensities as measured from bcd E1 embryos (2.19 ± 0.20 and 2.13 ± 0.22 for the dorsal and ventral sides, respectively; p = 0.27). In our analysis, the measured fluorescence intensities for both Hb and Otd were normalized by setting the mean of the highest values to 1 and the mean at non-expression regions to 0. See additional file 6 for Hb intensity profiles. Normalization of Hb intensities is required for our analysis, as well as in all published work in the field. This is due to the high variability of absolute Hb intensities among different embryos. However, since Hb expression has a clear boundary, unlike the Bcd gradient, such normalization is expected to cause problems that are relatively minor compared to normalization of the Bcd intensities. All calculations were performed using MatLab as described previously .
Plotting iso-intensity lines of Bcd onto an average embryo frame
We first compute the average shape of embryos by plotting the mean of normalized D-V positions along the perimeters of the embryos against the corresponding x/L. To compute the iso-intensity contour lines for Bcd, we choose a series of intensity values within the range of measured raw Bcd intensities as thresholds (t k ) to identify positions (x k , y k ) where BD = BV = t k by interpolation. Then we plot these positions onto the average embryo frame and connect the positions on the dorsal and ventral sides to visualize the distribution of iso-intensity contour lines (Fig. 4B). The iso-intensity lines are drawn for the anterior half of the embryo, where Bcd function is primarily relevant. In order to spread the iso-intensity lines for viewing quality, we use different Bcd intensity value intervals (tk+ 1 - t k ) for 1×-bcd, wt and 3×-bcd embryos (0.5, 1 and 1.5, respectively), whose average frames are calculated separately (N = 24, 28 and 29, respectively). For visual comparison, we also plot lines connecting the target boundary positions on two sides of the embryos using normalized average frames of 1×-bcd, wt and 3×-bcd embryos in the same presentation (Fig. 4D).
We use a 3-D simulation model, which is based on our recently developed 2-D model , to analyze the effects of both embryo geometry and cortex-yolk interaction on Bcd gradient formation. Briefly, our model is based on the equilibrium binding of Bcd molecules to non-specific, low affinity DNA binding sites throughout the Drosophila genome. In our model, Bcd molecules are produced at a constant rate within a bcd mRNA sphere located in the anterior. Unbound Bcd molecules at one location are allowed to diffuse to and from its six immediate neighboring locations thus truly in 3-D. Our model also considers the doubling of the nuclear number after each nuclear cycle and their relative locations (i.e., they are evenly distributed within the embryo prior to nuclear cycle 10, but migrate to the cortex to form a cortical layer at the onset of nuclear cycle 10). In our current work, an embryo was selected to resemble the average shape of a wt embryo with its dorsal and ventral curvatures on the midsagittal plane directly extracted from experimental measurements (Fig. 4B). This simulated embryo, which is shown in Fig. 5, is 560 μm in length and has a lateral diameter of 220 μm. We also simulated another embryo using an asymmetric approximation for its shape [45, 46] and obtained comparable results (see additional file 7 and the legend for further details).
In our model, the anterior pole of the simulated embryo has a coordinate of (0, 0, 0), and the bcd mRNA sphere from which Bcd is produced has a radius of 45 μm and a center coordinate of (75 μm, 0, 0). We tested our model by making reasonable changes in the shape or size of the mRNA distribution and obtained consistent Bcd protein gradient properties (not shown). The parameter value for the aggregate Bcd synthesis rate J for the entire embryo in our current model does not affect the experimentally observed three properties discussed in the main text (not shown), and we chose J = 1000 molecules/s so that at cycle 14, the maximum nuclear Bcd concentration on the dorsal side is 37.1 molecules/μm3, a value consistent with experimental estimates . The cortical layer has a depth of 10 μm in our current model and was kept unchanged between nuclear cycles 10 to 14; our simulations using thicker cortical layer similarly yielded the three experimentally observed properties (not shown). The simulated developmental time is based on experimental observations [18, 19], and the simulated results shown in Fig. 5 and additional file 7 are at 14 minutes into nuclear cycle 14. Other parameters used in our model are: the diffusion constant for free Bcd molecules, D = 2 μm2s-1; the spatially-uniform degradation rate of Bcd, ω = 0.00005 s-1, and the ratio of DNA-bound to free local Bcd concentrations at nuclear cycle 10, Bbound/Bfree = 0.15. As further detailed in , these parameter values were chosen through a systematic sampling such that the simulated Bcd profiles exhibit properties consistent with those observed experimentally, including proper λ values of the simulated Bcd profiles [6, 8], an enrichment of the Bcd molecules in the cortical layer , and stable nuclear Bcd concentrations at nuclear cycle 10-14 . These chosen values are biologically reasonable and, as shown in , are adjustable within our identified parameter space. In our simulations we tested grid size of 13, 23 and 53 μm3 and Δt = 0.1, 0.2, and 0.5 s, and obtained consistent results of Bcd gradient behaviors, indicating that our solutions are numerically stable. Results in Fig. 5 and additional file 7 were obtained with a grid size of 23 μm3 and Δt = 0.2 s. We also conducted simulations to evaluate the effects of other parameters on simulated Bcd gradient properties. As shown in additional file 8 (see legend for details), movements of the bcd mRNA sphere's center location along the A-P or D-V axis (Rx and Ry, respectively) or the Bbound/Bfree ratio at nuclear cycle 10, while affecting the specific values of the simulated profiles as expected, yielded results that remain generally consistent with the three experimentally observed properties. Together, these results show that our simulated Bcd properties that are in agreement with those observed experimentally are robust to reasonable variations (e.g., in bcd mRNA shape, size and location) that may naturally exist among real embryos.
We thank John Reinitz, Paul Macdonald, Tiffany Cook, Gary Struhl, and Christiane Nusslein-Volhard's lab for antibodies and fly strains, and Wei Wang and other lab members for discussions. This work was supported by grants from the NIH and NSF to JM, and the 973 programs (2009CB918702), the NSFC (90608029 and 30771217) to RJ.
- Kerszberg M, Wolpert L: Specifying positional information in the embryo: looking beyond morphogens. Cell. 2007, 130 (2): 205-209. 10.1016/j.cell.2007.06.038.View ArticlePubMedGoogle Scholar
- Lander AD: Morpheus unbound: reimagining the morphogen gradient. Cell. 2007, 128 (2): 245-256. 10.1016/j.cell.2007.01.004.View ArticlePubMedGoogle Scholar
- Arias AM, Hayward P: Filtering transcriptional noise during development: concepts and mechanisms. Nat Rev Genet. 2006, 7 (1): 34-44. 10.1038/nrg1750.View ArticlePubMedGoogle Scholar
- Lewis J: From signals to patterns: space, time, and mathematics in developmental biology. Science. 2008, 322 (5900): 399-403. 10.1126/science.1166154.View ArticlePubMedGoogle Scholar
- Jaeger J, Martinez-Arias A: Getting the measure of positional information. PLoS Biol. 2009, 7 (3): e81-10.1371/journal.pbio.1000081.View ArticlePubMedGoogle Scholar
- Gregor T, Tank DW, Wieschaus EF, Bialek W: Probing the limits to positional information. Cell. 2007, 130 (1): 153-164. 10.1016/j.cell.2007.05.025.PubMed CentralView ArticlePubMedGoogle Scholar
- Reinitz J: Developmental biology: a ten per cent solution. Nature. 2007, 448 (7152): 420-421. 10.1038/448420a.View ArticlePubMedGoogle Scholar
- He F, Wen Y, Deng J, Lin X, Lu LJ, Jiao R, Ma J: Probing intrinsic properties of a robust morphogen gradient in Drosophila. Dev Cell. 2008, 15 (4): 558-567. 10.1016/j.devcel.2008.09.004.PubMed CentralView ArticlePubMedGoogle Scholar
- He F, Saunders TE, Wen Y, Cheung D, Jiao R, ten Wolde PR, Howard M, Ma J: Shaping a morphogen gradient for positional precision. Biophys J. 2010Google Scholar
- Bollenbach T, Pantazis P, Kicheva A, Bokel C, Gonzalez-Gaitan M, Julicher F: Precision of the Dpp gradient. Development. 2008, 135 (6): 1137-1146. 10.1242/dev.012062.View ArticlePubMedGoogle Scholar
- Driever W, Nusslein-Volhard C: A gradient of bicoid protein in Drosophila embryos. Cell. 1988, 54 (1): 83-93. 10.1016/0092-8674(88)90182-1.View ArticlePubMedGoogle Scholar
- Ephrussi A, St Johnston D: Seeing is believing: the bicoid morphogen gradient matures. Cell. 2004, 116 (2): 143-152. 10.1016/S0092-8674(04)00037-6.View ArticlePubMedGoogle Scholar
- Struhl G, Struhl K, Macdonald PM: The gradient morphogen bicoid is a concentration-dependent transcriptional activator. Cell. 1989, 57 (7): 1259-1273. 10.1016/0092-8674(89)90062-7.View ArticlePubMedGoogle Scholar
- Driever W, Nusslein-Volhard C: The bicoid protein determines position in the Drosophila embryo in a concentration-dependent manner. Cell. 1988, 54 (1): 95-104. 10.1016/0092-8674(88)90183-3.View ArticlePubMedGoogle Scholar
- 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-154. 10.1038/342149a0.View ArticlePubMedGoogle Scholar
- Wolpert L: Positional information and the spatial pattern of cellular differentiation. J Theor Biol. 1969, 25 (1): 1-47. 10.1016/S0022-5193(69)80016-0.View ArticlePubMedGoogle Scholar
- 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-614. 10.1242/dev.031195.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- Gregor T, Wieschaus EF, McGregor AP, Bialek W, Tank DW: Stability and nuclear dynamics of the bicoid morphogen gradient. Cell. 2007, 130 (1): 141-152. 10.1016/j.cell.2007.05.026.PubMed CentralView ArticlePubMedGoogle Scholar
- Carroll SB, Scott MP: Zygotically active genes that affect the spatial expression of the fushi tarazu segmentation gene during early Drosophila embryogenesis. Cell. 1986, 45 (1): 113-126. 10.1016/0092-8674(86)90543-X.View ArticlePubMedGoogle Scholar
- Carroll SB, Winslow GM, Twombly VJ, Scott MP: Genes that control dorsoventral polarity affect gene expression along the anteroposterior axis of the Drosophila embryo. Development. 1987, 99 (3): 327-332.PubMedGoogle Scholar
- Gutjahr T, Frei E, Noll M: Complex regulation of early paired expression: initial activation by gap genes and pattern modulation by pair-rule genes. Development. 1993, 117 (2): 609-623.PubMedGoogle Scholar
- Yu Y, Pick L: Non-periodic cues generate seven ftz stripes in the Drosophila embryo. Mech Dev. 1995, 50 (2-3): 163-175. 10.1016/0925-4773(94)00333-I.View ArticlePubMedGoogle Scholar
- Kauffman SA: The Origins of Order: Self Organization and Selection in Evolution. 1993, New York: Oxford University PressGoogle Scholar
- Keranen SV, Fowlkes CC, Luengo Hendriks CL, Sudar D, Knowles DW, Malik J, Biggin MD: Three-dimensional morphology and gene expression in the Drosophila blastoderm at cellular resolution II: dynamics. Genome Biol. 2006, 7 (12): R124-10.1186/gb-2006-7-12-r124.PubMed CentralView ArticlePubMedGoogle Scholar
- Zeitlinger J, Zinzen RP, Stark A, Kellis M, Zhang H, Young RA, Levine M: Whole-genome ChIP-chip analysis of Dorsal, Twist, and Snail suggests integration of diverse patterning processes in the Drosophila embryo. Genes Dev. 2007, 21 (4): 385-390. 10.1101/gad.1509607.PubMed CentralView ArticlePubMedGoogle Scholar
- Li XY, MacArthur S, Bourgon R, Nix D, Pollard DA, Iyer VN, Hechmer A, Simirenko L, Stapleton M, Luengo Hendriks CL, et al: Transcription factors bind thousands of active and inactive regions in the Drosophila blastoderm. PLoS Biol. 2008, 6 (2): e27-10.1371/journal.pbio.0060027.PubMed CentralView ArticlePubMedGoogle Scholar
- Papatsenko D, Goltsev Y, Levine M: Organization of developmental enhancers in the Drosophila embryo. Nucleic Acids Res. 2009, 37 (17): 5665-5677. 10.1093/nar/gkp619.PubMed CentralView ArticlePubMedGoogle Scholar
- Ma X, Yuan D, Diepold K, Scarborough T, Ma J: The Drosophila morphogenetic protein Bicoid binds DNA cooperatively. Development. 1996, 122 (4): 1195-1206.PubMedGoogle Scholar
- Houchmandzadeh B, Wieschaus E, Leibler S: Establishment of developmental precision and proportions in the early Drosophila embryo. Nature. 2002, 415 (6873): 798-802.View ArticlePubMedGoogle Scholar
- Holloway DM, Harrison LG, Kosman D, Vanario-Alonso CE, Spirov AV: Analysis of pattern precision shows that Drosophila segmentation develops substantial independence from gradients of maternal gene products. Dev Dyn. 2006, 235 (11): 2949-2960. 10.1002/dvdy.20940.PubMed CentralView ArticlePubMedGoogle Scholar
- Spirov AV, Holloway DM: Making the body plan: precision in the genetic hierarchy of Drosophila embryo segmentation. Silico Biol. 2003, 3 (1-2): 89-100.Google Scholar
- Deng J, Wang W, Lu LJ, Ma J: A two-dimensional simulation model of the Bicoid gradient in Drosophila. PLoS ONE. 2010, 5 (4): e10275-10.1371/journal.pone.0010275.PubMed CentralView ArticlePubMedGoogle Scholar
- Crauk O, Dostatni N: Bicoid determines sharp and precise target gene expression in the Drosophila embryo. Curr Biol. 2005, 15 (21): 1888-1898. 10.1016/j.cub.2005.09.046.View ArticlePubMedGoogle Scholar
- Ochoa-Espinosa A, Yu D, Tsirigos A, Struffi P, Small S: Anterior-posterior positional information in the absence of a strong Bicoid gradient. Proc Natl Acad Sci USA. 2009, 106 (10): 3823-3828. 10.1073/pnas.0807878105.PubMed CentralView ArticlePubMedGoogle Scholar
- Surkova S, Spirov AV, Gursky VV, Janssens H, Kim AR, Radulescu O, Vanario-Alonso CE, Sharp DH, Samsonova M, et al: Canalization of gene expression in the Drosophila blastoderm by gap gene cross regulation. PLoS Biol. 2009, 7 (3): e1000049-PubMed CentralView ArticlePubMedGoogle Scholar
- Simpson-Brose M, Treisman J, Desplan C: Synergy between the hunchback and bicoid morphogens is required for anterior patterning in Drosophila. Cell. 1994, 78 (5): 855-865. 10.1016/S0092-8674(94)90622-X.View ArticlePubMedGoogle Scholar
- Lohr U, Chung HR, Beller M, Jackle H: Antagonistic action of Bicoid and the repressor Capicua determines the spatial limits of Drosophila head gene expression domains. Proc Natl Acad Sci USA. 2009, 106 (51): 21695-21700. 10.1073/pnas.0910225106.PubMed CentralView ArticlePubMedGoogle Scholar
- Jaeger J, Surkova S, Blagov M, Janssens H, Kosman D, Kozlov KN, Myasnikova E, Vanario-Alonso CE, Samsonova M, et al: Dynamic control of positional information in the early Drosophila embryo. Nature. 2004, 430 (6997): 368-371. 10.1038/nature02678.View ArticlePubMedGoogle Scholar
- Jaeger J, Blagov M, Kosman D, Kozlov KN, Myasnikova E, Surkova S, Vanario-Alonso CE, Samsonova M, Sharp DH, et al: Dynamical analysis of regulatory interactions in the gap gene system of Drosophila melanogaster. Genetics. 2004, 167 (4): 1721-1737. 10.1534/genetics.104.027334.PubMed CentralView ArticlePubMedGoogle Scholar
- Luengo Hendriks CL, Keranen SV, Fowlkes CC, Simirenko L, Weber GH, DePace AH, Henriquez C, Kaszuba DW, Hamann B, Eisen MB, et al: Three-dimensional morphology and gene expression in the Drosophila blastoderm at cellular resolution I: data acquisition pipeline. Genome Biol. 2006, 7 (12): R123-10.1186/gb-2006-7-12-r123.PubMed CentralView ArticlePubMedGoogle Scholar
- MacDonald PM: bicoid mRNA localization signal: phylogenetic conservation of function and RNA secondary structure. Development. 1990, 110 (1): 161-171.PubMedGoogle Scholar
- Kosman D, Small S, Reinitz J: Rapid preparation of a panel of polyclonal antibodies to Drosophila segmentation proteins. Dev Genes Evol. 1998, 208 (5): 290-294. 10.1007/s004270050184.View ArticlePubMedGoogle Scholar
- Xie B, Charlton-Perkins M, McDonald E, Gebelein B, Cook T: Senseless functions as a molecular switch for color photoreceptor differentiation in Drosophila. Development. 2007, 134 (23): 4243-4253. 10.1242/dev.012781.View ArticlePubMedGoogle Scholar
- Fowlkes CC, Hendriks CL, Keranen SV, Weber GH, Rubel O, Huang MY, Chatoor S, DePace AH, Simirenko L, Henriquez C, et al: A quantitative spatiotemporal atlas of gene expression in the Drosophila blastoderm. Cell. 2008, 133 (2): 364-374. 10.1016/j.cell.2008.01.053.View ArticlePubMedGoogle Scholar
- Umulis DM, Shimmi O, O'Connor MB, Othmer HG: Organism-scale modeling of early Drosophila patterning via bone morphogenetic proteins. Dev Cell. 2010, 18 (2): 260-274. 10.1016/j.devcel.2010.01.006.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.