Parasitoid wasps have developed different mechanisms to increase their parasitic success, usually at the expense of host survival (Fellowes and Godfray, 2000). Eggs of these insects are deposited inside the juvenile stages of their hosts, which in turn deploy several immune response strategies to eliminate or disable them (Yang et al., 2020). Drosophila melanogaster protects itself against parasitoid attacks through the production of specific elongated haemocytes called lamellocytes which form a capsule around the invading parasite (Lavine and Strand, 2002; Rizki and Rizki, 1992) and the subsequent activation of the phenol-oxidase cascade leading to the release of toxic radicals (Nappi et al., 1995). On the parasitoid side, robust responses have evolved to evade host immune defenses as for example the Drosophila-specific endoparasite Leptopilina boulardi, which releases venom during oviposition that modifies host behaviour (Varaldi et al., 2006) and inhibits encapsulation (Gueguen et al., 2011; Martinez et al., 2012).
Studies have shown that the wasp parasitic capacity is correlated to venom presence and its content (Colinet et al., 2009; Poirié et al., 2014), including that evolution of venom protein composition is driven by different levels of host susceptibility to infection (Cavigliasso et al., 2019). However, it had not been determined to this day, if and how parasitic range can affect venom protein composition and to which extent host specialization requires broad-spectrum factors or a plethora of specialized components.
These outstanding questions are now approached in a study by Cavigliasso and colleagues (Cavigliasso et al., 2021), where they perform experimental evolution of L. boulardi for 9 generations exposing it to different Drosophila host species and genetic backgrounds (two strains of D. melanogaster, D. simulans and D. yakuba). The authors tested whether the parasitic success of each selection regime was host-specific and how they influenced venom composition in parasitoids. For the first part, infection outcomes were assayed for each selection regime when cross-infecting different hosts. To get a finer measurement of the mechanisms under selection, the authors differentiated three phenotypes: overall parasitic success, encapsulation inhibition and escape from capsule. Throughout the course of experimental evolution, only encapsulation inhibition did not show an improved response upon selection on any host. Importantly, the cross-infection scenario revealed a clear specificity to the selected host for each evolved resistance.
As for venom composition, a trend of differential evolution was detected between host species, although a significant part of that was due to a larger differentiation in the D. yakuba regime, which showed a completely different directionality. Importantly, the authors could identify some of the specific proteins targeted by the several selection regimes, whether selected or counter-selected for. Interestingly, the D. yakuba regime is the only case where the key parasitoid protein LbSPNy (Colinet et al., 2009) was not counter-selected and the only regime in which the overall venom composition did not evolve towards the Ism strain, one of the two ancestral strains of L. boulardi used in the study. It is possible that these two results are correlated, since LbSPNy has been described to inhibit activation of the phenoloxidase cascade in D. yakuba and is one of the most abundant proteins in the ISy venom, making it a good target for selection (Colinet et al., 2013). The authors also discuss the possibility that this difference is related to the geographical distribution of the strains of L. boulardi, since each coincide with either D. melanogaster or D. yakuba.
This methodologically broad work by Cavigliasso and colleagues constitutes an important experimental contribution towards the understanding of how parasitoid adaptation to specific hosts is achieved at different phenotypic and mechanistic levels. It provides compelling evidence that venom composition evolves differently in response to specific parasitic ranges, particularly considering the evolutionary difference between the selective hosts. In line with this result, it is also concluded that the majority of venom proteins selected are lineage-specific, although a few broad-spectrum factors could also be detected. 
The question of whether parasitic range can affect venom composition and parasitic success is still open to more contributions. A potentially interesting long-term direction will be to use a similar setup of experimental evolution on the generalist L. heterotoma (Schlenke et al., 2007) . On a more immediate horizon, comparing the venom evolution of both L. heterotoma and L. boulardi under selection with different hosts and under cross-infection scenarios could reveal interesting patterns. The recent sequencing of the L. boulardi genome together with the vast number of studies addressing mechanisms of Drosophila resistance to parasitoid infection, will enable the thorough characterization of the genetic basis of host-parasitoid interactions and the deeper understanding of these ubiquitous and economically-relevant relationships.
 
*This recommendation text has been co-written with Tânia F. Paulo who is not a recommender of PCI Evol Biol

References

Cavigliasso, F., Mathé-Hubert, H., Gatti, J.-L., Colinet, D. and Poirié, M. (2021) Parasitic success and venom composition evolve upon specialization of parasitoid wasps to different host species. bioRxiv, 2020.10.24.353417, ver. 3 peer-reviewed and recommended by Peer Community in Evolutionary Biology. https://doi.org/10.1101/2020.10.24.353417

