Quantitative analyses reveal extracellular dynamics of Wnt ligands in Xenopus embryos
Abstract
The mechanism of intercellular transport of Wnt ligands is still a matter of debate. To better understand this issue, we examined the distribution and dynamics of Wnt8 in Xenopus embryos. While Venus-tagged Wnt8 was found on the surfaces of cells close to Wnt-producing cells, we also detected its dispersal over distances of 15 cell diameters. A combination of fluorescence correlation spectroscopy and quantitative imaging suggested that only a small proportion of Wnt8 ligands diffuses freely, whereas most Wnt8 molecules are bound to cell surfaces. Fluorescence decay after photoconversion showed that Wnt8 ligands bound on cell surfaces decrease exponentially, suggesting a dynamic exchange of bound forms of Wnt ligands. Mathematical modeling based on this exchange recapitulates a graded distribution of bound, but not free, Wnt ligands. Based on these results, we propose that Wnt distribution in tissues is controlled by a dynamic exchange of its abundant bound and rare free populations.
Introduction
The Wnt family of secreted signaling proteins has diverse roles in animal development, stem cell systems, and carcinogenesis (Clevers et al., 2014; Loh et al., 2016; Nusse and Clevers, 2017). It has been generally accepted that in the extracellular space, morphogenic Wnt ligands form a concentration gradient by dispersal (Clevers et al., 2014; Kiecker and Niehrs, 2001; Müller et al., 2013; Smith, 2009; Strigini and Cohen, 2000; Tabata and Takei, 2004; Yan and Lin, 2009; Zecca et al., 1996; Zhu and Scott, 2004). In contrast to this classical view, evidence also suggests dispersal-independent functions of Wnt ligands. For instance, a membrane-tethered form of Wingless (Wg) can recapitulate an almost normal pattern of Drosophila wings, suggesting that dispersal of Wg is dispensable for patterning (Alexandre et al., 2014). This dispersal-independent patterning can be explained by gradual attenuation of Wg expression in distally localized cells in which Wg was formerly expressed. However, it remains unclear to what extent dispersal-dependent and/or -independent mechanisms contribute to the graded distribution of Wnt proteins in tissue patterning.
Visualization of Wnt ligands is essential to understand their distributions. In the wing disc of Drosophila, Wg proteins are widely distributed from wing margin cells, where Wg is expressed (Strigini and Cohen, 2000; Zecca et al., 1996). Furthermore, long-range dispersal of Wg was evidenced by an experiment in which Wg was captured by distally expressed Frizzled2, a Wg receptor (Chaudhary et al., 2019). Similarly, endogenous Wnt ligands tagged with fluorescent proteins showed long-range distributions in C. elegans (Pani and Goldstein, 2018). In addition to these observations in invertebrates, we found that endogenous Wnt8 ligands disperse far from their source cells in Xenopus embryos (Mii et al., 2017). On the other hand, mouse Wnt3 accumulates within a few cell diameters of its source cells in the microenvironment of the intestine (Farin et al., 2016). These studies show that Wnt ligands apparently disperse in tissues and embryos, although the dispersal range varies. Importantly, in many of these studies, Wnt ligands accumulate locally on cell surfaces, showing punctate distributional patterns (Pani and Goldstein, 2018; Strigini and Cohen, 2000; Zecca et al., 1996). Furthermore, we demonstrated that Wnt8 and Frzb, a secreted Wnt inhibitor, accumulate separately and locally on cell surfaces in Xenopus embryos (Mii et al., 2017). However, these punctate accumulations on cell surfaces, largely ignored in the literature in the context of Wnt gradient formation, raise the question of whether such accumulations contribute to formation of concentration gradients in tissues and embryos.
Studies in Drosophila wing disc have shown that cell surface scaffolds, such as heparan sulfate (HS) proteoglycans (HSPGs), are required for both distribution and delivery of morphogens, including Wg, Hedgehog (Hh), and Decapentaplegic (Dpp) (Franch-Marro et al., 2005; Lin, 2004; Yan and Lin, 2009). From these studies, the ‘restricted diffusion’ model, in which morphogens are transferred extracellularly by interacting with cell surface scaffolds, has been proposed (Yan and Lin, 2009). In this model, the movement of each morphogen molecule is constrained in a ‘bucket brigade’ fashion by interactions with cell surface scaffolds. As a result of continuous interactions, morphogen molecules are slowly transferred (Han et al., 2005; Yan and Lin, 2009 #152; Kerszberg and Wolpert, 1998; Takei et al., 2004). However, it seems difficult to explain local accumulations of Wnt proteins by the restricted diffusion mechanism, because passive diffusion alone should result in smoothly decreasing gradients. On the other hand, we recently showed that HSPGs on cell surfaces are discretely distributed in a punctate manner, which varies with heparan sulfate (HS) modification, forming two different types of HS clusters, N-sulfo-rich and N-acetyl-rich forms (Mii et al., 2017). Notably, Wnt8 and Frzb, a secreted Frizzled-related protein (sFRP), accumulate separately on N-sulfo-rich and N-acetyl-rich HS clusters, respectively. Frzb expands the distribution and signaling range of Wnt8 by forming heterocomplexes (Mii and Taira, 2009), and Wnt8/Frzb complexes are colocalized with N-acetyl-rich HS clusters (Mii et al., 2017). N-sulfo-rich clusters are frequently internalized together with Wnt8, whereas N-acetyl-rich HS clusters tend to remain on the cell surface. This difference in stability on the cell surface may account for the short-range distribution of Wnt8 and the long-range distribution of Frzb (Mii and Taira, 2009; Mii et al., 2017) and suggests that the distribution of HS clusters should be considered in order to understand extracellular dynamics of Wnt ligands (Mii and Takada, 2020).
To explain the dynamics of Wnt ligands in tissues, quantitative analyses of Wnt ligands are required. Dynamics of secreted proteins have been investigated using fluorescence recovery after photobleaching (FRAP) (Sprague and McNally, 2005; Sprague et al., 2004) and fluorescence correlation spectroscopy (FCS), although optimal ranges for diffusion coefficients differ (Hess et al., 2002; Kicheva et al., 2012; Müller et al., 2013; Fradin, 2017). For example, FRAP measurements have shown that Dpp and Wg diffuse slowly in the Drosophila wing disc with diffusion coefficients ranging from 0.05 to 0.10 μm2/s, suggestive of the restricted diffusion model (Kicheva et al., 2007). In contrast, FCS measurements of FGF8 in zebrafish embryos showed fast, virtually free diffusion, with a diffusion coefficient of ~50 μm2/s (Yu et al., 2009). Furthermore, in contrast to the FRAP results, free diffusion of Dpp measured in the Drosophila wing disc using FCS yielded a diffusion coefficient of ~20 μm2/s (Zhou et al., 2012). FCS is based on fixed-point scanning within a confocal volume (typically sub-femtoliter) for several seconds, while FRAP evaluates considerably larger regions of photobleaching/photoconversion, containing tens or hundreds of cells (Rogers and Schier, 2011) and spanning long time windows (typically several hours). Under these experimental conditions for FRAP, it is proposed that diffusion of secreted proteins is affected by zigzag paths of the narrow intercellular space between polygonal epithelial cells, instead of an open, unobstructed space (hindered diffusion model) (Müller et al., 2013), and/or by endocytosis, which reduces the concentration of the diffusing species in the extracellular space. Thus, we need exercise caution when comparing data derived from FRAP and from FCS analyses.
In this study, we examined extracellular dynamics of Wnt8 and Frzb, both of which are involved in anteroposterior patterning of vertebrate embryos (Clevers and Nusse, 2012; Kiecker and Niehrs, 2001; MacDonald et al., 2009; Mii et al., 2017). First, we visualized their localization in Xenopus embryos by fusing them with fluorescent proteins and we examined their dispersion by capturing them in distant cells. We also examined their dispersal dynamics using FCS and fluorescence decay after photoconversion (FDAP) measurements in embryonic tissue. In particular, we refined FDAP-based analysis by focusing on a limited area at the cell boundary, which enabled us to quantify dynamics comparable to those measured by FCS. Based on these results and our previous findings, we propose a basic mathematical model to explain distribution and dispersion of secreted proteins.
Results
Extracellular distributions of secreted proteins depend on interactions with cell-surface molecules
As we have previously shown (Mii et al., 2017), Wnt8 and Frzb fused with monomeric Venus (mV) were visualized along cell boundaries when expressed in Xenopus embryos (Figure 1A). We note that biological activities of these proteins were not severely impaired by the fusion of mV and that the reduced activity of mV-tagged Wnt8 compared to untagged Wnt8 could possibly be due, at least in part, to differences in translation (Figure 1—figure supplement 1). In contrast, we found that only the secreted form of mV (sec-mV), which was not expected to bind specifically to the cell surface, was hardly visible along the cell boundary under the same conditions (Figure 1A, right). Since Wnt8 and Frzb colocalize with heparan sulfate clusters on cell surfaces, we speculated that binding to cell surface proteins, like heparan sulfate proteoglycans (HSPGs), affects the distribution of Wnt8 and Frzb. To examine this possibility, we added heparin-binding (HB) peptides, consisting of 16 (ARKKAAKA)2 (HB2) or 32 amino acids (ARKKAAKA)4 (HB4) (Verrecchio et al., 2000; Figure 1C) to sec-mV. Addition of HB peptides significantly increased the intensity of mVenus fluorescence in the intercellular region compared to that of sec-mV. This suggests that the intercellular distribution of secreted proteins depends on interactions with docking molecules on cell surfaces.
To directly examine this idea, we constructed a reconstitution system, consisting of HA-epitope-tagged secreted mVenus (sec-mV-2HA) and a membrane-tethered anti-HA antibody (‘tethered-anti-HA Ab’) (Figure 2A, see Figure 2—figure supplement 1 for cDNA cloning and validation of anti-HA antibody). This artificial protein and tethered-anti-HA Ab were expressed in separated areas in the animal cap region of Xenopus gastrulae. As with sec-mV, sec-mV-2HA was hardly visible in the intercellular space, even close to the source cells (Figure 2B). In contrast, sec-mV-2HA was observed around tethered-anti-HA Ab-expressing cells that were traced with memRFP, even though these cells were distantly located from the source cells (Figure 2B). Thus, interaction with cell surface proteins can affect distributions of secreted proteins.
This result also indicates that diffusing proteins are not readily visible using standard confocal microscopy, unless they are trapped by cell surface proteins. In fact, quantitative analysis of artificial secreted proteins revealed a slight, but significant increase of photon counts in the intercellular region by injection of mRNA for sec-mV, compared to uninjected embryos, indicating that sec-mV actually exists in the intercellular region (Figure 1D, Figure 1—figure supplement 2).
Populations of secreted Wnt8 and Frzb proteins disperse long distance
We next examined dispersal of molecules of mV-Wnt8 and mV-Frzb. Both mV-Wnt8 and mV-Frzb accumulated locally along the cell boundary at the subapical level (Figure 1A and B), consistent with previous observations (Mii et al., 2017), which indicated that populations of Wnt8 and Frzb in the intercellular space were bound to the cell surface at HS clusters. On the other hand, given that some molecules of mV-Wnt8 or mV-Frzb may drift away from the cell surface, these proteins would be almost undetectable with standard confocal microscopy, as exemplified by sec-mV (Figure 1A) and tethered-anti-HA Ab (Figure 2B). To examine such mobile proteins, we tried to capture them using ‘morphotrap’ located distantly from the source cells (Figure 2C). Morphotrap is a membrane-tethered form of anti-GFP nanobody, originally devised to block dispersal of Dpp-GFP from source cells (Harmansa et al., 2015). We supposed that morphotraps could be utilized to detect or visualize diffusible proteins, similar to tethered-anti-HA Ab. As expected, sec-mV accumulated on the surface of morphotrap-expressing cells remote from source cells (Figure 2C). Similarly, mV-Wnt8 and mV-Frzb were trapped (Figure 2C), evidencing the long-distance dispersal (over 15 cells/200 μm) of some of secreted mV-Wnt8 and mV-Frzb molecules. These proteins are not likely to be transferred by cell-movement-based mechanisms, including distant migration of source and morphotrap-expressing cells, because cells in the animal cap region form an epithelial sheet and are tightly packed. In addition to this dispersing population, mV-Wnt8 and mV-Frzb were also detectable in gradients from producing cells to morphotrap-expressing cells, unlike the case of sec-mV (Figure 2C and D). These results suggest that populations of mV-tagged Wnt8 and Frzb do not associate tightly with cell surfaces, thereby potentially dispersing far from source cells.
FCS analyses combined with quantitative imaging reveal cell-surface-bound and diffusing Wnt8 and Frzb proteins in the extracellular space
We next attempted to quantify the populations of Wnt8 or Frzb proteins associated with cell surfaces and diffusing in the extracellular space. For this purpose, we employed fluorescence correlation spectroscopy (FCS). FCS analyzes fluctuation of fluorescence by Brownian motion of fluorescent molecules in a sub-femtoliter confocal detection volume (Figure 3, A and B). By autocorrelation analysis (Figure 3C), FCS can measure diffusion coefficients (D) of mobile molecules and the number of particles in the detection volume, but inference of diffusion coefficients depends on mobile molecules (Hess et al., 2002; Fradin, 2017). FCS analyses were performed by injecting the same doses of mV-Wnt8 and sec-mV that were used in the experiments shown in Figure 1 (250 pg mRNA/embryo) to consider the relationship between photon counting from live-imaging and NoP from FCS. Furthermore, to measure the dynamics of mV-Wnt8 and mV-Frzb at a concentration equivalent to the endogenous concentration, a decreased amount of RNA was also injected (20 pg mRNA/embryo, see Figure 1—figure supplement 1A).
To analyze the data obtained by FCS measurements, we compared the suitability of one-component and two-component diffusion models using the Akaike information criterion (AIC) (Tsutsumi et al., 2016). AIC supported fitting with the two-component model, comprising fast and slow diffusing components (Figure 3—figure supplement 1A). Consistent with predictions, the result indicated that the number of particles (NoP) of mV-Wnt8 (250 pg/embryo) was significantly higher than that of mV-Wnt8 (20 pg/embryo, endogenous-equivalent level; Figure 3D and Figure 3—figure supplement 1B). Mean values of Dfast indicate that fast-diffusing components in all groups examined can be regarded as free diffusion (Figure 3D), because theoretical, as well as reported D values of freely diffusing proteins of similar size, range from 10 to 100 μm2/s (Pack et al., 2006; Yu et al., 2009; Zhou et al., 2012). Importantly, even at the endogenous-equivalent level, mV-Wnt8 and mV-Frzb show freely diffusing populations (Dfast >10 μm2/s, Figure 3—figure supplement 1B). We also note that the diffusion coefficient of the fast component of mV-Wnt8 (20 pg/embryo, endogenous-equivalent level) was significantly lower than that of mV-Wnt8 (250 pg/embryo), suggesting stronger constraints with the endogenous-equivalent expression level. Thus, we conclude that within the small volume of FCS measurements, a population of mV-Wnt8 and mV-Frzb molecules diffuses freely under physiological conditions.
As shown in Figure 1D, photon counts of mV-Wnt8 were much higher than those of sec-mV. However, under the same conditions as in Figure 1 (250 pg mRNA/embryo), NoP of mV-Wnt8 was similar to that of sec-mV (Figure 3—figure supplement 1B) or even smaller in another set of measurements (Figure 3—figure supplement 2). Thus, molecules detected in FCS appear not to contribute to the photon counts in the confocal imaging under these conditions. We speculate that FCS measurements might be biased to choose positions where HS-bound molecules are not abundant. Otherwise, HS-bound, immobile molecules cause strong photobleaching, which results in large drift of the fluorescence intensity. In general, such a data is not suitable for analysis.
Interestingly, slow components were observed not only in mV-Wnt8, but also in sec-mV (Figure 3D). To characterize these slow components, we examined the effects of HS-chain digestion with FCS. These analyses were performed with embryos injected with 250 pg/embryo RNA for mV-Wnt8 or sec-mV, because at the endogenous-equivalent level, measured values showed a large variance, possibly reflecting heterogeneity of the extracellular space, and also signal detection was difficult for sec-mV. For this purpose, we made a membrane-tethered form of Heparinase III (HepIII-HA-GPI, also known as heparitinase I) (Hashimoto et al., 2014). HepIII-HA-GPI enables us to digest HS chains in a region of interest (Figure 3—figure supplement 3), allowing us to examine the effects of HS-digestion in the same embryos. For mV-Wnt8, NoP and the fraction of fast components, Ffast was significantly increased by HepIII, suggesting release of mV-Wnt8 from HS chains (Figure 3E and F). Thus, we suggest that HS chains contribute to the slow components of mV-Wnt8. For sec-mV, although NoP was not significantly changed by HepIII, Ffast was slightly, but significantly increased by HepIII (Figure 3G and H). Furthermore, fluorescence cross-correlation spectroscopy (FCCS) analysis indicated that sec-mV did not interact with the cell membrane (Figure 3—figure supplement 1C and D). Thus, HS-chains showed some contribution to the slow components of sec-mV even without interaction. We speculate that such a slow population of sec-mV could be explained by hindered diffusion, in which torturous diffusion results from HS chains and other ECMs, because HS chains are highly hydrophilic and well hydrated (Figure 3—figure supplement 1E).
FDAP analyses suggest exchange of cell-surface-bound and unbound states of Wnt8 and Frzb proteins
Although FCS analysis is suitable for measuring diffusing molecules, it cannot directly analyze molecules with extremely low mobility (Hess et al., 2002). To directly analyze dynamics of such molecules, we next employed fluorescence decay after photostimulation/photoconversion (FDAP) assays (Matsuda et al., 2008; Müller et al., 2012) in the intercellular space of Xenopus embryos.
Since FRAP/FDAP measurements usually examine considerably larger regions (typically containing tens or hundreds of cells) than with FCS (Rogers and Schier, 2011), direct comparisons of dynamics between FRAP/FDAP and FCS may need careful consideration. Therefore, we restricted the area of photoconversion to a diameter of 1.66 μm and reduced the measurement time (16 s), allowing us to obtain dynamic data in the intercellular region under conditions comparable to those for FCS (Figure 4A). We refer to this FDAP mode as ‘cell-boundary FDAP.’ In this analysis, we fused a photoconvertible fluorescent protein, mKikGR (mK) (Habuchi et al., 2008) to Wnt8 and Frzb (mK-Wnt8 and mK-Frzb). These fusion proteins showed distributions in embryos similar to mV-tagged proteins (Figure 4A), and retained biological activities (Figure 4—figure supplement 1). Importantly, observed distribution patterns of mK-Wnt8 and mK-Frzb were stable for up to tens of minutes (Figure 4A). Therefore, we assumed a steady state during the FDAP analysis (16 s).
After photoconversion, red fluorescent intensity of mK-tagged proteins was measured in the same rectangular area as photoconversion (Figure 4B). Because puncta of Wnt8 are often internalized with HS clusters (Mii et al., 2017), we excluded data in which vesicular incorporation was observed during measurement of mK-Wnt8. As an immobilized control, mK-Frzb in formaldehyde-fixed embryos was similarly photoconverted. Its intensity within the rectangular region was fitted to the photobleaching model using repeated laser scanning, confirming that it actually was immobilized by formaldehyde fixation (Figure 4—figure supplement 2B). Compared with this fixed control, mK-Wnt8 and mK-Frzb showed faster decline of fluorescent intensities (Figure 4C and Figure 4—figure supplement 2A), confirming that a population of mK-Wnt8 and mK-Frzb moved away from the photoconverted area.
Given that the punctate distribution of mK-Wnt8 and -Frzb results from their binding to HS clusters (Mii et al., 2017), we considered whether a simple dissociation model (Equation 1) is suitable for curve-fitting of FDAP data. Indeed, bleaching-corrected FDAP curves of mK-Wnt8 and mK-Frzb were well fitted to this model (Figure 4D; residuals were mostly within 5% and all within 10%) with the indicated parameters (Figure 4E, the off-rate constant koff, and the rate of the constantly bound component C; for individual data plot, see Figure 4—figure supplement 2D). As a result, both mK-Wnt8 and mK-Frzb show large C values, indicating that the majority of these proteins can be considered immobile on the timescale of the measurement (Figure 4E). In addition, koff of mK-Wnt8 was significantly lower than that of mK-Frzb (Figure 4—figure supplement 2D), suggesting relatively rapid dissociation of mK-Frzb from the binding site. This difference appears to be consistent with FDAP spatial intensity profiles, in which photoconverted mK-Frzb, but not mK-Wnt8, accumulated in adjacent areas (Figure 4—figure supplement 2C, see also Videos 1, 2 and 3). Thus, we conclude that most mK-Wnt8 and mK-Frzb molecules are bound, but can be exchanged with unbound molecules, and also dissociation rate values of mK-Wnt8 and mK-Frzb differ significantly.
Mathematically modeling diffusion and distribution of secreted proteins
Based on our quantitative imaging (Figure 1) of tethered-anti-HA Ab and morphotrap (Figure 2) FCS (Figure 3) and FDAP (Figure 4), we conclude that most Wnt8 and Frzb molecules are bound to cell surfaces, while small numbers of freely diffusing molecules exist in the extracellular space. Furthermore, we have already shown that Wnt8 and Frzb utilize different types of HS clusters, N-sulfo-rich and N-acetyl-rich, as cell-surface scaffolds, respectively (Mii et al., 2017). Thus, we examined whether free diffusion and binding to HS clusters can explain the extracellular distribution or gradient formation of secreted proteins, using mathematical modeling.
Here, we consider two states of ligands: free and bound. The free state corresponds to the fast diffusing component in FCS, and we consider the bound component as immobile molecules. This model includes five dynamic processes: (i) ligand production, (ii) diffusion of free molecules in intercellular space, (iii) binding of ligands to HS clusters on cell surfaces, (iv) release of bound molecules from HS clusters and (v) internalization of bound molecules into cells. In one-dimensional space, the model is written as:
where u and v represent the concentration of free molecules and numbers of bound molecules, respectively, of a secreted protein. The symbols x and t are position and time, respectively. The symbols a(x), b, c, and g(x) represent binding, release, internalization, and production rates, respectively (Figure 5A); a(x) is equivalent to the amount of HS in HS clusters (for details, see Materials and methods). D (= 20 μm2/s) represents the diffusion coefficient of the free component in the extracellular space, which corresponds to the fast diffusing component measured by FCS (Figure 3D).
Under a wide range of appropriate parameter values, distributions of u and v converged to steady states within a few minutes. Compared to the fast diffusing component (Figure 5B), the contribution of the slow component (D = 0.50 μm2/s) to the distribution range is much smaller (Figure 5—figure supplement 1A). Hense, we mainly consider the fast component observed in FCS (Figure 3D), as the diffusing population in the model. The free component, u, quickly decreases, displaying a shallow continuous distribution pattern due to diffusion. In contrast, the bound component, v, shows a discrete distribution following a(x), and the level of v is much higher than u at any position, reflecting our conclusion that the majority of Wnt or sFRPs molecules in the extracellular space are bound. Given that activation of Wnt signaling requires internalization of the ligands (Kikuchi et al., 2009; Yamamoto et al., 2006) of the bound component, corresponding to cv in Equation 4, the distribution of the bound component in this model could be equivalent to the ‘actual’ gradient of Wnt signaling, even though it is not diffusing. We demonstrated that consistently some portion of Wnt8 ligands accumulated on N-sulfo-rich HS clusters initiate canonical Wnt signaling by forming the signaling complex ‘signalosome’ (Mii et al., 2017).
We have shown that N-sulfo-rich HS, but not N-acetyl-rich clusters, are frequently endocytosed (Mii et al., 2017). In this model, different internalization rates of N-acetyl-rich and N-sulfo-rich HS clusters can be reflected by varying the internalization rate of the docking sites, c. A smaller value of c results in long-range distribution (compare Figure 5B and C), explaining why Frzb shows a long-range distribution by binding to N-acetyl-rich HS clusters (Mii et al., 2017). We can evaluate the distribution by the decay length, λ. λ represents a distribution range when the steady state gradient is written as
(Kicheva et al., 2012; Kicheva et al., 2007). We calculated λ by curve-fitting the peak values of v to Equation 5. The value of λ with c = 0.1 or 0.01 is 6.346 or 10.79 μm (Figure 5B,C and Figure 5—figure supplement 1B,C for normalized plots), respectively, showing that internalization rates of HS clusters can affect distribution ranges, as observed between Wnt and sFRPs (Mii and Taira, 2009). In addition, we examined the contribution of dissociation from the bound to the diffusing state, suggested by our cell-boundary FDAP (Figure 4). Without dissociation, a shorter range distribution of the bound component was obtained (λ = 4.504 μm, Figure 5—figure supplement 1D) than with dissociation (λ = 6.346 μm, same data as in Figure 5B). Furthermore, in Xenopus embryos, Wnt8 in the intercellular space exhibited local accumulations (Figure 1B). In our model, when the binding rate an,max (Equation 7 in the Materials and methods) fluctuates randomly (i.e. the amount of HS at position x), the bound ligand component also fluctuates (Figure 5D, blue), reproducing the local accumulation of Wnt8 and Frzb in Xenopus embryos. Even under these conditions, the free component shows a continuously decreasing gradient (Figure 5D, red), which probably corresponds to the FCS-measured, diffusing component of the FGF8 gradient in zebrafish embryos (Yu et al., 2009; measuring concentrations by FCS in a wide field is technically difficult in larger Xenopus embryos). Thus, our mathematical model can generalize protein distributions in the extracellular space.
Discussion
As one of the major secreted signaling molecules, mechanisms of Wnt dispersal are crucial when we consider embryonic patterning and various other systems involving Wnt signaling (Routledge and Scholpp, 2019). Among many Wnt proteins, Wg distribution in the Drosophila wing disc has long been investigated as a morphogen gradient (Strigini and Cohen, 2000; Zecca et al., 1996). Various genetic studies show that the extracellular distribution of Wg largely depends on HSPGs, such as Dally and Dally-like glypicans (Baeg et al., 2004; Franch-Marro et al., 2005; Han et al., 2005). Furthermore, FRAP-based analysis suggests that the effective diffusion coefficient of Wg is much slower (0.05 μm2/s) than free diffusion (>10 μm2/s) (Kicheva et al., 2007). However, such dynamics of secreted signaling proteins still remain a matter of debate (Rogers and Schier, 2011). On the other hand, recently we found that HS chains on the cell surface are organized in clusters with varying degrees of N-sulfo modification in Xenopus embryos and HeLa cells. Furthermore, we demonstrated that endogenous Wnt8 protein visualized by immunostaining shows a punctate distribution, specifically associated with N-sulfo-rich HS clusters (Mii et al., 2017). Similar punctate distributions have also been observed with Wg in Drosophila (Strigini and Cohen, 2000; van den Heuvel et al., 1989), but the significance of these distributions has not yet been explained. Therefore, to gain insight into the mechanism of Wnt distribution, we examined Wnt8 protein dynamics.
Based on quantitative live-imaging techniques, we propose that most Wnt8 molecules distributed among cells are mostly cell-surface-bound, while a small portion of them are diffusing. Similarly, Wnt/EGL-20 shows that puncta mostly overlap with Frizzled and a small population of mobile/diffusing molecules is also suggested in C. elegans (Pani and Goldstein, 2018). In Xenopus embryos, Frizzled may also contribute to bind Wnt ligands because some Wnt8 puncta overlapped with Frizzled8 (Mii et al., 2017). Furthermore, a small population of diffusing Dpp has been shown in Drosophila wing disc (Zhou et al., 2012). Importantly, it has been suggested that these populations disperse over long distances, similar to our observation of mV-Wnt8 trapped using morphotrap (Figure 2C,D), generalizing the existence of long-dispersing populations in various model systems.
It is plausible that cell-surface-bound Wnt8 is mostly associated with HS clusters (Mii et al., 2017). The function of HSPGs in Wnt dispersal has been examined by genetic studies of Drosophila. These studies show that HSPGs are required for accumulation and transfer of Wnt ligands. Based on these results, it has been proposed that Wnt disperses by restricted diffusion, in which HSPGs transfer Wnt ligands in a bucket brigade manner (Yan and Lin, 2009). In our FDAP assay, most photoconverted mK-Wnt8 does not diffuse laterally, even when other puncta of Wnt8 exist near the site of photoconversion (Figure 4—figure supplement 2C, left panel, Video 1). We further considered this observation with modeling (Figure 5—figure supplement 1E). Unlike the experiment, modeling shows lateral dispersal of photoconverted molecules in neighboring regions, over time. For this difference, we mainly consider two possibilities: (i) Our imaging system may not be sufficiently sensitive to detect such a small increase, or the increase of the ligands may be obliterated by photobleaching. (ii) If binding dynamics of mK-Wnt8 are slower than in our model, ligands may diffuse away before re-binding in neighboring regions. On the other hand, mK-Frzb showed some lateral dispersal similar to the model. As previously discussed, these behaviors in FRAP experiments can be classified into some cases including ‘reaction dominant’ and ‘effective diffusion’ by a balance among the on-rate, the off-rate, and the diffusion coefficient (Sprague and McNally, 2005; Sprague et al., 2004). Restricted diffusion can be understood as a kind of effective diffusion in which dynamics of ligand binding/dissociation to HSPGs are similar to those of free diffusion. Although we did not derive binding constants, at least superficially, mK-Frzb showed an effective diffusion-like behavior, whereas mK-Wnt8 showed a reaction dominant-like behavior, in which free diffusion is much faster than binding/dissociation. In order to compare our data with those previously reported (Kicheva et al., 2012), we also performed curve-fitting with an effective/apparent diffusion model (Figure 4—figure supplement 3, Equation 2). As a result, the apparent diffusion coefficient Da (μm2/s) was calculated as 0.042 and 0.059 for mKikGR-Wnt8 and mKikGR-Frzb, respectively. These values are very close to a previously reported FRAP value for Wg (0.05 μm2/s) (Kicheva et al., 2007). Thus, such small values of Da relative to free diffusion could be interpreted as the result of interaction with cell surfaces, regardless of whether the protein of interest actually shows lateral diffusion in bucket brigade fashion.
We found that sec-mV is almost invisible with standard confocal microscopy (Figure 1B). Furthermore, binding to cell surface molecules such as HSPGs and membrane-tethered antibody was sufficient for visible distribution for artificial secreted proteins (Figures 1C and 2B). These findings are similar to recent demonstrations that secreted GFP can be synthetic morphogens with specific scaffold molecules in the Drosophila wing disc (Stapornwongkul et al., 2020) and in cultured cells (Toda et al., 2020). On the other hand, secreted GFP appears visible in some tissues, such as deep cells in early zebrafish embryos (Yu et al., 2009) and developing zebrafish brain (Veerapathiran et al., 2020). In the zebrafish brain, secreted EGFP did not show slow components (Veerapathiran et al., 2020), which is different from our observation in the Xenopus animal cap region (Figure 3D). Together, considering detection of diffusing molecules (Appendix), we speculate that these differences may reflect narrowness of the extracellular space in the tissues.
Our cell-boundary FDAP suggests that cell-surface-bound and diffusing populations are probably exchangeable. Although this result can be explained by dissociation of molecules from the bound state as described above, it also seems possible that endocytosis reduces the number of photoconverted molecules (Figure 4C). However, we consider this less likely. Endocytosis of Wnt8 is possibly mediated by caveolin (Mii et al., 2017; Yamamoto et al., 2006), and we have already shown that internalized Wnt8 was detected as puncta in the cell (Mii et al., 2017). However, in FDAP analyses in this study, we excluded observations with internalization of Wnt puncta from curve-fitting analysis. In our mathematical model, when dissociation from the cell-surface does not occur (b = 0 in Equation 3 and 4, Figure 5A), the range of the gradient (decay length, λ) was shortened from 6.35 to 4.50 μm (Figure 5—figure supplement 1D). Thus, at least in cases we analyzed, dissociation from the bound state seems to contribute to the long-range distribution and rapid formation of the gradient.
A goal of this study is to link quantitative measurements of local protein dynamics to larger spatiotemporal patterns of extracellular protein dispersal in embryos. We hypothesized that local dynamics of diffusion and interaction with HS chains measured by FCS and FDAP could be extrapolated to explain mechanisms for gradient formation across many cells. We mainly consider protein dispersal within a single plain, and this is exemplified in the animal cap region since mV-Wnt8 and mV-Frzb accumulated on the proximal side (to the source) of morphotrap-expressing cells (Figure 2C). But when we consider dispersal of secreted proteins in embryos, other routes can be involved. For example, a BMP antagonist, Chordin exhibits dispersal within the Brachet cleft, which is a fibronectin-rich ECM (Plouhinec et al., 2013). In addition, several other mechanisms, such as cell lineage-based dilution (Farin et al., 2016) and cytonemes/signaling filopodia (Roy et al., 2011; Stanganello et al., 2015) may contribute to dispersal of a morphogen. We emphasize that immobilization of morphogen molecules is a prerequisite for cytoneme/filopodium-mediated transfer of signaling. Gradient formation over long ranges has not been examined experimentally in this study. However, we attempted to understand the outcome of diffusion and binding, basic properties of morphogens. Thus, we propose a mathematical model consisting of free and bound components of Wnt based on observed local dynamics (Figure 5A). This model can be widely applied to secreted proteins that bind to cell surfaces, including sFRPs and other peptide growth factors. Notably, in our mathematical model, distributions of both free and bound components converged to steady states within a few minutes, showing rather fast dynamics in the context of embryonic patterning. This characteristic could solve perceived weaknesses of diffusion-base models (Müller et al., 2013), especially dilemmas related to the speed and stability of gradient formation. From this point of view, the combination of abundant cell-surface-bound and minimal diffusing populations would be beneficial for signaling stability and speed of pattern formation, respectively. Like tethered-anti-HA Ab (Figure 2B), atypical distributions of FGF (Shimokawa et al., 2011) and Nodal (Marjoram and Wright, 2011), in which ligands accumulate in locations distant from their sources, have been reported, although a theoretical explanation of these atypical distributions has proven elusive. In our model, atypical distributions can be reproduced if specific scaffolds for ligands (ligand binding proteins) are anchored on the surfaces of cells (Figure 5E). Furthermore, our model explains the puzzling localization of ligands in tissues. In mosaic analyses of the wing discs of Drosophila mutants, Hh and Dpp ligands accumulate at the edges of clones defective in HS synthesis (Takei et al., 2004; Yan and Lin, 2009). Distributional patterns of these ligands are explained by our model, which accounts for accumulations of ligand in regions lacking HS (Figure 5F). Thus, our model provides a basic framework to understand of the extracellular behavior of secreted proteins.
Materials and methods
All experiments using Xenopus laevis were approved by the Institutional Animal Care and Use Committee, National Institutes of Natural Sciences (Permit Number 18A038, 19A062, 20A053), or the Office for Life Science Research Ethics and Safety, University of Tokyo.
Xenopus embryo manipulation and microinjection
Request a detailed protocolUnfertilized eggs of Xenopus laevis were obtained by injection of gonadotropin (ASKA Pharmaceutical). These eggs were artificially fertilized with testis homogenates and dejellied using 4% L-cysteine (adjusted to pH 7.8 with NaOH). Embryos were incubated in 1/10x Steinberg’s solution at 14–17°C and were staged according to Nieuwkoop and Faber, 1967. Synthesized mRNAs were microinjected into early (2–16 cell) embryos. Amounts of injected mRNAs are described in figure legends.
Fluorescent image acquisition
Request a detailed protocolImage acquisition was performed using confocal microscopes (TSC SP8 system with HC PL APO ×10/NA0.40 dry objective or HC PL APO2 ×40/NA1.10 W CORR water immersion objective, Leica or LSM710 system with C-Apochromat 40x/1.2 W Corr M27 water immersion objective, Zeiss). Photon counting images were acquired with a HyD detector (Leica). Detailed conditions for imaging are available upon request. mV was constructed by introducing an A206K mutation to prevent protein aggregation (Zacharias et al., 2002). For FDAP and FCS measurements, gastrula embryos were embedded on 35 mm glass-based dishes (Iwaki) with 1.5% LMP agarose (#16520–050; Invitrogen) gel, which was made of 1/10x Steinberg’s solution. For other types of live-imaging, embryos were mounted in a silicone chamber made in-house with holes 1.8 mm in diameter. Fluorescent intensity was measured using Fiji, Image J (NIH) or Zen2009 (Zeiss).
Cell lines
Request a detailed protocolHybridomas (derived from mouse, 12CA5, anti-HA; 9E10, anti-Myc) were used to obtain their total RNA and subsequent cloning of immunoglobulin genes. These hybridomas have been neither authenticated nor tested for mycoplasma because no assays were performed with these hybridomas themselves. Instead, we confirmed generation of functional anti-HA or anti-Myc IgG from the cloned genes by co-IP assay.
Immunostaining of Xenopus embryos
Request a detailed protocolImmunostaining of Xenopus embryos was carried out according to a previous report (Mii et al., 2017). Briefly, embryos were fixed with MEMFA (0.1 M MOPS pH 7.4, 4 mM EGTA, 2 mM MgSO4, 3.7% formaldehyde) 2 hr at room temperature. Fixed embryos were dehydrated with EtOH (EtOH treatment improves staining with anti-Wnt8 and anti-HS antibodies). After rehydration, embryos were washed with TBT (1x TBS, 0.2% BSA, 0.1% Triton X-100) and blocked with TBTS (TBT supplemented with 10% heat-treated [70°C, 40 min] fetal bovine serum). The following procedures are similar for primary and secondary antibodies. Antibody was diluted with TBTS and was centrifuged 15 min at 15,000 rpm before use. Embryos were incubated with the supernatant of antibody solution overnight at 4°C. Then embryos were washed five times with TBT.
cDNA cloning of IgG from cultured hybridomas
Request a detailed protocolCultured hybridomas were harvested by centrifugation and total RNAs were prepared using ISOGEN (Nippon Gene), according to the manufacturer’s protocol. First strand cDNA pools were synthesized using SuperScript II reverse transcriptase (Invitrogen) and random hexamer oligo DNA. These cDNA pools were used as templates for PCR to isolate cDNAs for heavy chains and light chains of anti-HA and anti-Myc IgGs. See Supplementary file 1 for all primers used for PCR cloning. Full-length cDNAs were cloned into the pCSf107mT vector (Mii and Taira, 2009).
Cultured hybridoma cells were harvested by centrifugation and total RNAs were prepared using ISOGEN (Nippon Gene), according to the manufacturer’s protocol. First strand cDNA pools were synthesized using SuperScript II reverse transcriptase (Invitrogen) and random hexamer oligo DNA. These cDNA pools were used as templates for PCR to isolate cDNAs for heavy and light chains of anti-HA and anti-Myc IgGs.
Procedures for PCR cloning were as follows. IgG cDNAs for 3' regions of CDSs were obtained by PCR with degenerate primers (5' γ-F and 5' κ-F) and primers corresponding to constant regions of Ig genes (γ2b-const-R, γ1-const-R, 3' κ-R) (Wang et al., 2000). To obtain the complete CDSs, 5'RACE was carried out to obtain the first codons of Ig genes, using a modified protocol in which inosines are introduced into the G-stretch of the HSPst-G10 anchor (personal communications from Dr. Min K. Park). cDNAs were synthesized with gene-specific primers (HA-heavy-R1, Myc-heavy-R1, HA-light-R1, and Myc-light-R1), and tailed with poly-(C) by terminal deoxynucleotidyl transferase, and subsequently double-stranded cDNAs were synthesized with the HSPst-G10 anchor. 5' ends of cDNAs were amplified by PCR between the HSPst adaptor and gene-specific primers (HA-heavy-R2, Myc-heavy-R2, HA-light-R2 and Myc-light-R2) using the double-stranded cDNAs as templates. Full length CDSs were amplified using primers designed for both ends of the CDSs (HH-Bam-F, MH-Bam-F, HL-bam-F, and ML-Bgl-F for 5’ ends; 3' γ2b-R, 3' γ1-R and 3' κ-R for the 3’ end) and the first cDNA pools. See Supplementary file 1 for all primers used for PCR cloning. Full-length cDNAs were cloned into the pCSf107mT vector (Mii and Taira, 2009). Sequence data for anti-HA IgG genes have been deposited in Genbank/DDBJ under accession codes LC522514 and LC522515.
FDAP measurements
Request a detailed protocolFor expression in the animal cap region of Xenopus embryos, four-cell-stage Xenopus laevis embryos were microinjected with mRNAs for mK-Wnt8 and mK-Frzb (4.0 ng/embryo) at a ventral blastomere. Injected embryos were incubated at 14°C until the gastrula stage (st. 10.25–11.5) for subsequent confocal analysis. FDAP measurements were performed using the LSM710 system (Zeiss) with a C-Apochromat ×40, NA1.2 water immersion objective. Time-lapse image acquisition was carried out for 20 s each at 25 frames/s, and after 4 s (100 frames) from the start, intercellular mK-fusion proteins were photoconverted at a small rectangular region (1.66 × 2.49 μm) with 405 nm laser irradiation. After photoconversion, images were acquired for 16 s (400 frames). Red fluorescent intensities within the rectangular region where photoconversion was performed, were analysed by curve-fitting to Equation 1 (Figure 4D) or Equation 2 (Figure 4—figure supplement 3A), using the Curve Fitting Toolbox of MATLAB (Mathworks).
FCS measurements
Request a detailed protocolFCS measurements were carried out using a ConfoCor2 system (objective: C-Apochromat ×40, NA1.2 water immersion) (Zeiss; Figure 3—figure supplement 2 only) according to a previous report (Pack et al., 2006) or a TSC SP8 equipped with FCS (objective: HC PL APO 63x/1.20 W motCORR CS2) (Leica). mRNAs for mV-Wnt8 and sec-mV were microinjected into four- or eight-cell stage Xenopus embryos. Injected embryos were measured at gastrula stage (st. 10.5–11.5). Rhodamine 6G (Sigma-Aldrich) was used to calibrate detection volume, with a reported value of its diffusion coefficient (280 μm2/s) (Pack et al., 2006). PyCorrFit software (Müller et al., 2014) was used for curve-fitting analyses of FCS data from the Leica system. Models considering three-dimensional free diffusion with a Gaussian laser profile, including a triplet component (‘T + 3D’, a one-component model or ‘T + 3D + 3D’, a two-component model) were used for fitting. Akaike information criterion (AIC) was used to compare fitting with the one-component and two-component models according to a previous report (Tsutsumi et al., 2016).
Plasmid construction
Request a detailed protocolpCSf107mT (Mii and Taira, 2009) was used to make most plasmid constructs for mRNA synthesis. pCSf107SPw-mT and pCSf107SPf-mT were constructed, which have the original signal peptides of Wnt8 and Frzb, respectively. The coding sequence (CDS) for mVenus (mV) or mKikGR (mK) was inserted into the BamHI site of pCSf107SPw-mT or pCSf107SPf-mT to construct pCSf107SPw-mV-mT, pCSf107SPf-mV-mT, pCSf107SPw-mK-mT, and pCSf107SPf-mK-mT. Constructs for SP-mV, SP-mV-HB, and SP-mV-2HA were made with pCSf107SPf-mT. pCS2 +HA-IgH-TM-2FT (the heavy chain for anti-HA IgG with the transmembrane domain of a membrane-bound form of IgG heavy chain) was made by inserting the full length CDS of heavy chain of anti-HA IgG without the stop codon (using the EcoRI and BglII sites) and a partial CDS fragment corresponding to the IgG transmembrane domain (using the BglII and XbaI sites) into the EcoRI/XbaI sites of pCS2 +mcs-2FT-T. To construct pCSf107-SP-HepIII-HA-GPI, HepIII CDS was inserted into pCSf107-SP-mcs-4xHA-GPI.
Luciferase reporter assays
Request a detailed protocolLuciferase reporter assays were carried out as previously described (Mii and Taira, 2009). Multiple comparisons were carried out with pairwise Wilcoxon rank sum test (two-sided) in which significance levels (p-values) were adjusted by the Holm method, using R.
Mathematical modeling
Request a detailed protocolTwo ligand components were considered: free and bound. This model includes five dynamic processes: (i) production of ligands, (ii) diffusion of the free component in intercellular space, (iii) binding of ligands to dotted structures (‘docking sites’) such as HS clusters on the surface of cells, (iv) release of the bound component from ‘docking sites,’ and (v) internalization of the bound component into cells. The model in one-dimensional space is written as:
where u and v represent the amounts of free and bound components of morphogen molecules, respectively. The symbols a(x), b, c, and g(x) represent binding, release, internalization, and production rates, respectively. D represents the diffusion coefficient of the free component in the extracellular space. The ligand is assumed to be produced in a limited region using the following function:
We assumed that the binding rate a(x) depends on the position x, following heterogeneous distribution of HS clusters on the cell surface. The following function was used for a(x):
where p1 and p2 are the interval and width, respectively, of docking sites. We used the no-flux (Neumann) boundary condition at x = 0 and L. We calculated the model by numerical simulation. The initial distributions of u and v were set at 0 throughout the entire space. In the one-dimensional space, distributions of free (u) and bound (v) components of secreted proteins were obtained by computer simulation, where the spatial length L = 1000 (μm). Distributions of u (red) and v (blue) are presented at time t = 100 (sec), which almost reached steady states. We used the forward difference method with the spatial step Δx = 0.1 and temporal step Δt = 0.0001 in numerical calculations. In Figure 5B, parameter values are: D = 20.0, an,max = 10.0, b = 0.1, c = 0.1, gmax = 0.2, R = L/1000, p1 = 2, and p2 = 0.2. In other panels, distributions in specific conditions are shown (see figure legends).
Appendix 1
Detecting freely diffusing molecules with confocal microscopy
When we consider detection of diffusing molecules with confocal microscopy, diffusing molecules might not be detected because of their large displacement. In our imaging conditions, pixel dwell time was 24.4 μs and observed maximum Dfast in FCS was about 80 μm2/s (Figure S3E). During pixel dwell time t, mean square displacement (MSD) of molecules of diffusion coefficient D in three-dimensional space is 6Dt (Crank, 1975). Therefore, maximum MSD is estimated as
Accordingly, mean displacement of the molecules at the maximum is
On the other hand, when the numerical aperture (NA) of an objective lens is given, the Reyligh diffraction limit in our conditions is calculated as
This value of the diffraction limit is larger than the so that even the most rapidly diffusing molecules observed in FCS measurements, loss of fluorescence due to diffusion is not likely to occur. Indeed, we detected slight but significant increase of photon counts with sec-mV (Figure 1D and Figure 1—figure supplement 2) with a long pixel dwell time. We speculate that the difficulty of visualizing sec-mV is mainly due to narrowness of extracellular space in Xenopus embryos
Data availability
Sequence data for anti-HA IgG genes have been deposited in Genbank/DDBJ under accession codes LC522514 and LC522515.
-
NCBI GenBankID LC522514.1. Mus musculus mRNA for immunoglobulin gamma 2b, complete cds.
-
NCBI GenBankID LC522515.1. Mus musculus mRNA for immunoglobulin kappa, complete cds.
References
-
Isolation of monoclonal antibodies specific for human c-myc proto-oncogene productMolecular and Cellular Biology 5:3610–3616.https://doi.org/10.1128/MCB.5.12.3610
-
Purification of a RAS-responsive adenylyl cyclase complex from Saccharomyces cerevisiae by use of an epitope addition methodMolecular and Cellular Biology 8:2159–2165.https://doi.org/10.1128/MCB.8.5.2159
-
On the importance of protein diffusion in biological systems: the example of the bicoid morphogen gradientBiochimica Et Biophysica Acta (BBA) - Proteins and Proteomics 1865:1676–1686.https://doi.org/10.1016/j.bbapap.2017.09.002
-
Mechanisms for positional signalling by morphogen transport: a theoretical studyJournal of Theoretical Biology 191:103–114.https://doi.org/10.1006/jtbi.1997.0575
-
Investigating the principles of morphogen gradient formation: from tissues to cellsCurrent Opinion in Genetics & Development 22:527–532.https://doi.org/10.1016/j.gde.2012.08.004
-
Selective activation mechanisms of wnt signaling pathwaysTrends in Cell Biology 19:119–129.https://doi.org/10.1016/j.tcb.2009.01.003
-
A syngeneic monoclonal antibody to murine Meth-A sarcoma (HepSS-1) recognizes heparan sulfate glycosaminoglycan (HS-GAG): cell density and transformation dependent alteration in cell surface HS-GAG defined by HepSS-1Journal of Immunology 137:3900–3908.
-
Heparan sulfate proteoglycan clustering in wnt signaling and dispersalFrontiers in Cell and Developmental Biology 8:631.https://doi.org/10.3389/fcell.2020.00631
-
Morphogen gradients: from generation to interpretationAnnual Review of Cell and Developmental Biology 27:377–407.https://doi.org/10.1146/annurev-cellbio-092910-154148
-
Mechanisms of intercellular wnt transportDevelopment 146:dev176073.https://doi.org/10.1242/dev.176073
-
Fiji: an open-source platform for biological-image analysisNature Methods 9:676–682.https://doi.org/10.1038/nmeth.2019
-
Forming and interpreting gradients in the early Xenopus embryoCold Spring Harbor Perspectives in Biology 1:a002477.https://doi.org/10.1101/cshperspect.a002477
-
Analysis of binding reactions by fluorescence recovery after photobleachingBiophysical Journal 86:3473–3495.https://doi.org/10.1529/biophysj.103.026765
-
FRAP analysis of binding: proper and fittingTrends in Cell Biology 15:84–91.https://doi.org/10.1016/j.tcb.2004.12.001
-
Filopodia-based wnt transport during vertebrate tissue patterningNature Communications 6:5846.https://doi.org/10.1038/ncomms6846
-
Wingless gradient formation in the Drosophila wingCurrent Biology 10:293–300.https://doi.org/10.1016/S0960-9822(00)00378-X
-
Morphogens, their identification and regulationDevelopment 131:703–712.https://doi.org/10.1242/dev.01043
-
Assembly of protein complexes restricts diffusion of Wnt3a proteinsCommunications Biology 1:165.https://doi.org/10.1038/s42003-018-0172-x
-
Design of peptides with high affinities for heparin and endothelial cell proteoglycansJournal of Biological Chemistry 275:7701–7707.https://doi.org/10.1074/jbc.275.11.7701
-
Shaping morphogen gradients by proteoglycansCold Spring Harbor Perspectives in Biology 1:a002493.https://doi.org/10.1101/cshperspect.a002493
-
Incredible journey: how do developmental signals travel through tissue?Genes & Development 18:2985–2997.https://doi.org/10.1101/gad.1233104
Article and author information
Author details
Funding
Japan Science and Technology Agency (JPMJPR194B)
- Yusuke Mii
Japan Society for the Promotion of Science (24870031)
- Yusuke Mii
Japan Society for the Promotion of Science (15K14532)
- Yusuke Mii
Japan Society for the Promotion of Science (18K14720)
- Yusuke Mii
National Institutes of Natural Sciences (1311608)
- Yusuke Mii
National Institutes of Natural Sciences (01311801)
- Yusuke Mii
Japan Society for the Promotion of Science (24657147)
- Yusuke Mii
- Masanori Taira
Japan Society for the Promotion of Science (18H02447)
- Masanori Taira
Japan Society for the Promotion of Science (25251026)
- Masanori Taira
Japan Science and Technology Agency (JPMJCR13W6)
- Atsushi Mochizuki
Japan Science and Technology Agency (JPMJCR1922)
- Atsushi Mochizuki
Japan Society for the Promotion of Science (19H05670)
- Atsushi Mochizuki
Japan Society for the Promotion of Science (17K19418)
- Shinji Takada
Japan Society for the Promotion of Science (18H02454)
- Shinji Takada
Japan Society for the Promotion of Science (19H04797)
- Shinji Takada
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank Dr. Min Kyun Park for hybridoma culture and 5’RACE for IgG cDNAs. We also thank Dr. A Miyawaki for Venus cDNA; Dr. R Moon for pCS2+Xwnt8; Dr. M Affolter for morphotrap cDNA; Dr. W Hashimoto for the HepIII plasmid (pET21b-Phep_3797); Dr. Shinya Matsuda for critical reading; Drs. Hiroshi Koyama, Yohei Kondo and Motosuke Tsutsumi for discussion and Dr. Steven D Aird for editing and proofreading. We thank the Confocal Microscope core facility at the ConveRgence mEDIcine research cenTer (CREDIT), Asan Medical Center. This work was supported in part by following programs: KAKENHI (24870031, 15K14532, 18K14720 to YM; 19H05670 to AM; 24657147, 18H02447, 25251026 to MT; 17K19418, 18H02454, 19H04797 to ST), the NINS program for cross-disciplinary study (1311608, 01311801 to YM), Joint Usage/Research Center program of Institute for Frontier Life and Medical Sciences Kyoto University, JST-PRESTO (JPMJPR194B to YM), JST-CREST (JPMJCR13W6, JPMJCR1922 to AM).
Ethics
Animal experimentation: All experiments using Xenopus laevis were approved by the Institutional Animal Care and Use Committee, National Institutes of Natural Sciences (Permit Number 18A038, 19A062, 20A053), or the Office for Life Science Research Ethics and Safety, University of Tokyo.
Copyright
© 2021, Mii et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 2,047
- views
-
- 268
- downloads
-
- 14
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Developmental Biology
Neuronal stem cells generate a limited and consistent number of neuronal progenies, each possessing distinct morphologies and functions, which are crucial for optimal brain function. Our study focused on a neuroblast (NB) lineage in Drosophila known as Lin A/15, which generates motoneurons (MNs) and glia. Intriguingly, Lin A/15 NB dedicates 40% of its time to producing immature MNs (iMNs) that are subsequently eliminated through apoptosis. Two RNA-binding proteins, Imp and Syp, play crucial roles in this process. Imp+ MNs survive, while Imp−, Syp+ MNs undergo apoptosis. Genetic experiments show that Imp promotes survival, whereas Syp promotes cell death in iMNs. Late-born MNs, which fail to express a functional code of transcription factors (mTFs) that control their morphological fate, are subject to elimination. Manipulating the expression of Imp and Syp in Lin A/15 NB and progeny leads to a shift of TF code in late-born MNs toward that of early-born MNs, and their survival. Additionally, introducing the TF code of early-born MNs into late-born MNs also promoted their survival. These findings demonstrate that the differential expression of Imp and Syp in iMNs links precise neuronal generation and distinct identities through the regulation of mTFs. Both Imp and Syp are conserved in vertebrates, suggesting that they play a fundamental role in precise neurogenesis across species.
-
- Computational and Systems Biology
- Developmental Biology
The initially homogeneous epithelium of the early Drosophila embryo differentiates into regional subpopulations with different behaviours and physical properties that are needed for morphogenesis. The factors at top of the genetic hierarchy that control these behaviours are known, but many of their targets are not. To understand how proteins work together to mediate differential cellular activities, we studied in an unbiased manner the proteomes and phosphoproteomes of the three main cell populations along the dorso-ventral axis during gastrulation using mutant embryos that represent the different populations. We detected 6111 protein groups and 6259 phosphosites of which 3398 and 3433 were differentially regulated, respectively. The changes in phosphosite abundance did not correlate with changes in host protein abundance, showing phosphorylation to be a regulatory step during gastrulation. Hierarchical clustering of protein groups and phosphosites identified clusters that contain known fate determinants such as Doc1, Sog, Snail, and Twist. The recovery of the appropriate known marker proteins in each of the different mutants we used validated the approach, but also revealed that two mutations that both interfere with the dorsal fate pathway, Toll10B and serpin27aex do this in very different manners. Diffused network analyses within each cluster point to microtubule components as one of the main groups of regulated proteins. Functional studies on the role of microtubules provide the proof of principle that microtubules have different functions in different domains along the DV axis of the embryo.