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

Perhaps because of its deadliness, the 2013-2016 Ebola Virus (EBOV) epidemics in West-Africa has led to unprecedented publication and sharing of full virus genome sequences. This was both rapid (90 full genomes were shared within weeks [1]) and important (more than 1500 full genomes have been released overall [2]). Furthermore, the availability of the metadata (especially GPS location) has led to depth analyses of the geographical spread of the epidemics [3].
In their work, Dellicour et al. [4] pursue earlier phylogeographical investigations in an original and yet simple approach to address questions of key public health importance. The originality of the approach is dual. First, from a technical standpoint, they capture the spread of infectious diseases in a continuous framework using a novel model that allows for rare long-distance dispersal events. Second, in a more classical discrete meta-population framework, they simulate the effect of public health interventions by pruning the phylogenetic tree and assessing how this affects key parameters. For instance, to simulate the effect of closing borders they remove subsets of the phylogeny that involved dispersal between countries and to simulate the effect of protecting a region by quarantine they remove all the leaves (i.e. the infections sampled) from this region. This phylogeny pruning is both original and simple. It is however limited because it currently assumes that policies are 100% effective and earlier modelling work on human influenza showed that long distance travel bans had to be implemented with >99% efficiency in order to slow epidemic growth from a time scale of days to weeks [5].
From a biological standpoint, Dellicour et al. [4] corroborate earlier findings that highly populated locations (>1,000,000 inhabitants) were crucial in explaining the magnitude of the epidemics but also show the importance of the transmission between the three capital cities. They also show that rare long-distance dispersing events of the virus are not key to explaining the magnitude of the epidemics (even though they assume 100% efficiency of suppressing long-distance event). Finally, thanks to their continuous model they estimate the speed of spread of the epidemics and are able to detect the effect of border closing on this speed.
Overall, this study [4], which involves state-of-the-art Bayesian inference methods of infection phylogenies using MCMC, stands out because of its effort to simulate public health interventions. It stands as an encouragement for the development of intervention models with increased realism and for even faster and larger virus sequence data sharing.

References

[1] Gire et al. 2014. Genomic surveillance elucidates Ebola virus origin and transmission during the 2014 outbreak. Science 345: 1369–1372. doi: 10.1126/science.1259657.
[2] Holmes EC, Dudas G, Rambaut A and Andersen KG. 2016. The evolution of Ebola virus: insights from the 2013-2016 epidemic. Nature 538: 193–200. doi: 10.1038/nature19790.
[3] Dudas et al. 2017. Virus genomes reveal factors that spread and sustained the Ebola epidemic. Nature 544: 309–315 (2017). doi: 10.1038/nature22040.
[4] Dellicour S, Baele G, Dudas G, Faria NR, Pybus OG, Suchard MA, Rambaud A and Lemey P. 2018. Phylodynamic assessment of intervention strategies for the West African Ebola virus outbreak. bioRxiv, 163691, ver. 3 peer-reviewed by Peer Community In Evolutionary Biology. doi: 10.1101/163691.
[5] Hollingsworth TD, Ferguson NM and Anderson RM. 2006. Will travel restrictions control the international spread of pandemic influenza? Nature Medicine 12, 497–499. doi: 10.1038/nm0506-497.

A classic textbook example in evolutionary ecology for phenotypic plasticity—the expression of different phenotypes by a genotype under different environmental conditions—is on Daphnia (Dodson 1989). If various species of this small crustacean are exposed to predation risk or cues thereof, their offspring show induced defense phenotypes including helmets, neck teeth, or head- and tail-spines. These induced morphological changes lower the risk of being eaten by a predator. As in Daphnia, induction can span over generations, while other induced phenotypic plastic changes are almost instantaneous, including many responses in physiology and behaviour (Gabriel et al. 2005). 

Larvae male bulb mites also show plasticity in morphologies throughout development. They can develop into a costly adult fighter morphology or a less costly, but vulnerable, scrambler type. The question Deere and Smallegange (Deere & Smallegange 2023) address is whether male bulb mite larvae can anticipate which type will likely be adaptive once they become adult, or alternatively, whether the resource availability or population density they experience during their larvae phase determines frequencies of adult scrambler and fighter types. They explore this question through experimental evolution, by removing different fractions of developing intermediate larva types. They thereby manipulate the stage structure of populations and alter selective forces on these stages. The potential shift of fixed genetics, imposed by the experimental selection regimes, is evaluated by fitness assays in green garden experiments.

