The model muddle: in search of tumour growth laws

In this article we shall trace the historical development of tumour growth laws, which in a quantitative fashion describe the increase in tumour mass/volume over time. These models are usually formulated in terms of differential equations that relate the growth rate of the tumour to its current state, and range from the simple one-parameter exponential growth model, to more advanced models that contain a large number of parameters. Understanding the assumptions and consequences of such models is important, since they often underpin more complex models of tumour growth. The conclusion of this brief survey is that although much improvement has occurred over the last century, more effort and new models are required if we are to understand the intricacies of tumour growth

The impact of anticipation in dynamical systems

Collective motion in biology is often modelled as a dynamical system, in which individuals are represented as particles whose interactions are determined by the current state of the system. Many animals, however, including humans, have predictive capabilities, and presumably base their behavioural decisions---at least partially---upon an anticipated state of their environment. We explore a minimal version of this idea in the context of particles that interact according to a pairwise potential. Anticipation enters the picture by calculating the interparticle forces from linear extrapolations of the particle positions some time $\tau$ into the future. Simulations show that for intermediate values of $\tau$, compared to a transient time scale defined by the potential and the initial conditions, the particles form rotating clusters in which the particles are arranged in a hexagonal pattern. Analysis of the system shows that anticipation induces energy dissipation and we show that the kinetic energy asymptotically decays as $1/t$. Furthermore, we show that the angular momentum is not necessarily conserved for $\tau >0$, and that asymmetries in the initial condition therefore can cause rotational movement. These results suggest that anticipation could play an important role in collective behaviour, since it induces pattern formation and stabilises the dynamics of the system.Comment: Major revision compared to previous version. All figures replaced. Only introduction and discussion remain intac

Time scales and wave formation in non-linear spatial public goods games

The co-evolutionary dynamics of competing populations can be strongly affected by frequency-dependent selection and spatial population structure. As co-evolving populations grow into a spatial domain, their initial spatial arrangement and their growth rate differences are important factors that determine the long-term outcome. We here model producer and free-rider co-evolution in the context of a diffusive public good (PG) that is produced by the producers at a cost but evokes local concentration-dependent growth benefits to all. The benefit of the PG can be non-linearly dependent on public good concentration. We consider the spatial growth dynamics of producers and free-riders in one, two and three dimensions by modeling producer cell, free-rider cell and public good densities in space, driven by the processes of birth, death and diffusion (cell movement and public good distribution). Typically, one population goes extinct, but the time-scale of this process varies with initial conditions and the growth rate functions. We establish that spatial variation is transient regardless of dimensionality, and that structured initial conditions lead to increasing times to get close to an extinction state, called ε-extinction time. Further, we find that uncorrelated initial spatial structures do not influence this ε-extinction time in comparison to a corresponding well-mixed (non-spatial) system. In order to estimate the ε-extinction time of either free-riders or producers we derive a slow manifold solution. For invading populations, i.e. for populations that are initially highly segregated, we observe a traveling wave, whose speed can be calculated. Our results provide quantitative predictions for the transient spatial dynamics of cooperative traits under pressure of extinction

Stability Analysis of a Hybrid Cellular Automaton Model of Cell Colony Growth

Cell colonies of bacteria, tumour cells and fungi, under nutrient limited growth conditions, exhibit complex branched growth patterns. In order to investigate this phenomenon we present a simple hybrid cellular automaton model of cell colony growth. In the model the growth of the colony is limited by a nutrient that is consumed by the cells and which inhibits cell division if it falls below a certain threshold. Using this model we have investigated how the nutrient consumption rate of the cells affects the growth dynamics of the colony. We found that for low consumption rates the colony takes on a Eden-like morphology, while for higher consumption rates the morphology of the colony is branched with a fractal geometry. These findings are in agreement with previous results, but the simplicity of the model presented here allows for a linear stability analysis of the system. By observing that the local growth of the colony is proportional to the flux of the nutrient we derive an approximate dispersion relation for the growth of the colony interface. This dispersion relation shows that the stability of the growth depends on how far the nutrient penetrates into the colony. For low nutrient consumption rates the penetration distance is large, which stabilises the growth, while for high consumption rates the penetration distance is small, which leads to unstable branched growth. When the penetration distance vanishes the dispersion relation is reduced to the one describing Laplacian growth without ultra-violet regularisation. The dispersion relation was verified by measuring how the average branch width depends on the consumption rate of the cells and shows good agreement between theory and simulations.Comment: 8 pages, 6 figure

Neighborhood size-effects shape growing population dynamics in evolutionary public goods games

An evolutionary game emerges when a subset of individuals incur costs to provide benefits to all individuals. Public goods games (PGG) cover the essence of such dilemmas in which cooperators are prone to exploitation by defectors. We model the population dynamics of a non-linear\ua0PGG and consider density-dependence on the global level, while the game occurs within local neighborhoods. At low cooperation, increases in the public good provide increasing returns. At high cooperation, increases provide diminishing returns. This mechanism leads to diverse evolutionarily stable strategies, including monomorphic and polymorphic populations, and neighborhood-size-driven state changes, resulting in hysteresis between equilibria. Stochastic or strategy-dependent variations in neighborhood sizes favor coexistence by destabilizing monomorphic states. We integrate our model with experiments of cancer cell growth and confirm that our framework describes PGG dynamics observed in cellular populations. Our findings advance the understanding of how neighborhood-size effects in PGG shape the dynamics of growing populations. \ua9 2019, The Author(s)

