As plants cannot run from their enemies, natural selection has favoured the evolution of diverse chemical compounds (phytochemicals) to protect them against herbivores and pathogens. This provides an opportunity for plant feeders to exploit these compounds to combat their own enemies. Indeed, it is widely known that herbivores use such compounds as protection against predators [1]. Recently, this reasoning has been extended to pathogens, and elegant studies have shown that some herbivores feed on phytochemical-containing plants reducing both parasite abundance within hosts [2] and their virulence [3].
Suffering less from parasites is clearly beneficial for infected herbivores. Why then, is this behaviour not fixed in nature? There are, at least, two possible explanations. First, feeding on ‘medicinal’, often toxic, plants may impose costs upon uninfected individuals. Second, parasites may evolve resistance to such chemicals. Whereas the first possibility has been explored, and is taken as evidence for ‘self-medication’ [3], the second hypothesis requires investigation. A recent study by Palmer-Young et al. [4] fills this gap. This article reports evolution of resistance to two different phytochemicals, alone and in combination, in the trypanosome Crithidia bombi, a bumble bee (Bombus impatiens) parasite. To show this, the authors experimentally evolved a strain of C. bombi in the presence of thymol, eugenol or both simultaneously. These phytochemicals are commonly found in the nectar of several plant species, in particular those of the Lamiaceae, a family containing several aromatic herbs. Prior to selection both phytochemicals reduced C. bombi growth by about 50%. However, C. bombi rapidly evolved resistance in both single and the double phytochemical treatments. Moreover, no cost of resistance was detected when evolved parasites were placed in the ancestral, phytochemical-free environment. Therefore, resistance to phytochemicals would be expected to spread rapidly in natural populations of C. bombi. Clearly, thus, the herbivore strategy of sequestering plant secondary chemical compounds as a defence against their pathogens should fail. The question then is ‘why do they still do it’? Can self-medication work in the longer-term for bumblebees?
Well, yes. The very fact that resistance evolved shows that resistance is not fixed in natural C. bombi populations. This is surprising considering that resistance is not costly. This might be due to a number of reasons. Firstly, there may be costs of resistance that were not detected in this experiment. Second, it may not be possible to evolve universal resistance to the heterogeneity present in the phytochemical environment. Indeed, phytochemical environments are highly varied in time and space and bumblebees will visit different plants presenting different phytochemical cocktails throughout the season. Furthermore, migration of bees from populations exposed to different phytochemicals may prevent the fixation of one resistance type.
Or, it may be self-medication behaviour itself that prevents the evolution of resistance? Indeed, in the same way that infected bees seek cooler temperatures to slow growth of a parasitoid fly [5], they may also seek a more varied diet with diverse phytochemicals to which the parasite cannot evolve, but which reduces parasite growth. Further understanding of arthropod self-medication may thus be instrumental to prevent the observed worldwide decline of pollinators [6]. Furthermore, it may inform on mechanisms that prevent rapid evolution of drug resistance in other systems.

References

[1] Duffey SS. 1980. Sequestration of plant natural products by insects. Annual Review of Entomology 25: 447-477. doi: 10.1146/annurev.en.25.010180.002311

[2] Richardson LL et al. 2015. Secondary metabolites in floral nectar reduce parasite infections in bumblebees. Proceedings of the Royal Society of London B 282: 20142471. doi: 10.1098/rspb.2014.2471

[3] Lefèvre T et al. 2010. Evidence for trans-generational medication in nature. Ecology Letters 13: 1485-93. doi: 10.1111/j.1461-0248.2010.01537.x

[4] Palmer-Young EC, Sadd BM, Adler LS. 2017. Evolution of resistance to single and combined floral phytochemicals by a bumble bee parasite. Journal of Evolutionary Biology 30: 300-312. doi: 10.1111/jeb.13002

[5] Müller CB, Schmid-Hempel P. 1993. Exploitation of cold temperature as defence against parasitoids in bumblebees. Nature 363: 65-67. doi: 10.1038/363065a0

[6] Potts SG et al. 2010. Global pollinator declines: trends, impacts and drivers. Trends in Ecology and Evolution 25: 345-353. doi: 10.1016/j.tree.2010.01.007

