The evolution of cells within organisms can be an important determinant of disease. This is especially clear in the emergence of tumors and cancers from the underlying healthy tissue. In the healthy state, homeostasis is maintained through complex regulatory processes that ensure a relatively constant population size of cells, which is required for tissue function. Tumor cells escape this homeostasis, resulting in uncontrolled growth and consequent disease. Disease progression is driven by further evolutionary processes within the tumor, and so is the response of tumors to therapies. Therefore, evolutionary biology is an important component required for a better understanding of carcinogenesis and the treatment of cancers. In particular, evolutionary theory helps define the principles of mutant evolution and thus to obtain a clearer picture of the determinants of tumor emergence and therapy responses.     

The study by Raatz and Traulsen [1] makes an important contribution in this respect. They use mathematical and computational models to investigate trait evolution in the context of evolutionary rescue, motivated by the dynamics of cancer, and also bacterial infections. This study views the establishment of tumors as cell dynamics in harsh environments, where the population is prone to extinction unless mutants emerge that increase evolutionary fitness, allowing them to expand (evolutionary rescue). The core processes of the model include growth, death, and mutations. Random mutations are assumed to give rise to cell lineages with different trait combinations, where the birth and death rates of cells can change.  The resulting evolutionary trajectories are investigated in the models, and interesting new results were obtained. For example, the turnover of the population was identified as an important determinant of trait evolution. Turnover is defined as the balance between birth and death, with large rates corresponding to fast turnover and small rates to slow turnover. It was found that for fast cell turnover, a given adaptive step in the trait space results in a smaller increase in survival probability than for cell populations with slower turnover. In other words, evolutionary rescue is more difficult to achieve for fast compared to slow turnover populations. While more mutants can be produced for faster cell turnover rates, the analysis showed that this is not sufficient to overcome the barrier to the evolutionary rescue. This result implies that aggressive tumors with fast cell birth and death rates are less likely to persist and progress than tumors with lower turnover rates. This work emphasizes the importance of measuring the turnover rate in different tumors to advance our understanding of the determinants of tumor initiation and progression. The authors discuss that the well-documented heterogeneity in tumors likely also applies to cellular turnover. If a tumor consists of sub-populations with faster and slower turnover, it is possible that a slower turnover cell clone (e.g. characterized by a degree of dormancy) would enjoy a selective advantage. Another source of heterogeneity in turnover could be given by the hierarchical organization of tumors. Similar to the underlying healthy tissue, many tumors are thought to be maintained by a population of cancer stem cells, while the tumor bulk is made up of more differentiated cells. Tissue stem cells tend to be characterized by a lower turnover than progenitor or transit-amplifying cells. Depending on the assumptions about the self-renewal capacity of these different cell populations, the potential for evolutionary rescue could be different depending on the cell compartment in which the mutant emerges. This might be interesting to explore in the future.

There are also implications for treatment. Two types of treatment were investigated: density-affecting treatments in which the density of cells is reduced without altering their trait parameters, and trait-affecting treatments in which the birth and/or death rates are altered. Both types of treatment were found to change the trajectories of trait adaptation, which has potentially important practical implications. Interestingly, it was found that competitive release during treatment can result in situations where after treatment cessation, the non-extinct populations recover to reach sizes that were higher than in the absence of treatment. This points towards the potential of adaptive therapy approaches, where sensitive cells are maintained to some extent to suppress resistant clones [2] competitively. In this context, it is interesting that the success of such approaches might also depend on the turnover of the tumor cell population, as shown by a recent mathematical modeling study [3]. In particular, it was found that adaptive therapy is less likely to work for slow compared to fast turnover tumors. Yet, the current study by Raatz and Traulsen [1] suggests that tumors are more likely to evolve in a slow turnover setting.

While there is strong relevance of this analysis for tumor evolution, the results generated in this study have more general relevance. Besides tumors, the paper discusses applications to bacterial disease dynamics in some detail, which is also interesting to compare and contrast to evolutionary processes in cancer. Overall, this study provides insights into the dynamics of evolutionary rescue that represent valuable additions to evolutionary theory.  

References

[1] Raatz M, Traulsen A (2023) Promoting extinction or minimizing growth? The impact of treatment on trait trajectories in evolving populations. bioRxiv, 2022.06.17.496570, ver. 2 peer-reviewed and recommended by Peer Community in Evolutionary Biology. https://doi.org/10.1101/2022.06.17.496570

