Epithelial tissue architecture protects against cancer

The role of spatial structures of tissues in cancer initiation dynamics

It is widely believed that biological tissues evolved to lower the risks of cancer development. One of the specific ways to minimize the chances of tumor formation comes from proper spatial organization of tissues. However, the microscopic mechanisms of underlying processes remain not fully understood. We present a theoretical investigation on the role of spatial structures in cancer initiation dynamics. In our approach, the dynamics of single mutation fixations are analyzed using analytical calculations and computer simulations by mapping them to Moran processes on graphs with different connectivity that mimic various spatial structures. It is found that while the fixation probability is not affected by modifying the spatial structures of the tissues, the fixation times can change dramatically. The slowest dynamics is observed in "quasi-one-dimensional" structures, while the fastest dynamics is observed in "quasi-three-dimensional" structures. Theoretical calculations also suggest that there is a critical value of the degree of graph connectivity, which mimics the spatial dimension of the tissue structure, above which the spatial structure of the tissue has no effect on the mutation fixation dynamics. An effective discrete-state stochastic model of cancer initiation is utilized to explain our theoretical results and predictions. Our theoretical analysis clarifies some important aspects on the role of the tissue spatial structures in the cancer initiation processes.

Bernoulli and binomial proliferation on evolutionary graphs

In this paper we introduce random proliferation models on graphs. We consider two types of particles: type-1/mutant/invader/red particles proliferates on a population of type-2/wild-type/resident/blue particles. Unlike the well-known Moran model on graphs -as introduced in Lieberman et al. (2005)-, type-1 particles can occupy in a single iteration several neighbouring sites previously occupied by type-2 particles. Two variants are considered, depending on the random distribution involving the proliferation mechanism: Bernoulli and binomial proliferation. By comparison with fixation probability of type-1 particles in the Moran process, critical parameters are introduced. Properties of proliferation are studied and some particular cases are analytically solved. Finally, by updating the parameters that drive the processes through a density-dependent mechanism, it is possible to capture additional relevant features as fluctuating waves of type-1 particles over long periods of time. In fact, the models can be adapted to tackle more general, complex and realistic situations.

The evolution of metapopulation dynamics and the number of stem cells in intestinal crypts and other tissue structures in multicellular bodies

Carcinogenesis is a process of somatic evolution. Previous models of stem and transient amplifying cells in epithelial proliferating units like colonic crypts showed that intermediate numbers of stem cells in a crypt should optimally prevent progression to cancer. If a stem cell population is too small, it is easy for a mutator mutation to drift to fixation. If it is too large, it is easy for selection to drive cell fitness enhancing carcinogenic mutations to fixation. Here, we show that a multiscale microsimulation, that captures both within-crypt and between-crypt evolutionary dynamics, leads to a different conclusion. Epithelial tissues are metapopulations of crypts. We measured time to initiation of a neoplasm, implemented as inactivation of both alleles of a tumor suppressor gene. In our model, time to initiation is dependent on the spread of mutator clones in the crypts. The proportion of selectively beneficial and deleterious mutations in somatic cells is unknown and so was explored with a parameter. When the majority of non-neutral mutations are deleterious, the fitness of mutator clones tends to decline. When crypts are maintained by few stem cells, intercrypt competition tends to remove crypts with fixed mutators. When there are many stem cells within a crypt, there is virtually no crypt turnover, but mutator clones are suppressed by within-crypt competition. If the majority of non-neutral mutations are beneficial to the clone, then these results are reversed and intermediate-sized crypts provide the most protection against initiation. These results highlight the need to understand the dynamics of turnover and the mechanisms that control homeostasis, both at the level of stem cells within proliferative units and at the tissue level of competing proliferative units. Determining the distribution of fitness effects of somatic mutations will also be crucial to understanding the dynamics of tumor initiation and progression.

Optimizing homeostatic cell renewal in hierarchical tissues