Estimating the genetic and environmental components of variation in reproductive success is crucial to understanding the adaptive potential of populations to environmental change. To date, the heritability of lifetime reproductive success (fitness) has been estimated in a handful of wild animal population, mostly in mammals and birds, but has never been estimated for a marine species. The primary reason that such estimates are lacking in marine species is that most marine organisms have a dispersive larval phase, making it extraordinarily difficult to track the fate of offspring from one generation to the next.
In this study, Salles et al. [1] use an unprecedented 10 year data set for a wild population of orange clownfish (Amphiprion percula) to estimate the environmental, maternal and additive genetic components of life time reproductive success for the self-recruiting portion of the local population. Previous studies show that over 50% of juvenile clownfish recruiting to the population of clownfish at Kimbe Island (Kimbe Bay, PNG) are natal to the population. In other words, >50% of the juveniles recruiting to the population at Kimbe Island are offspring of parents from Kimbe Island. The identity and location of every adult clownfish in the Kimbe Island population was tracked over 10 years. At the same time newly recruiting juveniles were collected at regular intervals (biennially) and their parentage assigned with high confidence by 22 polymorphic microsatellite loci. Salles et al. then used a pedigree comprising 1735 individuals from up to 5 generations of clownfish at Kimbe Island to assess the contribution of every breeding pair of clownfish to self-recruitment within the local population. Because clownfish are site attached and live in close association with a host sea anemone, it was also possible to examine the contribution of reef location and host anemones species (either Heteractis magnifica or Stichodactyla gigantea) to reproductive success within the local population.
The study found that breeders from the eastern side of Kimbe Island, and mostly inhabiting S. gigantea sea anemones, produced more juveniles that recruited to the local population than breeders from other location around the island, or inhabiting H. magnifica. In fact, host anemone species and geographic location explained about 97% of the variance in reproductive success within the local population (i.e. excluding successful recruitment to other populations). By contrast, maternal and additive genetic effects explained only 1.9% and 1.3% of the variance, respectively. In other words, reef location and the species of host anemone inhabited had an overwhelming influence on the long-term contribution of breeding pairs of clownfish to replenishment of the local population. This overwhelming effect of the local habitat on reproductive success means that the population is potentially susceptible to rapid environmental changes - for example if S. giganta sea anemones are disproportionately susceptible to global warming, or reef habitats on the eastern side of the island are more susceptible to disturbance. By contrast, the small component of additive genetic variance in local reproductive success translated into low heritability and evolvability of lifetime reproductive success within the local population, as predicted by theory [2] and observed in some terrestrial species. Consequently, fitness would evolve slowly to environmental change.
Establishing the components of variation in fitness in a wild population of marine fishes is an astonishing achievement, made possible by the unprecedented long-term individual-level monitoring of the entire population of clownfish at Kimbe Island. A next step in this research would be to include other clownfish populations that are demographically and genetically connected to the Kimbe Island population through larval dispersal. It would be intriguing to establish the environmental, maternal and additive genetic components of reproductive success in the dispersing part of the Kimbe Island population, to see if this potentially differs among breeders who contribute more or less to replenishment within the local population.

References

[1] Salles, O. C., Almany, G. R., Berumen, M.L., Jones, G. P., Saenz-Agudelo, P., Srinivasan, M., Thorrold, S. R., Pujol, B., Planes, S. (2019). Strong habitat and weak genetic effects shape the lifetime reproductive success in a wild clownfish population. Zenodo, 3476529, ver. 3 peer-reviewed and recommended by Peer Community In Evolutionary Biology. doi: 10.5281/zenodo.3476529
[2] Fisher, R.A. (1930). The genetical theory of natural selection. Clarendon Press, Oxford, U.K.