[2] Gatenby RA, Silva AS, Gillies RJ, Frieden BR (2009) Adaptive Therapy. Cancer Research, 69, 4894–4903. https://doi.org/10.1158/0008-5472.CAN-08-3658

[3] Strobl MAR, West J, Viossat Y, Damaghi M, Robertson-Tessi M, Brown JS, Gatenby RA, Maini PK, Anderson ARA (2021) Turnover Modulates the Need for a Cost of Resistance in Adaptive Therapy. Cancer Research, 81, 1135–1147. https://doi.org/10.1158/0008-5472.CAN-20-0806

Several examples of chromosomal inversions carrying genes affecting mate choice have been reported from various organisms. Furthermore, inversions are also frequently involved in genetic isolation between populations or species. Past work has shown that inversions can spread when they capture not only some loci involved in mate choice but also loci involved in incompatibilities between hybridizing populations [1]. In this new paper [2], the authors derive analytical approximations for the selection coefficient associated with an inversion suppressing recombination between a locus involved in mate choice and one (or several) locus involved in Dobzhansky-Muller incompatibilities. Two mechanisms for mate choice are considered: assortative mating based on the allele present at a single locus, or a trait-preference model where one locus codes for the trait and another for the preference. The results show that such an inversion is generally favoured, the selective advantage associated with the inversion being strongest when hybridization is sufficiently frequent. Assuming pairwise epistatic interactions between loci involved in incompatibilities, selection for the inversion increases approximately linearly with the number of such loci captured by the inversion.

This paper is a good read for several reasons. First, it presents the problem clearly (e.g. the introduction provides a clear and concise presentation of the issue and past work) and its crystal-clear writing facilitates the reader's understanding of theoretical approaches and results. Second, the analysis is competently done and adds to previous work by showing that very general conditions are expected to be favourable to the spread of the type of inversion considered here. And third, it provides food for thought about the role of inversions in the origin or the reinforcement of divergence between nascent species. One result of this work is that an inversion linked to pre-zygotic isolation "is favoured so long as there is viability selection against recombinant genotypes", suggesting that genetic incompatibilities must have evolved first and that inversions capturing mating preference loci may then enhance pre-existing reproductive isolation. However, the results also show that inversions are more likely to be favoured in hybridizing populations among which gene flow is still high, rather than in more strongly isolated populations. This matches the observation that inversions are more frequently observed between sympatric species than between allopatric ones.

References

[1] Trickett AJ, Butlin RK. 1994. Recombination Suppressors and the Evolution of New Species. Heredity 73:339-345. doi: 10.1038/hdy.1994.180

[2] Dagilis AJ, Kirkpatrick M. 2016. Prezygotic isolation, mating preferences, and the evolution of chromosomal inversions. Evolution 70: 1465–1472. doi: 10.1111/evo.12954

Evolution is usually studied via two distinct approaches: by inferring evolutionary processes from relatedness patterns among living species or by observing evolution in action in the laboratory or field. A recent study by Baym and colleagues in Science [1] has now combined these approaches by taking advantage of the pattern left behind by spatially evolving bacterial populations.

Evolution is often considered too slow to see, and can only be inferred by studying patterns of relatedness using phylogenetic trees. Increasingly, however, researchers are moving nature into the lab and watching as evolution unfolds under their noses. The field of experimental evolution follows evolutionary change in the laboratory over 10s to 1000s of generations, yielding insights into bacterial, viral, plant, or fly evolution (among many other species) that are simply not possible in the field. Yet, as powerful as experimental evolution is, it lacks a posterchild. There is no Galapagos finch radiation, nor a stunning series of cichlids to showcase to our students to pique their interests. Let’s face it, E. coli is no stickleback! And while practitioners of experimental evolution can explain the virtues of examining 60,000 generations of bacterial evolution in action, appreciating this nevertheless requires a level of insight and imagination that often eludes students, who need to see “it” to get it.

Enter MEGA, an idea and a film that could become the new face of experimental evolution. It replaces big numbers of generations or images of scientists, with an actual picture of the scientific result. MEGA, or Microbial Evolution and Growth Arena, is essentially an enormous petri dish and is the brainchild of Michael Baym, Tami Leiberman and their colleagues in Roy Kishony’s lab at Technion Israel Institute of Technology and Harvard Medical School. The idea of MEGA is to allow bacteria to swim over a spatially defined landscape while adapting to the local conditions, in this case antibiotics. When bacteria are inoculated onto one end of the plate they consume resources while swarming forward from the plate edge. In a few days, the bacteria grow into an area with antibiotics to which they are susceptible. This stops growth until a mutation arises that permits the bacteria to jump this hurdle, after which growth proceeds until the next hurdle of a 10-fold higher antibiotic concentration, and so on. By this simple approach, Baym et al. [1] evolved E. coli that were nearly 105-fold more resistant to two different antibiotics in just over 10 days. In addition, they identified the mutations that were required for these changes, showed that mutations conferring smaller benefits were required before bacteria could evolve maximal resistance, observed changes to the mutation rate, and demonstrated the importance of spatial structure in constraining adaptation.