In order to maintain homeostasis, mature cells removed from the top compartment of hierarchical tissues have to be replenished by means of differentiation and self-renewal events happening in the more primitive compartments. As each cell division is associated with a risk of mutation, cell division patterns have to be optimized, in order to minimize or delay the risk of malignancy generation. Here we study this optimization problem, focusing on the role of division tree length, that is, the number of layers of cells activated in response to the loss of terminally differentiated cells, which is related to the balance between differentiation and self-renewal events in the compartments. Using both analytical methods and stochastic simulations in a metapopulation-style model, we find that shorter division trees are advantageous if the objective is to minimize the total number of one-hit mutants in the cell population. Longer division trees on the other hand minimize the accumulation of two-hit mutants, which is a more likely evolutionary goal given the key role played by tumor suppressor genes in cancer initiation. While division tree length is the most important property determining mutant accumulation, we also find that increasing the size of primitive compartments helps to delay two-hit mutant generation.

Cells in multicellular organisms are organized hierarchically. A stem cell gives rise to a chain of dividing and progressively differentiating offspring. At the end of this chain (called a lineage) are terminally differentiated cells that perform their function and undergo programmed cell death, to be replaced by new divisions of less differentiated cells. Here we are interested in the design of such lineages. At one extreme, one can imagine that a loss of terminally differentiated cells only results in divisions of cells in close hierarchical proximity to them, giving rise to very short division trees. On the other hand, it is possible that a long chain of increasingly primitive cells gets activated in response to the loss of differentiated cells. We expect that an important type of selection pressure acting upon tissue design is the minimization of mutations that happen in the course of everyday tissue maintenance (homeostasis). For example, tumor suppressor gene inactivation (two consecutive mutations) is an early rate-limiting step in many cancers. Using mathematical and computational methods, we find that the length of division trees is anti-correlated with the likelihood of double mutations, and lengthening the trees may provide an evolutionary advantage to the organism by delaying the onset of cancer.

The role of cell location and spatial gradients in the evolutionary dynamics of colon and intestinal crypts

BACKGROUND

Colon and intestinal crypts serve as an important model system for adult stem cell proliferation and differentiation. We develop a spatial stochastic model to study the rate of somatic evolution in a normal crypt, focusing on the production of two-hit mutants that inactivate a tumor suppressor gene. We investigate the effect of cell division pattern along the crypt on mutant production, assuming that the division rate of each cell depends on its location.

RESULTS

We find that higher probability of division at the bottom of the crypt, where the stem cells are located, leads to a higher rate of double-hit mutant production. The optimal case for delaying mutations occurs when most of the cell divisions happen at the top of the crypt. We further consider an optimization problem where the "evolutionary" penalty for double-hit mutant generation is complemented with a "functional" penalty that assures that fully differentiated cells at the top of the crypt cannot divide.

CONCLUSION

The trade-off between the two types of objectives leads to the selection of an intermediate division pattern, where the cells in the middle of the crypt divide with the highest rate. This matches the pattern of cell divisions obtained experimentally in murine crypts.

REVIEWERS

This article was reviewed by David Axelrod (nominated by an Editorial Board member, Marek Kimmel), Yang Kuang and Anna Marciniak-Czochra. For the full reviews, please go to the Reviewers' comments section.

Should tissue structure suppress or amplify selection to minimize cancer risk?

Background

It has been frequently argued that tissues evolved to suppress the accumulation of growth enhancing cancer inducing mutations. A prominent example is the hierarchical structure of tissues with high cell turnover, where a small number of tissue specific stem cells produces a large number of specialized progeny during multiple differentiation steps. Another well known mechanism is the spatial organization of stem cell populations and it is thought that this organization suppresses fitness enhancing mutations. However, in small populations the suppression of advantageous mutations typically also implies an increased accumulation of deleterious mutations. Thus, it becomes an important question whether the suppression of potentially few advantageous mutations outweighs the combined effects of many deleterious mutations.

Results

