Spatio-temporal morphology changes in and quenching effects on the 2D spreading dynamics of cell colonies in both plain and methylcellulose-containing culture media

To deal with complex systems, microscopic and global approaches become of particular interest. Our previous results from the dynamics of large cell colonies indicated that their 2D front roughness dynamics is compatible with the standard Kardar–Parisi–Zhang (KPZ) or the quenched KPZ equations either in plain or methylcellulose (MC)-containing gel culture media, respectively. In both cases, the influence of a non-uniform distribution of the colony constituents was significant. These results encouraged us to investigate the overall dynamics of those systems considering the morphology and size, the duplication rate, and the motility of single cells. For this purpose, colonies with different cell populations (N) exhibiting quasi-circular and quasi-linear growth fronts in plain and MC-containing culture media are investigated. For small N, the average radial front velocity and its change with time depend on MC concentration. MC in the medium interferes with cell mitosis, contributes to the local enlargement of cells, and increases the distribution of spatio-temporal cell density heterogeneities. Colony spreading in MC-containing media proceeds under two main quenching effects, I and II; the former mainly depending on the culture medium composition and structure and the latter caused by the distribution of enlarged local cell domains. For large N, colony spreading occurs at constant velocity. The characteristics of cell motility, assessed by measuring their trajectories and the corresponding velocity field, reflect the effect of enlarged, slow-moving cells and the structure of the medium. Local average cell size distribution and individual cell motility data from plain and MC-containing media are qualitatively consistent with the predictions of both the extended cellular Potts models and the observed transition of the front roughness dynamics from a standard KPZ to a quenched KPZ. In this case, quenching effects I and II cooperate and give rise to the quenched-KPZ equation. Seemingly, these results show a possible way of linking the cellular Potts models and the 2D colony front roughness dynamics.     

Disentangling cadherin-mediated cell-cell interactions in collective cancer cell migration

Cell dispersion from a confined area is fundamental in a number of biological processes, including cancer metastasis. To date, a quantitative understanding of the interplay of single-cell motility, cell proliferation, and intercellular contacts remains elusive. In particular, the role of E- and N-cadherin junctions, central components of intercellular contacts, is still controversial. Combining theoretical modeling with in vitro observations, we investigate the collective spreading behavior of colonies of human cancer cells (T24). The spreading of these colonies is driven by stochastic single-cell migration with frequent transient cell-cell contacts. We find that inhibition of E- and N-cadherin junctions decreases colony spreading and average spreading velocities, without affecting the strength of correlations in spreading velocities of neighboring cells. Based on a biophysical simulation model for cell migration, we show that the behavioral changes upon disruption of these junctions can be explained by reduced repulsive excluded volume interactions between cells. This suggests that in cancer cell migration, cadherin-based intercellular contacts sharpen cell boundaries leading to repulsive rather than cohesive interactions between cells, thereby promoting efficient cell spreading during collective migration.     

Long-range alteration of the physical environment mediates cooperation between Pseudomonas aeruginosa swarming colonies

Pseudomonas aeruginosa makes and secretes massive amounts of rhamnolipid surfactants that enable swarming motility over biogel surfaces. But how these rhamnolipids interact with biogels to assist swarming remains unclear. Here, I use a combination of optical techniques across scales and genetically engineered strains to demonstrate that rhamnolipids can induce agar gel swelling over distances >10,000× the body size of an individual cell. The swelling front is on the micrometric scale and is easily visible using shadowgraphy. Rhamnolipid transport is not restricted to the surface of the gel but occurs through the whole thickness of the plate and, consequently, the spreading dynamics depend on the local thickness. Surprisingly, rhamnolipids can cross the whole gel and induce swelling on the opposite side of a two-face Petri dish. The swelling front delimits an area where the mechanical properties of the surface properties are modified: water wets the surface more easily, which increases the motility of individual bacteria and enables collective motility. A genetically engineered mutant unable to secrete rhamnolipids (ΔrhlA), and therefore unable to swarm, is rescued from afar with rhamnolipids produced by a remote colony. These results exemplify the remarkable capacity of bacteria to change the physical environment around them and its ecological consequences.     