The exciting extension to classical experiments on phenotypic plasticity is that the authors aim at exploring eco-evolutionary feedbacks experimentally in a system that is a little more complex than basic host-parasite or predator-prey systems. The latter involve, for instance, rotifer-algae dynamics (Yoshida et al. 2003; Becks et al. 2012) or similar simple lab systems for which eco-evolutionary feedbacks have been demonstrated. The challenge for the exploration of more complex systems is revealed in the study by Deere and Smallegange. Their findings suggest that the frequencies of adult male morphotype is triggered by the environmental condition (nutrient availability) during the larval phase, i.e. a simple environmental induced plastic response. No fixed genetic shift in determining adult morphotype frequencies occurs. The trigger at the larva phase remains also not perfectly determined in their experiments, as population density and resource (food) availability are partly confounded. Additional complexity and selective aspects come into play in this system, as the targeted developmental stages that develop into fighter male morphs are also dispersal morphs. If selection on dispersal to avoid residing in food-limited environments is strong, triggering genetic shifts in fighter morphs by local population structure might be hard to experimentally achieve. Small sample sizes limit conclusions on complex interactions among duration of the experiment, population size, developmental stage types, and adult fighter frequencies. The presented study (Deere & Smallegange 2023) helps to bridge theoretical predictions and empirical evidence for eco-evolutionary feedbacks that goes beyond simple ecological-driven changes in population dynamics (Govaert et al. 2019).  

References

Becks, L., Ellner, S.P., Jones, L.E. & Hairston, N.G. (2012). The functional genomics of an eco-evolutionary feedback loop: linking gene expression, trait evolution, and community dynamics. Ecol. Lett., 15, 492-501.
https://doi.org/10.1111/j.1461-0248.2012.01763.x
 
Deere, J.A. & Smallegange, I.M. (2023). Individual differences in developmental trajectory leave a male polyphenic signature in bulb mite populations. bioRxiv, 2023.02.06.527265, ver. 3 peer-reviewed and recommended by Peer Community in Evolutionary Biology.
https://doi.org/10.1101/2023.02.06.527265
 
Dodson, S. (1989). Predator-induced Reaction Norms. Bioscience, 39, 447-452.
https://doi.org/10.2307/1311136
 
Gabriel, W., Luttbeg, B., Sih, A. & Tollrian, R. (2005). Environmental tolerance, heterogeneity, and the evolution of reversible plastic responses. Am. Nat., 166, 339-53.
https://doi.org/10.1086/432558
 
Govaert, L., Fronhofer, E.A., Lion, S., Eizaguirre, C., Bonte, D., Egas, M., et al. (2019). Eco-evolutionary feedbacks-Theoretical models and perspectives. Funct. Ecol., 33, 13-30.
https://doi.org/10.1111/1365-2435.13241
 
Yoshida, T., Jones, L.E., Ellner, S.P., Fussmann, G.F. & Hairston, N.G. (2003). Rapid evolution drives ecological dynamics in a predator-prey system. Nat. 2003 4246946, 424, 303-306.
https://doi.org/10.1038/nature01767

Understanding factors explaining both intra and interspecific variation in susceptibility to infection by parasites remains a key question in evolutionary biology. Within a species variation in susceptibility is often explained by differences in behaviour affecting exposure to infection and/or resistance affecting the degree by which parasite growth is controlled (Roy & Kirchner, 2000, Behringer et al., 2000). This can vary between the sexes (Kelly et al., 2018) and may be explained by the ability of a parasite to attack different organs or tissues (Brierley et al., 2019). However, what goes on within one species is not always relevant to another, making it unclear when patterns can be scaled up and generalised across species. This is also important to understand when parasites may jump hosts, or identify species that may be susceptible to a host jump (Longdon et al., 2015). Phylogenetic distance between hosts is often an important factor explaining susceptibility to a particular parasite in plant and animal hosts (Gilbert & Webb, 2007, Faria et al., 2013). 

In two separate experiments, Roberts and Longdon (Roberts & Longdon, 2022) investigated how sex and tissue tropism affected variation in the load of Drosophila C Virus (DCV) across multiple Drosophila species. DCV load has been shown to correlate positively with mortality (Longdon et al., 2015). Overall, they found that load did not vary between the sexes; within a species males and females had similar DCV loads for 31 different species. There was some variation in levels of DCV growth in different tissue types, but these too were consistent across males for 7 species of Drosophila. Instead, in both experiments, host phylogeny or interspecific variation, explained differences in DCV load with some species being more infected than others. 