For one thing, the rate of resistance evolution is impressive, and also quite scary given the mounting threat of antibiotic-resistant pathogens. However, MEGA also offers a uniquely visual insight into evolutionary change. By taking successive images of the MEGA plate, the group was able to watch the bacteria move, get trapped because of their susceptibility to the antibiotic, and then get past these traps as new mutations emerged that increased resistance. Each transition showcases evolution in real time. In addition, by leaving a spatial pattern of evolutionary steps behind, the MEGA plate offers unique opportunities to thoroughly investigate these steps when the experiment is finished. For instance, subsequent steps in mutational pathways can be characterized, but also their effects on fitness can be quantified in situ by measuring changes in survival and reproduction. This new method is undoubtedly a boon to the field of experimental evolution and offers endless opportunities for experimental elaboration. Perhaps of equal importance, MEGA is a tool that is portable to the classroom and to the public at large. Don’t believe in evolution? Watch this. You only have time for a short internship or lab practical? No problem. Don’t worry much about antibiotic resistance? Check this out. Like the best experimental tools, MEGA is simple but allows for complicated insights. And even if it is less charismatic than a finch, it still allows for the kinds of “gee-whiz” insights that will get students hooked on evolutionary biology.

Reference

[1] Baym M, Lieberman TD, Kelsic ED, Chait R, Gross R, Yelin I, Kishony R. 2016. Spatiotemporal microbial evolution on antibiotic landscapes. Science 353:1147-1151. doi: 10.1126/science.aag0822

In recent years many studies have documented shifts in phenology in response to climate change, be it in arrival times in migrating birds, budset in trees, adult emergence in butterflies, or flowering time in annual plants (Coen et al. 2018; Piao et al. 2019). While these changes are, in part, explained by phenotypic plasticity, more and more studies find that they involve also genetic changes, that is, they involve evolutionary change (e.g., Metz et al. 2020). Yet, evolutionary change may occur through genetic drift as well as selection. Therefore, in order to demonstrate adaptive evolutionary change in response to climate change, drift has to be excluded as an alternative explanation (Hansen et al. 2012). A new study by Gay et al. (2021) shows just how difficult this can be. 

The authors investigated a recent evolutionary shift in flowering time by in a population an annual plant that reproduces predominantly by self-fertilization. The population has recently been subjected to increased temperatures and reduced rainfalls both of which are believed to select for earlier flowering times. They used a “resurrection” approach (Orsini et al. 2013; Weider et al. 2018): Genotypes from the past (resurrected from seeds) were compared alongside more recent genotypes (from more recently collected seeds) under identical conditions in the greenhouse. Using an experimental design that replicated genotypes, eliminated maternal effects, and controlled for microenvironmental variation, they found said genetic change in flowering times: Genotypes obtained from recently collected seeds flowered significantly (about 2 days) earlier than those obtained 22 generations before. However, neutral markers (microsatellites) also showed strong changes in allele frequencies across the 22 generations, suggesting that effective population size, Ne, was low (i.e., genetic drift was strong), which is typical for highly self-fertilizing populations. In addition, several multilocus genotypes were present at high frequencies and persisted over the 22 generations, almost as in clonal populations (e.g., Schaffner et al. 2019). The challenge was thus to evaluate whether the observed evolutionary change was the result of an adaptive response to selection or may be explained by drift alone. 

Here, Gay et al. (2021) took a particularly careful and thorough approach. First, they carried out a selection gradient analysis, finding that earlier-flowering plants produced more seeds than later-flowering plants. This suggests that, under greenhouse conditions, there was indeed selection for earlier flowering times. Second, investigating other populations from the same region (all populations are located on the Mediterranean island of Corsica, France), they found that a concurrent shift to earlier flowering times occurred also in these populations. Under the hypothesis that the populations can be regarded as independent replicates of the evolutionary process, the observation of concurrent shifts rules out genetic drift (under drift, the direction of change is expected to be random). 

