The article “Drivers of genomic landscapes of differentiation across Populus divergence gradient” by Shang et al. describes an amazing dataset where genomic variations among 21 pairs of diverging poplar species are compared. Such comparisons are still quite rare and are needed to shed light on the processes shaping genomic divergence along the speciation gradient. Relying on two hundred whole-genome resequenced samples from 8 species that diverged from 1.3 to 4.8 million years ago, the authors aim at identifying the key factors involved in the genomic differentiation between species. They carried out a wide range of robust statistical tests aiming at characterizing the genomic differentiation along the genome of these species pairs. They highlight in particular the role of linked selection and gene flow in shaping the divergence along the genomes of species pairs. They also confirm the significance of introgression among species with a net divergence larger than the upper boundaries of the grey zone of speciation previously documented in animals (da from 0.005 to 0.02, Roux et al. 2016). Because these findings pave the way to research about the genomic mechanisms associated with speciation in species with allopatric and parapatric distributions, I warmingly recommend this article.

References

Roux C, Fraïsse C, Romiguier J, Anciaux Y, Galtier N, Bierne N (2016) Shedding Light on the Grey Zone of Speciation along a Continuum of Genomic Divergence. PLOS Biology, 14, e2000234. https://doi.org/10.1371/journal.pbio.2000234

Shang H, Rendón-Anaya M, Paun O, Field DL, Hess J, Vogl C, Liu J, Ingvarsson PK, Lexer C, Leroy T (2023) Drivers of genomic landscapes of differentiation across Populus divergence gradient. bioRxiv, 2021.08.26.457771, ver. 5 peer-reviewed and recommended by Peer Community in Evolutionary Biology. https://doi.org/10.1101/2021.08.26.457771

Evolutionary genomics, and population genetics in particular, aim to decipher the respective influence of neutral and selective forces shaping genetic polymorphism in a species/population. This is a much-needed requirement before scanning genome data for footprints of species adaptation to their biotic and abiotic environment (Johri et al. 2022). In general, we would like to quantify the proportion of the genome evolving neutrally and under selective (positive, balancing and negative) pressures (Kern and Hahn 2018, Johri et al. 2021). We thus need to understand patterns of linked selection along the genome, that is how the distribution of genetic polymorphisms is shaped by selected sites and the recombination landscape. The present contribution by Pavinato et al. (2022) provides an additional method in the population genomics toolbox to quantify the extent of linked positive and negative selection using temporal data.

The availability of genomics data for model and non-model species has led to improvement of the modeling framework for demography and selection (Johri et al. 2022), but also new inference methods making use of the full genome data based on the Sequential Markovian Coalescent (SMC, Li and Durbin 2011), Approximate Bayesian Computation (ABC, Jay et al. 2019), ABC and machine learning (Pudlo et al. 2016, Raynal et al. 2019) or Deep Learning (Sanchez et al. 2021). These methods are based on one sample in time and the use of the coalescent theory to reconstruct the past (demographic) history. However, it is also possible to obtain for many species temporal data sampled over several time points. For species with short generation time (in experimental evolution or monitored populations), one can sample a population every couple of generations as exemplified with Drosophila melanogaster (Bergland et al. 2010). For species with longer generation times that cannot be easily regularly sampled in time, it becomes possible to sequence available specimens from museums (e.g. Cridland et al. 2018) or ancient DNA samples. Methods using temporal data are based on the classical population genomics assumption that demography (migration, population subdivision, population size changes) leaves a genome-wide signal, while selection leaves a localized signal in the close vicinity of the causal mutation. Several methods do assess the demography of a population (change in effective population size, Ne, in time) using temporal data (e.g. Jorde and Ryman 2007) which can be used to calibrate the detection of loci under strong positive selection (Foll et al. 2014). Recently Buffalo and Coop (2020) used genome-wide covariance between allele frequency changes across time samples (and across replicates) to quantify the effects of linked selection over short timescales. 