Cavigliasso, F., Mathé-Hubert, H., Kremmer, L., Rebuf, C., Gatti, J.-L., Malausa, T., Colinet, D., Poiré, M. and  Léne. (2019). Rapid and Differential Evolution of the Venom Composition of a Parasitoid Wasp Depending on the Host Strain. Toxins, 11(629). https://doi.org/10.3390/toxins11110629

Colinet, D., Deleury, E., Anselme, C., Cazes, D., Poulain, J., Azema-Dossat, C., Belghazi, M., Gatti, J. L. and  Poirié, M. (2013). Extensive inter- and intraspecific venom variation in closely related parasites targeting the same host: The case of Leptopilina parasitoids of Drosophila. Insect Biochemistry and Molecular Biology, 43(7), 601–611. https://doi.org/10.1016/j.ibmb.2013.03.010

Colinet, D., Dubuffet, A., Cazes, D., Moreau, S., Drezen, J. M. and  Poirié, M. (2009). A serpin from the parasitoid wasp Leptopilina boulardi targets the Drosophila phenoloxidase cascade. Developmental and Comparative Immunology, 33(5), 681–689. https://doi.org/10.1016/j.dci.2008.11.013

Fellowes, M. D. E. and  Godfray, H. C. J. (2000). The evolutionary ecology of resistance to parasitoids by Drosophila. Heredity, 84(1), 1–8. https://doi.org/10.1046/j.1365-2540.2000.00685.x

Gueguen, G., Rajwani, R., Paddibhatla, I., Morales, J. and  Govind, S. (2011). VLPs of Leptopilina boulardi share biogenesis and overall stellate morphology with VLPs of the heterotoma clade. Virus Research, 160(1–2), 159–165. https://doi.org/10.1016/j.virusres.2011.06.005

Lavine, M. D. and  Strand, M. R. (2002). Insect hemocytes and their role in immunity. Insect Biochemistry and Molecular Biology, 32(10), 1295–1309. https://doi.org/10.1016/S0965-1748(02)00092-9

Martinez, J., Duplouy, A., Woolfit, M., Vavre, F., O’Neill, S. L. and  Varaldi, J. (2012). Influence of the virus LbFV and of Wolbachia in a host-parasitoid interaction. PloS One, 7(4), e35081. https://doi.org/10.1371/journal.pone.0035081

Nappi, A. J., Vass, E., Frey, F. and  Carton, Y. (1995). Superoxide anion generation in Drosophila during melanotic encapsulation of parasites. European Journal of Cell Biology, 68(4), 450–456.

Poirié, M., Colinet, D. and  Gatti, J. L. (2014). Insights into function and evolution of parasitoid wasp venoms. Current Opinion in Insect Science, 6, 52–60. https://doi.org/10.1016/j.cois.2014.10.004

Rizki, T. M. and  Rizki, R. M. (1992). Lamellocyte differentiation in Drosophila larvae parasitized by Leptopilina. Developmental and Comparative Immunology, 16(2–3), 103–110. https://doi.org/10.1016/0145-305X(92)90011-Z

Schlenke, T. A., Morales, J., Govind, S. and  Clark, A. G. (2007). Contrasting infection strategies in generalist and specialist wasp parasitoids of Drosophila melanogaster. PLoS Pathogens, 3(10), 1486–1501. https://doi.org/10.1371/journal.ppat.0030158

Varaldi, J., Petit, S., Boulétreau, M. and  Fleury, F. (2006). The virus infecting the parasitoid Leptopilina boulardi exerts a specific action on superparasitism behaviour. Parasitology, 132(Pt 6), 747–756. https://doi.org/10.1017/S0031182006009930

Yang, L., Qiu, L., Fang, Q., Stanley, D. W. and  Gong‐Yin, Y. (2020). Cellular and humoral immune interactions between Drosophila and its parasitoids. Insect Science. https://doi.org/10.1111/1744-7917.12863

Genome-wide data obtained from even a small number of individuals can provide unprecedented levels of detail about the evolutionary history of populations and species [1], determinants of genetic diversity [2], species boundaries and the process of speciation itself [3]. Loire and Galtier [4] present a clear example, using the emblematic Galápagos giant tortoise (Chelonoidis nigra), of how multi-species comparative population genomic approaches can provide valuable insights about population structure and species delimitation even when sample sizes are limited but the number of loci is large and distributed across the genome.

