The preprint by Bertels et al. [1] reports an interesting application of the well-accepted idea that positively selected traits (here variants) can appear several times independently; think about the textbook examples of flight capacity. Hence, the authors assume that reciprocally convergence implies positive selection. The methodology becomes then, in principle, straightforward as one can simply count variants in independent datasets to detect convergent mutations.
In this preprint, the authors have applied this counting strategy on 95 available sequence alignments of the env gene of HIV-1 [2,3] that corresponds to samples taken in different patients during the early phase of infection, at the very beginning of the onset of the immune system. They have compared the number and nature of the convergent mutations to a "neutral" model that assumes (a) a uniform distribution of mutations and (b) a substitution matrix estimated from the data. They show that there is an excess of convergent mutations when compared to the “neutral” expectations, especially for mutations that have arisen in 4+ patients. They also show that the gp41 gene is enriched in these convergent mutations. The authors then discuss in length the potential artifacts that could have given rise to the observed pattern.
I think that this preprint is remarkable in the proposed methodology. Samples are taken in different individuals, whose viral populations were founded by a single particle. Thus, there is no need for phylogenetic reconstruction of ancestral states that is the typical first step of trait convergent analyses. It simply becomes counting variants. This simple counting procedure needs nonetheless to be compared to a “neutral” expectation (a reference model), which includes the mutational process. In this article, the poor predictions of a specifically designed reference model is interpreted as an evidence for positive selection.
Whether the few mutations that are convergent in 4-7 samples out of 95 were selected or not is hard to assess with certainty. The authors have provided good evidence that they are, but only experimental validations will strengthen the claim. Nonetheless, beyond a definitive clue to the implication of selection on these particular mutations, I found the methodological strategy and the discussions on the potential biases highly stimulating. This article is an excellent starting point for further methodological developments that could be then followed by large-scale analyses of convergence in many different organisms and case studies.

References

[1] Bertels, F., Metzner, K. J., & Regoes R. R. (2018). Convergent evolution as an indicator for selection during acute HIV-1 infection. BioRxiv, 168260, ver. 4 peer-reviewed and recommended by PCI Evol Biol. doi: 10.1101/168260
[2] Keele, B. F., Giorgi, E. E., Salazar-Gonzalez, J. F., Decker, J. M., Pham, K.T., Salazar, M. G., Sun, C., Grayson, T., Wang, S., Li, H. et al. (2008). Identification and characterization of transmitted and early founder virus envelopes in primary HIV-1 infection. Proc Natl Acad Sci USA 105: 7552–7557. doi: 10.1073/pnas.0802203105
[3] Li, H., Bar, K. J., Wang, S., Decker, J. M., Chen, Y., Sun, C., Salazar-Gonzalez, J.F., Salazar, M.G., Learn, G.H., Morgan, C. J. et al. (2010). High multiplicity infection by HIV-1 in men who have sex with men. PLoS Pathogens 6:e1000890. doi: 10.1371/journal.ppat.1000890

Sexual selection occurs in flowering plants too. However it tends to be understudied in comparison to animal sexual selection, in part because the minuscule size and long dispersal distances of the individuals producing male gametes (pollen grains) seriously complicate the estimation of male siring success and thereby the measurement of sexual selection. Dorken and Perry [1] introduce a novel and clever approach to estimate sexual selection in plants, which bypasses the need for a direct quantification of absolute male mating success. This approach builds on the fact that the strength of sexual selection is directly related to the ability of individuals to monopolize mates [2]. In plants, mate monopolization can be assessed by examining the proportion of seeds produced by a given plant that are full-sibs, i.e. that share the same father. A nice feature of this proportion of full-sib seeds per maternal parent is it equals the coefficient of correlated paternity of Ritland [3], which can be readily obtained from the hundreds of plant mating system studies using genetic markers. A less desirable feature of the proportion of full sibs per maternal plant is that it is inversely related to population size, an effect that should be corrected for. The resulting index of mate monopolization is a simple product: (coefficient of correlated paternity)x(population size – 1). The authors test whether their index of mate monopolization is a good correlate of sexual selection, measured more traditionally as the selection differential on a trait influencing mating success, using a combination of theoretical and experimental approaches. Both approaches confirm that the two quantities are positively correlated, which suggests that the index of mate monopolization could be a convenient way to estimate the relative strength of sexual selection in flowering plants. These results call for further investigation, e.g. to verify that the effect of population size is well controlled for, or to assess the effects of non-random mating and inbreeding depression; however, this work paves the way for an expansion of sexual selection studies in flowering plants.

