Diffusional inhomogeneity in cell cultures
Zhang Jia-Zheng, Li Na, Chen Wei
State Key Laboratory of Surface Physics and Department of Physics, Fudan University, Shanghai 200433, China

 

† Corresponding author. E-mail: phchenwei@fudan.edu.cn

Abstract

Cell migrations in the cell cultures are found to follow non-Gaussian statistics. We recorded long-term cell migration patterns with more than six hundred cells located in 28 mm2. Our experimental data support the claim that an individual cell migration follows Gaussian statistics. Because the cell culture is inhomogeneous, the statistics of the cell culture exhibit a non-Gaussian distribution. We find that the normalized histogram of the diffusion velocity follows an exponential tail. A simple model is proposed based on the diffusional inhomogeneity to explain the exponential distribution of locomotion activity in this work. Using numerical calculation, we prove that our model is in great agreement with the experimental data.

1. Introduction

Bacteria, eukaryotic cells and animals like birds and insects[14] are believed to adopt perplexing random walk strategies to explore their environment. Many phenomenological models have been made to classify these random walk strategies by the characteristics of stochastic fluctuation. We notice that these systems have been treated in the same framework of self-propellers[5] with collective motion. The ensembles are always homogenous, if not, they will fall apart in the collective motion. Hence, the ensemble behavior is dependent on the diffusion features of a unit and the interaction relationships between the units.

We address that eukaryotic cell motion may be more complicated than the homogenous collective random motion. As a model system in this field, Dictyostelium discoideum (DD) cell is used in this paper. When individual DD cells are starving, they tend to aggregate to form a multicellular migrating slug.[6] Before the aggregation, cells move randomly to search the space for food in the local environment. Interestingly, many experiments showed that the aggregation is triggered by a specific chemical signal and then the movements of all the cells are synchronized and accelerated.[7] The aggregation process is much like the homogenous collective motion, while before synchronization and acceleration, the assumption of homogeneity may not be suitable for the searching period. In recent research, the inhomogeneity of the diffusion for acetylcholine receptors is reported,[8] suggesting the assumption of homogeneity is not straightforward. Here, we focus on the locomotion activity in the searching period. Anomalous dynamics have been put forward in this region.[912] Meanwhile, the search strategy adopted by cells attracted much attention in recent decades.[1315]

To understand the statistical strategy adapted by DD cells, we analyzed the cell trajectories in the searching process mentioned above. We wonder whether the assumption of homogeneity is still suitable for this system. Our experimental results suggest the cell culture is inhomogeneous. Cells with different diffusion constant exist simultaneously in the cell culture, despite the fact that all cells are divided from the same spore. Starting from the inhomogeneity, we proposed a simple model to describe the searching strategy. In this model, cells have their own diffusion velocity and exhibit different distributions of velocity fluctuations. Meanwhile, the distribution of cell velocities follows an exponential tail. This model gives us a new perspective and method to find the optimal strategy in order to explore the space. Furthermore, we present a straightforward test of our model by numerical calculation.

2. Experimental setup

KAx-3 cells of Dictyostelium Discoideum (DD) were used in our experiments and were obtained by the following procedures, similar to the procedures given in Ref. [16]. Cells were plated on the top of a solid medium with nutrition (20-g agar, 10-g glucose, 10-g trytone, 1-g yeast extract, 1.9-g KH2PO4, 0.78-g K2HPO4 3H2O, and 1-g MgSO4 per 1-L distilled water), grown in an incubator kept at 23 °C for about 39 hours, and then isolated by centrifuging at 100 g and re-suspended in phosphate buffer solution (PBS). The isolated cell of KAx-3 was respectively starved for 4 hours in a refrigerator kept at 4 °C before they were plated with the desired density on top of a 1.5% agar plate without any nutrition. The cell area density in the sample studied was kept 23 cells per mm2. All the procedures mentioned above were carried out under a clean hood and the environmental temperature was controlled at 25 °C within ±1 °C. An objective with 4×magnification was used because it can record as many cells as possible to achieve good statistics. Images were digitized and stored in a computer equipped with a digital camera (Canon 700D, spatial resolution 5138×3456). Four samples were recorded at the same time by four cameras for more than 50 hours, the time interval between images was 10 seconds. We used homemade software to recognize and track cells, and finally obtained the trajectories and other information of the cells.