We argue that the distribution of mutant fitness effects, e.g. the probability to hit a strong driver compared to many deleterious mutations, is crucial for the optimal organization of a cancer suppressing tissue architecture and should be taken into account in arguments for the evolution of such tissues.

Conclusion

We show that for systems that are composed of few cells reflecting the typical organization of a stem cell niche, amplification or suppression of selection can arise from subtle changes in the architecture. Moreover, we discuss special tissue structures that can suppress most types of non-neutral mutations simultaneously.

Reviewers

This article was reviewed by Benjamin Allen, Andreas Deutsch and Ignacio Rodriguez-Brenes. For the full reviews, please go to the Reviewers’ comments section.

A Hierarchical Probability Model of Colon Cancer

We consider a model of fixed sizeN= 2lin which there arelgenerations of daughter cells and a stem cell. In each generationithere are 2i−1daughter cells. At each integral time unit the cells split so that the stem cell splits into a stem cell and generation 1 daughter cell and the generationidaughter cells become two cells of generationi+1. The last generation is removed from the population. A stem cell acquires first and second mutations at ratesu1andu2, and a daughter cell acquires first and second mutations at ratesv1andv2. We find the distribution for the time it takes to acquire two mutations asNgoes to ∞ and the mutation rates go to 0. The mutation rates may tend to 0 at different speeds. We also find the distribution for the locations of the mutations. In particular, we determine whether or not the mutations occur on a stem cell and if not, at what generation in the daughter cells they occur. Several outcomes are possible, depending on how fast the rates go to 0. The model considered has been proposed by Komarova (2007) as a model for colon cancer.

A Hierarchical Probability Model of Colon Cancer

We consider a model of fixed sizeN= 2lin which there arelgenerations of daughter cells and a stem cell. In each generationithere are 2i−1daughter cells. At each integral time unit the cells split so that the stem cell splits into a stem cell and generation 1 daughter cell and the generationidaughter cells become two cells of generationi+1. The last generation is removed from the population. A stem cell acquires first and second mutations at ratesu1andu2, and a daughter cell acquires first and second mutations at ratesv1andv2. We find the distribution for the time it takes to acquire two mutations asNgoes to ∞ and the mutation rates go to 0. The mutation rates may tend to 0 at different speeds. We also find the distribution for the locations of the mutations. In particular, we determine whether or not the mutations occur on a stem cell and if not, at what generation in the daughter cells they occur. Several outcomes are possible, depending on how fast the rates go to 0. The model considered has been proposed by Komarova (2007) as a model for colon cancer.

Counterintuitive properties of the fixation time in network-structured populations

Evolutionary dynamics on graphs can lead to many interesting and counterintuitive findings. We study the Moran process, a discrete time birth-death process, that describes the invasion of a mutant type into a population of wild-type individuals. Remarkably, the fixation probability of a single mutant is the same on all regular networks. But non-regular networks can increase or decrease the fixation probability. While the time until fixation formally depends on the same transition probabilities as the fixation probabilities, there is no obvious relation between them. For example, an amplifier of selection, which increases the fixation probability and thus decreases the number of mutations needed until one of them is successful, can at the same time slow down the process of fixation. Based on small networks, we show analytically that (i) the time to fixation can decrease when links are removed from the network and (ii) the node providing the best starting conditions in terms of the shortest fixation time depends on the fitness of the mutant. Our results are obtained analytically on small networks, but numerical simulations show that they are qualitatively valid even in much larger populations.

Stem cell regulation: Implications when differentiated cells regulate symmetric stem cell division

We use a mathematical model to show that if symmetric stem cell division is regulated by differentiated cells, then changes in the population dynamics of the differentiated cells can lead to changes in the population dynamics of the stem cells. More precisely, the relative fitness of the stem cells can be affected by modifying the death rate of the differentiated cells. This result is interesting because stem cells are less sensitive than differentiated cells to environmental factors, such as medical therapy. Our result implies that stem cells can be manipulated indirectly by medical treatments that target the differentiated cells.

