The emergence and spread of insecticide resistance compromises the efficiency of insecticides as prevention tool against the transmission of insect-transmitted diseases (Moyes et al. 2017). In this context, the understanding of the genetic mechanisms of resistance and the way resistant alleles spread in insect populations is necessary and important to envision resistance management policies. A common and important mechanism of insecticide resistance is gene amplification and in particular amplification of insecticide detoxification genes, which leads to the overexpression of these genes (Bass & Field, 2011). Cattel and coauthors (2020) adopt a combination of experimental approaches to study the role of gene amplification in resistance to organophosphate insecticides in the mosquito Aedes aegypti and its occurrence in populations of South East Asia and to develop a molecular test to track resistance alleles.
Their first approach consists in performing an artificial selection on laboratory Ae. Aegypti populations started with individuals collected in Laos. In the selected population, an initial 90% mortality by adult exposure to the organophosphate insecticide malathion is imposed. This population shows a steep increase in resistance to malathion and other organophosphate insecticides, which is absent in the paired control population. The transcriptomic patterns of the control and the evolved populations as well as of a reference sensitive population reveals, among other differences, the over-expression of five carboxy/choline esterase (CCE) genes in the insecticide selected population. These five genes happen to be clustered in the Ae. aegypti genome and whole genome sequencing of a highly resistant population combined to qPCR test on genomic DNA showed that the overexpression of these genes is due to gene amplification. Although it would have been more elegant to have replicate selected and control populations and to perform the transcriptomic and the genomic analyses directly on the experimental populations, the authors gather a set of experimental evidence which combined to previous knowledge on the function of the amplified and over-expressed genes and on their implication in organophosphate insecticide resistance in other species allow to discard the possibility that this gene amplification spread by drift in the selected population.
In a second part of the paper, copy number variation for CCE genes is checked in field sample populations. This test reveals the presence of resistance alleles in half of the fourteen South East Asia populations sampled. Very interestingly, it also reveals a high level of complexity and diversity among the resistance alleles: it shows first the existence, both in the experimental and the field populations, of at least two amplified alleles (differing by the number of genes amplified) and second a high variation in the copy number of amplified genes. This indicates that gene amplification as a molecular resistance mechanism has actually lead to a high diversity of resistance alleles. These alleles are likely to differ both by the level of resistance conferred and the fitness cost imposed in the absence of the insecticide and these two values are affecting the evolution of their frequency in the field and ultimately the spread of resistance.
The last part of the paper is devoted to the development of a high-throughput Taqman assay which allows to determine rapidly the copy number of one of the esterase genes amplified in the resistance alleles described earlier. This assay is nicely validated and will definitely be a useful tool to determine the occurrence of these resistance alleles in field population. The fact that it gives access to the copy number will also allow to follow its copy number across time and get insight into the complexity of resistance evolution by gene amplification.
To sum up, this paper studies the implication of carboxy/choline esterase genes amplification in organophosphate resistance evolution in Ae. aegypti, reveals the diversity among individuals and populations of this resistance mechanism, because of variation both in the identity of the genes amplified and in their copy number and sets up a fast and efficient tool to detect and follow the spread of these resistant alleles in the field. Additionally, the different experimental approaches adopted have generated genomic and transcriptomic data, of which only the part related to CCE gene amplification has been exploited. These data are very likely to reveal other genomic and expression determinants of resistance that will give access to an extra degree of complexity in organophosphate insecticide resistance determinism and evolution.

References

Bass C, Field LM (2011) Gene amplification and insecticide resistance. Pest Management Science, 67, 886–890. https://doi.org/10.1002/ps.2189
Cattel J, Haberkorn C, Laporte F, Gaude T, Cumer T, Renaud J, Sutherland IW, Hertz JC, Bonneville J-M, Arnaud V, Nous C, Fustec B, Boyer S, Marcombe S, David J-P (2020) A genomic amplification affecting a carboxylesterase gene cluster confers organophosphate resistance in the mosquito Aedes aegypti: from genomic characterization to high-throughput field detection. bioRxiv, 2020.06.08.139741, ver. 4 peer-reviewed and recommended by PCI Evolutionary Biology. https://doi.org/10.1101/2020.06.08.139741
Moyes CL, Vontas J, Martins AJ, Ng LC, Koou SY, Dusfour I, Raghavendra K, Pinto J, Corbel V, David J-P, Weetman D (2017) Contemporary status of insecticide resistance in the major Aedes vectors of arboviruses infecting humans. PLOS Neglected Tropical Diseases, 11, e0005625. https://doi.org/10.1371/journal.pntd.0005625