In spite of its name, the effective population size, Ne, has a complex and often distant relationship to census population size, as we usually understand it. In truth, it is primarily an abstract concept aimed at measuring the amount of genetic drift occurring in a population at any given time. The standard way to model random genetic drift in population genetics is the Wright-Fisher model and, with a few exceptions, definitions of the effective population size stems from it: “a certain model has effective population size, Ne, if some characteristic of the model has the same value as the corresponding characteristic for the simple Wright-Fisher model whose actual size is Ne” (Ewens 2004). Since Sewall Wright introduced the concept of effective population size in 1931 (Wright 1931), it has flourished and there are today numerous definitions of it depending on the process being examined (genetic diversity, loss of alleles, efficacy of selection) and the characteristic of the model that is considered. These different definitions of the effective population size were generally introduced to address specific aspects of the evolutionary process. One aspect that has been hotly debated since the first estimates of genetic diversity in natural populations were published is the so-called Lewontin’s paradox (1974). Lewontin noted that the observed variation in heterozygosity across species was much smaller than one would expect from the neutral expectations calculated with the actual size of the species.
In essence, what Galtier and Rousselle propose in their clever paper is to introduce a new approach to compare effective population sizes across species and thereby a new way to address Lewontin’s paradox. Classically, the effective population size in this type of comparative genomic studies is simply estimated from nucleotide diversity at putatively neutral sites using the equation relating levels of diversity (θ) to mutation rate per generation (μ) and effective population size, Ne, θ = 4Neμ. As Galtier and Rousselle point out there are many issues with this approach. In particular, although we can now estimate θ very precisely, we generally do not have a reliable estimate of the mutation rate, and the method rests on many, unwarranted, assumptions; for example that the population is at mutation-drift equilibrium. Instead they propose to estimate the effective population size from the load of segregating deleterious mutations which can be summarized by the ratio of nonsynonymous to synonymous mutations, πN/πS: small-Ne species are expected to accumulate more deleterious mutations and carry a higher load than large-Ne ones at selection/drift equilibrium (Ohta et al. 1973; Welch et al. 2008). At first glance, this suggestion seems counterintuitive since considering sites under selection undoubtedly adds a new layer of complexity to an already intricate situation. Indeed, one is now bringing to the brew another elusive object, namely, the Distribution of Fitness Effect of mutations (DFE). However, estimating Ne from the load of segregating deleterious mutations may actually simplify the situation in two important ways. First, using πN/πS does not require assumption about μ (as the mutation rate will cancel out in the ratio). Second, the ratio of nonsynonymous to synonymous nucleotide diversity, reaches equilibrium faster than the nucleotide diversity at synonymous site after a change in population size, so we can hope for less sensitivity to (often unknown) recent demographic history (Brandvain and Wright 2016).
Extending recent developments in the estimation of the DFE, Galtier and Rousselle eventually obtain estimates of the average deleterious effect, , where Ne is the effective population size and is the mean fitness effect of non-synonymous mutations. Assuming further that distinct species share a common DFE and therefore a common , they obtain estimates of the between species ratio of Ne from the between species ratio of . Applying their newly developed approach to various datasets they conclude that the power of drift varies by a factor of at least 500 between large-Ne (Drosophila) and small-Ne species (H. sapiens). This is an order of magnitude larger than what would be obtained by comparing estimates of the variation in neutral diversity. Hence the proposed approach seems to have gone some way in making Lewontin’s paradox less paradoxical. But, perhaps more importantly, as the authors tersely point out at the end of the abstract their results further questions the meaning of Ne parameters in population genetics. And arguably this could well be the most important contribution of their study and something that is badly needed.

References

Brandvain Y, Wright SI (2016) The Limits of Natural Selection in a Nonequilibrium World. Trends in Genetics, 32, 201–210. doi: 10.1016/j.tig.2016.01.004
Ewens WJ (2010) Mathematical Population Genetics: Theoretical Introduction. Springer-Verlag New York Inc., New York, NY. doi: 10.1007/978-0-387-21822-9
Galtier N, Rousselle M (2020) How much does Ne vary among species? bioRxiv, 861849, ver. 4 peer-reviewed and recommended by PCI Evolutionary Biology. doi: 10.1101/861849
Lewontin RC (1974) The genetic basis of evolutionary change. Columbia University Press, New York. Ohta T (1973) Slightly Deleterious Mutant Substitutions in Evolution. Nature, 246, 96–98. doi: 10.1038/246096a0
Welch JJ, Eyre-Walker A, Waxman D (2008) Divergence and Polymorphism Under the Nearly Neutral Theory of Molecular Evolution. Journal of Molecular Evolution, 67, 418–426. doi: 10.1007/s00239-008-9146-9
Wright S (1931) Evolution in Mendelian Populations. Genetics, 16, 97–159. url: https://www.genetics.org/content/16/2/97

