Morphology and dynamic scaling analysis of cell colonies with linear growth fronts

Dimensional crossover in Kardar-Parisi-Zhang growth

Two-dimensional (2D) Kardar-Parisi-Zhang (KPZ) growth is usually investigated on substrates of lateral sizes L_{x}=L_{y}, so that L_{x} and the correlation length (ξ) are the only relevant lengths determining the scaling behavior. However, in cylindrical geometry, as well as in flat rectangular substrates L_{x}≠L_{y} and, thus, the surfaces can become correlated in a single direction, when ξ∼L_{x}≪L_{y}. From extensive simulations of several KPZ models, we demonstrate that this yields a dimensional crossover in their dynamics, with the roughness scaling as W∼t^{β_{2D}} for t≪t_{c} and W∼t^{β_{1D}} for t≫t_{c}, where t_{c}∼L_{x}^{1/z_{2D}}. The height distributions (HDs) also cross over from the 2D flat (cylindrical) HD to the asymptotic Tracy-Widom Gaussian orthogonal ensemble (Gaussian unitary ensemble) distribution. Moreover, 2D to one-dimensional (1D) crossovers are found also in the asymptotic growth velocity and in the steady-state regime of flat systems, where a family of universal HDs exists, interpolating between the 2D and 1D ones as L_{y}/L_{x} increases. Importantly, the crossover scalings are fully determined and indicate a possible way to solve 2D KPZ models.

Effect of substrate heterogeneity and topology on epithelial tissue growth dynamics

Tissue growth kinetics and interface dynamics depend on the properties of the tissue environment and cell-cell interactions. In cellular environments, substrate heterogeneity and geometry arise from a variety factors, such as the structure of the extracellular matrix and nutrient concentration. We used the CellSim3D model, a kinetic cell division simulator, to investigate the growth kinetics and interface roughness dynamics of epithelial tissue growth on heterogeneous substrates with varying topologies. The results show that the presence of quenched disorder has a clear effect on the colony morphology and the roughness scaling of the interface in the moving interface regime. In a medium with quenched disorder, the tissue interface has a smaller interface roughness exponent, α, and a larger growth exponent, β. The scaling exponents also depend on the topology of the substrate and cannot be categorized by well-known universality classes.

In silico testing of the universality of epithelial tissue growth

The universality of interfacial roughness in growing epithelial tissue has remained a controversial issue. Kardar-Parisi-Zhang (KPZ) and molecular beam epitaxy (MBE) universality classes have been reported among other behaviors including a total lack of universality. Here, we simulate tissues using the cellsim3d kinetic division model for deformable cells to investigate cell-colony scaling. With seemingly minor model changes, it can reproduce both KPZ- and MBE-like scaling in configurations that mimic the respective experiments. Tissue growth with strong cell-cell adhesion in a linear geometry is KPZ like, while weakly adhesive tissues in a radial geometry are MBE like. This result neutralizes the apparent scaling controversy.

Morphologies and dynamics of free surfaces of crystals composed of active particles

Active matter exhibits various collective motions and nonequilibrium phases, such as crystals; however, their surface properties have been poorly explored. Here, we use Brownian dynamics simulations to investigate the surface morphology and dynamics of two-dimensional active crystals during and after growth. For crystal growth on a substrate, the position and roughness of the crystal surface reach steady states at different times. In the steady state, the surface exhibits superdiffusive behaviour at the short time, and the roughness is insensitive to the roughening process and particle activity. We observe two-stage and three-stage surface roughening at different Péclet numbers. The result of dynamic scaling analysis shows that the surface is similar to anomalous roughening, which is distinct from the normal roughening typically found in conventional passive systems. Capillary wave theory for a thermal equilibrium system can describe the active surface fluctuations only in the long-wavelength regime, indicating that active particles mainly drive the surface out of equilibrium locally. These similarities and differences between the active and passive crystal surfaces are essential for understanding active crystals and interfaces.

Morphologies and dynamics of the interfaces between active and passive phases