Direct upstream motility in Escherichia coli

We provide an experimental demonstration of positive rheotaxis (rapid and continuous upstream motility) in wild-type Escherichia coli freely swimming over a surface. This hydrodynamic phenomenon is dominant below a critical shear rate and robust against Brownian motion and cell tumbling. We deduce that individual bacteria entering a flow system can rapidly migrate upstream (>20 μm/s) much faster than a gradually advancing biofilm. Given a bacterial population with a distribution of sizes and swim speeds, local shear rate near the surface determines the dominant hydrodynamic mode for motility, i.e., circular or random trajectories for low shear rates, positive rheotaxis for moderate flow, and sideways swimming at higher shear rates. Faster swimmers can move upstream more rapidly and at higher shear rates, as expected. Interestingly, we also find on average that both swim speed and upstream motility are independent of cell aspect ratio.     

The amoebal cell of Physarum polycephalum: colony formation and growth.

Two stages of colony growth were observed during microscopic studies of Physarum polycephalum amoebae. During the first stage, “spreading growth,” the colony is composed of dispersed single cells. During the second stage, “aggregate growth,” most of the active cells in a colony are aggregated in a ring at the colony boundary. Measurements of cell movement as a function of bacterial concentration indicate that, during both spreading and aggregate growth, cell movements are not affected by changes in bacterial concentration but that the transition from spreading to aggregate growth occurs earlier on plates with lower bacterial concentrations. These results indicate that autonomous characteristics of the amoebae are more important for the determination of colony form than local variations in the concentrations of nutrients.The genetic determination of colony form is demonstrated by the existence of mutants that display specific alterations in colony morphology. Because the aggregate rings of these mutants move at an increased rate, mutant clones appear as variant sectors of wild-type colonies. The increased rate of mutant ring movement suggests that this selection method may be a useful technique for isolating mutant myxamoebae with defects in movement and behavior.     

Computer-aided cell colony counting

Counting cell colonies is a tedious task when performed with the light microscope. Moreover, unless strict double-blind protocols are adhered to, biased counts are difficult to avoid. Presented here is a computer software application that performs accurate, reproducible cell colony counts with a minimum of user generated bias. The application is based upon the Apple IICX computer system with Image software and AppleScan. Colonies are grown on 24-well plates and prepared in such a way as to permit good quality scanning. The scans are then transferred to Image and the individual colonies in each well are counted. Good correlation with counts done by light microscopy is achieved.     

Automated time-lapse microscopy and high-resolution tracking of cell migration

We describe a novel fully automated high-throughput time-lapse microscopy system and evaluate its performance for precisely tracking the motility of several glioma and osteoblastic cell lines. Use of this system revealed cell motility behavior not discernable with conventional techniques by collecting data (1) from closely spaced time points (minutes), (2) over long periods (hours to days), (3) from multiple areas of interest, (4) in parallel under several different experimental conditions. Quantitation of true individual and average cell velocity and path length was obtained with high spatial and temporal resolution in “scratch” or “wound healing” assays. This revealed unique motility dynamics of drug-treated and adhesion molecule-transfected cells and, thus, this is a considerable improvement over current methods of measurement and analysis. Several fluorescent vital labeling methods commonly used for end-point analyses (GFP expression, DiO lipophilic dye, and Qtracker nanocrystals) were found to be useful for time-lapse studies under specific conditions that are described. To illustrate one application, fluorescently labeled tumor cells were seeded onto cell monolayers expressing ectopic adhesion molecules, and this resulted in consistently reduced tumor cell migration velocities. These highly quantitative time-lapse analysis methods will promote the creation of new cell motility assays and increase the resolution and accuracy of existing assays.Joseph S. Fotos and Vivek P. Patel contributed equally to this work     

An integrated approach to the characterization of cell movement.