In a landmark paper, Flores et al. (2011) showed that the interactions between bacteria and their viruses could be nicely described using a bipartite infection networks.  Two quantitative properties of these networks were of particular interest, namely modularity and nestedness.  Modularity emerges when groups of host species (or genotypes) shared groups of viruses.  Nestedness provided a view of the degree of specialization of both partners: high nestedness suggests that hosts differ in their susceptibility to infection, with some highly susceptible host genotypes selecting for very specialized viruses while strongly resistant host genotypes select for generalist viruses.  Translated to the plant pathology parlance, this extreme case would be equivalent to a gene-for-gene infection model (Flor 1956): new mutations confer hosts with resistance to recently evolved viruses while maintaining resistance to past viruses.  Likewise, virus mutations for expanding host range evolve without losing the ability to infect ancestral host genotypes.  By contrast, a non-nested network would represent a matching-allele infection model (Frank 2000) in which each interacting organism evolves by losing its capacity to resist/infect its ancestral partners, resembling a Red Queen dynamic.  Obviously, the reality is more complex and may lie anywhere between these two extreme situations.

Recently, Valverde et al. (2020) developed a model to explain the emergence of nestedness and modularity in plant-virus infection networks across diverse habitats.  They found that local modularity could coexist with global nestedness and that intraspecific competition was the main driver of the evolution of ecosystems in a continuum between nested-modular and nested networks.  These predictions were tested with field data showing the association between plant host species and different viruses in different agroecosystems (Valverde et al. 2020).  The effect of interspecific competition in the structure of empirical plant host-virus infection networks was also tested by McLeish et al. (2019).  Besides data from agroecosystems, evolution experiments have also shown the pervasive emergence of nestedness during the diversification of independently-evolved lineages of potyviruses in Arabidopsis thaliana genotypes that differ in their susceptibility to infection (Hillung et al. 2014; González et al. 2019; Navarro et al. 2020).

In their study, Moury et al. (2021) have expanded all these previous observations to a diverse set of pathosystems that range from viruses, bacteria, oomycetes, fungi, nematodes to insects.  While modularity was barely seen in only a few of the systems, nestedness was a common trend (observed in ~94% of all systems).  This nestedness, as seen in previous studies and as predicted by theory, emerged as a consequence of the existence of generalist and specialist strains of the parasites that differed in their capacity to infect more or less resistant plant genotypes.

As pointed out by Moury et al. (2021) in their conclusions, the ubiquity of nestedness in plant-parasite infection matrices has strong implications for the evolution and management of infectious diseases.

References

Flor, H. H. (1956). The complementary genic systems in flax and flax rust. In Advances in genetics, 8, 29-54. https://doi.org/10.1016/S0065-2660(08)60498-8

Flores, C. O., Meyer, J. R., Valverde, S., Farr, L., and Weitz, J. S. (2011). Statistical structure of host–phage interactions. Proceedings of the National Academy of Sciences, 108, E288-E297. https://doi.org/10.1073/pnas.1101595108

Frank, S. A. (2000). Specific and non-specific defense against parasitic attack. Journal of Theoretical Biology, 202, 283-304. https://doi.org/10.1006/jtbi.1999.1054

González, R., Butković, A., and Elena, S. F. (2019). Role of host genetic diversity for susceptibility-to-infection in the evolution of virulence of a plant virus. Virus evolution, 5(2), vez024. https://doi.org/10.1093/ve/vez052

Hillung, J., Cuevas, J. M., Valverde, S., and Elena, S. F. (2014). Experimental evolution of an emerging plant virus in host genotypes that differ in their susceptibility to infection. Evolution, 68, 2467-2480. https://doi.org/10.1111/evo.12458

