The evolutionary trade-off between stem cell niche size, aging, and tumorigenesis

Mutations, substitutions, and selection: Linking mutagenic processes to cancer using evolutionary theory

Cancers are the product of evolutionary events, where molecular variation occurs and accumulates in tissues and tumors. Sequencing of this molecular variation informs not only which variants are driving tumorigenesis, but also the mechanisms behind what is fueling mutagenesis. Both of these details are crucial for preventing premature deaths due to cancer, whether it is by targeting the variants driving the cancer phenotype or by measures to prevent exogenous mutations from contributing to somatic evolution. Here, we review tools to determine both molecular signatures and cancer drivers, and avenues by which these metrics may be linked.

Adult stem cell activity in naked mole rats for long-term tissue maintenance

The naked mole rat (NMR), Heterocephalus glaber, the longest-living rodent, provides a unique opportunity to explore how evolution has shaped adult stem cell (ASC) activity and tissue function with increasing lifespan. Using cumulative BrdU labelling and a quantitative imaging approach to track intestinal ASCs (Lgr5+) in their native in vivo state, we find an expanded pool of Lgr5+ cells in NMRs, and these cells specifically at the crypt base (Lgr5+CBC) exhibit slower division rates compared to those in short-lived mice but have a similar turnover as human LGR5+CBC cells. Instead of entering quiescence (G0), NMR Lgr5+CBC cells reduce their division rates by prolonging arrest in the G1 and/or G2 phases of the cell cycle. Moreover, we also observe a higher proportion of differentiated cells in NMRs that confer enhanced protection and function to the intestinal mucosa which is able to detect any chemical imbalance in the luminal environment efficiently, triggering a robust pro-apoptotic, anti-proliferative response within the stem/progenitor cell zone.

Minimizing cell number fluctuations in self-renewing tissues with a stem-cell niche

Self-renewing tissues require that a constant number of proliferating cells is maintained over time. This maintenance can be ensured at the single-cell level or the population level. Maintenance at the population level leads to fluctuations in the number of proliferating cells over time. Often, it is assumed that those fluctuations can be reduced by increasing the number of asymmetric divisions, i.e., divisions where only one of the daughter cells remains proliferative. Here, we study a model of cell proliferation that incorporates a stem-cell niche of fixed size, and explicitly model the cells inside and outside the niche. We find that in this model, fluctuations are minimized when the difference in growth rate between the niche and the rest of the tissue is maximized and all divisions are symmetric divisions, producing either two proliferating or two nonproliferating daughters. We show that this optimal state leaves visible signatures in clone size distributions and could thus be detected experimentally.

Reuniting philosophy and science to advance cancer research

Cancers rely on multiple, heterogeneous processes at different scales, pertaining to many biomedical fields. Therefore, understanding cancer is necessarily an interdisciplinary task that requires placing specialised experimental and clinical research into a broader conceptual, theoretical, and methodological framework. Without such a framework, oncology will collect piecemeal results, with scant dialogue between the different scientific communities studying cancer. We argue that one important way forward in service of a more successful dialogue is through greater integration of applied sciences (experimental and clinical) with conceptual and theoretical approaches, informed by philosophical methods. By way of illustration, we explore six central themes: (i) the role of mutations in cancer; (ii) the clonal evolution of cancer cells; (iii) the relationship between cancer and multicellularity; (iv) the tumour microenvironment; (v) the immune system; and (vi) stem cells. In each case, we examine open questions in the scientific literature through a philosophical methodology and show the benefit of such a synergy for the scientific and medical understanding of cancer.

Are genetic drift and stem cell adherence in laboratory culture issues for cultivated meat production?

Mesenchymal stem cell-based cultivated meat is a promising solution to the ecological and ethical problems posed by traditional meat production, since it exhibits a protein content and composition that is more comparable to original meat proteins than any other source of cultivated meat products, including plants, bacteria, and fungi. Nonetheless, the nature and laboratory behavior of mesenchymal stem cells pose two significant challenges for large-scale production: genetic drift and adherent growth in culture. Culture conditions used in the laboratory expose the cells to a selective pressure that causes genetic drift, which may give rise to oncogene activation and the loss of "stemness." This is why genetic and functional analysis of the cells during culture is required to determine the maximum number of passages within the laboratory where no significant mutations or loss of function are detected. Moreover, the adherent growth of mesenchymal stem cells can be an obstacle for their large-scale production since volume to surface ratio is limited for high volume containers. Multi-tray systems, roller bottles, and microcarriers have been proposed as potential solutions to scale-up the production of adherent cells required for cultivated meat. The most promising solutions for the safety problems and large-scale obstacles for cultivated meat production are the determination of a limit number of passages based on a genetic analysis and the use of microcarriers from edible materials to maximize the volume to surface proportion and decrease the downstream operations needed for cultivated meat production.