This study is neat in that it incorporates and explores simultaneously both intra and interspecific variation in infection-related life-history traits which is not often done (but see (Longdon et al., 2015, Imrie et al., 2021, Longdon et al., 2011, Johnson et al., 2012). Indeed, most studies to date explore either inter-specific differences in susceptibility to a parasite (it can or can’t infect a given species) (Davies & Pedersen, 2008, Pfenning-Butterworth et al., 2021) or intra-specific variability in infection-related traits (infectivity, resistance etc.) due to factors such as sex, genotype and environment (Vale et al., 2008, Lambrechts et al., 2006). This work thus advances on previous studies, while at the same time showing that sex differences in parasite load are not necessarily pervasive. 

References

Behringer DC, Butler MJ, Shields JD (2006) Avoidance of disease by social lobsters. Nature, 441, 421–421. https://doi.org/10.1038/441421a

Brierley L, Pedersen AB, Woolhouse MEJ (2019) Tissue tropism and transmission ecology predict virulence of human RNA viruses. PLOS Biology, 17, e3000206. https://doi.org/10.1371/journal.pbio.3000206

Davies TJ, Pedersen AB (2008) Phylogeny and geography predict pathogen community similarity in wild primates and humans. Proceedings of the Royal Society B: Biological Sciences, 275, 1695–1701. https://doi.org/10.1098/rspb.2008.0284

Faria NR, Suchard MA, Rambaut A, Streicker DG, Lemey P (2013) Simultaneously reconstructing viral cross-species transmission history and identifying the underlying constraints. Philosophical Transactions of the Royal Society B: Biological Sciences, 368, 20120196. https://doi.org/10.1098/rstb.2012.0196

Gilbert GS, Webb CO (2007) Phylogenetic signal in plant pathogen–host range. Proceedings of the National Academy of Sciences, 104, 4979–4983. https://doi.org/10.1073/pnas.0607968104

Imrie RM, Roberts KE, Longdon B (2021) Between virus correlations in the outcome of infection across host species: Evidence of virus by host species interactions. Evolution Letters, 5, 472–483. https://doi.org/10.1002/evl3.247

Johnson PTJ, Rohr JR, Hoverman JT, Kellermanns E, Bowerman J, Lunde KB (2012) Living fast and dying of infection: host life history drives interspecific variation in infection and disease risk. Ecology Letters, 15, 235–242. https://doi.org/10.1111/j.1461-0248.2011.01730.x

Kelly CD, Stoehr AM, Nunn C, Smyth KN, Prokop ZM (2018) Sexual dimorphism in immunity across animals: a meta-analysis. Ecology Letters, 21, 1885–1894. https://doi.org/10.1111/ele.13164

Lambrechts L, Chavatte J-M, Snounou G, Koella JC (2006) Environmental influence on the genetic basis of mosquito resistance to malaria parasites. Proceedings of the Royal Society B: Biological Sciences, 273, 1501–1506. https://doi.org/10.1098/rspb.2006.3483

Longdon B, Hadfield JD, Day JP, Smith SCL, McGonigle JE, Cogni R, Cao C, Jiggins FM (2015) The Causes and Consequences of Changes in Virulence following Pathogen Host Shifts. PLOS Pathogens, 11, e1004728. https://doi.org/10.1371/journal.ppat.1004728

Longdon B, Hadfield JD, Webster CL, Obbard DJ, Jiggins FM (2011) Host Phylogeny Determines Viral Persistence and Replication in Novel Hosts. PLOS Pathogens, 7, e1002260. https://doi.org/10.1371/journal.ppat.1002260

Pfenning-Butterworth AC, Davies TJ, Cressler CE (2021) Identifying co-phylogenetic hotspots for zoonotic disease. Philosophical Transactions of the Royal Society B: Biological Sciences, 376, 20200363. https://doi.org/10.1098/rstb.2020.0363

Roberts KE, Longdon B (2023) Heterogeneities in infection outcomes across species: examining sex and tissue differences in virus susceptibility. bioRxiv 2022.11.01.514663, ver. 2 peer-reviewed and recommended by Peer Community in Evolutionary Biology. https://doi.org/10.1101/2022.11.01.514663 

Roy BA, Kirchner JW (2000) Evolutionary Dynamics of Pathogen Resistance and Tolerance. Evolution, 54, 51–63. https://doi.org/10.1111/j.0014-3820.2000.tb00007.x

Vale PF, Stjernman M, Little TJ (2008) Temperature-dependent costs of parasitism and maintenance of polymorphism under genotype-by-environment interactions. Journal of Evolutionary Biology, 21, 1418–1427. https://doi.org/10.1111/j.1420-9101.2008.01555.x

Awareness has been growing that structural variants in the genome of species play a fundamental role in adaptive evolution and diversification [1]. Here, Le Moan and co-authors [2] report empirical genomic-wide SNP data on the European plaice (Pleuronectes platessa) across a major environmental transmission zone, ranging from the North Sea to the Baltic Sea. Regions of high linkage disequilibrium suggest the presence of two structural variants that appear to have evolved 220 kya. These two putative structural variants show weak signatures of isolation by distance when contrasted against the rest of the genome, but the frequency of the different putative structural variants appears to co-vary in some parts of the studied range with the environment, indicating the involvement of both selective and neutral processes. This study adds to the mounting body of evidence that structural genomic variants harbour significant information that allows species to respond and adapt to the local environmental context.

References

[1] Wellenreuther, M., Mérot, C., Berdan, E., & Bernatchez, L. (2019). Going beyond SNPs: the role of structural genomic variants in adaptive evolution and species diversification. Molecular ecology, 28(6), 1203-1209. doi: 10.1111/mec.15066
[2] Le Moan, A. Bekkevold, D. & Hemmer-Hansen J. (2020). Evolution at two time-frames: ancient and common origin of two structural variants involved in local adaptation of the European plaice (Pleuronectes platessa). bioRxiv, 662577, ver. 5 peer-reviewed and recommended by PCI Evol Biol. doi: 10.1101/662577