In the present contribution, Pavinato et al. (2022) make use of temporal data to draw the joint estimation of demographic and selective parameters using a simulation-based method (ABC-Random Forests). This study by Pavinato et al. (2022) builds a framework allowing to infer the census size of the population in time (N) separately from the effect of genetic drift, which is determined by change in effective population size (Ne) in time, as well estimates of genome-wide parameters of selection. In a nutshell, the authors use a forward simulator and summarize genome data by genomic windows using classic statistics (nucleotide diversity, Tajima’s D, FST, heterozygosity) between time samples and for each sample. They specifically use the distributions (higher moments) of these statistics among all windows. The authors combine as input for the ABC-RF, vectors of summary statistics, model parameters and five latent variables: Ne, the ratio Ne/N, the number of beneficial mutations under strong selection, the average selection coefficient of strongly selected mutations, and the average substitution load. Indeed, the authors are interested in three different types of selection components: 1) the adaptive potential of a population which is estimated as the population mutation rate of beneficial mutations (θb), 2) the number of mutations under strong selection (irrespective of whether they reached fixation or not), and 3) the overall population fitness which is a function of the genetic load. In other words, the novelty of this method is not to focus on the detection of loci under selection, but to infer key parameters/distributions summarizing the genome-wide signal of demography and (positive and negative) selection. As a proof of principle, the authors then apply their method to a dataset of feral populations of honey bees (Apis mellifera) collected in California across many years and recovered from Museum samples (Cridland et al. 2018). The approach yields estimates of Ne which are on the same order of magnitude of previous estimates in hymenopterans, and the authors discuss why the different populations show various values of Ne and N which can be explained by different history of admixture with wild but also domesticated lineages of bees.

This study focuses on quantifying the genome-wide joint footprints of demography, and strong positive and negative selection to determine which proportion of the genome evolves neutrally or not. Further application of this method can be anticipated, for example, to study species with ecological and life-history traits which generate discrepancies between census size and Ne, for example for plants with selfing or seed banking (Sellinger et al. 2020), and for which the genome-wide effect of linked selection is not fully understood.

References

Johri P, Aquadro CF, Beaumont M, Charlesworth B, Excoffier L, Eyre-Walker A, Keightley PD, Lynch M, McVean G, Payseur BA, Pfeifer SP, Stephan W, Jensen JD (2022) Recommendations for improving statistical inference in population genomics. PLOS Biology, 20, e3001669. https://doi.org/10.1371/journal.pbio.3001669

Kern AD, Hahn MW (2018) The Neutral Theory in Light of Natural Selection. Molecular Biology and Evolution, 35, 1366–1371. https://doi.org/10.1093/molbev/msy092

Johri P, Riall K, Becher H, Excoffier L, Charlesworth B, Jensen JD (2021) The Impact of Purifying and Background Selection on the Inference of Population History: Problems and Prospects. Molecular Biology and Evolution, 38, 2986–3003. https://doi.org/10.1093/molbev/msab050

Pavinato VAC, Mita SD, Marin J-M, Navascués M de (2022) Joint inference of adaptive and demographic history from temporal population genomic data. bioRxiv, 2021.03.12.435133, ver. 6 peer-reviewed and recommended by Peer Community in Evolutionary Biology. https://doi.org/10.1101/2021.03.12.435133

Li H, Durbin R (2011) Inference of human population history from individual whole-genome sequences. Nature, 475, 493–496. https://doi.org/10.1038/nature10231

Jay F, Boitard S, Austerlitz F (2019) An ABC Method for Whole-Genome Sequence Data: Inferring Paleolithic and Neolithic Human Expansions. Molecular Biology and Evolution, 36, 1565–1579. https://doi.org/10.1093/molbev/msz038

Pudlo P, Marin J-M, Estoup A, Cornuet J-M, Gautier M, Robert CP (2016) Reliable ABC model choice via random forests. Bioinformatics, 32, 859–866. https://doi.org/10.1093/bioinformatics/btv684