3. Results
3.1. Definition of the velocity and the mean velocity of a trajectory

In the experiment, we obtained the DD cell trajectories from consecutive images of the cell culture with spatial resolution of and captured 633 cells located in 28 mm2 in each frame. More than frames were used. The trajectories are recorded between collisions, because cells are able to run into another cell and their trajectories break. We only analyze the trajectories of individual cells in this work. Using a homemade analysis program, we determined their position at time t and the two-dimensional displacement vector over the delay time τ from 7000 trajectories.

Figure 1 shows the measured mean square displacement (MSD) as a function of delay time τ for the individual cell trajectories. We used the Ornstein–Uhlenbeck process to fit our data as an approximation with the resistant time s and D is the diffusion coefficient. The fitted resistant time s indicates cells lose their direction after 580 s.[17] Hence, we can define the velocity from the displacements , when τ is much shorter than τr. Without further loss of generality, we define with τ = 50 s. For every trajectory, we define , averaged on all the velocities in each trajectory. represents the mean velocity of a trajectory.

Fig. 1. (color online) Measured MSD as a function of delay time τ for the individual cell trajectories. Three data points are skipped between every two points in the figure for clarity. The red line is a non-linear fit to the data points using OU process equation, , with resistant time .

Figure 2 shows the normalized histograms of v and . The exponential tail in the histogram of v was hard to explain under the framework of OU process. Meanwhile, the histogram of shown in Fig. 2(b) indicates the mean velocities of trajectories are widely dispersed, and the distribution also exhibits an exponential tail.

Fig. 2. (color online) (a) Normalized histogram of the cell velocity v measured by every displacement in 50 s. (b) Normalized histogram of the mean velocity measured by every trajectory. The blue lines in panels (a) and (b) are exponential fit to the curve tail, and , with , .
3.2. Stochastic fluctuation of each diffusion group and their coupling with the group velocity

We classified trajectories into five bins with equal width from the slowest velocity to the fastest velocity. Details are in the caption of Fig. 3. In each bin, trajectories belong to the same velocity group. Stochastic fluctuations of each velocity group are shown in Fig. 3(a). In each group, the fluctuation seems to be Gaussian. We fitted the distributions with Gaussian functions, and got the fitted parameters of the mean value and standard deviation . We shifted the x coordinate, and plotted versus in Fig. 3(b), where . The data of five groups overlap each other after shift and the linearity of both tails indicate the Gaussian fits of the fluctuations in five groups are accurate.

Fig. 3. (color online) (a) Normalized histogram of velocity measured by all displacements from the trajectories in each bin. The red lines are Gaussian fit of the data in each group. Trajectories are classified into five bins, (black squares), (red circles), and (green up triangles), (blue down triangles), and (light blue rhombuses). (b) The plot for log( ) versus , with , where and σ are the fitting mean value and the standard deviation in the Gaussian fits in panel (a), respectively. (c) The standard deviation as a function of . The blue line is a linear fit of the data, , with and κ = 0.088.

Figure 3(c) shows the standard deviation as a function of . We fit the data linearly and the fitting empirical equation reads where the fitting parameters and κ = 0.088.

4. Discussion

With the help of empirical equation (1), we proposed a straightforward model to describe this inhomogeneity. First, we tried to rebuild the distribution in Fig. 2(a) with the distributions in Fig. 2(b) where is the probability distribution function (PDF) of v before classification, is the distributions of velocity fluctuations of each group, and is the weight of each group, which directly related to the number of cells in the group. Second, the bin-width of the cell groups were decreased continuously and the discrete equation (2) would change into an integral form as the following equation (4). Here was replaced by the mean velocity of a trajectory and was now replaced by the histogram number of trajectories with a mean velocity of . Hence, equations (1) and (2) become as where is the histogram as shown in Fig. 2(b).

