The increase in the prevalence of (antibiotic) resistance has become a major global health concern and is an excellent example of the impact of real-time evolution on human society. This has led to a boom of studies that investigate the mechanisms and factors involved in the evolution of resistance, and to the spread of the concept of "costs of resistance". This concept refers to the relative fitness disadvantage of a drug-resistant genotype compared to a non-resistant reference genotype in the ancestral (untreated) environment.

In their paper, Lenormand et al. [1] discuss the history of this concept and highlight its caveats and limitations. The authors address both practical and theoretical problems that arise from the simplistic view of "costly resistance" and argue that they can be prejudicial for antibiotic resistance studies. For a better understanding, they visualize their points of critique by means of Fisher's Geometric model.

The authors give an interesting historical overview of how the concept arose and speculate that it emerged (during the 1980s) in an attempt by ecologists to spread awareness that fitness can be environment-dependent, and because of the concept's parallels to trade-offs in life-history evolution. They then identify several problems that arise from the concept, which, besides the conceptual misunderstandings that they can cause, are important to keep in mind when designing experimental studies.

The authors highlight and explain the following points:
1. Costs of resistance do not necessarily imply pleiotropic effects of a resistance mutation, and pleiotropy is not necessarily the cause of fitness trade-offs.
2. Any non-treated environment and any treatment dose can result in a different cost.
3. Different reference genotypes may result in different costs. Specifically, the reference genotype has to be "optimally" adapted to the reference environment to provide an accurate measurement of costs.

Lenormand et al.'s paper [1] is a timely perspective piece in light of the ever-increasing efforts to understand and tackle resistance evolution [2]. Although some readers may shy away from the rather theoretical presentation of the different points of concern, it will be useful for both theoretical and empirical readers by illustrating the misconceptions that can arise from the concept of the cost of resistance. Ultimately, the main lesson to be learned from this paper may not be to ban the term "cost of resistance" from one's vocabulary, but rather to realize that the successful fight against drug resistance requires more differential information than the measurement of fitness effects in a drug-treated vs. non-treated environment in the lab [3-4]. Specifically, a better integration of the ecological aspects of drug resistance evolution and maintenance is needed [5], and we are far from a general understanding of how environmental factors interact and influence an organism's (absolute and relative) fitness and the effect of resistance mutations.

References

[1] Lenormand T, Harmand N, Gallet R. 2018. Cost of resistance: an unreasonably expensive concept. bioRxiv 276675, ver. 3 peer-reviewed by Peer Community In Evolutionary Biology. doi: 10.1101/276675
[2] Andersson DI and Hughes D. Persistence of antibiotic resistance in bacterial populations. 2011. FEMS Microbiology Reviews, 35: 901-911. doi: 10.1111/j.1574-6976.2011.00289.x
[3] Chevereau G, Dravecká M, Batur T, Guvenek A, Ayhan DH, Toprak E, Bollenbach T. 2015. Quantifying the determinants of evolutionary dynamics leading to drug resistance. PLoS biology 13, e1002299. doi: 10.1371/journal.pbio.1002299
[4] Bengtsson-Palme J, Kristiansson E, Larsson DGJ. 2018. Environmental factors influencing the development and spread of antibiotic resistance. FEMS Microbiology Reviews 42: 68–80. doi: 10.1093/femsre/fux053
[5] Hiltunen T, Virta M, Laine AL. 2017. Antibiotic resistance in the wild: an eco-evolutionary perspective. Philosophical Transactions of the Royal Society B: Biological Sciences 372: 20160039. doi: 10.1098/rstb.2016.0039

