Insecticide resistance in mosquitoes represents a notable challenge to public health efforts aimed at controlling vector-borne diseases. Among mosquito species, Aedes aegypti is particularly significant due to its extensive geographic spread and its ability to transmit arboviruses causing diseases such as dengue, yellow fever, Zika, and chikungunya (Brown et al., 2014). Insecticide resistance typically develops through two main mechanisms: target-site mutations, which affect the insecticide's interaction with its target, and metabolic resistance, in which insecticide detoxification is enhanced in mosquitoes. While target-site mutations are well characterized, the mechanisms underlying metabolic resistance are understudied. 

The study by Bacot and colleagues (2024) contributes to our understanding of the genetic and evolutionary mechanisms driving insecticide resistance, focusing on a case of metabolic resistance in Aedes aegypti from French Guiana. Following the recent identification of a copy number variant region on chromosome 1, potentially linked to overexpression of detoxification enzymes (Cattel et al., 2020), this study explores the region’s genomic architecture, its likely origin and provides compelling evidence for its role in insecticide resistance.

Through RNA sequencing and whole-genome pool sequencing, the authors reveal that this 220 kilobase duplication increases the expression level of several clustered P450 genes. Cytochrome P450s are known to play a role in breaking down pyrethroids like deltamethrin, a commonly used insecticide. The role of P450 enzymes in detoxification was demonstrated by treating mosquitoes with piperonyl butoxide, a P450 enzyme inhibitor, and observing reduction in deltamethrin resistance, further confirmed by RNA interference experiments. Despite the clear advantages of this genomic duplication in conferring resistance, the study also uncovers a fitness cost associated with carrying the duplication. Through experimental evolution, the authors find that mosquitoes with the duplication experience reduced fitness in the absence of insecticide pressure. Given the regions structural complexity, the authors could not completely disassociate the effect of the duplicated region and that of a target-site mutation. However, they developed an assay that can accurately track the presence of this resistance allele in mosquito populations. 

Altogether, the study by Bacot et al. (2024) highlights the challenges of characterizing the phenotypic effect of copy number variant regions in complex genomes, such as that of Aedes aegypti. It emphasizes the need for further studies on the origin and spread of this duplication to better understand how similar resistance mechanisms might evolve and disseminate. Overall, the completeness and coherence of the narrative, the detailed and thorough analysis, and the insightful discussion, make this work not only a significant contribution to the field of insecticide resistance but an interesting read for the general evolutionary biology community.   

References

Brown, J. E., Evans, B. R., Zheng, W., Obas, V., Barrera-Martinez, L., Egizi, A., Zhao, H., Caccone, A., & Powell, J. R. (2014). Human impacts have shaped historical and recent evolution in Aedes aegypti, the dengue and yellow fever mosquito. Evolution, 68(2), 514–525. https://doi.org/10.1111/evo.12281

Cattel, J., Faucon, F., Le Péron, B., Sherpa, S., Monchal, M., Grillet, L., Gaude, T., Laporte, F., Dusfour, I., Reynaud, S., & David, J. P. (2019). Combining genetic crosses and pool targeted DNA-seq for untangling genomic variations associated with resistance to multiple insecticides in the mosquito Aedes aegypti. Evolutionary applications, 13(2), 303–317. https://doi.org/10.1111/eva.12867

Tiphaine Bacot, Chloe Haberkorn, Joseph Guilliet, Julien Cattel, Mary Kefi, Louis Nadalin, Jonathan Filee, Frederic Boyer, Thierry Gaude, Frederic Laporte, Jordan Tutagata, John Vontas, Isabelle Dusfour, Jean-Marc Bonneville, Jean-Philippe David (2024) A genomic duplication spanning multiple P450s contributes to insecticide resistance in the dengue mosquito Aedes aegypti. bioRxiv, ver.5 peer-reviewed and recommended by PCI Evol Biol https://doi.org/10.1101/2024.04.03.587871