Recent advances in next-generation sequencing are allowing us to uncover the evolution of sex chromosomes in non-model organisms. This study [1] represents an example of this application to birds of two Sylvioidea species from the genus Zosterops (commonly known as white-eyes). The study is exemplary in the amount and types of data generated and in the thoroughness of the analysis applied. Both male and female genomes were sequenced to allow the authors to identify sex-chromosome specific scaffolds. These data were augmented by generating the transcriptome (RNA-seq) data set. The findings after the analysis of these extensive data are intriguing: neoZ and neoW chromosome scaffolds and their breakpoints were identified. Novel sex chromosome formation appears to be accompanied by translocation events. The timing of formation of novel sex chromosomes was identified using molecular dating and appears to be relatively recent. Yet first signatures of distinct evolutionary patterns of sex chromosomes vs. autosomes could be already identified. These include the accumulation of transposable elements and changes in GC content. The changes in GC content could be explained by biased gene conversion and altered recombination landscape of the neo sex chromosomes. The authors also study divergence and diversity of genes located on the neo sex chromosomes. Here their findings appear to be surprising and need further exploration. The neoW chromosome already shows unique patterns of divergence and diversity at protein-coding genes as compared with genes on either neoZ or autosomes. In contrast, the genes on the neoZ chromosome do not display divergence or diversity patterns different from those for autosomes. This last observation is puzzling and I believe should be explored in further studies. Overall, this study significantly advances our knowledge of the early stages of sex chromosome evolution in vertebrates, provides an example of how such a study could be conducted in other non-model organisms, and provides several avenues for future work.

References

[1] Leroy T., Anselmetti A., Tilak M.K., Bérard S., Csukonyi L., Gabrielli M., Scornavacca C., Milá B., Thébaud C. and Nabholz B. (2019). A bird’s white-eye view on neo-sex chromosome evolution. bioRxiv, 505610, ver. 4 peer-reviewed and recommended by PCI Evolutionary Biology. doi: 10.1101/505610

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

The present-day distribution of extant species is the result of the interplay between their past population demography (e.g., expansion, contraction, isolation, and migration) and adaptation to the environment. Shedding light on the timing and magnitude of key demographic events helps identify potential drivers of such events and interaction of those drivers, such as life history traits and past episodes of environmental shifts.

The understanding of the key factors driving species evolution gives important insights into how the species may respond to changing conditions, which can be particularly relevant for the management of harmful species, such as agricultural pests (e.g. [1]). Meaningful demographic inferences present major challenges. These include formulating evolutionary scenarios fitting species biology and the eco-geographical context and choosing informative molecular markers and accurate quantitative approaches to statistically compare multiple demographic scenarios and estimate the parameters of interest. A further issue comes with result interpretation. Accurately dating the inferred events is far from straightforward since reliable calibration points are necessary to translate the molecular estimates of the evolutionary time into absolute time units (i.e. years). This can be attempted in different ways, such as by using fossil and archaeological records, heterochronous samples (e.g. ancient DNA), and/or mutation rate estimated from independent data (e.g. [2], [3] for review). Nonetheless, most experimental systems rarely meet these conditions, hindering the comprehensive interpretation of results.

The contribution of Chapuis et al. [4] addresses these issues to investigate the recent history of the African insect pest Schistocerca gregaria (desert locust). They apply Approximate Bayesian Computation-Random Forest (ABC-RF) approaches to microsatellite markers. Owing to their fast mutation rate microsatellite markers offer at least two advantages: i) suitability for analyzing recently diverged populations, and ii) direct estimate of the germline mutation rate in pedigree samples. The work of Chapuis et al. [4] benefits of both these advantages, since they have estimates of mutation rate and allele size constraints derived from germline mutations in the species [5].

The main aim of the study is to infer the history of divergence of the two subspecies of the desert locust, which have spatially disjoint distribution corresponding to the dry regions of North and West-South Africa. They first use paleo-vegetation maps to formulate hypotheses about changes in species range since the last glacial maximum. Based on them, they generate 12 divergence models. For the selection of the demographic model and parameter estimation, they apply the recently developed ABC-RF approach, a powerful inferential tool that allows optimizing the use of summary statistics information content, among other advantages [6]. Some methodological novelties are also introduced in this work, such as the computation of the error associated with the posterior parameter estimates under the best scenario. The accuracy of timing estimate is assured in two ways: i) by the use of microsatellite markers with known evolutionary dynamics, as underlined above, and ii) by assessing the divergence time threshold above which posterior estimates are likely to be biased by size homoplasy and limits in allele size range [7]. The best-supported model suggests a recent divergence event of the subspecies of S. gregaria (around 2.6 kya) and a reduction of populations size in one of the subspecies (S. g. flaviventris) that colonized the southern distribution area. As such, results did not support the hypothesis that the southward colonization was driven by the expansion of African dry environments associated with the last glacial maximum, as it has been postulated for other arid-adapted species with similar African disjoint distributions [8]. The estimated time of divergence points at a much more recent origin for the two subspecies, during the late Holocene, in a period corresponding to fairly stable arid conditions similar to current ones [9,10].