Active matters exhibit interesting collective behaviors and novel phases, which provide an important platform for the study of nonequilibrium physics. Mixtures of active and passive particles have been intensively studied in motility-induced phase separation, but the morphology of the active-passive interface has been poorly explored. In this work, we investigate the interface morphology in two-dimensional mixtures of active and passive particles using Brownian dynamics simulations. By systematically changing the Péclet number (Pe) and area fraction (ρ), we obtain the phase diagram of the active-passive interface, including rough sharp, rough invasive and flat interdiffusive interfaces. For a sharp interface, dynamic scaling analysis in the propagation stage shows that the roughness exponent α, the growth exponent β, the time exponent κ, and the dynamic exponent z satisfy z = α/(β - κ). Such anomalous scaling indicates that the roughening behavior does not belong to the conventional universality classes with Family-Vicsek scaling for the growth of passive interfaces. On the other hand, the interface in the middle-wavelength regime during the morphology relaxation stage can be described by capillary wave theory. The mean interface position propagates with time as t^1/2, which is robust at different ρ and Pe values in the propagation stage and exhibits superdiffusion in the morphology relaxation stage. These similarities and differences between the active-inactive interfaces and passive interfaces cast light on the interfacial growth of active matter.

Roughness and dynamics of proliferating cell fronts as a probe of cell-cell interactions

Juxtacellular interactions play an essential but still not fully understood role in both normal tissue development and tumour invasion. Using proliferating cell fronts as a model system, we explore the effects of cell-cell interactions on the geometry and dynamics of these one-dimensional biological interfaces. We observe two distinct scaling regimes of the steady state roughness of in-vitro propagating Rat1 fibroblast cell fronts, suggesting different hierarchies of interactions at sub-cell lengthscales and at a lengthscale of 2-10 cells. Pharmacological modulation significantly affects the proliferation speed of the cell fronts, and those modulators that promote cell mobility or division also lead to the most rapid evolution of cell front roughness. By comparing our experimental observations to numerical simulations of elastic cell fronts with purely short-range interactions, we demonstrate that the interactions at few-cell lengthscales play a key role. Our methodology provides a simple framework to measure and characterise the biological effects of such interactions, and could be useful in tumour phenotyping.

Kardar-Parisi-Zhang Interfaces with Curved Initial Shapes and Variational Formula

We study fluctuations of interfaces in the Kardar-Parisi-Zhang (KPZ) universality class with curved initial conditions. By simulations of a cluster growth model and experiments with liquid-crystal turbulence, we determine the universal scaling functions that describe the height distribution and the spatial correlation of the interfaces growing outward from a ring. The scaling functions, controlled by a single dimensionless time parameter, show crossover from the statistical properties of the flat interfaces to those of the circular interfaces. Moreover, employing the KPZ variational formula to describe the case of the ring initial condition, we find that the formula, which we numerically evaluate, reproduces the numerical and experimental results precisely without adjustable parameters. This demonstrates that precise numerical evaluation of the variational formula is possible at all, and underlines the practical importance of the formula, which is able to predict the one-point distribution of KPZ interfaces for general initial conditions.

Effect of heterogeneity and spatial correlations on the structure of a tumor invasion front in cellular environments

Analysis of invasion front has been widely used to decipher biological properties, as well as the growth dynamics of the corresponding populations. Likewise, the invasion front of tumors has been investigated, from which insights into the biological mechanisms of tumor growth have been gained. We develop a model to study how tumors' invasion front depends on the relevant properties of a cellular environment. To do so, we develop a model based on a nonlinear reaction-diffusion equation, the Fisher-Kolmogorov-Petrovsky-Piskunov equation, to model tumor growth. Our study aims to understand how heterogeneity in the cellular environment's stiffness, as well as spatial correlations in its morphology, the existence of both of which has been demonstrated by experiments, affects the properties of tumor invasion front. It is demonstrated that three important factors affect the properties of the front, namely the spatial distribution of the local diffusion coefficients, the spatial correlations between them, and the ratio of the cells' duplication rate and their average diffusion coefficient. Analyzing the scaling properties of tumor invasion front computed by solving the governing equation, we show that, contrary to several previous claims, the invasion front of tumors and cancerous cell colonies cannot be described by the well-known models of kinetic growth, such as the Kardar-Parisi-Zhang equation.

Geometry dependence in linear interface growth

The effect of geometry in the statistics of nonlinear universality classes for interface growth has been widely investigated in recent years, and it is well known to yield a split of them into subclasses. In this work, we investigate this for the linear classes of Edwards-Wilkinson and of Mullins-Herring in one and two dimensions. From comparison of analytical results with extensive numerical simulations of several discrete models belonging to these classes, as well as numerical integrations of the growth equations on substrates of fixed size (flat geometry) or expanding linearly in time (radial geometry), we verify that the height distributions (HDs) and the spatial and the temporal covariances are universal but geometry-dependent. In fact, the HDs are always Gaussian, and, when defined in terms of the so-called "KPZ ansatz" [h≃v_{∞}t+(Γt)^{β}χ], their probability density functions P(χ) have mean null, so that all their cumulants are null, except by their variances, which assume different values in the flat and radial cases. The shape of the (rescaled) covariance curves is analyzed in detail and compared with some existing analytical results for them. Overall, these results demonstrate that the splitting of such university classes is quite general, being not restricted to the nonlinear ones.