Spectral dynamics of guided edge removals and identifying transient amplifiers for death-Birth updating

The paper deals with two interrelated topics: (1) identifying transient amplifiers in an iterative process, and (2) analyzing the process by its spectral dynamics, which is the change in the graph spectra by edge manipulation. Transient amplifiers are networks representing population structures which shift the balance between natural selection and random drift. Thus, amplifiers are highly relevant for understanding the relationships between spatial structures and evolutionary dynamics. We study an iterative procedure to identify transient amplifiers for death-Birth updating. The algorithm starts with a regular input graph and iteratively removes edges until desired structures are achieved. Thus, a sequence of candidate graphs is obtained. The edge removals are guided by quantities derived from the sequence of candidate graphs. Moreover, we are interested in the Laplacian spectra of the candidate graphs and analyze the iterative process by its spectral dynamics. The results show that although transient amplifiers for death-Birth updating are generally rare, a substantial number of them can be obtained by the proposed procedure. The graphs identified share structural properties and have some similarity to dumbbell and barbell graphs. We analyze amplification properties of these graphs and also two more families of bell-like graphs and show that further transient amplifiers for death-Birth updating can be found. Finally, it is demonstrated that the spectral dynamics possesses characteristic features useful for deducing links between structural and spectral properties. These feature can also be taken for distinguishing transient amplifiers among evolutionary graphs in general.

BRAVO self-confined expression through WOX5 in the Arabidopsis root stem-cell niche

In animals and plants, stem-cell niches are local microenvironments that are tightly regulated to preserve their unique identity while communicating with adjacent cells that will give rise to specialized cell types. In the primary root of Arabidopsis thaliana, two transcription factors, BRAVO and WOX5, among others, are expressed in the stem-cell niche. Intriguingly, BRAVO, a repressor of quiescent center divisions, confines its own gene expression to the stem-cell niche, as evidenced in a bravo mutant background. Here, we propose through mathematical modeling that BRAVO confines its own expression domain to the stem-cell niche by attenuating a WOX5-dependent diffusible activator of BRAVO. This negative feedback drives WOX5 activity to be spatially restricted as well. The results show that WOX5 diffusion and sequestration by binding to BRAVO are sufficient to drive the experimentally observed confined BRAVO expression at the stem-cell niche. We propose that the attenuation of a diffusible activator can be a general mechanism acting at other stem-cell niches to spatially confine genetic activity to a small region while maintaining signaling within them and with the surrounding cells.

Cancer evolution: Darwin and beyond

Clinical and laboratory studies over recent decades have established branched evolution as a feature of cancer. However, while grounded in somatic selection, several lines of evidence suggest a Darwinian model alone is insufficient to fully explain cancer evolution. First, the role of macroevolutionary events in tumour initiation and progression contradicts Darwin's central thesis of gradualism. Whole-genome doubling, chromosomal chromoplexy and chromothripsis represent examples of single catastrophic events which can drive tumour evolution. Second, neutral evolution can play a role in some tumours, indicating that selection is not always driving evolution. Third, increasing appreciation of the role of the ageing soma has led to recent generalised theories of age-dependent carcinogenesis. Here, we review these concepts and others, which collectively argue for a model of cancer evolution which extends beyond Darwin. We also highlight clinical opportunities which can be grasped through targeting cancer vulnerabilities arising from non-Darwinian patterns of evolution.

Spectral analysis of transient amplifiers for death-birth updating constructed from regular graphs

A central question of evolutionary dynamics on graphs is whether or not a mutation introduced in a population of residents survives and eventually even spreads to the whole population, or becomes extinct. The outcome naturally depends on the fitness of the mutant and the rules by which mutants and residents may propagate on the network, but arguably the most determining factor is the network structure. Some structured networks are transient amplifiers. They increase for a certain fitness range the fixation probability of beneficial mutations as compared to a well-mixed population. We study a perturbation method for identifying transient amplifiers for death–birth updating. The method involves calculating the coalescence times of random walks on graphs and finding the vertex with the largest remeeting time. If the graph is perturbed by removing an edge from this vertex, there is a certain likelihood that the resulting perturbed graph is a transient amplifier. We test all pairwise nonisomorphic regular graphs up to a certain order and thus cover the whole structural range expressible by these graphs. For cubic and quartic regular graphs we find a sufficiently large number of transient amplifiers. For these networks we carry out a spectral analysis and show that the graphs from which transient amplifiers can be constructed share certain structural properties. Identifying spectral and structural properties may promote finding and designing such networks.