Galápagos giant tortoises are endemic to the Galápagos Islands and are currently recognized as an endangered, multi-species complex including both extant and extinct taxa. Taxonomic definitions are based on morphology, geographic isolation and population genetic evidence based on short DNA sequences of the mitochondrial genome (mtDNA) and/or a dozen or so nuclear microsatellite loci [5-8]. The species complex enjoys maximal protection. Population recoveries have been quite successful and spectacular conservation programs based on mitochondrial genes and microsatellites are ongoing. This includes for example individual translocations, breeding program, “hybrid” sterilization or removal, and resurrection of extinct lineages).

In 2013, Loire et al. [9] provided the first population genomic analyses based on genome scale data (~1000 coding loci derived from blood-transcriptomes) from five individuals, encompassing three putative “species”: Chelonnoidis becki, C. porteri and C. vandenburghi. Their results raised doubts about the validity/accuracy of the currently accepted designations of “genetic distinctiveness”. However, the implications for conservation and management have remained unnoticed.

In 2017, Loire and Galtier [4] have re-appraised this issue using an original multi-species comparative population genomic analysis of their previous data set [9]. Based on a comparison of 53 animal species, they show that both the level of genome-wide neutral diversity (πS) and level of population structure estimated using the inbreeding coefficient (F) are much lower than would be expected from a sample covering multiple species. The observed values are more comparable to those typically reported at the “among population” level within a single species such as human (Homo sapiens). The authors go to great length to assess the sensitivity of their method to detect population structure (or lack thereof) and show that their results are robust to potential issues, such as contamination and sequencing errors that can occur with Next Generation Sequencing techniques; and biases related to the small sample size and sub-sampling of individuals. They conclude that published mtDNA and microsatellite-based assessment of population structure and species designations may be biased towards over-splitting.

This manuscript is a very good read as it shows the potential of the now relatively affordable genome-wide data for helping to both resolve and clarify population and species boundaries, illuminate demographic trends, evolutionary trajectories of isolated groups, patterns of connectivity among them, and test for evidence of local adaptation and even reproductive isolation. The comprehensive information provided by genome-wide data can critically inform and assist the development of the best strategies to preserve endangered populations and species. Loire and Galtier [4] make a strong case for applying genomic data to the Galápagos giant tortoises, which is likely to redirect conservation efforts more effectively and at lower cost. The case of the Galápagos giant tortoises is certainly a very emblematic example, which will find an echo in many other endangered species conservation programs.

References

[1] Li H and Durbin R. 2011. Inference of human population history from individual whole-genome sequences. Nature, 475: 493–496. doi: 10.1038/nature10231

[2] Romiguier J, Gayral P, Ballenghien M, Bernard A, Cahais V, Chenuil A, Chiari Y, Dernat R, Duret L, Faivre N, Loire E, Lourenco JM, Nabholz B, Roux C, Tsagkogeorga G, Weber AA-T, Weinert LA, Belkhir K, Bierne N, Glémin S and Galtier N. 2014. Comparative population genomics in animals uncovers the determinants of genetic diversity. Nature, 515: 261–263. doi: 10.1038/nature13685

[3] Roux C, Fraïsse C, Romiguier J, Anciaux Y, Galtier N and Bierne N. 2016. Shedding light on the grey zone of speciation along a continuum of genomic divergence. PLoS Biology, 14: e2000234. doi: 10.1371/journal.pbio.2000234

[4] Loire E and Galtier N. 2017. Lacking conservation genomics in the giant Galápagos tortoise. bioRxiv 101980, ver. 4 of September 26, 2017. doi: 10.1101/101980

[5] Beheregaray LB, Ciofi C, Caccone A, Gibbs JP and Powell JR. 2003. Genetic divergence, phylogeography and conservation units of giant tortoises from Santa Cruz and Pinzón, Galápagos Islands. Conservation Genetics, 4: 31–46. doi: 10.1023/A:1021864214375

[6] Ciofi C, Milinkovitch MC, Gibbs JP, Caccone A and Powell JR. 2002. Microsatellite analysis of genetic divergence among populations of giant Galápagos tortoises. Molecular Ecology, 11: 2265–2283. doi: 10.1046/j.1365-294X.2002.01617.x

[7] Garrick RC, Kajdacsi B, Russello MA, Benavides E, Hyseni C, Gibbs JP, Tapia W and Caccone A. 2015. Naturally rare versus newly rare: demographic inferences on two timescales inform conservation of Galápagos giant tortoises. Ecology and Evolution, 5: 676–694. doi: 10.1002/ece3.1388