Eden model with nonlocal growth rules and kinetic roughening in biological systems

We investigate an off-lattice Eden model where the growth of new cells is performed with a probability dependent on the availability of resources coming externally towards the growing aggregate. The concentration of nutrients necessary for replication is assumed to be proportional to the voids connecting the replicating cells to the outer region, introducing therefore a nonlocal dependence on the replication rule. Our simulations point out that the Kadar-Parisi-Zhang (KPZ) universality class is a transient that can last for long periods in plentiful environments. For conditions of nutrient scarcity, we observe a crossover from regular KPZ to unstable growth, passing by a transient consistent with the quenched KPZ class at the pinning transition. Our analysis sheds light on results reporting on the universality class of kinetic roughening in akin experiments of biological growth.

Kardar-Parisi-Zhang Interfaces with Inward Growth

We study the (1+1)-dimensional Kardar-Parisi-Zhang (KPZ) interfaces growing inward from ring-shaped initial conditions, experimentally and numerically, using growth of a turbulent state in liquid-crystal electroconvection and an off-lattice Eden model, respectively. To realize the ring initial condition experimentally, we introduce a holography-based technique that allows us to design the initial condition arbitrarily. Then, we find that fluctuation properties of ingrowing circular interfaces are distinct from those for the curved or circular KPZ subclass and, instead, are characterized by the flat subclass. More precisely, we find an asymptotic approach to the Tracy-Widom distribution for the Gaussian orthogonal ensemble and the Airy_{1} spatial correlation, as long as time is much shorter than the characteristic time determined by the initial curvature. Near this characteristic time, deviation from the flat KPZ subclass is found, which can be explained in terms of the correlation length and the circumference. Our results indicate that the sign of the initial curvature has a crucial role in determining the universal distribution and correlation functions of the KPZ class.

Pinning-depinning transition in a stochastic growth model for the evolution of cell colony fronts in a disordered medium

We study a stochastic lattice model for cell colony growth, which takes into account proliferation, diffusion, and rotation of cells, in a culture medium with quenched disorder. The medium is composed of sites that inhibit any possible change in the internal state of the cells, representing the disorder, as well as by active medium sites that do not interfere with the cell dynamics. By means of Monte Carlo simulations we find that the velocity of the growing interface, which is taken as the order parameter of the model, strongly depends on the density of active medium sites (ρ_{A}). In fact, the model presents a (continuous) second-order pinning-depinning transition at a certain critical value of ρ_{A}^{crit}, such as, for ρ_{A}>ρ_{A}^{crit}, the interface moves freely across the disordered medium, but for ρ_{A}<ρ_{A}^{crit} the interface becomes irreversible pinned by the disorder. By determining the relevant critical exponents, our study reveals that within the depinned phase the interface can be rationalized in terms of the Kardar-Parisi-Zhang universality class, but when approaching the critical threshold, the nonlinear term of the Kardar-Parisi-Zhang equation tends to vanish and then the pinned interface belongs to the quenched Edwards-Wilkinson universality class.

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.

The universal growth rate behavior and regime transition in adherent cell colonies

In this work, we used five cell lineages, cultivated in vitro, to show they follow a common functional form to the growth rate: a sigmoidal curve, suggesting that competition and cooperation (usual mechanisms for systems with this behavior) might be present. Both theoretical and experimental investigations, on the causes of this behavior, are challenging for the research field; since the sigmoidal form to the growth rate seems to absorb important properties of such systems, e.g., cell deformation and statistical interactions. We shed some light on this subject by showing how cell spreading affects the radius behavior of the growing colonies. Doing numerical time derivatives of the experimental data, we obtained the growth rates. Using reduced variables for the time and rates, we obtained the collapse of all colonies growth rates onto one curve with sigmoidal shape. This suggests a universal-type behavior, with regime transition related to a morphological transition of adherent cell colonies.

Simulation of avascular tumor growth by agent-based game model involving phenotype-phenotype interactions