Autocrine signaling can explain the emergence of Allee effects in cancer cell populations

In many human cancers, the rate of cell growth depends crucially on the size of the tumour cell population. Low, zero, or negative growth at low population densities is known as the Allee effect; this effect has been studied extensively in ecology, but so far lacks a good explanation in the cancer setting. Here, we formulate and analyze an individual-based model of cancer, in which cell division rates are increased by the local concentration of an autocrine growth factor produced by the cancer cells themselves. We show, analytically and by simulation, that autocrine signaling suffices to cause both strong and weak Allee effects. Whether low cell densities lead to negative (strong effect) or reduced (weak effect) growth rate depends directly on the ratio of cell death to proliferation, and indirectly on cellular dispersal. Our model is consistent with experimental observations from three patient-derived brain tumor cell lines grown at different densities. We propose that further studying and quantifying population-wide feedback, impacting cell growth, will be central for advancing our understanding of cancer dynamics and treatment, potentially exploiting Allee effects for therapy

Structural correlations in bacterial metabolic networks

<p>Abstract</p> <p>Background</p> <p>Evolution of metabolism occurs through the acquisition and loss of genes whose products acts as enzymes in metabolic reactions, and from a presumably simple primordial metabolism the organisms living today have evolved complex and highly variable metabolisms. We have studied this phenomenon by comparing the metabolic networks of 134 bacterial species with known phylogenetic relationships, and by studying a neutral model of metabolic network evolution.</p> <p>Results</p> <p>We consider the 'union-network' of 134 bacterial metabolisms, and also the union of two smaller subsets of closely related species. Each reaction-node is tagged with the number of organisms it belongs to, which we denote organism degree (OD), a key concept in our study. Network analysis shows that common reactions are found at the centre of the network and that the average OD decreases as we move to the periphery. Nodes of the same OD are also more likely to be connected to each other compared to a random OD relabelling based on their occurrence in the real data. This trend persists up to a distance of around five reactions. A simple growth model of metabolic networks is used to investigate the biochemical constraints put on metabolic-network evolution. Despite this seemingly drastic simplification, a 'union-network' of a collection of unrelated model networks, free of any selective pressure, still exhibit similar structural features as their bacterial counterpart.</p> <p>Conclusions</p> <p>The OD distribution quantifies topological properties of the evolutionary history of bacterial metabolic networks, and lends additional support to the importance of horizontal gene transfer during bacterial metabolic evolution where new reactions are attached at the periphery of the network. The neutral model of metabolic network growth can reproduce the main features of real networks, but we observe that the real networks contain a smaller common core, while they are more similar at the periphery of the network. This suggests that natural selection and biochemical correlations can act both to diversify and to narrow down metabolic evolution.</p

The contribution of age structure to cell population responses to targeted therapeutics

Cells grown in culture act as a model system for analyzing the effects of anticancer compounds, which may affect cell behavior in a cell cycle position-dependent manner. Cell synchronization techniques have been generally employed to minimize the variation in cell cycle position. However, synchronization techniques are cumbersome and imprecise and the agents used to synchronize the cells potentially have other unknown effects on the cells. An alternative approach is to determine the age structure in the population and account for the cell cycle positional effects post hoc. Here we provide a formalism to use quantifiable age distributions from live cell microscopy experiments to parameterize an age-structured model of cell population response

Using DNA Methylation Patterns to Infer Tumor Ancestry

Background: Exactly how human tumors grow is uncertain because serial observations are impractical. One approach to reconstruct the histories of individual human cancers is to analyze the current genomic variation between its cells. The greater the variations, on average, the greater the time since the last clonal evolution cycle (‘‘a molecular clock hypothesis’’). Here we analyze passenger DNA methylation patterns from opposite sides of 12 primary human colorectal cancers (CRCs) to evaluate whether the variation (pairwise distances between epialleles) is consistent with a single clonal expansion after transformation. Methodology/Principal Findings: Data from 12 primary CRCs are compared to epigenomic data simulated under a single clonal expansion for a variety of possible growth scenarios. We find that for many different growth rates, a single clonal expansion can explain the population variation in 11 out of 12 CRCs. In eight CRCs, the cells from different glands are all equally distantly related, and cells sampled from the same tumor half appear no more closely related than cells sampled from opposite tumor halves. In these tumors, growth appears consistent with a single ‘‘symmetric’ ’ clonal expansion. In three CRCs, the variation in epigenetic distances was different between sides, but this asymmetry could be explained by a single clonal expansion with one region of a tumor having undergone more cell division than the other. The variation in one CRC was complex and inconsistent with a simple single clonal expansion