BACKGROUND: Most phenomena in developmental biology involve or depend upon cell migration. This article describes a comprehensive framework for the characterization and analysis of trajectories defined by cell movement. The following two perspectives are considered: (a) the behavior of each individual cell and (b) interactions between neighboring pairs of cells. METHODS: The measurements considered for individual trajectories include the velocity magnitude and orientation, maximum spatial dispersion, displacement effectiveness, and displacement entropies. Interactions between two trajectories are characterized by comparing the respective velocities. RESULTS: The potential of the overall framework is illustrated using data of moving cells in different biological environments. The work shows that it is possible to use the new algorithm presented here to characterize cell motility. CONCLUSIONS: The features of the algorithm were successful in determining the motility changes under different experimental conditions.     

Effects of long-term serial cell passaging on cell spreading,migration, and cell-surface ultrastructures of cultured vascular endothelial cells

The effects of serial cell passaging on cell spreading, migration, and cell-surface ultrastructures have been less investigated directly. This study evaluated the effects of long-term serial cell passaging (totally 35 passages) on cultured human umbilical vein endothelial cells which were pre-stored at −80 °C as usual. Percentage- and spread area-based spreading assays, measurements of fluorescently labeled actin filaments, migration assay, and measurements of cell-surface roughness were performed and quantitatively analyzed by confocal microscopy or atomic force microscopy. We found that the abilities of cell spreading and migration first increased at early passages and then decreased after passage 15, in agreement with the changes in average length of actin filaments. Recovery from cold storage and effects of cell passaging were potentially responsible for the increases and decreases of the values, respectively. In contrast, the average roughness of cell surfaces (particularly the nucleus-surrounding region) first dropped at early passages and then rose after passage 15, which might be caused by cold storage- and cell passaging-induced endothelial microparticles. Our data will provide important information for understanding serial cell passaging and implies that for pre-stored adherent cells at −80 °C cell passages 5–10 are optimal for in vitro studies.     

A hydrodynamics approach to the evolution of multicellularity: flagellar motility and germ-soma differentiation in volvocalean green algae

During the unicellular-multicellular transition, there are opportunities and costs associated with larger size. We argue that germ-soma separation evolved to counteract the increasing costs and requirements of larger multicellular colonies. Volvocalean green algae are uniquely suited for studying this transition because they range from unicells to multicellular individuals with germ-soma separation. Because Volvocales need flagellar beating for movement and to avoid sinking, their motility is modeled and analyzed experimentally using standard hydrodynamics. We provide comparative hydrodynamic data of an algal lineage composed of organisms of different sizes and degrees of complexity. In agreement with and extending the insights of Koufopanou, we show that the increase in cell specialization as colony size increases can be explained in terms of increased motility requirements. First, as colony size increases, soma must evolve, the somatic-to-reproductive cell ratio increasing to keep colonies buoyant and motile. Second, increased germ-soma specialization in larger colonies increases motility capabilities because internalization of nonflagellated germ cells decreases colony drag. Third, our analysis yields a limiting maximum size of the volvocalean spheroid that agrees with the sizes of the largest species known. Finally, the different colony designs in Volvocales reflect the trade-offs between reproduction, colony size, and motility.     

'Dicty dynamics': Dictyostelium motility as persistent random motion

We model the motility of Dictyostelium cells in a systematic data-driven manner. We deduce a minimal dynamical model that reproduces the statistical features of experimental trajectories. These are trajectories of the centroid of the cell perimeter, which is more sensitive to pseudopod activity than the usual tracking by centroid or nucleus. Our data account for cell individuality and dictate a model that extends the cell-type specific models recently derived for mammalian cells. Two generalized Langevin equations model stochastic periodic pseudopod motion parallel and orthogonal to the amoeba's direction of motion. This motion propels the amoeba with a random periodic left-right waddle in a direction that has a long persistence time. The model fully accounts for the statistics of the experimental trajectories, including velocity power spectra and auto-correlations, non-Gaussian velocity distributions, and multiplicative noise. Thus, we find neither need nor place in our data for an interpretation in terms of anomalous diffusion. The model faithfully captures cell individuality as different parameter values in the model, and serves as a basis for integrating the local mechanics of cell motion with our observed long-term behavior.     

Quantitative image analysis of connective tissue progenitors