A calibrated agent-based computer model of stochastic cell dynamics in normal human colon crypts useful for in silico experiments

BACKGROUND

Normal colon crypts consist of stem cells, proliferating cells, and differentiated cells. Abnormal rates of proliferation and differentiation can initiate colon cancer. We have measured the variation in the number of each of these cell types in multiple crypts in normal human biopsy specimens. This has provided the opportunity to produce a calibrated computational model that simulates cell dynamics in normal human crypts, and by changing model parameter values, to simulate the initiation and treatment of colon cancer.

RESULTS

An agent-based model of stochastic cell dynamics in human colon crypts was developed in the multi-platform open-source application NetLogo. It was assumed that each cell's probability of proliferation and probability of death is determined by its position in two gradients along the crypt axis, a divide gradient and in a die gradient. A cell's type is not intrinsic, but rather is determined by its position in the divide gradient. Cell types are dynamic, plastic, and inter-convertible. Parameter values were determined for the shape of each of the gradients, and for a cell's response to the gradients. This was done by parameter sweeps that indicated the values that reproduced the measured number and variation of each cell type, and produced quasi-stationary stochastic dynamics. The behavior of the model was verified by its ability to reproduce the experimentally observed monocolonal conversion by neutral drift, the formation of adenomas resulting from mutations either at the top or bottom of the crypt, and by the robust ability of crypts to recover from perturbation by cytotoxic agents. One use of the virtual crypt model was demonstrated by evaluating different cancer chemotherapy and radiation scheduling protocols.

CONCLUSIONS

A virtual crypt has been developed that simulates the quasi-stationary stochastic cell dynamics of normal human colon crypts. It is unique in that it has been calibrated with measurements of human biopsy specimens, and it can simulate the variation of cell types in addition to the average number of each cell type. The utility of the model was demonstrated with in silico experiments that evaluated cancer therapy protocols. The model is available for others to conduct additional experiments.

Stem cell control, oscillations, and tissue regeneration in spatial and non-spatial models

Normal human tissue is organized into cell lineages, in which the highly differentiated mature cells that perform tissue functions are the end product of an orderly tissue-specific sequence of divisions that start with stem cells or progenitor cells. Tissue homeostasis and effective regeneration after injuries requires tight regulation of these cell lineages and feedback loops play a fundamental role in this regard. In particular, signals secreted from differentiated cells that inhibit stem cell division and stem cell self-renewal are important in establishing control. In this article we study in detail the cell dynamics that arise from this control mechanism. These dynamics are fundamental to our understanding of cancer, given that tumor initiation requires an escape from tissue regulation. Knowledge on the processes of cellular control can provide insights into the pathways that lead to deregulation and consequently cancer development.

Dynamics of mutant cells in hierarchical organized tissues

Most tissues in multicellular organisms are maintained by continuous cell renewal processes. However, high turnover of many cells implies a large number of error-prone cell divisions. Hierarchical organized tissue structures with stem cell driven cell differentiation provide one way to prevent the accumulation of mutations, because only few stem cells are long lived. We investigate the deterministic dynamics of cells in such a hierarchical multi compartment model, where each compartment represents a certain stage of cell differentiation. The dynamics of the interacting system is described by ordinary differential equations coupled across compartments. We present analytical solutions for these equations, calculate the corresponding extinction times and compare our results to individual based stochastic simulations. Our general compartment structure can be applied to different tissues, as for example hematopoiesis, the epidermis, or colonic crypts. The solutions provide a description of the average time development of stem cell and non stem cell driven mutants and can be used to illustrate general and specific features of the dynamics of mutant cells in such hierarchically structured populations. We illustrate one possible application of this approach by discussing the origin and dynamics of PIG-A mutant clones that are found in the bloodstream of virtually every healthy adult human. From this it is apparent, that not only the occurrence of a mutant but also the compartment of origin is of importance.