McLeish, M., Sacristán, S., Fraile, A., and García-Arenal, F. (2019). Coinfection organizes epidemiological networks of viruses and hosts and reveals hubs of transmission. Phytopathology, 109, 1003-1010. https://doi.org/10.1094/PHYTO-08-18-0293-R

Moury B, Audergon J-M, Baudracco-Arnas S, Krima SB, Bertrand F, Boissot N, Buisson M, Caffier V, Cantet M, Chanéac S, Constant C, Delmotte F, Dogimont C, Doumayrou J, Fabre F, Fournet S, Grimault V, Jaunet T, Justafré I, Lefebvre V, Losdat D, Marcel TC, Montarry J, Morris CE, Omrani M, Paineau M, Perrot S, Pilet-Nayel M-L and Ruellan Y (2021) The quasi-universality of nestedness in the structure of quantitative plant-parasite interactions. bioRxiv, 2021.03.03.433745, ver. 4 recommended and peer-reviewed by PCI Evolutionary Biology. https://doi.org/10.1101/2021.03.03.433745

Navarro, R., Ambros, S., Martinez, F., Wu, B., Carrasco, J. L., and Elena, S. F. (2020). Defects in plant immunity modulate the rates and patterns of RNA virus evolution. bioRxiv. doi: https://doi.org/10.1101/2020.10.13.337402

Valverde, S., Vidiella, B., Montañez, R., Fraile, A., Sacristán, S., and García-Arenal, F. (2020). Coexistence of nestedness and modularity in host–pathogen infection networks. Nature ecology & evolution, 4, 568-577. https://doi.org/10.1038/s41559-020-1130-9

The observed levels of genomic diversity in contemporary populations are the result of changes imposed by several evolutionary processes. Among them, natural selection is known to dramatically shape the genetic diversity of loci associated with phenotypes which affect the fitness of carriers. As such, many efforts have been dedicated towards developing methods to detect signatures of natural selection from genomes of contemporary samples [1].
Recent technological advances made the generation of large-scale genomic data from temporal samples, either from experimental populations or historical or ancient samples, accessible to a wide scientific community [2]. Notably, temporal population genomic data allow for a direct observation and study of how, for instance, allele frequencies change through time in response to evolutionary stimuli. Such information can be exploited to detect loci under natural selection, either via mathematical modelling or by investigating empirical distributions [3].
However, most of current methods to detect selection from temporal genomic data have largely ignored selfing populations, despite the latter comprising a significant proportion of species with social and economic importance. Selfing changes genomic patterns by reducing the effective recombination rate, which makes distinguishing between neutral evolution and natural selection even more challenging than for the case of outcrossing populations [4]. Nevertheless, an outlier-approach based on temporal genomic data for the selfing Arabidopsis thaliana population revealed loci under selection [5].
This study suggested the promise of detecting selection for selfing populations and encouraged further investigations to test the power of selection scans under different mating systems.
To address this question, Navascués et al. [6] extended a previously proposed approach for temporal genome scan [7] to incorporate partial self-fertilization. In the original implementation [7], it is assumed that, under neutrality, all loci provide levels of genetic differentiation drawn from the same distribution. If some of the loci are under selection, such distribution should show heterogeneity. Navascués et al. [6] proposed a test for the homogeneity between loci-specific and genome-wide differentiation by deriving a null distribution of FST via simulations using SLiM [8]. After filtering for low-frequency variants and correct for multiple tests, authors derived a statistical test for selection and assess its power under a wide range of scenarios of selfing rate, selection coefficient, duration and type of selection [6].
The newly proposed test achieved good performance to distinguish between neutral and selected loci in most tested scenarios.
As expected, the test's performance significantly drops for scenarios of high selfing rates and selection from standing variation. Additionally, the probability to correctly detect selection decreases with increasing distance from the causal variant. Intriguingly, the test showed high power when the selected ancestral allele had an initial low frequency, and when the selected derived allele had a high initial frequency. When applied to a data set of around 1,000 SNPs from the highly selfing Medicago truncatula population, an annual plant of the legume family [9], the test did not provide any candidate loci under selection [6].
In summary, the detection of loci under selection in selfing populations is and largely remains a challenging task even when explictly account for the different mating system. However, recombination events that occurred before the selective pressure allow ancestral beneficial alleles to exhibit a detectable pattern of non-neutrality. As such, in partially selfing populations, the strength of the footprint of selection depends on several factors, mostly on the selfing rate, the time of onset and type of selection.
One major assumption of this study is that the model implies unstructured population and continuity between samples obtained from the same geographical location over time. As such assumptions are typically violated in real populations, further research into the effect of more complex demographic scenarios is desired to fully understand the power to detect selection in selfing populations. Furthermore, more power could be gained by including additional genomic information at each time point. In this context, recent approaches that make full use of genomic data based on deep learning [10] may contribute significantly towards this goal. Similarly, the effect of data filtering on the power to detect selection should be further explored, especially in the context of DNA resequencing experiments. These analyses will help elucidate the power offered by selection scans from temporal genomic data in selfing populations.