[8] Poulakakis N, Edwards DL, Chiari Y, Garrick RC, Russello MA, Benavides E, Watkins-Colwell GJ, Glaberman S, Tapia W, Gibbs JP, Cayot LJ and Caccone A. 2015. Description of a new Galápagos giant tortoise species (Chelonoidis; Testudines: Testudinidae) from Cerro Fatal on Santa Cruz island. PLoS ONE, 10: e0138779. doi: 10.1371/journal.pone.0138779

[9] Loire E, Chiari Y, Bernard A, Cahais V, Romiguier J, Nabholz B, Lourenço JM and Galtier N. 2013. Population genomics of the endangered giant Galápagos tortoise. Genome Biology, 14: R136. doi: 10.1186/gb-2013-14-12-r136

Genetic diversity in temperate and boreal forests tree species has been strongly affected by late Pleistocene climate oscillations [2,3,5], but also by anthropogenic forces. Particularly in Europe, where a long history of human intervention has re-distributed species and populations, it can be difficult to know if a given forest arose through natural regeneration and gene flow or through some combination of natural and human-mediated processes. This uncertainty can confound inferences of the causes and consequences of standing genetic variation, which may impact our interpretation of demographic events that shaped species before humans became dominant on the landscape. In their paper entitled "Genomic data provides new insights on the demographic history and the extent of recent material transfers in Norway spruce", Chen et al. [1] used a genome-wide dataset of 400k SNPs to infer the demographic history of Picea abies (Norway spruce), the most widespread and abundant spruce species in Europe, and to understand its evolutionary relationship with two other spruces (Picea obovata [Siberian spruce] and P. omorika [Serbian spruce]). Three major Norway spruce clusters were identified, corresponding to central Europe, Russia and the Baltics, and Scandinavia, which agrees with previous studies. The density of the SNP data in the present paper enabled inference of previously uncharacterized admixture between these groups, which corresponds to the timing of postglacial recolonization following the last glacial maximum (LGM). This suggests that multiple migration routes gave rise to the extant distribution of the species, and may explain why Chen et al.'s estimates of divergence times among these major Norway spruce groups were older (15mya) than those of previous studies (5-6mya) – those previous studies may have unknowingly included admixed material [4]. Treemix analysis also revealed extensive admixture between Norway and Siberian spruce over the last ~100k years, while the geographically-restricted Serbian spruce was both isolated from introgression and had a dramatically smaller effective population size (Ne) than either of the other two species. This small Ne resulted from a bottleneck associated with the onset of the iron age ~3000 years ago, which suggests that anthropogenic depletion of forest resources has severely impacted this species. Finally, ancestry of Norway spruce samples collected in Sweden and Denmark suggest their recent transfer from more southern areas of the species range. This northward movement of genotypes likely occurred because the trees performed well relative to local provenances, which is a common observation when trees from the south are planted in more northern locations (although at the potential cost of frost damage due to inappropriate phenology). While not the reason for the transfer, the incorporation of southern seed sources into the Swedish breeding and reforestation program may lead to more resilient forests under climate change. Taken together, the data and analysis presented in this paper allowed inference of the intra- and interspecific demographic histories of a tree species group at a very high resolution, and suggest caveats regarding sampling and interpretation of data from areas with a long history of occupancy by humans.

References

[1] Chen, J., Milesi, P., Jansson, G., Berlin, M., Karlsson, B., Aleksić, J. M., Vendramin, G. G., Lascoux, M. (2018). Genomic data provides new insights on the demographic history and the extent of recent material transfers in Norway spruce. BioRxiv, 402016. ver. 3 peer-reviewed and recommended by PCI Evol Biol. doi: 10.1101/402016
[2] Holliday, J. A., Yuen, M., Ritland, K., & Aitken, S. N. (2010). Postglacial history of a widespread conifer produces inverse clines in selective neutrality tests. Molecular Ecology, 19(18), 3857–3864. doi: 10.1111/j.1365-294X.2010.04767.x
[3] Ingvarsson, P. K. (2008). Multilocus patterns of nucleotide polymorphism and the demographic history of Populus tremula. Genetics, 180, 329-340. doi: 10.1534/genetics.108.090431
[4] Lockwood, J. D., Aleksić, J. M., Zou, J., Wang, J., Liu, J., & Renner, S. S. (2013). A new phylogeny for the genus Picea from plastid, mitochondrial, and nuclear sequences. Molecular Phylogenetics and Evolution, 69(3), 717–727. doi: 10.1016/j.ympev.2013.07.004
[5] Pyhäjärvi, T., Garcia-Gil, M. R., Knürr, T., Mikkonen, M., Wachowiak, W., & Savolainen, O. (2007). Demographic history has influenced nucleotide diversity in European Pinus sylvestris populations. Genetics, 177(3), 1713–1724. doi: 10.1534/genetics.107.077099 "

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