All tumors, both benign and metastatic, undergo an avascular growth stage with nutrients supplied by the surrounding tissue. This avascular growth process is much easier to carry out in more qualitative and quantitative experiments starting from tumor spheroids in vitro with reliable reproducibility. Essentially, this tumor progression would be described as a sequence of phenotypes. Using agent-based simulation in a two-dimensional spatial lattice, we constructed a composite growth model in which the phenotypic behavior of tumor cells depends on not only the local nutrient concentration and cell count but also the game among cells. Our simulation results demonstrated that in silico tumors are qualitatively similar to those observed in tumor spheroid experiments. We also found that the payoffs in the game between two living cell phenotypes can influence the growth velocity and surface roughness of tumors at the same time. Finally, this current model is flexible and can be easily extended to discuss other situations, such as environmental heterogeneity and mutation.

Kinetic interfaces of patchy particles

We study the irreversible adsorption of patchy particles on substrates in the limit of advective mass transport. Recent numerical results show that the interface roughening depends strongly on the particle attributes, such as, patch-patch correlations, bond flexibility and strength of the interactions, uncovering new absorbing phase transitions. Here, we revisit these results and discuss in detail the transitions. In particular, we present new evidence that the tricritical point, observed in systems of particles with flexible patches, is in the tricritical directed percolation universality class. A scaling analysis of the time evolution of the correlation length for the aggregation of patchy particles with distinct bonding energies confirms that the critical regime is in the Kardar-Parisi-Zhang with quenched disorder universality class.

Skewness in (1+1)-dimensional Kardar-Parisi-Zhang-type growth

We use the (1+1)-dimensional Kardar-Parisi-Zhang equation driven by a Gaussian white noise and employ the dynamic renormalization-group of Yakhot and Orszag without rescaling [J. Sci. Comput. 1, 3 (1986)]. Hence we calculate the second- and third-order moments of height distribution using the diagrammatic method in the large-scale and long-time limits. The moments so calculated lead to the value S=0.3237 for the skewness. This value is comparable with numerical and experimental estimates.

Dynamic scaling analysis of two-dimensional cell colony fronts in a gel medium: a biological system approaching a quenched Kardar-Parisi-Zhang universality

The interfacial two-dimensional spreading dynamics of quasilinear Vero cell colony fronts in methylcellulose (MC)-containing culture medium, under a constant average front displacement velocity regime, was investigated. Under comparable experimental conditions, the average colony front displacement velocity becomes lower than that reported for a standard culture medium. Initially, the presence of MC in the medium hinders both the colony spreading, due to a gradual change in the average size and shape of cells and their distribution in the colony, and the cell motility in the gelled medium. Furthermore, at longer culture times enlarged cells appear at random in the border region of the colony. These cells behave as obstacles (pinning sites) for the displacement of smaller cells towards the colony front. The dynamic scaling analysis of rough fronts yields the set of exponents α=0.63±0.04,β=0.75±0.05, and z=0.84±0.05, which is close to that expected for a quenched Kardar-Parisi-Zhang model.

Influence of individual cell motility on the 2D front roughness dynamics of tumour cell colonies

The dynamics of in situ 2D HeLa cell quasi-linear and quasi-radial colony fronts in a standard culture medium is investigated. For quasi-radial colonies, as the cell population increased, a kinetic transition from an exponential to a constant front average velocity regime was observed. Special attention was paid to individual cell motility evolution under constant average colony front velocity looking for its impact on the dynamics of the 2D colony front roughness. From the directionalities and velocity components of cell trajectories in colonies with different cell populations, the influence of both local cell density and cell crowding effects on individual cell motility was determined. The average dynamic behaviour of individual cells in the colony and its dependence on both local spatio-temporal heterogeneities and growth geometry suggested that cell motion undergoes under a concerted cell migration mechanism, in which both a limiting random walk-like and a limiting ballistic-like contribution were involved. These results were interesting to infer how biased cell trajectories influenced both the 2D colony spreading dynamics and the front roughness characteristics by local biased contributions to individual cell motion. These data are consistent with previous experimental and theoretical cell colony spreading data and provide additional evidence of the validity of the Kardar-Parisi-Zhang equation, within a certain range of time and colony front size, for describing the dynamics of 2D colony front roughness.

Crossover from growing to stationary interfaces in the Kardar-Parisi-Zhang class

This Letter reports on how the interfaces in the (1+1)-dimensional Kardar-Parisi-Zhang (KPZ) class undergo, in the course of time, a transition from the flat, growing regime to the stationary one. Simulations of the polynuclear growth model and experiments on turbulent liquid crystal reveal universal functions of the KPZ class governing this transition, which connect the distribution and correlation functions for the growing and stationary regimes. This in particular shows how interfaces realized in experiments and simulations actually approach the stationary regime, which is never attained unless a stationary interface is artificially given as an initial condition.