Physically, equation (4) supplies a clear insight into the dynamics of the cell culture. Generally, the instantaneous velocity will be forgotten by a cell after the typical persistent time. Because the persistent time was much shorter than the length of the trajectory, the velocities could be regarded as independent. These independent velocities falling into a small bin should obey Gaussian distribution under the Central-Limit Theorem. The width of this Gaussian distribution depends on the mean velocity, as shown in Fig. 3(c). Meanwhile, the velocities of the cell are not thesame. plays a determinant role in the system, which responds to the observed non-Gaussian distributions indeed.

Furthermore, to thoroughly understand the significance of the exponential tail of , we made a numerical simulation as follows. For simplicity, we construct the formula of to represent the form of as follows: where and is chosen to overlap the tail of , is chosen to keep the areas under two curves equal, and is the fitting parameter from Fig. 2(b). The experimental data and the constructed are shown in Fig. 4(a).

Fig. 4. (color online) (color online) Numerical calculation. (a) Normalized histogram of cell velocity (black squares) and numerical simulation data (red circles) from . (b) Simulation result of the velocity distribution (red circles), compared with experimental data (dark squares).

Finally, one obtains where C2 and C3 are constants. Equation (6) suggests that the tail of the distribution of velocity is only dependent on λ, which is the characteristic length of the exponential tail of .

Equation (6) is not easy to calculate analytically, instead we find out the answer by numerical calculation. We calculate of a cell culture with cells following Eqs. (3) and (5). The result of is shown in Fig. 4(b). It is well matched with experimental data.

Now we are able to review our understanding of the dynamics of cell motion. Previous studies believed that the cell motion was governed by random motion with an unclear form. We suggest understanding this system in the following steps. First, for an individual cell, the velocity fluctuates around a constant value and the fluctuation noise is multiplicative. The equation follows: where is a noise term with unit intensity. The values of and κ determine the difference between values of the β1 and β2 in Fig. 2. Second, the values of vary from cell to cell and they obey a probability distribution with an exponential tail.

5. Summary

To summarize, with the help of a digital camera having spatial resolution of 5138×3456, we analyzed migrating trajectories of cells, and found the exponential tail of the probability distribution function of the cell velocity, which stems from the inhomogeneous weight number of cells for different diffusion motion. We believe this simple model is significant in understanding the random motion of eukaryotic cells, and also in finding the optimal strategy to explore the space.

Acknowledgment

We thank Prof. Lian-Sheng Hou for the supply of KAx-3 cell samples.

Reference
[1] Berg H C Purcell E M 1977 Biophys. J. 20 193
[2] Czirok A Schlett K Madarasz E Vicsek T 1998 Phys. Rev. Lett. 81 3038
[3] Couzin I D Krause J 2003 Adv. Study Behav. 32 1
[4] Ni R Puckett J G DUfresne E R Ouellette N T 2015 Phys. Rev. Lett. 115 118104
[5] Nourhani A Crespi V H Lammert P E 2014 Phys. Rev. 89 062304
[6] Maree A F M Hogeweg P 2014 Proc. Natl. Acad Sci. USA 98 3879
[7] Varnum B Soll D R 2014 J. Cell Biol. 99 1151
[8] He W Song H Su Y Geng L Ackerson B J Peng H B Tong P 2016 Nat. Commun. 7 11701
[9] Li L Norrelykke S F Cox E C 2008 Plos One 3 e2093
[10] Dieterich P Klages R Preuss R Schwab A 2008 Proc. Natl. Acad. Sci. USA 105 459
[11] Takagi H Sato M J Yanagida T Ueda M 2008 Plos One 3 e2648
[12] Shenderove A D Sheetz M P 1997 Biolphys. J. 72 2382
[13] Viswanathan G M Afanasyev V Buldyrev S V Havlin S da Luz M G E Raposo E P Stanley H E 1999 Nature 401 911
[14] Viswanathan G M Bartumeus F Buldyrev S V Catalan J Fulco U L Havlin S da Luz M G E Lyra M L Raposo E P Stanley H E 2002 Physica 314 208
[15] Reynolds A 2002 J. Plrev. 14 59
[16] Chen S Li N Hsu S F Zhang J Z Lai P Y Chan C K Chen W 2014 Soft Matter 10 3421
[17] Uhlenbeck G E Ornstein L S 1930 Phys. Rev. 36 823