Raynal L, Marin J-M, Pudlo P, Ribatet M, Robert CP, Estoup A (2019) ABC random forests for Bayesian parameter inference. Bioinformatics, 35, 1720–1728. https://doi.org/10.1093/bioinformatics/bty867

Sanchez T, Cury J, Charpiat G, Jay F (2021) Deep learning for population size history inference: Design, comparison and combination with approximate Bayesian computation. Molecular Ecology Resources, 21, 2645–2660. https://doi.org/10.1111/1755-0998.13224

Bergland AO, Behrman EL, O’Brien KR, Schmidt PS, Petrov DA (2014) Genomic Evidence of Rapid and Stable Adaptive Oscillations over Seasonal Time Scales in Drosophila. PLOS Genetics, 10, e1004775. https://doi.org/10.1371/journal.pgen.1004775

Cridland JM, Ramirez SR, Dean CA, Sciligo A, Tsutsui ND (2018) Genome Sequencing of Museum Specimens Reveals Rapid Changes in the Genetic Composition of Honey Bees in California. Genome Biology and Evolution, 10, 458–472. https://doi.org/10.1093/gbe/evy007

Jorde PE, Ryman N (2007) Unbiased Estimator for Genetic Drift and Effective Population Size. Genetics, 177, 927–935. https://doi.org/10.1534/genetics.107.075481

Foll M, Shim H, Jensen JD (2015) WFABC: a Wright–Fisher ABC-based approach for inferring effective population sizes and selection coefficients from time-sampled data. Molecular Ecology Resources, 15, 87–98. https://doi.org/10.1111/1755-0998.12280

Buffalo V, Coop G (2020) Estimating the genome-wide contribution of selection to temporal allele frequency change. Proceedings of the National Academy of Sciences, 117, 20672–20680. https://doi.org/10.1073/pnas.1919039117

Sellinger TPP, Awad DA, Moest M, Tellier A (2020) Inference of past demography, dormancy and self-fertilization rates from whole genome sequence data. PLOS Genetics, 16, e1008698. https://doi.org/10.1371/journal.pgen.1008698

The unique challenges associated with living on an island often result in organisms displaying a specific suite of traits commonly referred to as “island syndrome” (Adler and Levins, 1994; Burns, 2019; Baeckens and Van Damme, 2020). Large phenotypic shifts such as changes in size (e.g. shifts to gigantism or dwarfism, Lomolino, 2005) or coloration (Doutrelant et al., 2016) abound in the literature. However, less obvious phenotypes may also play a key role in adaptation to islands.

One such trait, reduced immune function, has important implications for the future of island populations in the face of anthropogenic-induced changes. Due to lower parasite pressure caused by a less diverse and less virulent parasite population, island hosts may show a decrease in immune defenses (Beadell et al., 2006; Pérez‐Rodríguez et al., 2013). However, this hypothesis has been challenged, as many studies have found ambiguous or conflicting results (Matson, 2006; Illera et al., 2015).

While most previous work has examined various immunological parameters (e.g., antibody concentrations), here, Barthe et al. (2022) take the novel approach of examining molecular signatures of immune genes. Using comparative genomic data from 34 different species of birds the authors examine the ratio of synonymous substitutions (i.e., not changing an amino acid) to non-synonymous substitutions (i.e., changing an amino acid) in innate and acquired immune genes (Pn/Ps ratio). Because population sizes on islands are lower which will affect molecular evolution, they compare these results to data from 97 control genes.  Assuming relaxed selection on islands predicts that the difference between the Pn/Ps ratio of immune genes and of control genes (ΔPn/Ps) is greater in island species compared to mainland ones.

As with previous work the authors found that the results differ depending on the category of immune genes. Both forms of innate defense: beta-defensins and Toll-like receptors did not show higher ΔPn/Ps for island populations. As these genes still have a higher Pn/Ps than control genes, the authors argue these results are in line with these genes being under purifying selection but lacking an “island effect”. Instead, the authors argue that demographic effects (i.e., reductions in Ne) may lead to the decreased immunity documented in other studies. In contrast, there was a reduction in Pn/Ps in MHC II genes, known to be under balancing selection. This reduction was stronger in island species and thus the authors argue that this is the only class of genes where a role for relaxed selection can be invoked. 