References

[1] Stern AJ, Nielsen R (2019) Detecting Natural Selection. In: Handbook of Statistical Genomics , pp. 397–40. John Wiley and Sons, Ltd. https://doi.org/10.1002/9781119487845.ch14
[2] Leonardi M, Librado P, Der Sarkissian C, Schubert M, Alfarhan AH, Alquraishi SA, Al-Rasheid KAS, Gamba C, Willerslev E, Orlando L (2017) Evolutionary Patterns and Processes: Lessons from Ancient DNA. Systematic Biology, 66, e1–e29. https://doi.org/10.1093/sysbio/syw059
[3] Dehasque M, Ávila‐Arcos MC, Díez‐del‐Molino D, Fumagalli M, Guschanski K, Lorenzen ED, Malaspinas A-S, Marques‐Bonet T, Martin MD, Murray GGR, Papadopulos AST, Therkildsen NO, Wegmann D, Dalén L, Foote AD (2020) Inference of natural selection from ancient DNA. Evolution Letters, 4, 94–108. https://doi.org/10.1002/evl3.165
[4] Vitalis R, Couvet D (2001) Two-locus identity probabilities and identity disequilibrium in a partially selfing subdivided population. Genetics Research, 77, 67–81. https://doi.org/10.1017/S0016672300004833
[5] Frachon L, Libourel C, Villoutreix R, Carrère S, Glorieux C, Huard-Chauveau C, Navascués M, Gay L, Vitalis R, Baron E, Amsellem L, Bouchez O, Vidal M, Le Corre V, Roby D, Bergelson J, Roux F (2017) Intermediate degrees of synergistic pleiotropy drive adaptive evolution in ecological time. Nature Ecology and Evolution, 1, 1551–1561. https://doi.org/10.1038/s41559-017-0297-1
[6] Navascués M, Becheler A, Gay L, Ronfort J, Loridon K, Vitalis R (2020) Power and limits of selection genome scans on temporal data from a selfing population. bioRxiv, 2020.05.06.080895, ver. 4 peer-reviewed and recommended by PCI Evol Biol. https://doi.org/10.1101/2020.05.06.080895
[7] Goldringer I, Bataillon T (2004) On the Distribution of Temporal Variations in Allele Frequency: Consequences for the Estimation of Effective Population Size and the Detection of Loci Undergoing Selection. Genetics, 168, 563–568. https://doi.org/10.1534/genetics.103.025908
[8] Messer PW (2013) SLiM: Simulating Evolution with Selection and Linkage. Genetics, 194, 1037–1039. https://doi.org/10.1534/genetics.113.152181
[9] Siol M, Prosperi JM, Bonnin I, Ronfort J (2008) How multilocus genotypic pattern helps to understand the history of selfing populations: a case study in Medicago truncatula. Heredity, 100, 517–525. https://doi.org/10.1038/hdy.2008.5
[10] Sanchez T, Cury J, Charpiat G, Jay F Deep learning for population size history inference: Design, comparison and combination with approximate Bayesian computation. Molecular Ecology Resources, n/a. https://doi.org/10.1111/1755-0998.13224