A comprehensive study by Roux et al 2024 investigates the impact of eusociality on the efficacy of natural selection in termites, with and additional focus on whether higher levels of social complexity are associated with lower effective population size (Ne) and relaxed purifying selection.

Eusociality is characterized by a division of reproductive labor, cooperative care of offspring, and overlapping generations, and has evolved independently across various animal taxa with the most complex social systems found in Hymenoptera (bees, wasps, ants) and termites. Because reproduction is limited to a few individuals, this leads to a reduced effective population size (Ne), which impacts genome evolution. Smaller Ne increases the influence of genetic drift, weakening the efficiency of natural selection and allowing the accumulation of weakly deleterious mutations. This phenomenon, known as the "drift barrier," alters the mutation-selection balance in eusocial organisms.

Studies in a range of social arthropods including ants, termites, crustaceans and spiders have shown an elevated ratio of nonsynonymous to synonymous substitutions (dN/dS), indicating relaxed purifying selection due to small Ne. In termite species with complex social structures, such as those with large colonies and high caste specialization, there are reports of higher dN/dS ratios compared to simpler social species. This suggests that higher social complexity, reflected in traits like nesting strategies and developmental pathways, further reduces Ne and the effectiveness of natural selection. 

The authors address these hypotheses by exploring the genomic impact of eusociality in termites (Isoptera) taking two approaches: First, they analyze transcriptome data from 66 Blattodea species and calculate the ratio of non-synonymous to synonymous mutations (dN/dS) as an indicator of natural selection efficiency and effective population size. They analyses reveal an increased dN/dS ratio in termites compared to other Blattodea species, reinforcing the notion that convergent evolution toward eusociality significantly reduces effective population size and weakens natural selection efficiency across the genome. Additionally, a comparison of 68 termite transcriptomes shows that this effect is more pronounced in species with higher social complexity. This is exciting as it advances our understanding of how increasing complexity in social organization decreases Ne and the efficiency of natural selection. The study substantiates the notion that social transitions follow evolutionary trajectories where lower and Ne and increasing drift have negative consequences for genome evolution (Ma et al 2024).

References

Camille Roux, Alice Ha, Arthur Weyna, Morgan Lode, Jonathan Romiguier (2024) The impact of social complexity on the efficacy of natural selection in termites. bioRxiv, ver.2 peer-reviewed and recommended by PCI Evol Biol. https://doi.org/10.1101/2024.04.26.591327

Jilong Ma, Jesper Bechsgaard, Anne Aagaard, Palle Villesen, Trine Bilde, Mikkel Heide Schierup (2024) Sociality in spiders is an evolutionary dead-end. 
bioRxiv 2024.04.22.590577. https://doi.org/10.1101/2024.04.22.590577

The introduction of Ludwigia species as ornamental plants in both North America and Europe dates back almost two centuries, during which time they expanded as naturalized and later invasive species in these territories (Dandelot et al. 2005, Okada et al. 2009). Repeated deliberate or non-deliberate introductions over time of this species complex that can hybridize has given rise to an evolutionarily complex scenario, which is compounded by the difficulty in delimiting some of these species and by the diversity of their modes of reproduction.

Dandelot (2004) and Dandelot et al. (2005) determined the presence of two Ludwigia taxa in France, L. peploides subsp. montevidensis (Spreng.) P.H.Raven (here after Lpm), and L. grandiflora subsp. hexapetala (Hook. & Arn.) G.L.Nesom & Kartesz (here after Lgh), based on their cytotypes (2n = 16, and 2n = 80, respectively) and without evidence of hybridization between them. Furthermore, despite a predominantly vegetative reproduction observed for both species, they differed in their breeding systems. While Lpm is self-compatible and produce a high number of viable seeds in all populations, Lgh is self-incompatible and its populations may drastically differ in seed viability (Dandelot 2004). Several years later, Portillo-Lemus et al. (2021) determined that the differences in seed production between some populations of Lgh are due to the existence of a heteromorphic reproductive system in this taxon, involving a self-incompatible morph (long-style morph; hereafter L-morph), and a self-compatible morph (short-style morph: hereafter S-morph). Moreover, Portillo-Lemus et al. (2022) observed that self-pollen in the L-morph flowers stop growing lately (i.e., in the ovarian area) without fertilizing the ovules, concluding that a late-acting self-incompatible system (hereafter, LSI) is present in this morph.