Together these results demonstrate that the changes in immunity experienced by island species are complex and that different categories of immune genes can experience different selective pressures. By including control genes in their study, they particularly highlight the importance of accounting for shifts in Ne when examining patterns of island species evolution. Hopefully, this kind of framework will be applied to other taxa to determine if these results are widespread or more specific to birds. 

References

Adler GH, Levins R (1994) The Island Syndrome in Rodent Populations. The Quarterly Review of Biology, 69, 473–490. https://doi.org/10.1086/418744

Baeckens S, Van Damme R (2020) The island syndrome. Current Biology, 30, R338–R339. https://doi.org/10.1016/j.cub.2020.03.029

Barthe M, Doutrelant C, Covas R, Melo M, Illera JC, Tilak M-K, Colombier C, Leroy T, Loiseau C, Nabholz B (2022) Evolution of immune genes in island birds: reduction in population sizes can explain island syndrome. bioRxiv, 2021.11.21.469450, ver. 4 peer-reviewed and recommended by Peer Community in Evolutionary Biology. https://doi.org/10.1101/2021.11.21.469450

Beadell JS, Ishtiaq F, Covas R, Melo M, Warren BH, Atkinson CT, Bensch S, Graves GR, Jhala YV, Peirce MA, Rahmani AR, Fonseca DM, Fleischer RC (2006) Global phylogeographic limits of Hawaii’s avian malaria. Proceedings of the Royal Society B: Biological Sciences, 273, 2935–2944. https://doi.org/10.1098/rspb.2006.3671

Burns KC (2019) Evolution in Isolation: The Search for an Island Syndrome in Plants. Cambridge University Press, Cambridge. https://doi.org/10.1017/9781108379953

Doutrelant C, Paquet M, Renoult JP, Grégoire A, Crochet P-A, Covas R (2016) Worldwide patterns of bird colouration on islands. Ecology Letters, 19, 537–545. https://doi.org/10.1111/ele.12588

Illera JC, Fernández-Álvarez Á, Hernández-Flores CN, Foronda P (2015) Unforeseen biogeographical patterns in a multiple parasite system in Macaronesia. Journal of Biogeography, 42, 1858–1870. https://doi.org/10.1111/jbi.12548

Lomolino MV (2005) Body size evolution in insular vertebrates: generality of the island rule. Journal of Biogeography, 32, 1683–1699. https://doi.org/10.1111/j.1365-2699.2005.01314.x

Matson KD (2006) Are there differences in immune function between continental and insular birds? Proceedings of the Royal Society B: Biological Sciences, 273, 2267–2274. https://doi.org/10.1098/rspb.2006.3590

Pérez-Rodríguez A, Ramírez Á, Richardson DS, Pérez-Tris J (2013) Evolution of parasite island syndromes without long-term host population isolation: parasite dynamics in Macaronesian blackcaps Sylvia atricapilla. Global Ecology and Biogeography, 22, 1272–1281. https://doi.org/10.1111/geb.12084