Arraying cell cultures using PEG-DMA micromolding in standard culture dishes

A robust and effortless procedure is presented, which allows for the microstructuring of standard cell culture dishes. Cell adhesion and proliferation are controlled by three-dimensional poly(ethylene glycol)-dimethacrylate (PEG-DMA) microstructures. The spacing between microwells can be extended to millimeter size in order to enable the combination with robotic workstations. Cell arrays of microcolonies can be studied under boundary-free growth conditions by lift-off of the PEG-DMA layer in which the growth rate is accessible via the evolution of patch areas. Alternatively, PEG-DMA stencils can be used as templates for plasma-induced patterning.

Reproduction-time statistics and segregation patterns in growing populations

Pattern formation in microbial colonies of competing strains under purely space-limited population growth has recently attracted considerable research interest. We show that the reproduction time statistics of individuals has a significant impact on the sectoring patterns. Generalizing the standard Eden growth model, we introduce a simple one-parameter family of reproduction time distributions indexed by the variation coefficient δ∈[0,1], which includes deterministic (δ=0) and memory-less exponential distribution (δ=1) as extreme cases. We present convincing numerical evidence and heuristic arguments that the generalized model is still in the Kardar-Parisi-Zhang (KPZ) universality class, and the changes in patterns are due to changing prefactors in the scaling relations, which we are able to predict quantitatively. With the example of Saccharomyces cerevisiae, we show that our approach using the variation coefficient also works for more realistic reproduction time distributions.

Growth dynamics of cancer cell colonies and their comparison with noncancerous cells

The two-dimensional (2D) growth dynamics of HeLa (cervix cancer) cell colonies was studied following both their growth front and the pattern morphology evolutions utilizing large population colonies exhibiting linearly and radially spreading fronts. In both cases, the colony profile fractal dimension was d(f)=1.20±0.05 and the growth fronts displaced at the constant velocity 0.90±0.05 μm min(-1). Colonies showed changes in both cell morphology and average size. As time increased, the formation of large cells at the colony front was observed. Accordingly, the heterogeneity of the colony increased and local driving forces that set in began to influence the dynamics of the colony front. The dynamic scaling analysis of rough colony fronts resulted in a roughness exponent α = 0.50±0.05, a growth exponent β = 0.32±0.04, and a dynamic exponent z=1.5±0.2. The validity of this set of scaling exponents extended from a lower cutoff l(c)≈60 μm upward, and the exponents agreed with those predicted by the standard Kardar-Parisi-Zhang continuous equation. HeLa data were compared with those previously reported for Vero cell colonies. The value of d(f) and the Kardar-Parisi-Zhang-type 2D front growth dynamics were similar for colonies of both cell lines. This indicates that the cell colony growth dynamics is independent of the genetic background and the tumorigenic nature of the cells. However, one can distinguish some differences between both cell lines during the growth of colonies that may result from specific cooperative effects and the nature of each biosystem.

Dynamics and morphology characteristics of cell colonies with radially spreading growth fronts

The dynamics of two-dimensional (2D) radially spreading growth fronts of Vero cell colonies was investigated utilizing two types of colonies, namely type I starting from clusters with a small number of cells, which initially exhibited arbitrary-shaped rough growth fronts and progressively approached quasicircular ones as the cell population increased; and type II colonies, starting from a relatively large circular three-dimensional (3D) cell cluster. For large cell population colonies, the fractal dimension of the fronts was D(F) = 1.20±0.05. For low cell populations, the mean colony radius increased exponentially with time, but for large ones the constant radial front velocity 0.20±0.02 μm min(-1) was reached. Colony spreading was accompanied by changes in both cell morphology and average size, and by the formation of very large cells, some of them multinuclear. Therefore the heterogeneity of colonies increased and local driving forces that set in began to influence the 2D growth front kinetics. The retardation effect related to the exponential to constant radial front velocity transition was assigned to a number of possible interferences including the cell duplication and 3D growth in the bulk of the colony. The dynamic scaling analysis of overhang-corrected rough colony fronts, after arc-radius coordinate system transformation, resulted in roughness exponent α = 0.50±0.05 and growth exponent β = 0.32±0.04, for arc lengths greater than 100 μm. This set of scaling exponents agreed with that predicted by the Kardar, Parisi, and Zhang continuous equation. For arc lengths shorter than 2-3 cell diameters, the value α = 0.85±0.05 would be related to a cell front roughening caused by temporarily membrane deformations occasionally interfered by cell proliferation.