We investigate the average stem cell driven dynamics of cell counts in an abstract multi compartment model. Within this framework one can represent different tissue structures, as for example hematopoiesis, the skin or the colonic crypt. Our analysis is based on an individual cell model in which cells can differentiate, reproduce or die. We give closed solutions to the corresponding system of coupled differential equations, that describe the average dynamics of all cell types. There are three cases of interest: (i) Mutations at the stem cell level, (ii) Mutations in downstream compartments associated with more mature, non stem cell types, (iii) Mutations in downstream compartments with cells acquiring stem cell like properties. The average dynamics shows for (i) and (iii) an increase of mutants towards an equilibrium, in case (ii) the average mutant cell count goes through a maximum, but mutants die out in the long run. We calculate the corresponding extinction times for every compartment. We discuss applications to hematopoietic diseases such as, PIG-A mutant cells or the classic oncogene BCR-ABL. Although the abstract model is a simplified sketch of cell differentiation, it is capable of describing many aspects of a wide variety of such tissues and associated diseases.

Organ aging and susceptibility to cancer may be related to the geometry of the stem cell niche

Telomere loss at each cell replication limits the proliferative capacity of normal cells, including adult stem cells. Entering replicative senescence protects dividing cells from neoplastic transformation, but also contributes to aging of the tissue. Recent experiments have shown that intestinal mouse stem cells divide symmetrically, at random make decisions to remain stem cells or to differentiate, and gradually lose telomeric DNA. A cell's decision whether to differentiate or to remain a stem cell depends on the local cellular and chemical environment and thus tissue architecture is expected to play role in cell proliferation dynamics. To take into account the structure of the stem cell niche in determining its proliferative potential and susceptibility to cancer, a theoretical model is introduced and the niche proliferative potential is quantified for different architectures. The niche proliferative potential is quantitatively related to the proliferative potential of the individual stem cells for different structural classes of the stem cell niche. Stem cells at the periphery of a niche are under pressure to divide and to differentiate, as well as to maintain the stem cell niche boundary, and thus the geometry of the stem cell niche is expected to play a role in determining the stem cell division sequence and differentiation. Smaller surface-to-volume ratio is associated with higher susceptibility to cancer, higher tissue renewal capacity, and decreased aging rate. Several testable experimental predictions are discussed, as well the presence of stochastic effects.

Pathways to tumorigenesis--modeling mutation acquisition in stem cells and their progeny

Most adult tissues consist of stem cells, progenitors, and mature cells, and this hierarchical architecture may play an important role in the multistep process of carcinogenesis. Here, we develop and discuss the important predictions of a simple mathematical model of cancer initiation and early progression within a hierarchically structured tissue. This work presents a model that incorporates both the sequential acquisition of phenotype altering mutations and tissue hierarchy. The model simulates the progressive effect of accumulating mutations that lead to an increase in fitness or the induction of genetic instability. A novel aspect of the model is that symmetric self-renewal, asymmetric division, and differentiation are all incorporated, and this enables the quantitative study of the effect of mutations that deregulate the normal, homeostatic stem cell division pattern. The model is also capable of predicting changes in both tissue composition and in the progression of cells along their lineage at any given time and for various sequences of mutations. Simulations predict that the specific order in which mutations are acquired is crucial for determining the pace of cancer development. Interestingly, we find that the importance of genetic stability differs significantly depending on the physiological expression of mutations related to symmetric self-renewal and differentiation of stem and progenitor cells. In particular, mutations that lead to the alteration of the stem cell division pattern or the acquisition of some degree of immortality in committed progenitors lead to an early onset of cancer and diminish the impact of genetic instability.

Cancer as an evolutionary and ecological process

Neoplasms are microcosms of evolution. Within a neoplasm, a mosaic of mutant cells compete for space and resources, evade predation by the immune system and can even cooperate to disperse and colonize new organs. The evolution of neoplastic cells explains both why we get cancer and why it has been so difficult to cure. The tools of evolutionary biology and ecology are providing new insights into neoplastic progression and the clinical control of cancer.