Most mutations affecting fitness are deleterious and they have many evolutionary consequences. The dynamics and consequences of deleterious mutations are a long-standing question in evolutionary biology and a strong theoretical background has already been developed, for example, to predict the mutation load, inbreeding depression or background selection. One of the classical results is that inbreeding helps purge partially recessive deleterious mutations by exposing them to selection in homozygotes. However, this mainly results from single-locus considerations. When interactions among several, more or less linked, deleterious mutations are taken into account, peculiar dynamics can emerge. One of them, called pseudo-overdominance (POD), corresponds to the maintenance in a population of two (or more) haplotype blocks composed of several recessive deleterious mutations in repulsion that mimics overdominance. Indeed, homozygote individuals for one of the haplotype blocks expose many deleterious mutations to selection whereas they are reciprocally masked in heterozygotes, leading to higher fitness of heterozygotes compared to both homozygotes. A related process, called associative overdominance (AOD) is the effect of such deleterious alleles in repulsion on the linked neutral variation that can be increased by AOD. Although this possibility has been recognized for a long time (Otha and Kimura 1969), it has been mainly considered an anecdotal process. Recently, both theoretical (Zhao and Charlesworth 2016) and genomic analyses (Gilbert et al. 2020) have renewed interest in such a process, suggesting that it could be important in weakly recombining regions of a genome. Donald Waller (2021) - one of the co-authors of the current work - also recently proposed that POD could be quantitatively important with broad implications, and could resolve some unexplained observations such as the maintenance of inbreeding depression in highly selfing species. Yet, a proper theoretical framework analysing the effect of inbreeding on POD was lacking.

In this theoretical work, Diala Abu Awad and Donald Waller (2022) addressed this question through an elegant combination of analytical predictions and intensive multilocus simulations. They determined the conditions under which POD can be maintained and how long it could resist erosion by recombination, which removes the negative association between deleterious alleles (repulsion) at the core of the mechanism. They showed that under tight linkage, POD regions can persist for a long time and generate substantial segregating load and inbreeding depression, even under inbreeding, so opposing (for a while) to the purging effect. They also showed that background selection can affect the genomic structure of POD regions by rapidly erasing weak POD regions but maintaining strong POD regions (i.e with many tightly linked deleterious alleles).

These results have several implications. They can explain the maintenance of inbreeding depression despite inbreeding (as anticipated by Waller 2021), which has implications for the evolution of mating systems. If POD can hardly emerge under high selfing, it can persist from an outcrossing ancestor long after the transition towards a higher selfing rate and could explain the maintenance of mixed mating systems(which is possible with true overdominance, see Uyenoyama and Waller 1991). The results also have implications for genomic analyses, pointing to regions of low or no recombination where POD could be maintained, generating both higher diversity and heterozygosity than expected and variance in fitness. As structural variations are likely widespread in genomes with possible effects on suppressing recombination (Mérot et al. 2020), POD regions should be checked more carefully in genomic analyses (see also Gilbert et al. 2020).

Overall, this work should stimulate new theoretical and empirical studies, especially to assess how quantitatively strong and widespread POD can be. It also stresses the importance of properly considering genetic linkage genome-wide, and so the role of recombination landscapes in determining patterns of diversity and fitness effects.

Awad DA, Waller D (2022) Conditions for maintaining and eroding pseudo-overdominance and its contribution to inbreeding depression. bioRxiv, 2021.12.16.473022, ver. 3 peer-reviewed and recommended by Peer Community in Evolutionary Biology. https://doi.org/10.1101/2021.12.16.473022

Gilbert KJ, Pouyet F, Excoffier L, Peischl S (2020) Transition from Background Selection to Associative Overdominance Promotes Diversity in Regions of Low Recombination. Current Biology, 30, 101-107.e3. https://doi.org/10.1016/j.cub.2019.11.063

Mérot C, Oomen RA, Tigano A, Wellenreuther M (2020) A Roadmap for Understanding the Evolutionary Significance of Structural Genomic Variation. Trends in Ecology & Evolution, 35, 561–572. https://doi.org/10.1016/j.tree.2020.03.002

Ohta T, Kimura M (1969) Linkage disequilibrium at steady state determined by random genetic drift and recurrent mutation. Genetics, 63, 229–238. https://doi.org/10.1093/genetics/63.1.229

Uyenoyama MK, Waller DM (1991) Coevolution of self-fertilization and inbreeding depression II. Symmetric overdominance in viability. Theoretical Population Biology, 40, 47–77. https://doi.org/10.1016/0040-5809(91)90046-I