Although the authors cannot exclude that their microsatellite data bear limited information on older colonization events than the last one, they bring arguments in favour of alternative explanations. The hypothesis privileged does not involve climatic drivers, but the particularly efficient dispersal behaviour of the species, whose individuals are able to fly over long distances (up to thousands of kilometers) under favourable windy conditions. A single long-distance dispersal event by a few individuals would explain the genetic signature of the bottleneck. There is a growing number of studies in phylogeography in arid regions in the Southern hemisphere, but the impact of past climate changes on the species distribution in this region remains understudied relative to the Northern hemisphere [11,12].

The study presented by Chapuis et al. [4] offers several important insights into demographic changes and the evolutionary history of an agriculturally important pest species in Africa, which could also mirror the history of other organisms in the continent. As the authors point out, there are necessarily some uncertainties associated with the models of past ecosystems and climate, especially for Africa. Interestingly, the authors argue that the information on paleo-vegetation turnover was more informative than climatic niche modeling for the purpose of their study since it made them consider a wider range of bio-geographical changes and in turn a wider range of evolutionary scenarios (see discussion in Supplementary Material). Microsatellite markers have been offering a useful tool in population genetics and phylogeography for decades, but their popularity is perhaps being taken over by single nucleotide polymorphism (SNP) genotyping and whole-genome sequencing (WGS) (the peak year of the number of the publication with “microsatellite” is in 2012 according to PubMed).

This study reaffirms the usefulness of these classic molecular markers to estimate past demographic events, especially when species- and locus-specific microsatellite mutation features are available and a powerful inferential approach is adopted. Nonetheless, there are still hurdles to overcome, such as the limitations in scenario choice associated with the simulation software used (e.g. not allowing for continuous gene flow in this particular case), which calls for further improvement of simulation tools allowing for more flexible modeling of demographic events and mutation patterns. In sum, this work not only contributes to our understanding of the makeup of the African biodiversity but also offers a useful statistical framework, which can be applied to a wide array of species and molecular markers (microsatellites, SNPs, and WGS).

References

[1] Lehmann, P. et al. (2018). Complex responses of global insect pests to climate change. bioRxiv, 425488. doi: https://dx.doi.org/10.1101/425488

[2] Donoghue, P. C., & Benton, M. J. (2007). Rocks and clocks: calibrating the Tree of Life using fossils and molecules. Trends in Ecology & Evolution, 22(8), 424-431. doi: https://dx.doi.org/10.1016/j.tree.2007.05.005

[3] Ho, S. Y., Lanfear, R., Bromham, L., Phillips, M. J., Soubrier, J., Rodrigo, A. G., & Cooper, A. (2011). Time‐dependent rates of molecular evolution. Molecular ecology, 20(15), 3087-3101. doi: https://dx.doi.org/10.1111/j.1365-294X.2011.05178.x

[4] Chapuis, M.-P., Raynal, L., Plantamp, C., Meynard, C. N., Blondin, L., Marin, J.-M. and Estoup, A. (2020). A young age of subspecific divergence in the desert locust Schistocerca gregaria, inferred by ABC Random Forest. bioRxiv, 671867, ver. 4 peer-reviewed and recommended by PCI Evolutionary Biology. doi: https://dx.doi.org/10.1101/671867

5] Chapuis, M.-P., Plantamp, C., Streiff, R., Blondin, L., & Piou, C. (2015). Microsatellite evolutionary rate and pattern in Schistocerca gregaria inferred from direct observation of germline mutations. Molecular ecology, 24(24), 6107-6119. doi: https://dx.doi.org/10.1111/mec.13465

[6] Raynal, L., Marin, J. M., Pudlo, P., Ribatet, M., Robert, C. P., & Estoup, A. (2018). ABC random forests for Bayesian parameter inference. Bioinformatics, 35(10), 1720-1728. doi: https://dx.doi.org/10.1093/bioinformatics/bty867

[7] 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: https://dx.doi.org/10.1046/j.1365-294X.2002.01576.x

[8] Moodley, Y. et al. (2018). Contrasting evolutionary history, anthropogenic declines and genetic contact in the northern and southern white rhinoceros (Ceratotherium simum). Proceedings of the Royal Society B, 285(1890), 20181567. doi: https://dx.doi.org/10.1098/rspb.2018.1567