At this point, it is relevant to understand the possible interactions between populations of different morphs in Lgh, and the implications that they may have on their expansive success in non-native areas in order to develop more effective management plans. To achieve this goal, Stoeckel et al. (2024) analyzed the population genetics in 53 Lgh populations in western Europe (without finding any Lpm population in the sampling area), 40 of which exclusively presented the L-morph and 13 the S-morph. This fact offered the opportunity to compare and interpret the differences between populations of different morphs in Lgh. Other previous works on genetic diversity of Lgh in peripatric or non-native areas pointed to a high clonality and an extremely low genetic diversity (Okada et al. 2009, Armitage et al. 2013), concluding in a monoclonal or few ancestral original clones for these invasive populations.

However, the investigations of Stoeckel et al. (2024) found a high genetic diversity in all populations of Lgh studied despite their predominant clonal reproduction. Interestingly, they found that sexual reproduction is also present, not only in the S-morph by selfing, but also in the L-morph, although limited and preferably by allogamy. They discuss the advantages and drawbacks of the different modes of reproduction observed in Lgh populations, the interactions among them, and the implications that both, the scarcely documented LSI (Gibbs 2014) and selfing, have in the reproductive success and in the maintenance of the high genetic diversity observed in Lgh in western Europe.

The contrasting results with the previous ones (Okada et al. 2009, Armitage et al. 2013) stress the relevance of using appropriate markers and analyses to assess the genetic diversity in autoployploid species, as well as the necessity of knowing the modes of reproduction in the populations studied for an optimal interpretation of the genetic metrics. The approach of the study by Stoeckel et al. (2024) had the challenge of having found suitable markers to deal with a taxon of complex origin such as Lgh, whose genome is compound by a set of autotetraploid chromosomes shared with Lpm and traces of ancient hybridizations of other diploid lineages (Barloy et al. 2024). Using RAD-Seq, Stoeckel et al. (2024) generated an original set of 36 polymorphic SNPs shared between Lgh and Lpm ensuring that these SNPs belong to the tetraploid part of the Lgh genome derived from Lpm. Another interesting contribution of this work is the exhaustive analysis of several genetic descriptors (indexes) and the interpretative guide they provide for each of them in relation to the different modes of reproduction of the study system. Finally, they propose a pair of very useful synthetic indices (i.e., clonality index and selfing index), since they allow to classify populations according to their levels of clonality and selfing.

Stoeckel et al. (2024) conclude the relevance that selfing and LSI populations, and the hybridization between them may have on the expansion and success of invasive plant species, and the necessity to know the modes of reproduction of these populations jointly with their genetic diversity in order to develop appropriate management plans. This study raises new questions such as the modes of reproduction and genetic diversity and structure have other Lgh populations, both invasive and native, and the dynamics of these populations under different future scenarios.

References

Armitage, J. D., Könyves, K., Bailey, J. P., David, J. C., & Culham, A. (2013). A molecular, morphological and cytological investigation of the identity of non-native Ludwigia (Onagraceae) populations in Britain. New Journal of Botany, 3(2), 88–95. https://doi.org/10.1179/2042349713Y.0000000023

Barloy, D., Lemus, L. P.-, Krueger-Hadfield, S. A., Huteau, V., & Coriton, O. (2024). Genomic relationships among diploid and polyploid species of the genus Ludwigia L. section Jussiaea using a combination of cytogenetic, morphological, and crossing investigations. ver. 4 peer-reviewed and recommended by Peer Community in Evolutionary Biology. https://doi.org/10.1101/2023.01.02.522458