The study may well have stopped here, concluding that there is good evidence for an adaptive response to selection for earlier flowering times in these self-fertilizing plants, at least under the hypothesis that selection gradients estimated in the greenhouse are relevant to field conditions. However, the authors went one step further. They used the change in the frequencies of the multilocus genotypes across the 22 generations as an estimate of realized fitness in the field and compared them to the phenotypic assays from the greenhouse. The results showed a tendency for high-fitness genotypes (positive frequency changes) to flower earlier and to produce more seeds than low-fitness genotypes. However, a simulation model showed that the observed correlations could be explained by drift alone, as long as Ne is lower than ca. 150 individuals. The findings were thus consistent with an adaptive evolutionary change in response to selection, but drift could only be excluded as the sole explanation if the effective population size was large enough. 

The study did provide two estimates of Ne (19 and 136 individuals, based on individual microsatellite loci or multilocus genotypes, respectively), but both are problematic. First, frequency changes over time may be influenced by the presence of a seed bank or by immigration from a genetically dissimilar population, which may lead to an underestimation of Ne (Wang and Whitlock 2003). Indeed, the low effective size inferred from the allele frequency changes at microsatellite loci appears to be inconsistent with levels of genetic diversity found in the population. Moreover, high self-fertilization reduces effective recombination and therefore leads to non-independence among loci. This lowers the precision of the Ne estimates (due to a higher sampling variance) and may also violate the assumption of neutrality due to the possibility of selection (e.g., due to inbreeding depression) at linked loci, which may be anywhere in the genome in case of high degrees of self-fertilization. 

There is thus no definite answer to the question of whether or not the observed changes in flowering time in this population were driven by selection. The study sets high standards for other, similar ones, in terms of thoroughness of the analyses and care in interpreting the findings. It also serves as a very instructive reminder to carefully check the assumptions when estimating neutral expectations, especially when working on species with complicated demographies or non-standard life cycles. Indeed the issues encountered here, in particular the difficulty of establishing neutral expectations in species with low effective recombination, may apply to many other species, including partially or fully asexual ones (Hartfield 2016). Furthermore, they may not be limited to estimating Ne but may also apply, for instance, to the establishment of neutral baselines for outlier analyses in genome scans (see e.g, Orsini et al. 2012). 

References

Cohen JM, Lajeunesse MJ, Rohr JR (2018) A global synthesis of animal phenological responses to climate change. Nature Climate Change, 8, 224–228. https://doi.org/10.1038/s41558-018-0067-3

Gay L, Dhinaut J, Jullien M, Vitalis R, Navascués M, Ranwez V, Ronfort J (2021) Evolution of flowering time in a selfing annual plant: Roles of adaptation and genetic drift. bioRxiv, 2020.08.21.261230, ver. 4 recommended and peer-reviewed by Peer Community in Evolutionary Biology. https://doi.org/10.1101/2020.08.21.261230

Hansen MM, Olivieri I, Waller DM, Nielsen EE (2012) Monitoring adaptive genetic responses to environmental change. Molecular Ecology, 21, 1311–1329. https://doi.org/10.1111/j.1365-294X.2011.05463.x

Hartfield M (2016) Evolutionary genetic consequences of facultative sex and outcrossing. Journal of Evolutionary Biology, 29, 5–22. https://doi.org/10.1111/jeb.12770

Metz J, Lampei C, Bäumler L, Bocherens H, Dittberner H, Henneberg L, Meaux J de, Tielbörger K (2020) Rapid adaptive evolution to drought in a subset of plant traits in a large-scale climate change experiment. Ecology Letters, 23, 1643–1653. https://doi.org/10.1111/ele.13596

Orsini L, Schwenk K, De Meester L, Colbourne JK, Pfrender ME, Weider LJ (2013) The evolutionary time machine: using dormant propagules to forecast how populations can adapt to changing environments. Trends in Ecology & Evolution, 28, 274–282. https://doi.org/10.1016/j.tree.2013.01.009

Orsini L, Spanier KI, Meester LD (2012) Genomic signature of natural and anthropogenic stress in wild populations of the waterflea Daphnia magna: validation in space, time and experimental evolution. Molecular Ecology, 21, 2160–2175. https://doi.org/10.1111/j.1365-294X.2011.05429.x

Piao S, Liu Q, Chen A, Janssens IA, Fu Y, Dai J, Liu L, Lian X, Shen M, Zhu X (2019) Plant phenology and global climate change: Current progresses and challenges. Global Change Biology, 25, 1922–1940. https://doi.org/10.1111/gcb.14619