References

[1] Dorken ME and Perry LE. 2017. Correlated paternity measures mate monopolization and scales with the magnitude of sexual selection. Journal of Evolutionary Biology 30: 377-387 doi: 10.1111/jeb.13013

[2] Klug H, Heuschele J, Jennions M and Kokko H. 2010. The mismeasurement of sexual selection. Journal of Evolutionary Biology 23:447-462. doi: 10.1111/j.1420-9101.2009.01921.x

[3] Ritland K. 1989. Correlated matings in the partial selfer Mimulus guttatus. Evolution 43:848-859. doi: 10.2307/2409312

Understanding how populations evolve under thermal stress and how this process shapes the response of other stress responses is an important research topic in the context of thermal adaptation and climate change. In a thermal experimental evolution study in Drosophila subobscura, Castañeda (2024) addressed the correlated responses to selection for increasing knockdown temperature in different resistance traits, either directly related to thermal stress (e.g. knockdown time at different temperatures and CTmax) or not (e.g. desiccation and starvation resistance). 

The author found that the evolution of higher knockdown temperature did in fact lead to correlated responses in other stress traits. While such correlations might be expected for the thermal stress traits measured (knockdown time and CTmax), it was perhaps less expectable for desiccation and starvation resistance. However, the general occurrence of correlated evolutionary responses between stressors has been previously described, namely in Drosophila (e.g. see Bubliy and Loeschcke 2005), pointing to a possible genetic link between distinct (thermal) stress traits. 

There are however some features that make the findings of this study rather appealing. First, the evidence that the correlated stress responses depend on the intensity of thermal selection (i.e. the warming rate) and on the sex of the organisms. Second, correlated patterns of both desiccation and starvation resistance highlight the possibility of the evolution of a cross-tolerance response, which might positively impact on population ability to evolve under sustained stressful environments (Rodgers and Gomez Izasa 2023). However, it is important to point out that the correlated patterns between these two resistance traits (desiccation and starvation) were not exactly consistent. In fact, the negative correlated response observed for female starvation resistance is thought provoking and argues again a general scenario of cross-tolerance. 

While these findings are a step forward for a more multifaceted understanding of thermal adaptation in the context of stressful environments, they also highlight the need for further studies of thermal adaptation namely 1) addressing the underlying physiological and genomic mechanisms that link male and female heat tolerance and the response to other stress resistance traits (namely starvation resistance); 2) testing the extent to which cross-resistance patterns can be generalized to different thermal selection contexts and populations. 

In addition, this study also opens new questions considering the scope of correlated evolution to other stress traits, that might be relevant in diverse ecological scenarios. For instance, does selection towards higher heat resistance lead to correlated evolution of cold resistance? And under which circumstances (e.g. different heat selection intensities)?  In fact, the occurrence of a positive (or negative) correlation cold and heat stress responses is a topic of high interest, with relevant ecological implications particularly considering the increased thermal fluctuations in natural environments because of climate warming. Cross-tolerance between cold and heat stress responses has been described (Singh 2022, Rodgers and Gomez Izasa 2023). On the other hand, negative correlations (i.e. trade-offs) between these stress traits (Stazione et al. 2020; Schou et al 2022) can impact negatively on populations’ ability to withstand thermal variability. 

As climatic changes proceed leading to increasing environmental variability, empirical studies such as that of Castañeda (2024) are critical in the pursue for a multivariate perspective on trait evolution in scenarios of climate change adaptation. Understanding how tolerance to different environmental stressors may evolve and which factors can act as drivers of that variation will ultimately enable better forecasts of climate change effects on biodiversity in nature.

References

Castañeda, LE. Cross-tolerance evolution is driven by selection on heat tolerance in Drosophila subobscura. Biorxiv, ver. 4 peer-reviewed and recommended by Peer Community in Evolutionary Biology (2024). https://www.doi.org/10.1101/2023.09.05.556367

Bubliy, OA, Loeschcke, V. Correlated responses to selection for stress resistance and longevity in a laboratory population of Drosophila melanogaster. J Evol Biol. 18(4):789-803 (2005). https://www.doi.org/10.1111/j.1420-9101.2005.00928.x. 

Rodgers, EM, Gomez Isaza, DF. The mechanistic basis and adaptive significance of cross-tolerance: a 'pre-adaptation' to a changing world? J Exp Biol. 226(11):jeb245644 (2023). https://www.doi.org/10.1242/jeb.245644. 