Waller DM (2021) Addressing Darwin’s dilemma: Can pseudo-overdominance explain persistent inbreeding depression and load? Evolution, 75, 779–793. https://doi.org/10.1111/evo.14189

Zhao L, Charlesworth B (2016) Resolving the Conflict Between Associative Overdominance and Background Selection. Genetics, 203, 1315–1334. https://doi.org/10.1534/genetics.116.188912

Antimalarial resistance is a major hurdle to malaria eradication efforts. The spread of drug resistance follows basic evolutionary principles, with competitive interactions between resistant and susceptible malaria strains being central to the fitness of resistant parasites. These competitive interactions can be used to design resistance management strategies, whereby a fitness cost of resistant parasites can be exploited through maintaining competitive suppression of the more fit drug-susceptible parasites. This can potentially be achieved using lower drug dosages or lower frequency of drug treatments. This approach has been demonstrated to work empirically in a rodent malaria model [1,2] and has been demonstrated to have clinical success in cancer treatments [3]. However, these resistance management approaches assume a fitness cost of the resistant pathogen, and, in the case of malaria parasites in general, and for artemisinin resistant parasites in particular, there is limited information on the presence of such fitness cost. The best suggestive evidence for the presence of fitness costs comes from the discontinuation of the use of the drug, which, in the case of chloroquine, was followed by a gradual drop in resistance frequency over the following decade [see e.g. 4,5]. However, with artemisinin derivative drugs still in use, alternative ways to study the presence of fitness costs need to be undertaken. 
There are several good in vitro studies demonstrating artemisinin resistant parasites being competitively suppressed by wildtype parasites [see e.g. 6–9], however, these have the limitation that they will only be able to detect the fitness cost during the blood stage of the infection and in an artificial environment. So far, there have not been animal models that have thoroughly studied the presence of resistance fitness costs for artemisinin resistant parasites. Moreover, in these types of studies, the focus is mostly on the fitness cost as detected in the vertebrate host. However, malaria parasites spent a significant portion of their life cycle in the mosquito host, where fitness costs could also be expressed. Overall, it is the fitness over the entire life cycle of the parasite that would determine if and to what extent a reduction in resistance frequency would be observed when the use of a drug is stopped. 
Here, Villa and colleagues present a study to quantify such fitness cost of artesunate-resistant parasites, not only in a vertebrate host, but also in the mosquito vector [10]. They used the underutilized model system of avian malaria species Plasmodium relictum in canaries. Villa and colleagues selected for several different resistance strains, which had a similar delayed clearance phenotype as observed in the field. Interestingly, they did not find evidence of a fitness cost in the vertebrate host. In fact, the resistant strains reached greater parasitaemia than the susceptible strains. From this set of experiments it is unclear whether this is an anomaly or a relevant result. Future work should establish this, though fitness benefits associated with drug resistance have been seen before in Leishmania parasites [11]. An important caveat to the present study is that the parasites were grown in the absence of competition and it is feasible that a cost is not detected when growing by themselves, but would become apparent when in competition. However, these types of experiments are technologically more challenging to perform as it would require an accurate quantification methodology able to distinguish based on one SNP. This problem has been circumvented by either using relative peak height in sanger sequencing [12], or via the likely more accurate route of pyrosequencing [7–9], though these methodologies only give relative frequencies rather than absolute densities. 
 
The most interesting observation in the study by Villa et al is that the authors detected a fitness cost being played out in the mosquito vector, where the resistant strains had a decreased infectivity compared to the susceptible strain. This finding is important because 1) it demonstrates that the whole life cycle needs to be taken into account when understanding fitness costs, 2) resistance management strategies that are based on treatment within the vertebrate host may not have the intended effect if the cost does not play out in this host, and 3) it opens new research avenues to explore the possibility of exploiting fitness costs in mosquito vector. Future research should focus on incorporating these assays on fitness costs in mosquitoes, particularly for P. falciparum parasites. Additionally, it would be interesting to expand this work in a competitive environment, both in the vertebrate host as in the mosquito host. Finally, it would be important to establish the generalizability of such fitness cost in mosquitoes. If it is a significant factor, mathematical models could incorporate this effect in predictions on the spread of resistance.