References

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

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

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

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

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

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

Each individual multicellular organism, each of our bodies, is a small universe. Every living surface -skin, cuticle, bark, mucosa- is the home place to milliards of bacteria, fungi and viruses. They constitute our microbiota. Some of them are essential for certain organisms. Other could not live without their hosts. For many species, the relationship between host and microbiota is so close that their histories are inseparable. The recognition of this biological inextricability has led to the notion of holobiont as the organism ensemble of host and microbiota. When individuals of a particular animal or plant species expand their geographical range, it is the holobiont that expands. And these processes of migration, expansion and colonization are often accompanied by evolutionary and ecological innovations in the interspecies relationships, at the macroscopic level (e.g. novel predator-prey or host-parasite interactions) and at the microscopic level (e.g. changes in the microbiota composition). From the human point of view, these novel interactions can be economically disastrous if they involve and threaten important crop or cattle species. And this is especially worrying in the present context of genetic standardization and intensification for mass-production on the one hand, and of climate change on the other.

With this perspective, the international team led by Amandine Cornille presents a study aiming at understanding the evolutionary history of the rosy apple aphid Dysaphis plantaginea Passerini, a major pest of the cultivated apple tree Malus domestica Borkh (1). The apple tree was probably domesticated in Central Asia, and later disseminated by humans over the world in different waves, and it was probably introduced in Europe by the Greeks. It is however unclear when and where D. plantaginea started parasitizing the cultivated apple tree. The ancestral D. plantaginea could have already infected the wild ancestor of current cultivated apple trees, but the aphid is not common in Central Asia. Alternatively, it may have gained access only later to the plant, possibly via a host jump, from Pyrus to Malus that may have occurred in Asia Minor or in the Caucasus. In the present preprint, Olvera-Vázquez and coworkers have analysed over 650 D. plantaginea colonies from 52 orchards in 13 countries, in Western, Central and Eastern Europe as well as in Morocco and the USA. The authors have analysed the genetic diversity in the sampled aphids, and have characterized as well the composition of the associated endosymbiont bacteria. The analyses detect substantial recent admixture, but allow to identify aphid subpopulations slightly but significantly differentiated and isolated by distance, especially those in Morocco and the USA, as well as to determine the presence of significant gene flow. This process of colonization associated to gene flow is most likely indirectly driven by human interactions. Very interestingly, the data show that this genetic diversity in the aphids is not reflected by a corresponding diversity in the associated microbiota, largely dominated by a few Buchnera aphidicola variants. In order to determine polarity in the evolutionary history of the aphid-tree association, the authors have applied approximate Bayesian computing and machine learning approaches. Albeit promising, the results are not sufficiently robust to assess directionality nor to confidently assess the origin of the crop pest. Despite the large effort here communicated, the authors point to the lack of sufficient data (in terms of aphid isolates), especially originating from Central Asia. Such increased sampling will need to be implemented in the future in order to elucidate not only the origin and the demographic history of the interaction between the cultivated apple tree and the rosy apple aphid. This knowledge is needed to understand how this crop pest struggles with the different seasonal and geographical selection pressures while maintaining high genetic diversity, conspicuous gene flow, differentiated populations and low endosymbiontic diversity.

References

1. Olvera-Vazquez SG, Remoué C, Venon A, Rousselet A, Grandcolas O, Azrine M, Momont L, Galan M, Benoit L, David GM, Alhmedi A, Beliën T, Alins G, Franck P, Haddioui A, Jacobsen SK, Andreev R, Simon S, Sigsgaard L, Guibert E, Tournant L, Gazel F, Mody K, Khachtib Y, Roman A, Ursu TM, Zakharov IA, Belcram H, Harry M, Roth M, Simon JC, Oram S, Ricard JM, Agnello A, Beers EH, Engelman J, Balti I, Salhi-Hannachi A, Zhang H, Tu H, Mottet C, Barrès B, Degrave A, Razmjou J, Giraud T, Falque M, Dapena E, Miñarro M, Jardillier L, Deschamps P, Jousselin E, Cornille A (2021) Large-scale geographic survey provides insights into the colonization history of a major aphid pest on its cultivated apple host in Europe, North America and North Africa. bioRxiv, 2020.12.11.421644, ver. 3 peer-reviewed and recommended by Peer Community in Evolutionary Biology. https://doi.org/10.1101/2020.12.11.421644