[9] Kröpelin, S. et al. (2008). Climate-driven ecosystem succession in the Sahara: the past 6000 years. science, 320(5877), 765-768. doi: https://dx.doi.org/10.1126/science.1154913

[10] Maley, J. et al. (2018). Late Holocene forest contraction and fragmentation in central Africa. Quaternary Research, 89(1), 43-59. doi: https://dx.doi.org/10.1017/qua.2017.97

[11] Beheregaray, L. B. (2008). Twenty years of phylogeography: the state of the field and the challenges for the Southern Hemisphere. Molecular Ecology, 17(17), 3754-3774. doi: https://dx.doi.org/10.1111/j.1365-294X.2008.03857.x

[12] Dubey, S., & Shine, R. (2012). Are reptile and amphibian species younger in the Northern Hemisphere than in the Southern Hemisphere?. Journal of evolutionary biology, 25(1), 220-226. doi: https://dx.doi.org/10.1111/j.1420-9101.2011.02417.x

*****

A video about this preprint is available here:

Differentiation-based genome scans lie at the core of speciation and adaptation genomics research. Dating back to Lewontin & Krakauer (1973), they have become very popular with the advent of genomics to identify genome regions of enhanced differentiation relative to neutral expectations. These regions may represent genetic barriers between divergent lineages and are key for studying reproductive isolation. However, genome scan methods can generate a high rate of false positives, primarily if the neutral population structure is not accounted for (Bierne et al. 2013). Moreover, interpreting genome scans can be challenging in the context of secondary contacts between diverging lineages (Bierne et al. 2011), because the coupling between different components of reproductive isolation (local adaptation, intrinsic incompatibilities, mating preferences, etc.) can occur readily, thus preventing the causes of differentiation from being determined.

Smadja and collaborators (2022) applied a sophisticated genome scan for trait association (BAYPASS, Gautier 2015) to underlie the genetic basis of a polygenetic behaviour: assortative mating in hybridizing mice. My interest in this neat study mainly relies on two reasons. First, the authors used an ingenious geographical setting (replicate pairs of “Choosy” versus “Non-Choosy” populations) with multi-way comparisons to narrow down the list of candidate regions resulting from BAYPASS. The latter corrects for population structure, handles cost-effective pool-seq data and allows for gene-based analyses that aggregate SNP signals within a gene. These features reinforce the set of outlier genes detected; however, not all are expected to be associated with mating preference. 

The second reason why this study is valuable to me is that Smadja et al. (2022) complemented the population genomic approach with functional predictions to validate the genetic signal. In line with previous behavioural and chemical assays on the proximal mechanisms of mating preferences, they identified multiple olfactory and vomeronasal receptor genes as highly significant candidates. Therefore, combining genomic signals with functional analyses is a clever way to provide insights into the causes of reproductive isolation, especially when multiple barriers are involved. This is typically true for reinforcement (Butlin & Smadja 2018), suspected to occur in these mice because, in that case, assortative mating (a prezygotic barrier) evolves in response to the cost of hybridization (for example, due to hybrid inviability). 

As advocated by the authors, their study paves the way for future work addressing the genetic basis of reinforcement, a trait of major evolutionary importance for which we lack empirical data. They also make a compelling case using complementary approaches that olfactory and vomeronasal receptors have a central role in mammal speciation.

References:

Bierne N, Welch J, Loire E, Bonhomme F, David P (2011) The coupling hypothesis: why genome scans may fail to map local adaptation genes. Mol Ecol 20: 2044–2072. https://doi.org/10.1111/j.1365-294X.2011.05080.x

Bierne N, Roze D, Welch JJ (2013) Pervasive selection or is it…? why are FST outliers sometimes so frequent? Mol Ecol 22: 2061–2064. https://doi.org/10.1111/mec.12241 

Butlin RK, Smadja CM (2018) Coupling, Reinforcement, and Speciation. Am Nat 191:155–172. https://doi.org/10.1086/695136 

Gautier M (2015) Genome-Wide Scan for Adaptive Divergence and Association with Population-Specific Covariates. Genetics 201:1555–1579. https://doi.org/10.1534/genetics.115.181453 

Lewontin RC, Krakauer J (1973) Distribution of gene frequency as a test of the theory of selective neutrality of polymorphisms. Genetics 74: 175–195. https://doi.org/10.1093/genetics/74.1.175 

Smadja CM, Loire E, Caminade P, Severac D, Gautier M, Ganem G (2022) Divergence of olfactory receptors associated with the evolution of assortative mating and reproductive isolation in mice. bioRxiv, 2022.07.21.500634, ver. 3 peer-reviewed and recommended by Peer Community in Evolutionary Biology. https://doi.org/10.1101/2022.07.21.500634