Fixation probabilities in graph-structured populations under weak selection

A population's spatial structure affects the rate of genetic change and the outcome of natural selection. These effects can be modeled mathematically using the Birth-death process on graphs. Individuals occupy the vertices of a weighted graph, and reproduce into neighboring vertices based on fitness. A key quantity is the probability that a mutant type will sweep to fixation, as a function of the mutant's fitness. Graphs that increase the fixation probability of beneficial mutations, and decrease that of deleterious mutations, are said to amplify selection. However, fixation probabilities are difficult to compute for an arbitrary graph. Here we derive an expression for the fixation probability, of a weakly-selected mutation, in terms of the time for two lineages to coalesce. This expression enables weak-selection fixation probabilities to be computed, for an arbitrary weighted graph, in polynomial time. Applying this method, we explore the range of possible effects of graph structure on natural selection, genetic drift, and the balance between the two. Using exhaustive analysis of small graphs and a genetic search algorithm, we identify families of graphs with striking effects on fixation probability, and we analyze these families mathematically. Our work reveals the nuanced effects of graph structure on natural selection and neutral drift. In particular, we show how these notions depend critically on the process by which mutations arise.

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.

Tumorigenesis as the Paradigm of Quasi-neutral Molecular Evolution

In the absence of both positive and negative selections, coding sequences evolve at a neutral rate (R = 1). Such a high genomic rate is generally not achievable due to the prevalence of negative selection against codon substitutions. Remarkably, somatic evolution exhibits the seemingly neutral rate R ∼ 1 across normal and cancerous tissues. Nevertheless, R ∼ 1 may also mean that positive and negative selections are both strong, but equal in intensity. We refer to this regime as quasi-neutral. Indeed, individual genes in cancer cells often evolve at a much higher, or lower, rate than R ∼ 1. Here, we show that 1) quasi-neutrality is much more likely when populations are small (N < 50); 2) stem-cell populations in single normal tissue niches, from which tumors likely emerge, have a small N (usually <50) but selection at this stage is measurable and strong; 3) when N dips below 50, selection efficacy decreases precipitously; and 4) notably, N is smaller in the stem-cell niche of the small intestine than in the colon. Hence, the ∼70-fold higher rate of phenotypic evolution (observed as cancer risk) in the latter can be explained by the greater efficacy of selection, which then leads to the fixation of more advantageous and fewer deleterious mutations in colon cancers. In conclusion, quasi-neutral evolution sheds a new light on a general evolutionary principle that helps to explain aspects of cancer evolution.

Neutral Theory and the Somatic Evolution of Cancer

Kimura's neutral theory argued that positive selection was not responsible for an appreciable fraction of molecular substitutions. Correspondingly, quantitative analysis reveals that the vast majority of substitutions in cancer genomes are not detectably under selection. Insights from the somatic evolution of cancer reveal that beneficial substitutions in cancer constitute a small but important fraction of the molecular variants. The molecular evolution of cancer community will benefit by incorporating the neutral theory of molecular evolution into their understanding and analysis of cancer evolution-and accepting the use of tractable, predictive models, even when there is some evidence that they are not perfect.

The evolutionary trade-off between stem cell niche size, aging, and tumorigenesis

Many epithelial tissues within multicellular organisms are continually replenished by small independent populations of stem cells largely responsible for maintaining tissue homeostasis. These continually dividing populations are subject to mutations that can lead to tumorigenesis but also contribute to aging. Mutations accumulate in stem cell niches and change the rate of cell division and differentiation; the pace of this process and the fate of specific mutations depend strongly on niche population size. Here, we create a mathematical model of the intestinal stem cell niche, crypt system, and epithelium. We calculate the expected effect of fixed mutations in stem cell niches and their effect on tissue homeostasis throughout the intestinal epithelium over organismal lifetime. We find that, due to the small population size of stem cell niches, mutations predominantly fix via genetic drift and decrease stem cell fitness, leading to niche and tissue attrition, and contributing to organismal aging. We also explore mutation accumulation at various stem cell niche sizes and demonstrate that an evolutionary trade-off exists between niche size, tissue aging, and the risk of tumorigenesis. Further, mouse and human niches exist at a size that minimizes the probability of tumorigenesis, at the expense of accumulating deleterious mutations due to genetic drift. Finally, we show that the trade-off between the probability of tumorigenesis and the extent of aging depends on whether or not mutational effects confer a selective advantage in the stem cell niche.