The influence of environmental change over geological time on the tempo and mode of biological diversification is a hot topic in biogeography. Of central interest are questions about where, when, and how fast lineages proliferated, suffered extinction, and migrated in response to tectonic events, the waxing and waning of dominant biomes, etc. In this context, the dynamic conditions of the Miocene have received much attention, from studies of many clades and biogeographic regions. Here, Chazot et al. [1] present an exemplary analysis of butterflies (tribe Ithomiini) in the Neotropics, examining their diversification across the Andes and Amazon. They infer sharp contrasts between these regions in the late Miocene: accelerated diversification during orogeny of the Andes, and greater extinction in the Amazon associated during the Pebas system, with interchange and local diversification increasing following the Pebas during the Pliocene. Two features of this study stand out. First is the impressive taxon sampling (340 out of 393 extant species). Second is the use of ancestral range reconstructions to compute per-lineage rates of colonization between regions, and rates of speciation within regions, through time. The latter allows for relatively fine-grained comparisons across the 2 fundamental dimensions of historical biogeography, space and time, and is key to the main results. The method resonated with me because I performed a similar analysis in a study showing evidence for uplift-driven diversification in the Hengduan Mountains of China [2]. This analysis is complemented by a variety of other comparative methods for inferring variable diversification across clades, through time, and in response to external factors. Overall, it represents a very nice contribution to our understanding of the effects of Miocene/Pliocene environmental change on the evolution of Neotropical biodiversity. **References** [1] Chazot N, Willmott KR, Lamas G, Freitas AVL, Piron-Prunier F, Arias CF, Mallet J, De-Silva DL and Elias M. 2017. Renewed diversification following Miocene landscape turnover in a Neotropical butterfly radiation. BioRxiv 148189, ver 4 of 19th December 2017. doi: [10.1101/148189](https://doi.org/10.1101/148189) [2] Xing Y, and Ree RH. 2017. Uplift-driven diversification in the Hengduan Mountains, a temperate biodiversity hotspot. Proceedings of the National Academy of Sciences of the United States of America, 114: E3444-E3451. doi: [10.1073/pnas.1616063114](https://doi.org/10.1073/pnas.1616063114)

Many viruses are transmitted by biological vectors, i.e. organisms that transfer the virus from one host to another. Dengue virus (DENV) is one of them. Dengue is a mosquito-borne viral disease that has rapidly spread around the world since the 1940s. One recent estimate indicates 390 million dengue infections per year [1]. As many arthropod-borne vertebrate viruses, DENV has to cross several anatomical barriers in the vector, to multiply in its body and to invade its salivary glands before getting transmissible. As a consequence, vectors are not passive carriers but genuine hosts of the viruses that potentially have important effects on the composition of virus populations and, ultimately, on virus epidemiology and virulence. Within infected vectors, virus populations are expected to acquire new mutations and to undergo genetic drift and selection effects. However, the intensity of these evolutionary forces and the way they shape virus genetic diversity are poorly known.

In their study, Lequime et al. [2] finely disentangled the effects of genetic drift and selection on DENV populations during their infectious cycle within mosquito (Aedes aegypti) vectors. They evidenced that the genetic diversity of viruses within their vectors is shaped by genetic drift, selection and vector genotype. The experimental design consisted in artificial acquisition of purified virus by mosquitoes during a blood meal. The authors monitored the diversity of DENV populations in Ae. aegypti individuals at different time points by high-throughput sequencing (HTS). They estimated the intensity of genetic drift and selection effects exerted on virus populations by comparing the DENV diversity at these sampling time points with the diversity in the purified virus stock (inoculum).

Disentangling the effects of genetic drift and selection remains a methodological challenge because both evolutionary forces operate concomitantly and both reduce genetic diversity. However, selection reduces diversity in a reproducible manner among experimental replicates (here, mosquito individuals): the fittest variants are favoured at the expense of the weakest ones. In contrast, genetic drift reduces diversity in a stochastic manner among replicates. Genetic drift acts equally on all variants irrespectively of their fitness. The strength of genetic drift is frequently evaluated with the effective population size Ne: the lower Ne, the stronger the genetic drift [3]. The estimation of the effective population size of DENV populations by Lequime et al. [2] was based on single-nucleotide polymorphisms (SNPs) that were (i) present both in the inoculum and in the virus populations sampled at the different time points and (ii) that were neutral (or nearly-neutral) and therefore subjected to genetic drift only and insensitive to selection. As expected for viruses that possess small and constrained genomes, such neutral SNPs are extremely rare. Starting from a set of >1800 SNPs across the DENV genome, only three SNPs complied with the neutrality criteria and were enough represented in the sequence dataset for a precise Ne estimation. Using the method described by Monsion et al. [4], Lequime et al. [2] estimated Ne values ranging from 5 to 42 viral genomes (95% confidence intervals ranged from 2 to 161 founding viral genomes). Consequently, narrow bottlenecks occurred at the virus acquisition step, since the blood meal had allowed the ingestion of ca. 3000 infectious virus particles, on average. Interestingly, bottleneck sizes did not differ between mosquito genotypes. Monsion et al.’s [4] formula provides only an approximation of Ne. A corrected formula has been recently published [5]. We applied this exact Ne formula to the means and variances of the frequencies of the three neutral markers estimated before and after the bottlenecks (Table 1 in [2]), and nearly identical Ne estimates were obtained with both formulas.

Selection intensity was estimated from the dN/dS ratio between the nonsynonymous and synonymous substitution rates using the HTS data on DENV populations. DENV genetic diversity increased following initial infection but was restricted by strong purifying selection during virus expansion in the midgut. Again, no differences were detected between mosquito genotypes. However and importantly, significant differences in DENV genetic diversity were detected among mosquito genotypes. As they could not be related to differences in initial genetic drift or to selection intensity, the authors raise interesting alternative hypotheses, including varying rates of de novo mutations due to differences in replicase fidelity or differences in the balancing selection regime. Interestingly, they also suggest that this observation could simply result from a methodological issue linked to the detection threshold of low-frequency SNPs.
 

References

[1] Bhatt S, Gething PW, Brady OJ, Messina JP, Farlow AW, Moyes CL, Drake JM, et al. 2013. The global distribution and burden of dengue. Nature 496: 504–7 doi: 10.1038/nature12060

[2] Lequime S, Fontaine A, Gouilh MA, Moltini-Conclois I and Lambrechts L. 2016. Genetic drift, purifying selection and vector genotype shape dengue virus intra-host genetic diversity in mosquitoes. PloS Genetics 12: e1006111 doi: 10.1371/journal.pgen.1006111

[3] Charlesworth B. 2009. Effective population size and patterns of molecular evolution and variation. Nature Reviews Genetics 10: 195-205 doi: 10.1038/nrg2526

[4] Monsion B, Froissart R, Michalakis Y and Blanc S. 2008. Large bottleneck size in cauliflower mosaic virus populations during host plant colonization. PloS Pathogens 4: e1000174 doi: 10.1371/journal.ppat.1000174

[5] Thébaud G and Michalakis Y. 2016. Comment on ‘Large bottleneck size in cauliflower mosaic virus populations during host plant colonization’ by Monsion et al. (2008). PloS Pathogens 12: e1005512 doi: 10.1371/journal.ppat.1005512

Theoretical work, initiated by Levene (1953) [1] and Dempster (1955) [2], suggests that within a given environment, the way populations are regulated and contribute to the next generation is a key factor for the maintenance of local adaptation polymorphism. In this theoretical context, hard selection describes the situation where the genetic composition of each population affects its contribution to the next generation whereas soft selection describes the case where the contribution of each population is fixed, whatever its genetic composition. Soft selection is able to maintain polymorphism, whereas hard selection invariably leads to the fixation of one of the alleles. Although the specific conditions (e.g. of migration between populations or drift level) in which this prediction holds have been studied in details by theoreticians, experimental tests have mainly failed, usually leading to the conclusion that the allele frequency dynamics was driven by other mechanisms in the experimental systems and conditions used. Gallet, Froissart and Ravigné [3] have set up a bacterial experimental system which allowed them to convincingly demonstrate that soft selection generates the conditions for polymorphism maintenance when hard selection does not, everything else being equal. The key ingredients of their experimental system are (1) the possibility to accurately produce hard and soft selection regimes when daily transferring the populations and (2) the ability to establish artificial well-characterized reproducible trade-offs. To do so, they used two genotypes resisting each one to one antibiotic and combined, across habitats, low antibiotic doses and difference in medium productivity. The experimental approach contains two complementary parts: the first one is looking at changes in the frequencies of two genotypes, initially introduced at around 50% each, over a small number of generations (ca 40) in different environments and selection regimes (soft/hard) and the second one is convincingly showing polymorphism protection by establishing that in soft selection regimes, the lowest fitness genotype is not eliminated even when introduced at low frequency. In this manuscript, a key point is the dialog between theoretical and experimental approaches. The experiments have been thought and designed to be as close as possible to the situations analysed in theoretical work. For example, the experimental polymorphism protection test (experiment 2) closely matches the equivalent analysis classically performed in theoretical approaches. This close fit between theory and experiment is clearly a strength of this study. This said, the experimental system allowing them to realise this close match also has some limitations. For example, changes in allele frequencies could only be monitored over a quite low number of generations because a longer time-scale would have allowed the contribution of de novo mutations and the likely emergence of a generalist genotype resisting to both antibiotics used to generate the local adaptation trade-offs. These limitations, as well as the actual significance of the experimental tests, are discussed in deep details in the manuscript.

References

[1] Levene H. 1953. Genetic equilibrium when more than one niche is available. American Naturalist 87: 331–333. doi: 10.1086/281792

[2] Dempster ER. 1955. Maintenance of genetic heterogeneity. Cold Spring Harbor Symposia on Quantitative Biology. 20: 25–32. doi: 10.1101/SQB.1955.020.01.005

[3] Gallet R, Froissart R, Ravigné V. 2017. Things softly attained are long retained: dissecting the impacts of selection regimes on polymorphism maintenance in experimental spatially heterogeneous environments. bioRxiv 100743; doi: 10.1101/100743

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