A Mathematical Model of Collective Cell Migration in a Three-Dimensional, Heterogeneous Environment

<div><p>Cell migration is essential in animal development, homeostasis, and disease progression, but many questions remain unanswered about how this process is controlled. While many kinds of individual cell movements have been characterized, less effort has been directed towards understanding how clusters of cells migrate collectively through heterogeneous, cellular environments. To explore this, we have focused on the migration of the border cells during Drosophila egg development. In this case, a cluster of different cell types coalesce and traverse as a group between large cells, called nurse cells, in the center of the egg chamber. We have developed a new model for this collective cell migration based on the forces of adhesion, repulsion, migration and stochastic fluctuation to generate the movement of discrete cells. We implement the model using Identical Math Cells, or IMCs. IMCs can each represent one biological cell of the system, or can be aggregated using increased adhesion forces to model the dynamics of larger biological cells. The domain of interest is filled with IMCs, each assigned specific biophysical properties to mimic a diversity of cell types. Using this system, we have successfully simulated the migration of the border cell cluster through an environment filled with larger cells, which represent nurse cells. Interestingly, our simulations suggest that the forces utilized in this model are sufficient to produce behaviors of the cluster that are observed <i>in vivo</i>, such as rotation. Our framework was developed to capture a heterogeneous cell population, and our implementation strategy allows for diverse, but precise, initial position specification over a three- dimensional domain. Therefore, we believe that this model will be useful for not only examining aspects of <i>Drosophila</i> oogenesis, but also for modeling other two or three-dimensional systems that have multiple cell types and where investigating the forces between cells is of interest.</p></div

Simulating the three dimensional model results in collective migration.

<p>A simulation showing six border cells (green), two polar cells (red), the epithelium (transparent green), and the surface of the oocyte (black, right) at three time points during the migration. Fifteen nurse cells are situated inside the egg chamber, but are not plotted so as to maintain clarity of this three dimensional structure. Polar cells are surrounded by border cells, making them hard to distinguish. (A) At 2 minutes, cells are beginning to invade between nurse cells. (B) At 2.4 hours, the cluster is about halfway to its destination. (C) At 5.6 hours, the border cell cluster has reached the edge of the oocyte. See also Supplemental <a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0122799#pone.0122799.s003" target="_blank">S3 Movie</a>.</p

Simulations with four, six, and eight border cells at the same time point (<i>t</i> = 1.8 hours) during migration.

<p>The cluster with four border cells (A) has moved significantly less distance than the cluster with six (B) or eight (C) border cells.</p

The total time and relative time taken for computational clusters of 4, 6 and 8 border cells to complete migration.

<p>The total time and relative time taken for computational clusters of 4, 6 and 8 border cells to complete migration.</p

Composite non-dimensional parameter values used for simulation.

<p>The parameter <i>α</i> is the time scale (an hour) and <i>D</i> is the diameter of an IMC (7<i>μ</i>m). We can consider the <i>α</i> and <i>D</i> scaling the damping viscosity coefficient, <i>μ</i>, then the parameters in the table represent a ratio of force exerted on the IMC by movement through the heterogeneous medium and the force exerted by the adhesion, repulsion, migration and stochastic forces.</p

Polar cell positions along main axis of migration.

<p>(A) The distance of the polar cells from the anterior of the egg chamber versus time. (B) The relative positions of the two polar cells to one another, along the axis that runs from anterior to posterior through the egg chamber. Each line corresponds to one of the polar cells. As the cluster moves forward, we observe that the polar cells are changing position with respect to one another along this axis, including a complete switch at 0.8 hours. This simulation modeled six border cells.</p

Border cell migration in <i>Drosophila melanogaster</i> egg development.

<p>(a) At the beginning of stage eight, the polar cells (yellow) and border cells (green) lie in the follicular epithelium of the developing egg chamber. In stage nine, these cells coalesce to form a cluster that detaches from the epithelium. The cluster then translocates between large nurse cells through the egg chamber to reach the developing oocyte (gray) by stage ten. The border cell cluster migrates about 150<i>μ</i>m over approximately 4–6 hours. (b) Still images from a time-lapse movie of wild-type border cell migration. The motile cells are marked in green by expression of Slbo-life-Act-GFP. The oocyte, which autofluoresces, is indicated by the dashed line. The nuclei of all cells, including the large, polyploid nurse cells, are seen in blue. In the image at time 0, border cells have already clustered (arrow) and begun moving towards the oocyte. In this example, the border cells reached the oocyte border by 3.5h (arrow on right-most panel). (c) Still images at a higher magnification from a time-lapse movie of a different egg chamber. Images differ by 30 minute intervals. The border cells are marked by a membrane-tethered GFP, and show wild-type behavior. Cells can be observed to change relative positions with respect to the front of the cluster as they move toward the right (arrow and arrowhead indicate the same cell over time). See also Supplemental <a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0122799#pone.0122799.s001" target="_blank">S1</a> and <a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0122799#pone.0122799.s002" target="_blank">S2</a> Movies.</p