Schou, MF, Engelbrecht, A, Brand, Z, Svensson, EI, Cloete, S, Cornwallis, CK. Evolutionary trade-offs between heat and cold tolerance limit responses to fluctuating climates. Sci Adv. 8(21):eabn9580 (2022). https://www.doi.org/10.1126/sciadv.abn9580. 

Singh, K, Arun Samant, M, Prasad, NG. Evolution of cross-tolerance in Drosophila melanogaster as a result of increased resistance to cold stress. Sci Rep. 12(1):19536 (2022). https://www.doi.org/10.1038/s41598-022-23674-z.

Stazione, L, Norry, FM, Gomez, FH, Sambucetti, P. Heat knockdown resistance and chill-coma recovery as correlated responses to selection on mating success at high temperature in Drosophila buzzatii. Ecol Evol. 10(4):1998-2006 (2020). https://www.doi.org/10.1002/ece3.6032.

Genetic markers are used for in modern population genetics/genomics to uncover the past neutral and selective history of population and species. Besides Single Nucleotide Polymorphisms (SNPs) obtained from whole genome data, microsatellites (or Short Tandem Repeats, SSR) have been common markers of choice in numerous population genetics studies of non-model species with large sample sizes [1]. Microsatellites can be used to uncover and draw inference of the past population demography (e.g. expansion, decline, bottlenecks…), population split, population structure and gene flow, but also life history traits and modes of reproduction (e.g. [2,3]). These markers are widely used in conservation genetics [4] or to study parasites or disease vectors [5]. Microsatellites do show higher mutation rate than SNPs increasing, on the one hand, the statistical power to infer recent events (for example crop domestication, [2,3]), while, on the other hand, decreasing their statistical power over longer time scales due to homoplasy [6].
To perform such analyses, however, an excellent and reliable quality of data is required. As emphasized in the article by De Meeûs et al. [7] three main issues do bias the observed heterozygosity at microsatellites: null alleles, short allele dominance (SAD) and stuttering. These originates from poor PCR amplification. As a result, an excess of homozygosity is observed at the microsatellite loci leading to overestimation of the variation statistics FIS and FST as well as increased linage disequilibrium (LD). For null alleles, several methods and software do help to reduce the bias, and in the present study, De Meeûs et al. [7] propose a way to tackle issues with SAD and stuttering.
The authors study a dataset consisting of 387 samples from 61 subsamples genotyped at nine loci of the species Ixodes scapularis, i.e. ticks transmitting the Lyme disease. Based on correlation methods and FST, FIS they can uncover null alleles and SAD. Stuttering is detected by evaluating the heterozygote deficit between alleles displaying a single repeat difference. Without correction, six loci are affected by one of these amplification problems generating a large deficit of heterozygotes (measured by significant FIS and FST) remaining so after correction for the false discovery rate (FDR). These results would be classically interpreted as a strong Wahlund effect and/or selection at several loci.
After correcting for null alleles, the authors apply two novel corrections: 1) a re-examination of the chromatograms reveals previously disregarded larger alleles thus decreasing SAD, and 2) pooling alleles close in size decreasing stuttering. The corrected dataset shows then a significant excess of heterozygotes as could be expected in a dioecious species with strong population structure. The FDR correction removes then the significant excess of homozygotes and LD between pairs of loci. FST on the cured dataset is used to demonstrate the strong population structure and small effective subpopulation sizes. This is confirmed by a clustering analysis using discriminant analysis of principal components (DAPC).
While based on a specific dataset of ticks from different populations sampled across the USA, the generality of the authors’ approach is presented in Figure 6 in which they provide a step by step flowchart to cure microsatellite datasets from null alleles, SAD and stuttering. Several criteria based on FIS, FST and LD between loci are used as decision keys in the flowchart. An excel file is also provided as help for the curation steps. This study and the proposed methodology are thus extremely useful for all population geneticists working on non-model species with large number of samples genotyped at microsatellite markers. The method not only allows more accurate estimates of heterozygosity but also prevents the thinning of datasets due to the removal of problematic loci. As a follow-up and extension of this work, an exhaustive simulation study could investigate the influence of these data quality issues on past demographic and population structure inference under a wide range of scenarios. This would allow to quantify the current biases in the literature and the robustness of the methodology devised by De Meeûs et al. [7].

References