Dandelot, S. (2004). Les Ludwigia spp. invasives du Sud de la France: Historique, Biosystématique, Biologie et Ecologie [PhD thesis, Aix-Marseille 3]. https://www.theses.fr/2004AIX30052

Dandelot, S., Verlaque, R., Dutartre, A., & Cazaubon, A. (2005). Ecological, dynamic and taxonomic problems due to Ludwigia (Onagraceae) in France. Hydrobiologia, 551(1), 131–136. https://doi.org/10.1007/s10750-005-4455-0

Gibbs, P. E. (2014). Late-acting self-incompatibility – the pariah breeding system in flowering plants. New Phytologist, 203(3), 717–734. https://doi.org/10.1111/nph.12874

Okada, M., Grewell, B. J., & Jasieniuk, M. (2009). Clonal spread of invasive Ludwigia hexapetala and L. grandiflora in freshwater wetlands of California. Aquatic Botany, 91(3), 123–129. https://doi.org/10.1016/j.aquabot.2009.03.006

Portillo Lemus, L. O., Bozec, M., Harang, M., Coudreuse, J., Haury, J., Stoeckel, S., & Barloy, D. (2021). Self-incompatibility limits sexual reproduction rather than environmental conditions in an invasive water primrose. Plant-Environment Interactions, 2(2), 74–86. https://doi.org/10.1002/pei3.10042

Portillo Lemus, L. O., Harang, M., Bozec, M., Haury, J., Stoeckel, S., & Barloy, D. (2022). Late-acting self-incompatible system, preferential allogamy and delayed selfing in the heteromorphic invasive populations of Ludwigia grandiflora subsp. hexapetala. Peer Community Journal, 2. https://doi.org/10.24072/pcjournal.108

Stoeckel, S., Becheler, R., Portillo-Lemus, L., Harang, M., Besnard, A.-L., Lassalle, G., Causse-Védrines, R., Michon-Coudouel, S., Park D. J., Pope, B. J., Petit, E. J. & Barloy, D. (2024) Reproductive modes in populations of late-acting self-incompatible and self-compatible polyploid Ludwigia grandiflora subsp. hexapetala in western Europe. biorxiv, ver.4 peer-reviewed and recommended by PCI Evol Biol https://doi.org/10.1101/2024.03.21.586104

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.

Spatial genetic structure, i.e. the non-random spatial distribution of genotypes, arises in populations because of different processes including spatially limited dispersal and selection. Knowledge on the spatial genetic structure of plant populations is important to assess biological parameters such as gene dispersal distances and the potential for local adaptations, as well as for applications in conservation management and breeding. In their work, Niskanen and colleagues demonstrate a multifaceted approach to characterise the spatial genetic structure in two replicate sites of a continuously distributed Scots pine population in South-Eastern Finland. They mapped and assessed the ages of 469 naturally regenerated adults and genotyped them using a SNP array which resulted in 157 325 filtered polymorphic SNPs. Their dataset is remarkably powerful because of the large numbers of both individuals and SNPs genotyped. This made it possible to characterise precisely the decay of genetic relatedness between individuals with spatial distance despite the extensive dispersal capacity of Scots pine through pollen, and ensuing expectations of an almost panmictic population.

The authors’ data analysis was particularly thorough. They demonstrated that two metrics of pairwise relatedness, the genomic relationship matrix (GRM, Yang et al. 2011) and the kinship coefficient (Loiselle et al. 1995) were strongly correlated and produced very similar inference of family relationships: >99% of pairs of individuals were unrelated, and the remainder exhibited 2nd (e.g., half-siblings) to 4th degree relatedness. Pairwise relatedness decayed with spatial distance which resulted in extremely weak but statistically significant spatial genetic structure in both sites, quantified as Sp=0.0005 and Sp=0.0008. These estimates are at least an order of magnitude lower than estimates in the literature obtained in more fragmented populations of the same species or in other conifers. Estimates of the neighbourhood size, the effective number of potentially mating individuals belonging to a within-population neighbourhood (Wright 1946), were relatively large with Nb=1680-3210 despite relatively short gene dispersal distances, σg = 36.5–71.3m, which illustrates the high effective density of the population. 