References

[1] Huijben S, Bell AS, Sim DG, Tomasello D, Mideo N, Day T, et al. 2013. Aggressive chemotherapy and the selection of drug resistant pathogens. PLoS Pathog. 9(9): e1003578. https://doi.org/10.1371/journal.ppat.1003578
 
[2] Huijben S, Nelson WA, Wargo AR, Sim DG, Drew DR, Read AF. 2010. Chemotherapy, within-host ecology and the fitness of drug-resistant malaria parasites. Evolution (N Y). 64(10): 2952-68. https://doi.org/10.1111/j.1558-5646.2010.01068.x
 
[3] Zhang J, Cunningham JJ, Brown JS, Gatenby RA. 2017. Integrating evolutionary dynamics into treatment of metastatic castrate-resistant prostate cancer. Nat Commun. 8(1). https://doi.org/10.1038/s41467-017-01968-5
 
[4] Laufer MK, Takala-Harrison S, Dzinjalamala FK, Stine OC, Taylor TE, Plowe C v. 2010. Return of chloroquine-susceptible falciparum malaria in Malawi was a reexpansion of diverse susceptible parasites. J Infect Dis. 202(5): 801-808. https://doi.org/10.1086/655659 

[5] Mharakurwa S, Matsena-Zingoni Z, Mudare N, Matimba C, Gara TX, Makuwaza A, et al. 2021. Steep rebound of chloroquine-sensitive Plasmodium falciparum in Zimbabwe. J Infect Dis. 223(2): 306-9. https://doi.org/10.1093/infdis/jiaa368
 
[6] Tirrell AR, Vendrely KM, Checkley LA, Davis SZ, McDew-White M, Cheeseman IH, et al. 2019. Pairwise growth competitions identify relative fitness relationships among artemisinin resistant Plasmodium falciparum field isolates. Malar J. 18: 295. https://doi.org/10.1186/s12936-019-2934-4
 
[7] Hott A, Tucker MS, Casandra D, Sparks K, Kyle DE. 2015. Fitness of artemisinin-resistant Plasmodium falciparum in vitro. J Antimicrob Chemother. 70(10): 2787-2796. https://doi.org/10.1093/jac/dkv199
 
[8] Straimer J, Gnädig NF, Stokes BH, Ehrenberger M, Crane AA, Fidock DA. 2017. Plasmodium falciparum K13 mutations differentially impact ozonide susceptibility and parasite fitness in vitro. mBio. 8(2): e00172-17. https://doi.org/10.1128/mBio.00172-17
 
[9] Nair S, Li X, Arya GA, McDew-White M, Ferrari M, Nosten F, et al. 2018. Fitness costs and the rapid spread of kelch13-C580Y substitutions conferring artemisinin resistance. Antimicrob Agents Chemother. 62(9). https://doi.org/10.1128/AAC.00605-18
 
[10] Villa M, Berthomieu A, Rivero A. Fitness costs and benefits in response to artificial artesunate selection in Plasmodium. 2022. bioRxiv, 20220128478164, ver 3 peer-reviewed and recommended by Peer Community in Evolutionary Biology. https://doi.org/10.1101/2022.01.28.478164
 
[11] Vanaerschot M, Decuypere S, Berg M, Roy S, Dujardin JC. 2013. Drug-resistant microorganisms with a higher fitness--can medicines boost pathogens? Crit Rev Microbiol. 39(4): 384-394. https://doi.org/10.3109/1040841X.2012.716818
 
[12] Hassett MR, Roepe PD. In vitro growth competition experiments that suggest consequences of the substandard artemisinin epidemic that may be accelerating drug resistance in P. falciparum malaria. 2021. PLoS One. 16(3): e0248057. https://doi.org/10.1371/journal.pone.0248057