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.