OBJECTIVE: To develop an image analysis system to automatically identify colony-forming units (CFUs) in in vitro cell cultures of connective tissue progenitors. This system was designed to quantitatively assess colony morphology and number of colonies in 4-cm(2) culture wells. STUDY DESIGN: Large field-of-view high-resolution fluorescence images of 4',6-diamidino-2-phenylindole (DAPI)- and alkaline phosphatase (AP)-stained bone marrow cell cultures were obtained using an epi-fluorescence microscope and automated scanning stage. Cell nuclei were identified in the DAPI-stained images after removal of fluorescent debris from the image. An Euclidean distance map (EDM) of the segmented cell nuclei was used to cluster cell nuclei into colonies. The automated system was evaluated using 40 tissue culture wells of bone marrow aspirate samples. The results of the automated analysis were compared to the manual tracings of colonies by 3 reviewers. RESULTS: The automated method agreed with all 3 reviewers on average 87.5% of the time. Additionally, reviewers identified other colonies not outlined by the reviewers on average 2.7 times more than the automated method. CONCLUSION: The automated method is a less biased method for identifying CFUs than individual reviewers, it provides more quantitative information about colony morphology than can be obtained manually and it is less time consuming.     

Nonlinearity in the growth of bacterial colonies: conditions and causes

The universally recognized kinetic model of colony growth, introduced by Pirt, predicts a linear increase of colony size. The linearity follows from the assumption that the colony expands through the growth of only such cells that are located immediately behind the moving colony front, in the so-called peripheral zone of constant width and density. In this work, Pirt's model was tested on two bacteria--Alcaligenes sp. and Pseudomonas fluorescens--having markedly distinct cultural properties and grown on agarized medium with pyruvate. The colony size dynamics was followed for different densities of the inoculum, ranging from a single cell to a microdroplet of bacterial suspension (10(5)-10(6) cells), and for different depths of the agar layer, determining the amount of available substrate. A linear growth mode was observed only with P. fluorescens and only in the case of growth from a microdroplet. When originating from a single cell, colonies of both organisms displayed nonlinear growth with a distinct peak of Kr (the rate of colony radius increase) occurring after 2-3 days of growth. The growth of P. fluorescens colonies showed virtually no dependence on the depth of the agarized medium, whereas the rate of colony size increase of Alcaligenes sp. turned out to be directly related to the medium layer thickness. The departure from linearity is consistently explained by a new kinetic chart stipulating a possible contribution to the colony growth not only of peripheral cells but also (much more distinct in Alcaligenes) of cells at the colony center. The colony growth dynamics is determined not only by the concentration of the limiting substrate but also by the amount of autoinhibitor, the synthesis of which is governed by age of cells. The distinctions of growth from a single cell and microdroplet could also originate as a result of dissociation into the R- and S-forms and competition between the corresponding subpopulations for oxygen and the common substrate.     

Two continuum models for the spreading of myxobacteria swarms

We analyze the phenomenon of spreading of a Myxococcus xanthus bacterial colony on plates coated with nutrient. The bacteria spread by gliding on the surface. In the first few hours, cell growth is irrelevant to colony spread. In this case, bacteria spread through peninsular protrusions from the edge of the initial colony. We analyze the diffusion through the narrowing reticulum of cells on the surface mathematically and derive formulae for the spreading rates. On the time scale of tens of hours, effective diffusion of the bacteria, combined with cell division and growth, causes a constant linear increase in the colony's radius. Mathematical analysis and numerical solution of reaction-diffusion equations describing the bacterial and nutrient dynamics demonstrate that, in this regime, the spreading rate is proportional to the square root of both the effective diffusion coefficient and the nutrient concentration. The model predictions agree with the data on spreading rate dependence on the type of gliding motility.     

Reproductive allocation within a polygyne, polydomous colony of the ant Myrmica rubra

