Modeling and analysis of collective cell migration in an in vivo three-dimensional environment
- aMolecular, Cellular and Developmental Biology, University of California Santa Barbara, CA 93106;
- bDepartment of Biological Chemistry, The Johns Hopkins School of Medicine, Baltimore, MD 21205;
- cDepartment of Biological Sciences, Indian Institute of Science Education and Research Kolkata, West Bengal 741252, India;
- dDepartment of Chemical Physics, The Weizmann Institute of Science, Rehovot 76100, Israel
See allHide authors and affiliations
Edited by Peter N. Devreotes, The Johns Hopkins University School of Medicine, Baltimore, MD, and approved March 1, 2016 (received for review November 17, 2015)

Significance
Many cells travel within an organism to perform important duties in development, homeostasis, and disease. Questions remain as to why cells migrate in diverse forms: some travel as single cells, whereas others move collectively in groups of different sizes. What is also unclear is how extracellular conditions influence migration. In our study, which combines live imaging with theoretical modeling, we address these questions and reveal how both physical and chemical features of the microenvironment influence the optimum size of a migrating cell group.
Abstract
A long-standing question in collective cell migration has been what might be the relative advantage of forming a cluster over migrating individually. Does an increase in the size of a collectively migrating group of cells enable them to sample the chemical gradient over a greater distance because the difference between front and rear of a cluster would be greater than for single cells? We combined theoretical modeling with experiments to study collective migration of the border cells in-between nurse cells in the Drosophila egg chamber. We discovered that cluster size is positively correlated with migration speed, up to a particular point above which speed plummets. This may be due to the effect of viscous drag from surrounding nurse cells together with confinement of all of the cells within a stiff extracellular matrix. The model predicts no relationship between cluster size and velocity for cells moving on a flat surface, in contrast to movement within a 3D environment. Our analyses also suggest that the overall chemoattractant profile in the egg chamber is likely to be exponential, with the highest concentration in the oocyte. These findings provide insights into collective chemotaxis by combining theoretical modeling with experimentation.
The ability to sense and follow directional signals is essential for migrating cells. Gradients of chemical signals are believed to guide moving cells to their targets (1, 2). Classically, the formation of such a gradient involves a “source” that continuously produces a signal that freely diffuses through extracellular space and a “sink” that actively eliminates it (3⇓⇓–6). To study cellular responses, investigators have artificially created chemical gradients in vitro using gels (7, 8), micropipettes (9), and microfluidic devices (10⇓–12). However, endogenous gradients are challenging to measure directly. Recently, endogenous promoter-driven fluorescence fusion proteins have been used to visualize a chemical gradient in vivo (13); however, such reporters usually give weak signals and the fluorescent tag can alter the molecule’s properties. Even when a signal can be detected, it reveals the bulk of the chemical, whereas the protein that is available for the migrating cell to sense may be a tiny, invisible fraction of the total (13). Moreover, in vivo, cells probably integrate information from multiple signals (6, 14), both biochemical and physical. Both the signal and the response may be quite different in 3D in vivo environments compared with 2D experimental paradigms. For all these reasons, it is important to analyze migrating cells in their native environments.
It is becoming increasingly apparent that, even though some cells migrate as individuals, many cells move in interconnected clusters, strands, or sheets (15). Clusters of cells are larger than individual cells and at least theoretically should be able to sample a chemical gradient over a greater distance, so that the difference between front and rear of a cluster would be more pronounced than for single cells. Thus, a large cluster of cells could in principle be more sensitive to weaker signals, more polarized, and thus migrate more effectively. However, existing experimental data from neural crest explants of different sizes showed that migration speed across a flat surface was unaffected by cluster size (16). In vivo, collectively migrating cells move through a 3D microenvironment, raising the possibility that a size–speed correlation, if it exists, might be missed in artificial cell culture environments. Moreover, if any size–speed correlation does exist, it remains to be seen whether it is positive or negative.
Theoretical modeling has been instrumental in probing basic questions in developmental biology for decades (17⇓–19). It can reveal insights that are counterintuitive or difficult to ascertain experimentally. In this study, we combined live imaging of a simple and well-studied model of collective cell migration in the Drosophila ovary—the border cells—with theoretical modeling to investigate the chemical and physical features that influence their migration speed. We found experimentally, in vivo, a size–speed correlation that supports the theoretical prediction that larger clusters move faster in 3D but not on 2D surfaces. Moreover, our analysis suggests that the sum total of all chemoattractants that the border cells respond to, which includes at least four known ligands for receptor tyrosine kinases (RTKs), likely forms an exponential gradient. Finally, we report that the external physical constraints of the environment ultimately limit the increase in speed conferred by increasing cluster size. Thus, multiple features of the in vivo environment, both physical and chemical, exert effects on the speed and optimal size of a migrating group of cells.
The Drosophila ovary is made up of 12–16 ovarioles, each of which serves as an “assembly line” of developing egg chambers (Fig. 1A). An egg chamber is composed of 16 germ cells surrounded by a single layer of somatic epithelial cells. A pair of specialized cells called polar cells develops at each end of each egg chamber and secretes Unpaired (Upd), a cytokine that activates JAK/STAT signaling in the surrounding epithelial follicle cells. At the anterior end, during developmental stage 9, these four to eight so-called border cells that receive high levels of JAK/STAT become migratory and carry the polar cells, which are not mobile, in-between the germ-line cells called nurse cells, until they reach the oocyte (20⇓–22) (Fig. 1B and Movie S1). During migration, the cells respond to at least four secreted factors produced in the germ line: PDGF- and VEGF-related factor 1 (PVF1), which binds to and activates its receptor (PVR), as well as Spitz, Keren, and Gurken (Grk), which bind and activate epidermal growth factor receptor (EGFR), which like PVR is a receptor tyrosine kinase. Inhibiting both receptors in the border cells at the same time causes the cells to migrate poorly and sometimes off-track (23⇓⇓–26). It would be interesting to directly measure the spatial profiles of the ligands; however, antibodies are only available against PVF1 and Grk, and the extracellular protein is not abundant enough to detect. Therefore, we have no direct information about the shapes of the extracellular gradients or even if there are gradients at all.
Drosophila ovary and border cell migration simulations. (A) Drosophila ovariole containing egg chambers of increasing maturity. (B) Stage 9 egg chamber showing border cells at onset of migration. Border cells are in blue. (C) Schematic illustration of our model of chemotaxis of cellular clusters: due to the chemoattractant gradient, the cells at the cluster surface protrude and induce a force that is directed at the local outward normal, and is increasing with the local concentration of the chemical signal (black arrows). (A–C) Polar cells are in purple. (D) Schematic illustration of the drag sources on the border cell cluster (black circle) as it moves along the axis of the egg. One mechanism of drag is due to direct molecular binding with the surrounding cells (red rim of the cluster), and the second mechanism comes from the cytoplasmic flow around the cluster (streamlines with arrows).
Here, we developed a theoretical model that predicted that larger clusters would migrate faster than smaller ones in 3D but not in 2D. We then generated clusters of different sizes, both smaller and larger than those found in wild type (WT), and evaluated their migration speeds. Our results reveal that migration speed increases with increasing cluster size up to a point beyond which speed decreases precipitously. The drop in speed suggested that, in this confined 3D tissue, cluster speed is likely greatly affected by viscous drag and led us to refine the model. We also compared the observed changes in speed along the migration path to those predicted by simple mathematical models of responses to either linear or exponential gradients. Together, the combination of empirical observations and theoretical modeling emphasize the importance of analyzing cells in realistic microenvironments and show that both chemical and mechanical features influence the optimal size of a collectively migrating cluster in vivo.
Theoretical Model
During chemotaxis, cells migrate directionally in response to chemoattractants, moving up or down spatial gradients. However, because the functional chemoattractant molecules that are actually available to the migrating cell usually constitute only a small fraction of the total ligands, and multiple ligands may impinge simultaneously on a given cell or cell group, the shape of chemoattractant gradient seen by the cell, in most systems is unknown. In the border cell system, PVF1 and three ligands for the EGF receptor—Spitz, Keren, and Gurken—may form either exponential or linear gradients, or some combination, to guide border cell migration. An exponential profile would most likely form by continuous ligand secretion by the oocyte and concomitant removal of ligands by surrounding follicle cells, all of which also express PVR and EGFR.
To model border cell migration, we assume for simplicity that each cluster has a spherical shape, whose surface is made up of a single layer of border cells. Each cell at the cluster surface produces a pulling force that points along the outward normal to the sphere, and is an increasing function of the local concentration (c) of the chemoattractant (Fig. 1C). This pulling force is exerted by the border cells on the surrounding nurse cells, through the extension of lamellipodia and filopodia, which adhere to the surrounding matrix and then contract, thereby pulling the border cells forward. The relation is linear for low concentrations, whereas at high concentrations we expect this relationship to saturate at a maximal cellular pulling force. The force per unit area of cell membrane in contact with the chemoattractant, depends on the response function S(c), which is the fraction of activated receptors at the cell surface, and for a simple process (Supporting Information, Eqs. S1–S6, and see Table S1 for full list of variables and parameters used in the model), we get as follows:
List of variables
Behavior of the resultant force (Eq. 4) for a linear chemical profile. The y axes are force in log scale. (A) As a function of the maximal concentration cm. We find a linear increase of the force for shallow gradients, and a saturation-driven slowing down for steep gradients [calculated at a fixed position (
This force (Eq. 4) can be approximated by a simpler expression, for example by expanding Eq. 4 for large
An exponential concentration profile has the following form:
(A) Calculated velocity (purple line, from Eq. S9) of a cluster moving up an exponential chemical gradient (blue line). The response function S(x) is shown by the yellow line, and we see that the peak in the cluster velocity corresponds to the steepest part of the response (using: c0 = 0, cm = 10,
The center-of-mass velocity of the cluster is determined by the following Langevin equation (in the limit of vanishing inertia), where the effective friction force balances the pulling forces:
The total drag is the combination of these two dissipation mechanisms, so that we have the following: λ = λad + λvisc. Due to the finite lateral confinement of the hard wall of the egg chamber, the viscous drag increases above the simple Stokes’ law when the ratio R/Regg approaches unity. This effect is usually described in terms of a polynomial expansion of this ratio (27):
Note that our system is different from the scenario of a 2D, in vitro layer of cells where the viscous drag from the infinite overlying buffer is indeed negligible compared with the effective contact friction with the substrate (28). In the confined 3D volume of the egg, it is not clear a priori whether the viscous contribution is important, which is why we write the model in a general form containing both friction contributions. It turns out (see below) that the observed scaling strongly supports the conclusion that the drag from viscous flow is dominant in this system.
From Eqs. 3 and 8, we see that the average velocity in the direction of the chemical gradient increases with the cluster radius: linearly when surface contacts dominate ⇒
We next want to compute the velocity fluctuations of the clusters. The data on the fluctuations of the cluster velocity in the direction orthogonal to the chemical gradient (y), so presumably unaffected by it, indicate a rather constant background noise level that is inherent to the cluster motion. This noise may be described as an additive noise source, inherent to the pulling forces that each cell at the cluster surface produces. There are two limiting cases in which the random forces induced by each cell are (i) uncorrelated or (ii) correlated with its neighbors. Cell–cell correlations are evident in collective cellular motion (30, 31). We will characterize the random active forces induced by the individual cells as having a typical force f0, a mean burst duration τnoise (persistence time of the cellular pulling force), and mean waiting time between bursts τ. We then solve the Langevin equation for the center-of-mass motion of a cluster (32), in the limit of vanishing inertial effects λτnoise ≫ 1. We therefore get the following:
From the previous section, we found that the drag has two contributions: For contact dominated friction λ = λ0Ns ∝ R2, we get the following:
Results
Dependence of Speed on Cluster Size.
We predict that the resultant pulling force due to the chemical gradient (
To test the prediction that migration speed is positively correlated with cluster size in a 3D environment, we measured the average diameter and average speed of border cell cluster in each movie, and plotted diameter–speed data from 18 different WT movies (Fig. 2A, square). To obtain even smaller clusters, we knocked down E-cadherin specifically in central nonmotile polar cells, which causes clusters to break apart (29). Compared with other methods that generate smaller border cell clusters (34, 35), this method does not directly alter the migrating cells themselves, thus presumably minimally changing their properties, and the smaller clusters generated are able to migrate. We measured average speed and average cluster size from 10 UpdGal4, UAS-EcadRNAi movies (Fig. 2A, circles), and combined them with 18 WT movies (Fig. 2A, squares). Our data revealed that cluster speed increases up to a diameter of 24 µm (R = 12 µm), beyond which the velocity saturates (Fig. 2A). The increase of speed with size is found to be consistent with the predicted quadratic behavior, thereby indicating that the dominant mode of friction is most likely due to the viscous (Stokes) flow of the nurse cells’ cytoplasm around the cluster.
Dependence of mean cluster velocity on the cluster radius. (A) Symbols show data points averaged from different experiments: Circles, E-cadherin k.d. (UpdGal4, UAS-EcadRNAi); squares, WT; diamonds, mutant extralarge clusters. Solid black line gives the fit to Eq. 16, using
The data also had an indication that the velocity saturates for the largest clusters in this dataset. This led us to reason that due to the confinement of the rigid extracellular matrix (ECM), if the clusters increase further in size they approach the limit of a “plug” where the flow around the cluster is highly inhibited (Fig. 1D), and the motion is blocked. This phenomenon is described by the following general expression that combines the two drag mechanisms (Eqs. 8 and 9):
Border cells move in a confined space delineated by nurse cells and the surrounding ECM. Moreover, the huge nurse cell nuclei further decrease the effective radius of the rigid confinement Regg, such that it is much smaller than the radius of the actual egg chamber (which is ∼50–100 µm depending on the position within the egg chamber). When the border cell cluster is small (i.e., diameter ≤ 24 µm), the ability to sample a larger distance within chemical gradient causes faster migration for larger clusters. However, as the cluster becomes larger (i.e., diameter > 24 µm), the confinement of the viscous flows in the surrounding nurse cells increases sharply. Therefore, further increases in cluster size have a negative impact on speed (Movie S2). Note that extra large border cell clusters tend to be more elongated than WT and polar cell Ecad k.d. clusters. However, average velocity does not correlate strongly with shape (Fig. 2 D and E). We used y-axis diameter in calculations because x-diameter elongation could be affected by thin extensions of cells between the surrounding nurse cells, and these extensions do not conform to the treatment of drag arising from bulk viscous flow. Indeed, for extralarge clusters, y-axis diameter (Fig. 2C) correlates better than x-axis diameter (Fig. 2B) with average cluster velocity.
Fluctuation in Migration Speed.
The model predicts that, if the individual noise in the cells’ pulling forces is uncorrelated, and friction is dominated by cluster surface contacts, then the velocity variance in the direction perpendicular to the chemical gradient should decrease as follows:
(A) Dependence of mean cluster velocity on cluster radius, with SD shown. (B) Cluster velocity in relation to cluster position along the migration path. (C) Variance in the y-direction velocity
Shape of the Chemoattractant Gradient.
To understand whether the overall chemical gradient in the egg chamber is linear or exponential, we extracted migration speed data from WT border cell migration movies, and plotted speed against distance traveled, expressed as a fraction of the total migration path (Fig. 3A and Fig. S4). Averaging 20 different WT movies, we obtained a percentage migration speed line shown in solid line. Each individual movie shows decrease of cluster migration speed to around zero at the end of migration, consistent with our hypothesis that, as border cells migrate up the chemical gradient, their guidance receptors saturate, causing them to lose front–back polarity. However, instead of a steady decrease in speed, border cells accelerate first after they start migrating, and then decelerate as they move toward the oocyte, an effect that resembles the theoretical response to an exponential chemoattractant gradient model, as shown in the dashed line (Fig. 3A). This initial increase in speed does not result from border cells escaping from the anterior end of the egg chamber, because in the majority of movies where the acceleration occurred at the onset of migration, detachment from the anterior end had already taken place before the acceleration (Fig. S5C). Therefore, the overall chemoattractant gradient likely increases exponentially with highest concentration at the oocyte. We also plot a fit to the motion in a linear chemical profile (Fig. 3D, red dashed line), which monotonously decreases and does not reproduce the peak in velocity observed at 20–30% of the trajectory. However, the noise in the measured motion, and in particular the large oscillations in the motion due to several slow phases (“pauses”) (Fig. 3 and Fig. S4), does not allow us to definitely discriminate between the exponential and linear profiles.
Migration profiles of WT and receptor-deficient clusters. Migration speeds at different positions along migration path in WT (A), slboGal4; EGFRDN (B), and slboGal4; PVRDN (C) egg chambers. Symbols of same color and shape indicate data points from the same movie. Solid lines give the average of all of the experiments. Dashed line gives the fit to Eq. 17. (D) Theoretical fits to WT (red), EGFRDN (green), and PVRDN (black) migration speed using exponential (solid lines) and linear (dashed lines) chemoattractant concentration profile.
Velocity up the gradient, for several experiments. Colored lines are individual WT (A), EGFRDN (B), and PVRDN (C) experiments, averaged every five time points for smoothening. The thick black line gives the average of all of the experiments, and the thick dashed-dot black line gives the fit to Eq. 17. The peak value of the theoretical expression is normalized to velocity of 1.2 µm/min.
Nurse cell influence on border cell migration. (A) Illustration of border cell cluster untrapping from a nurse cell–nurse cell junction. NC, nurse cell. (B) Time spent untrapping from a junction in relation to cluster size. Linear regression was done and linear plot was shown in the graph. (C) Percentage of border cell migration movies that show an initial increase in speed that have detachment (blue bar) and no detachment (red bar) involved.
We next want to understand the relative contributions of different ligands to the overall chemoattractant profile. The migrating border cells express both PVR and EGFR. We expressed a dominant-negative (DN) form of EGFR (EGFRDN) (24) in border cells, so that they can only respond to PVF1 ligand. We plotted speed vs. percentage of migration from 10 different experiments and averaged all of the experimental data (solid line in Fig. 3B). A dashed black line gives the model calculation using an exponential gradient with these parameters:
Discussion
Border cell migration in the Drosophila ovary has been an excellent model to study collective cell migration in a native, 3D environment. Having identified appropriate culture and imaging techniques, we are now able to capture the entire 4- to 6-h process live and obtain a wealth of data. In this paper, we combined theoretical modeling with experimental measurements to answer an important question in the field of collective cell migration: which features of the microenvironment influence the speed and optimal size of a collectively migrating group of cells?
In this study, we discovered a correlation between size and speed of migrating clusters to address a long-standing question. Does the increase in size in collectively migrating cells help them sense the gradient better and achieve more efficient migration? A recent study of Xenopus neural crest cells reported no size–speed correlation in migrating explants observed in 2D culture (16). In contrast, we find that size matters for migrating border cell clusters. One possible explanation for the discrepancy is that the neural crest cells were analyzed as they moved on a 2D surface, although their normal migration occurs in a 3D environment. We analyzed border cell migration in their natural 3D environment as they move in-between nurse cells. Our model predicts just this difference between 2D and 3D. Other studies have reported key differences between cellular behaviors and molecular mechanisms in 2D vs. 3D (38). It will be interesting to determine whether, in other 3D chemotactic cell migration models, larger clusters move faster as predicted by our model. Another study recently undertook a more complex model of border cell migration, which also predicted increasing speed with increasing size within the normal range of four to eight cells (39). It is interesting that these complex simulations predict qualitatively the same behavior as our simple and thus analytically tractable model. This makes our experimental validation of the predictions satisfying. We go even further to predict, and then to show, that this relationship only holds for a particular range of cluster sizes.
The benefits of increasing cluster size do not increase without limit. For border cells, speed increases until cluster reaches 24 µm. Further increase in cluster size leads to a precipitous decrease in speed, attributable to stronger viscous drag due to confinement by the surrounding nurse cells. Potential caveats exist because we assume that the border cell cluster is a perfect sphere. The influence of viscous drag could be reduced for larger clusters if they become elongated rather than spheroid. This may explain the observed elongated morphology of, for example, the zebrafish lateral line primordium (40). Another solution to reducing the influence of viscous drag was recently described for cranial neural crest cells, which migrate collectively as a large group. Instead of splitting into small clusters, these cells reduce the strength of cell–cell adhesions so that they can move in a more fluid-like manner (41). Note that we find that chemotactic migration in 3D mimics the kinematic behavior of objects falling under gravity or bubbles rising through a confined column of viscous medium. Although the sources of the forces are different, the resulting behavior is strikingly similar.
We also used the measured speed along the migration path to deduce that the overall chemoattractant profile in the egg chamber is more likely to be exponential, with highest concentration in the oocyte, than linear. This is difficult to measure directly because (i) antibodies only reveal the bulk of the chemical, whereas the protein that is available for the migrating cell to sense may be a tiny, invisible fraction of the total; and (ii) border cells integrate responses to multiple chemical signals. Moreover, experimental manipulation using PVRDN and EGFRDN lines reveals the relative contributions of PVR and EGFR ligands to the overall chemical profile: PVF1 forms a sharp exponential gradient, whereas EGFR ligands Gurken, Keren, and Spitz exhibit a shallower exponential with a much longer decay length. One interpretation of these data is that they interact with lower affinity and thus are removed more slowly.
Other studies have compared linear vs. exponential chemoattractant gradients in cell migration. For example, using microfluidic devices, Herzmark et al. (42) generated linear and exponential f-Met-Leu-Phe (fMLP) gradients and watched the chemotactic responses of neutrophil-differentiated HL60 cells in these gradients. Similar to our modeling, they concluded that HL60 cells chemotax best at the low end of a linear concentration gradient, whereas in an exponential gradient, HL60 cells show highest chemotactic migration velocity in the middle of migration (42). However, it remains uncertain what chemoattractant profiles cells encounter in vivo. Our finding here in the Drosophila border cell migration system is an in vivo example that suggests the existence of multiple exponential chemoattractant gradients.
Our results raise the question as to why the sizes of WT clusters are poised close to the highest migration velocities (Fig. 2A). We wonder if there is an advantage to rapid movement. The complete process of oogenesis requires several days so it is not immediately clear why border cell migration speed would need to be optimized. However, migrating border cell clusters encounter nurse cell junctions along the migration path, and can get trapped in those junctions (see the dips in migration speed in Fig. S4). A simple analysis of the forces in such junctions (Supporting Information) suggests that it is necessary for a migrating cluster to achieve certain size so that chemotactic forces are able to overcome the osmotic pressure from nurse cells, and pull the cluster out of nurse cell junctions (Fig. S5). Indeed, we see that, when a border cell cluster gets larger, it takes less time to escape nurse cell junctions (Fig. S5). Evolutionarily, this might be another reason that larger cluster size is selected and suggests that multiple physical constraints of a tissue exert powerful influences on optimal cluster size and migration speed.
Examples of collective cell migrations abound in normal development, in tissue homeostasis, and in diseases, most notably cancer (reviewed in ref. 43). Collective migration may be the rule rather than the exception in metastasis, despite the physical challenges posed by squeezing groups of cells into blood and lymphatic vessels (44). Additional examples of collective migrations include blood vessel sprouting, epithelial closures, and innumerable embryonic morphogenetic movements. Cells move in groups of diverse sizes and varying arrangements through equally complex microenvironments. The border cells represent a clear example of how collective movement can enhance direction sensing (29) and how the combined chemical and physical properties of the microenvironment influences the optimal cluster size. The results presented here therefore lead us to propose that the details of the diverse physical and chemical microenvironments likely exert a strong influence on the observed sizes and arrangements of collectively migrating cells in vivo.
Materials and Methods
Fly Stocks and Egg Chamber Live Imaging.
SlboGal4, UAS-dsRedNLS; UAS-moesinGFP, and slbo-Lifeact-GFP flies were used for live imaging and were described before (26, 29). UAS-PVRDN, UAS-EGFRDN, and UAS-EcadRNAi (VDRC: KK103962) were described previously (23, 24, 26, 29). Fly fattening, dissection, and live imaging follow the established protocol (45).
Theoretical Modeling.
Please see Supporting Information for the details of theoretical modeling.
Registration of Egg Chambers.
To correct for egg chamber movement and/or growth during imaging, we measured border cell speed relative to the center of egg chamber by registering the outlines from frame to frame. We first used StackReg plugin (rigid body) in Fiji [ImageJ 1.49c; National Institutes of Health (46)] to align the egg chambers. Then we designed a Matlab code to further align the egg chambers (Supporting Information).
SI Theoretical Modeling
Calculated Response to the External Signal: Constant Receptor Recycling Rate.
Consider membrane receptors that are in an empty form (ρn), and those that got attached to the signaling molecule (ρa). We consider endocytic activity that removes both receptors indiscriminately, with constant rate Γ.
The equations of motion are as follows:
Exponential Concentration Profile.
First, how does an exponential concentration gradient arise naturally in the organism? When a chemical is produced at one end of a (roughly) 1D object, and then diffuses from that source (with diffusion coefficient D) while being degraded at a constant-rate γ everywhere, we end up with the following exponential concentration profile:
Let us explore the behavior of the velocity in the exponential gradient. For the location of the peak, we can find an approximate expression:
Similarly, we show the effect of changing the value of the parameter ε in Fig. S2C.
Cluster Untrapping While Migrating In-Between Nurse Cells.
Migrating border cells frequently (two to approximately three junctions per entire migration) encounter nurse cell–nurse cell junctions that can slow down the migration and trap border cells for several minutes up to an hour (Fig. S4). WT border cells usually untrap themselves (Fig. S4).
We can use a simple calculation to analyze this problem (Fig. S5B). The energy cost of the untrapping move can be estimated to be proportional to the work done pushing the surrounding cells:
The critical radius
SI Materials and Methods
To calculate the nurse cell untrapping time, SlboGal4, UAS-dsRedNLS; UAS-moesinGFP movies were taken using transmitted light detector (transmission–photomultiplier tube) to visualize nurse cell boundaries, and GFP/RFP GaAsP spectral detectors on a Zeiss LSM 780 microscope. Nurse cell junction is defined as a point along the border cell migration path where at least three nurse cell boundaries intersect, and is verified over multiple stacks in z direction. The process of border cells untrapping from nurse cell junction starts from when front of the cluster first touches the junction and ends when back of the cluster leaves it. We only looked at the junctions where border cell migration speed significantly decreased (26.3% of all junctions). We first masked the GFP channel with surface in Imaris, and then used BoundingboxAA (x direction) to measure cluster diameter. Cluster diameter is averaged over the time period when cluster is at nurse cell junction.
Acknowledgments
We thank Amy Peterson for assistance with movie tracking and Thomas Kenney for assistance with confocal microscopy. This work was supported by NIH Grant GM R01 46425 (to D.J.M.). N.S.G is the recipient of the Lee and William Abramowitz Professorial Chair of Biophysics, and acknowledges Israeli Science Foundation Grant 580/12 for support. We acknowledge the use of the Neuroscience Research Institute-Molecular, Cellular and Developmental Biology Microscopy Facility and the instrumentation supported by the Office of The Director, National Institutes of Health of the NIH under Award S10OD010610.
Footnotes
- ↵1To whom correspondence may be addressed. Email: denise.montell{at}lifesci.ucsb.edu or nir.gov{at}weizmann.ac.il.
Author contributions: D.C., N.S.G., and D.J.M. designed research; D.C., W.D., M.P., and N.S.G. performed research; J.L. and N.S.G. contributed new reagents/analytic tools; D.C., W.D., N.S.G., and D.J.M. analyzed data; and D.C., N.S.G., and D.J.M. wrote the paper.
The authors declare no conflict of interest.
This article is a PNAS Direct Submission.
This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10.1073/pnas.1522656113/-/DCSupplemental.
References
- ↵.
- Parent CA,
- Devreotes PN
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵.
- Heit B,
- Tavener S,
- Raharjo E,
- Kubes P
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵.
- Steinberg M
- ↵
- ↵.
- Grimm O,
- Coppey M,
- Wieschaus E
- ↵
- ↵
- ↵
- ↵
- ↵.
- Duchek P,
- Rørth P
- ↵
- ↵
- ↵.
- Bohlin T
- ↵.
- Vitorino P,
- Meyer T
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵.
- Liu Y,
- Montell DJ
- ↵.
- Bai J,
- Montell D
- ↵
- ↵
- ↵
- ↵.
- Kuriyama S, et al.
- ↵.
- Herzmark P, et al.
- ↵
- ↵.
- Cheung KJ, et al.
- ↵.
- Prasad M,
- Wang X,
- He L,
- Montell DJ
- ↵
Citation Manager Formats
Article Classifications
- Biological Sciences
- Cell Biology