The authors showed the implications of their findings for selection. The capacity for local adaptation depends on dispersal distances and the strength of the selection coefficient. In the study population, the authors inferred that local adaptation can only occur if environmental heterogeneity occurs over a distance larger than approximately one kilometre (or larger, if considering long-distance dispersal). Interestingly, in Scots pine, no local adaptation has been described on similar geographic scales, in contrast to some other European or Mediterranean conifers (Scotti et al. 2023).

The authors’ results are relevant for the management of conservation and breeding. They showed that related individuals occurred within sites only and that they shared a higher number of rare alleles than unrelated ones. Since rare alleles are enriched in new and recessive deleterious variants, selecting related individuals could have negative consequences in breeding programmes. The authors also showed, in their response to reviewers, that their powerful dataset was not suitable to obtain a robust estimate of effective population size, Ne, based on the linkage disequilibrium method (Do et al. 2014). This illustrated that the estimation of Ne used for genetic indicators supported in international conservation policy (Hoban et al. 2020, CBD 2022) remains challenging in large and continuous populations (see also Santo-del-Blanco et al. 2023, Gargiulo et al. 2024).

References

CBD (2022) Kunming-Montreal Global Biodiversity Framework. https://www.cbd.int/doc/decisions/cop-15/cop-15-dec-04-en.pdf

Do C, Waples RS, Peel D, Macbeth GM, Tillett BJ, Ovenden JR (2014). NeEstimator v2: re-implementation of software for the estimation of contemporary effective population size (Ne ) from genetic data. Molecular Ecology Resources 14: 209–214. https://doi.org/10.1111/1755-0998.12157

Gargiulo R, Decroocq V, González-Martínez SC, Paz-Vinas I, Aury JM, Kupin IL, Plomion C, Schmitt S, Scotti I, Heuertz M (2024) Estimation of contemporary effective population size in plant populations: limitations of genomic datasets. Evolutionary Applications, in press, https://doi.org/10.1101/2023.07.18.549323

Hoban S, Bruford M, D’Urban Jackson J, Lopes-Fernandes M, Heuertz M, Hohenlohe PA, Paz-Vinas I, et al. (2020) Genetic diversity targets and indicators in the CBD post-2020 Global Biodiversity Framework must be improved. Biological Conservation 248: 108654. https://doi.org/10.1016/j.biocon.2020.108654

Loiselle BA, Sork VL, Nason J & Graham C (1995) Spatial genetic structure of a tropical understorey shrub, Psychotria officinalis (Rubiaceae). American Journal of Botany 82: 1420–1425. https://doi.org/10.1002/j.1537-2197.1995.tb12679.x

Santos-del-Blanco L, Olsson S, Budde KB, Grivet D, González-Martínez SC, Alía R, Robledo-Arnuncio JJ (2022). On the feasibility of estimating contemporary effective population size (Ne) for genetic conservation and monitoring of forest trees. Biological Conservation 273: 109704. https://doi.org/10.1016/j.biocon.2022.109704

Scotti I, Lalagüe H, Oddou-Muratorio S, Scotti-Saintagne C, Ruiz Daniels R, Grivet D, et al. (2023) Common microgeographical selection patterns revealed in four European conifers. Molecular Ecology 32: 393-411. https://doi.org/10.1111/mec.16750

Wright S (1946) Isolation by distance under diverse systems of mating. Genetics 31: 39–59. https://doi.org/10.1093/genetics/31.1.39

Yang J, Lee SH, Goddard ME & Visscher PM (2011) GCTA: a tool for genome-wide complex trait analysis. The American Journal of Human Genetics 88: 76–82. https://www.cell.com/ajhg/pdf/S0002-9297(10)00598-7.pdf