1. Ant colonies commonly have multiple egg‐laying queens (secondary polygyny). Polygyny is frequently associated with polydomy (single colonies occupy multiple nest sites) and restricted dispersal of females. The production dynamics and reproductive allocation patterns within a population comprising one polygyne, polydomous colony of the red ant Myrmica rubra were studied. 2. Queen number per nest increased with nest density and the number of adult workers increased with the number of resident queens and with nest density. This suggests that nest site limitation promotes polygyny and that workers accumulate in nest units incapable of budding. 3. Nest productivity increased with the number of adult workers and production per queen was independent of queen number. Productivity increased with nest density, suggesting local resource enhancement. This shows that productivity can be a linear function of queen numbers and that the limiting factor is not the egg‐laying capacity of queens. 4. The total and per capita production of reproductives decreased towards the periphery of the colony, suggesting that the spatial location of nest units affects sexual production. Thus nests at the periphery of the colony invested more heavily in new workers. This is consistent with earlier observations in plants and could either represent investment in future budding or increased defence. 5. The colony produced only five new queens and 2071 males, hence the sex ratio was extremely male biased.     

Flow and Diffusion in Channel-Guided Cell Migration

Collective migration of mechanically coupled cell layers is a notable feature of wound healing, embryonic development, and cancer progression. In confluent epithelial sheets, the dynamics have been found to be highly heterogeneous, exhibiting spontaneous formation of swirls, long-range correlations, and glass-like dynamic arrest as a function of cell density. In contrast, the flow-like properties of one-sided cell-sheet expansion in confining geometries are not well understood. Here, we studied the short- and long-term flow of Madin-Darby canine kidney (MDCK) cells as they moved through microchannels. Using single-cell tracking and particle image velocimetry (PIV), we found that a defined averaged stationary cell current emerged that exhibited a velocity gradient in the direction of migration and a plug-flow-like profile across the advancing sheet. The observed flow velocity can be decomposed into a constant term of directed cell migration and a diffusion-like contribution that increases with density gradient. The diffusive component is consistent with the cell-density profile and front propagation speed predicted by the Fisher-Kolmogorov equation. To connect diffusion-mediated transport to underlying cellular motility, we studied single-cell trajectories and occurrence of vorticity. We discovered that the directed large-scale cell flow altered fluctuations in cellular motion at short length scales: vorticity maps showed a reduced frequency of swirl formation in channel flow compared with resting sheets of equal cell density. Furthermore, under flow, single-cell trajectories showed persistent long-range, random-walk behavior superimposed on drift, whereas cells in resting tissue did not show significant displacements with respect to neighboring cells. Our work thus suggests that active cell migration manifests itself in an underlying, spatially uniform drift as well as in randomized bursts of short-range correlated motion that lead to a diffusion-mediated transport.     

Flow and Diffusion in Channel-Guided Cell Migration

Collective migration of mechanically coupled cell layers is a notable feature of wound healing, embryonic development, and cancer progression. In confluent epithelial sheets, the dynamics have been found to be highly heterogeneous, exhibiting spontaneous formation of swirls, long-range correlations, and glass-like dynamic arrest as a function of cell density. In contrast, the flow-like properties of one-sided cell-sheet expansion in confining geometries are not well understood. Here, we studied the short- and long-term flow of Madin-Darby canine kidney (MDCK) cells as they moved through microchannels. Using single-cell tracking and particle image velocimetry (PIV), we found that a defined averaged stationary cell current emerged that exhibited a velocity gradient in the direction of migration and a plug-flow-like profile across the advancing sheet. The observed flow velocity can be decomposed into a constant term of directed cell migration and a diffusion-like contribution that increases with density gradient. The diffusive component is consistent with the cell-density profile and front propagation speed predicted by the Fisher-Kolmogorov equation. To connect diffusion-mediated transport to underlying cellular motility, we studied single-cell trajectories and occurrence of vorticity. We discovered that the directed large-scale cell flow altered fluctuations in cellular motion at short length scales: vorticity maps showed a reduced frequency of swirl formation in channel flow compared with resting sheets of equal cell density. Furthermore, under flow, single-cell trajectories showed persistent long-range, random-walk behavior superimposed on drift, whereas cells in resting tissue did not show significant displacements with respect to neighboring cells. Our work thus suggests that active cell migration manifests itself in an underlying, spatially uniform drift as well as in randomized bursts of short-range correlated motion that lead to a diffusion-mediated transport.