[1] Jarne, P., and Lagoda, P. J. (1996). Microsatellites, from molecules to populations and back. Trends in ecology & evolution, 11(10), 424-429. doi: 10.1016/0169-5347(96)10049-5
[2] Cornille, A., Giraud, T., Bellard, C., Tellier, A., Le Cam, B., Smulders, M. J. M., Kleinschmit, J., Roldan-Ruiz, I. and Gladieux, P. (2013). Postglacial recolonization history of the E uropean crabapple (Malus sylvestris M ill.), a wild contributor to the domesticated apple. Molecular Ecology, 22(8), 2249-2263. doi: 10.1111/mec.12231
[3] Parat, F., Schwertfirm, G., Rudolph, U., Miedaner, T., Korzun, V., Bauer, E., Schön C.-C. and Tellier, A. (2016). Geography and end use drive the diversification of worldwide winter rye populations. Molecular ecology, 25(2), 500-514. doi: 10.1111/mec.13495
[4] Broquet, T., Ménard, N., & Petit, E. (2007). Noninvasive population genetics: a review of sample source, diet, fragment length and microsatellite motif effects on amplification success and genotyping error rates. Conservation Genetics, 8(1), 249-260. doi: 10.1007/s10592-006-9146-5
[5] Koffi, M., De Meeûs, T., Séré, M., Bucheton, B., Simo, G., Njiokou, F., Salim, B., Kaboré, J., MacLeod, A., Camara, M., Solano, P., Belem, A. M. G. and Jamonneau, V. (2015). Population genetics and reproductive strategies of African trypanosomes: revisiting available published data. PLoS neglected tropical diseases, 9(10), e0003985. doi: 10.1371/journal.pntd.0003985
[6] Estoup, A., Jarne, P., & Cornuet, J. M. (2002). Homoplasy and mutation model at microsatellite loci and their consequences for population genetics analysis. Molecular ecology, 11(9), 1591-1604. doi: 10.1046/j.1365-294X.2002.01576.x
[7] De Meeûs, T., Chan, C. T., Ludwig, J. M., Tsao, J. I., Patel, J., Bhagatwala, J., and Beati, L. (2019). Deceptive combined effects of short allele dominance and stuttering: an example with Ixodes scapularis, the main vector of Lyme disease in the USA. bioRxiv, 622373, ver. 4 peer-reviewed and recommended by Peer Community In Evolutionary Biology. doi: 10.1101/622373

Interspecific matings are by definition rare events in nature, but when they occur they can be very important, and not only because they might condition gene flow between species. Even when such matings have no genetic consequence, for instance if they do not yield any fertile hybrid offspring, they can still have an impact on the population dynamics of the species involved [1]. Such atypical pairings between heterospecific partners are usually regarded as detrimental or undesired; as they interfere with the occurrence or success of intraspecific matings, they are expected to cause a decline in absolute fitness.
The story is not always so simple however, and it might all depend on the timing of events and on the identity of the partners. Using the herbivorous mite Tetranychus urticae as a model, Clemente et al. [2] experimentally arranged matings with two other Tetranychus species that commonly share the same host plants as T. urticae. They carefully controlled the history of events: heterospecific matings could occur just before, just after, 24h before, or 24h after, a conspecific mating. Interestingly, the oviposition rate (total fecundity) of females was increased when mating with a heterospecific individual. This suggests that heterospecic sperm can stimulate oogenesis just as conspecific sperm does. Such a positive effect was observed for matings involving T. ludeni females and T. urticae males, but a negative effect is found in the interaction with T. evansi. Sex-ratio (fertilization success in those species) could also be impacted but, unlike fertilization, this occurred when the mating events were distant in time. This is is at odds with what is observed in conspecific matings, where sperm displacement occurs only if mating events are temporally close. Overall, the effects of heterospecific mating were quite variable and it is challenging to predict a single, general, effect of interspecific matings. The net effect will likely be context-dependent, depending on the relative frequency of the difference mating sequences and on how fecundity and sex-ratio contribute to overall fitness, both aspect strongly influenced by the population dynamics and structure.

References

[1] Gröning J. & Hochkirch A. 2008. Reproductive interference between animal species. The Quarterly Review of Biology 83: 257-282. doi: 10.1086/590510

[2] Clemente SH, Santos I, Ponce AR, Rodrigues LR, Varela SAM & Magalhaes S. 2017 Despite reproductive interference, the net outcome of reproductive interactions among spider mite species is not necessarily costly. bioRxiv 113274, ver. 4 of the 30th of June 2017. doi: 10.1101/113274