Schaffner LR, Govaert L, De Meester L, Ellner SP, Fairchild E, Miner BE, Rudstam LG, Spaak P, Hairston NG (2019) Consumer-resource dynamics is an eco-evolutionary process in a natural plankton community. Nature Ecology & Evolution, 3, 1351–1358. https://doi.org/10.1038/s41559-019-0960-9

Wang J, Whitlock MC (2003) Estimating Effective Population Size and Migration Rates From Genetic Samples Over Space and Time. Genetics, 163, 429–446. PMID: 12586728

Weider LJ, Jeyasingh PD, Frisch D (2018) Evolutionary aspects of resurrection ecology: Progress, scope, and applications—An overview. Evolutionary Applications, 11, 3–10. https://doi.org/10.1111/eva.12563

Differentiation-based genome scans lie at the core of speciation and adaptation genomics research. Dating back to Lewontin & Krakauer (1973), they have become very popular with the advent of genomics to identify genome regions of enhanced differentiation relative to neutral expectations. These regions may represent genetic barriers between divergent lineages and are key for studying reproductive isolation. However, genome scan methods can generate a high rate of false positives, primarily if the neutral population structure is not accounted for (Bierne et al. 2013). Moreover, interpreting genome scans can be challenging in the context of secondary contacts between diverging lineages (Bierne et al. 2011), because the coupling between different components of reproductive isolation (local adaptation, intrinsic incompatibilities, mating preferences, etc.) can occur readily, thus preventing the causes of differentiation from being determined.

Smadja and collaborators (2022) applied a sophisticated genome scan for trait association (BAYPASS, Gautier 2015) to underlie the genetic basis of a polygenetic behaviour: assortative mating in hybridizing mice. My interest in this neat study mainly relies on two reasons. First, the authors used an ingenious geographical setting (replicate pairs of “Choosy” versus “Non-Choosy” populations) with multi-way comparisons to narrow down the list of candidate regions resulting from BAYPASS. The latter corrects for population structure, handles cost-effective pool-seq data and allows for gene-based analyses that aggregate SNP signals within a gene. These features reinforce the set of outlier genes detected; however, not all are expected to be associated with mating preference. 

The second reason why this study is valuable to me is that Smadja et al. (2022) complemented the population genomic approach with functional predictions to validate the genetic signal. In line with previous behavioural and chemical assays on the proximal mechanisms of mating preferences, they identified multiple olfactory and vomeronasal receptor genes as highly significant candidates. Therefore, combining genomic signals with functional analyses is a clever way to provide insights into the causes of reproductive isolation, especially when multiple barriers are involved. This is typically true for reinforcement (Butlin & Smadja 2018), suspected to occur in these mice because, in that case, assortative mating (a prezygotic barrier) evolves in response to the cost of hybridization (for example, due to hybrid inviability). 

As advocated by the authors, their study paves the way for future work addressing the genetic basis of reinforcement, a trait of major evolutionary importance for which we lack empirical data. They also make a compelling case using complementary approaches that olfactory and vomeronasal receptors have a central role in mammal speciation.

References:

Bierne N, Welch J, Loire E, Bonhomme F, David P (2011) The coupling hypothesis: why genome scans may fail to map local adaptation genes. Mol Ecol 20: 2044–2072. https://doi.org/10.1111/j.1365-294X.2011.05080.x

Bierne N, Roze D, Welch JJ (2013) Pervasive selection or is it…? why are FST outliers sometimes so frequent? Mol Ecol 22: 2061–2064. https://doi.org/10.1111/mec.12241 

Butlin RK, Smadja CM (2018) Coupling, Reinforcement, and Speciation. Am Nat 191:155–172. https://doi.org/10.1086/695136 

Gautier M (2015) Genome-Wide Scan for Adaptive Divergence and Association with Population-Specific Covariates. Genetics 201:1555–1579. https://doi.org/10.1534/genetics.115.181453 

Lewontin RC, Krakauer J (1973) Distribution of gene frequency as a test of the theory of selective neutrality of polymorphisms. Genetics 74: 175–195. https://doi.org/10.1093/genetics/74.1.175 

Smadja CM, Loire E, Caminade P, Severac D, Gautier M, Ganem G (2022) Divergence of olfactory receptors associated with the evolution of assortative mating and reproductive isolation in mice. bioRxiv, 2022.07.21.500634, ver. 3 peer-reviewed and recommended by Peer Community in Evolutionary Biology. https://doi.org/10.1101/2022.07.21.500634