Viruses establish intimate relationships with the cells they infect. The virocell is a novel entity, different from the original host cell and beyond the mere combination of viral and cellular genetic material. In these close encounters, viral and cellular genomes often hybridise, combine, recombine, merge and excise. Such chemical promiscuity leaves genomics scars that can be passed on to descent, in the form of deletions or duplications and, importantly, insertions and back and forth exchange of genetic material between viruses and their hosts.
In this preprint [1], Di Giovanni and coworkers report the identification of 13 genes present in the extant genomes of members of the Leptopilina wasp genus, bearing sound signatures of having been horizontally acquired from an ancestral virus. Importantly the authors identify Leptopilina boulardi filamentous virus (LbFV) as an extant relative of the ancestral virus that served as donor for the thirteen horizontally transferred genes. While pinpointing genes with a likely possible viral origin in eukaryotic genomes is only relatively rare, identifying an extant viral lineage related to the ancestral virus that continues to infect an extant relative of the ancestral host is remarkable. But the amazing evolutionary history of the Leptopilina hosts and these filamentous viruses goes beyond this shared genes. These wasps are endoparasitoids of Drosophila larvae, the female wasp laying the eggs inside the larvae and simultaneously injecting venom that hinders the immune response. The composition of the venoms is complex, varies between wasp species and also between individuals within a species, but a central component of all these venoms are spiked structures that vary in morphology, symmetry and size, often referred to as virus-like particles (VLPs).
In this preprint, the authors convincingly show that the expression pattern in the Leptopilina wasps of the thirteen genes identified to have been horizontally acquired from the LbFV ancestor coincides with that of the production of VLPs in the female wasp venom gland. Based on this spatio-temporal match, the authors propose that these VLPs have a viral origin. The data presented in this preprint will undoubtedly stimulate further research on the composition, function, origin, evolution and diversity of these VLP structures, which are highly debated (see for instance [2] and [3]).

References

[1] Di Giovanni, D., Lepetit, D., Boulesteix, M., Ravallec, M., & Varaldi, J. (2018). A behavior-manipulating virus relative as a source of adaptive genes for parasitoid wasps. bioRxiv, 342758, ver. 5 peer-reviewed and recommended by PCI Evol Biol. doi: 10.1101/342758
[2] Poirié, M., Colinet, D., & Gatti, J. L. (2014). Insights into function and evolution of parasitoid wasp venoms. Current Opinion in Insect Science, 6, 52-60. doi: 10.1016/j.cois.2014.10.004
[3] Heavner, M. E., Ramroop, J., Gueguen, G., Ramrattan, G., Dolios, G., Scarpati, M., ... & Govind, S. (2017). Novel organelles with elements of bacterial and eukaryotic secretion systems weaponize parasites of Drosophila. Current Biology, 27(18), 2869-2877. doi: 10.1016/j.cub.2017.08.019

Accurate information flow is central to living systems. The continuity of genomes through generations as well as the reproducible functioning and survival of the individual organisms require a faithful information transfer during replication, transcription and translation. The differential efficiency of natural selection against “mistakes” results in decreasing fidelity rates for replication, transcription and translation. At each level in the information flow chain (replication, transcription, translation), numerous complex molecular systems have evolved and been selected for preventing, identifying and, when possible, correcting or removing such “mistakes” arising during information transfer.

However, fidelity cannot be improved ad infinitum. First, because of the limits imposed by the physical nature of the processes of copying and recoding information over different molecular supports: all mechanisms ensuring fidelity during biological information transfer ultimately rely on chemical kinetics and thermodynamics. The more accurate a copying process is, the lower the synthesis rate and the higher the energetic cost of correcting errors. Second, because of the limits imposed by random genetic drift: natural selection cannot effectively act on an allele that contributes with a small differential advantage unless effective population size is large. If s <1/Ne (or s <1/(2Ne) in diploids) the allele frequency in the population is de facto subject to neutral drift processes.

In their preprint “Random genetic drift sets an upper limit on mRNA splicing accuracy in metazoans”, Bénitière, Necsulea and Duret explore the validity of this last mentioned “drift barrier” hypothesis for the case study of alternative splicing diversity in eukaryotes (Bénitière et al. 2022). Splicing refers to an ensemble of eukaryotic molecular processes mediated by a large number of proteins and ribonucleoproteins and involving nucleotide sequence recognition, that uses as a molecular substrate a precursor messenger RNA (mRNA), directly transcribed from the DNA, and produces a mature mRNA by removing introns and joining exons (Chow et al. 1977). Alternative splicing refers to the case in which different molecular species of mature mRNAs can be produced, either by cis-splicing processes acting on the same precursor mRNA, e.g. by varying the presence/absence of different exons or by varying the exon-exon boundaries, or by trans-splicing processes, joining exons from different precursor mRNA molecules.

The diversity of mRNA molecular species generated by alternative splicing enlarges the molecular phenotypic space that can be generated from the same genotype. In humans, alternative splicing occurs in around 95% of the ca. 20,000 genes, resulting in ca. 100,000 medium-to-high abundance transcripts (Pan et al. 2008). In multicellular organisms, the frequency of alternatively spliced mRNAs varies between tissues and across ontogeny, often in a switch-like pattern (Wang et al. 2008). In the molecular and cell biology community, it is commonly accepted that splice variants contribute with specific functions (Marasco and Kornblihtt 2023) although there exists a discussion around the functional nature of low-frequency splice variants (see for instance the debate between Tress et al. 2017 and Blencowe 2017). The origin, diversity, regulation and evolutionary advantage of alternative splicing constitutes thus a playground of the selectionist-neutralist debate, with one extreme considering that most splice variants are mere “mistakes” of the splicing process (Pickrell et al. 2010), and the other extreme considering that alternative splicing is at the core of complexity in multicellular organisms, as it increases the genome coding potential and allows for a large repertoire of cell types (Chen et al. 2014).

In their manuscript, Bénitière, Necsulea and Duret set the cursor towards the neutralist end of the gradient and test the hypothesis of whether the high alternative splice rate in “complex” organisms corresponds to a high rate of splicing “mistakes”, arising from the limit imposed by the drift barrier effect on the power of natural selection to increase accuracy (Bush et al. 2017). In their preprint, the authors convincingly show that in metazoans a fraction of the variation of alternative splicing rate is explained by variation in proxies of population size, so that species with smaller Ne display higher alternative splice rates. They communicate further that abundant splice variants tend to preserve the reading frame more often than low-frequency splice variants, and that the nucleotide splice signals in abundant splice variants display stronger evidence of purifying selection than those in low-frequency splice variants. From all the evidence presented in the manuscript, the authors interpret that “variation in alternative splicing rate is entirely driven by variation in the efficacy of selection against splicing errors”.

The authors honestly present some of the limitations of the data used for the analyses, regarding i) the quality of the proxies used for Ne (i.e. body length, longevity and dN/dS ratio); ii) the heterogeneous nature of the RNA sequencing datasets (full organisms, organs or tissues; different life stages, sexes or conditions); and iii) mostly short RNA reads that do not fully span individual introns. Further, data from bacteria do not verify the herein communicated trends, as it has been shown that bacterial species with low population sizes do not display higher transcription error rates (Traverse and Ochman 2016). Finally, it will be extremely interesting to introduce a larger evolutionary perspective on alternative splicing rates encompassing unicellular eukaryotes, in which an intriguing interplay between alternative splicing and gene duplication has been communicated (Hurtig et al. 2020).

The manuscript from Bénitière, Necsulea and Duret makes a significant advance to our understanding of the diversity, the origin and the physiology of post-transcriptional and post-translational mechanisms by emphasising the fundamental role of non-adaptive evolutionary processes and the upper limits to splicing accuracy set by genetic drift.

References

Bénitière F, Necsulea A, Duret L. 2023. Random genetic drift sets an upper limit on mRNA splicing accuracy in metazoans. bioRxiv, ver. 4 peer-reviewed and recommended by Peer Community in Evolutionary Biology. https://doi.org/10.1101/2022.12.09.519597 

Blencowe BJ. 2017. The Relationship between Alternative Splicing and Proteomic Complexity. Trends Biochem Sci 42:407–408. https://doi.org/10.1016/j.tibs.2017.04.001

Bush SJ, Chen L, Tovar-Corona JM, Urrutia AO. 2017. Alternative splicing and the evolution of phenotypic novelty. Philos Trans R Soc Lond B Biol Sci 372:20150474. https://doi.org/10.1098/rstb.2015.0474

Chen L, Bush SJ, Tovar-Corona JM, Castillo-Morales A, Urrutia AO. 2014. Correcting for differential transcript coverage reveals a strong relationship between alternative splicing and organism complexity. Mol Biol Evol 31:1402–1413. https://doi.org/10.1093/molbev/msu083

Chow LT, Gelinas RE, Broker TR, Roberts RJ. 1977. An amazing sequence arrangement at the 5’ ends of adenovirus 2 messenger RNA. Cell 12:1–8. https://doi.org/10.1016/0092-8674(77)90180-5

Hurtig JE, Kim M, Orlando-Coronel LJ, Ewan J, Foreman M, Notice L-A, Steiger MA, van Hoof A. 2020. Origin, conservation, and loss of alternative splicing events that diversify the proteome in Saccharomycotina budding yeasts. RNA 26:1464–1480. https://doi.org/10.1261/rna.075655.120

Marasco LE, Kornblihtt AR. 2023. The physiology of alternative splicing. Nat Rev Mol Cell Biol 24:242–254. https://doi.org/10.1038/s41580-022-00545-z

Pan Q, Shai O, Lee LJ, Frey BJ, Blencowe BJ. 2008. Deep surveying of alternative splicing complexity in the human transcriptome by high-throughput sequencing. Nat Genet 40:1413–1415. https://doi.org/10.1038/ng.259

Pickrell JK, Pai AA, Gilad Y, Pritchard JK. 2010. Noisy splicing drives mRNA isoform diversity in human cells. PLoS Genet 6:e1001236. https://doi.org/10.1371/journal.pgen.1001236

Traverse CC, Ochman H. 2016. Conserved rates and patterns of transcription errors across bacterial growth states and lifestyles. Proc Natl Acad Sci U S A 113:3311–3316. https://doi.org/10.1073/pnas.1525329113

Tress ML, Abascal F, Valencia A. 2017. Alternative Splicing May Not Be the Key to Proteome Complexity. Trends Biochem Sci 42:98–110. https://doi.org/10.1016/j.tibs.2016.08.008

Wang ET, Sandberg R, Luo S, Khrebtukova I, Zhang L, Mayr C, Kingsmore SF, Schroth GP, Burge CB. 2008. Alternative isoform regulation in human tissue transcriptomes. Nature 456:470–476. https://doi.org/10.1038/nature07509

Duplications, mutations and recombination may be considered the main sources of genomic variation and evolution. In addition, sexual recombination is essential in purging deleterious mutations and allowing advantageous allelic combinations to occur (Glémin et al. 2019). However, in parthenogenetic asexual organisms, variation cannot be explained by sexual recombination, and other mechanisms must account for it. Although it is known that transposable elements (TE) may influence on genome structure and gene expression patterns, their role as a primary source of genomic variation and rapid adaptability has received less attention. An important role of TE on adaptive genome evolution has been documented for fungal phytopathogens (Faino et al. 2016), suggesting that TE activity might explain the evolutionary dynamics of this type of organisms.
The phytopathogen nematode Meloidogyne incognita is one of the worst agricultural pests in warm climates (Savary et al. 2019). This species, as well as other root-knot nematodes (RKN), shows a wide geographical distribution range infecting diverse groups of plants. Although allopolyploidy may have played an important role on the wide adaptation of this phytopathogen, it may not explain by itself the rapid changes required to overcome plant resistance in a few generations. Paradoxically, M. incognita reproduces asexually via mitotic parthenogenesis (Trudgill and Blok 2001; Castagnone-Sereno and Danchin 2014) and only few single nucleotide variations were identified between different host races isolates (Koutsovoulos et al. 2020). Therefore, this is an interesting model to explore other sources of genomic variation such the TE activity and its role on the success and adaptability of this phytopathogen.
To address these questions, Kozlowski et al. (2020) estimated the TE mobility across 12 geographical isolates that presented phenotypic variations in Meloidogyne incognita, concluding that recent activity of TE in both genic and regulatory regions might have given rise to relevant functional differences between genomes. This was the first estimation of TE activity as a mechanism probably involved in genome plasticity of this root-knot nematode. This study also shed light on evolutionary mechanisms of asexual organisms with an allopolyploid origin. These authors re-annotated the 185 Mb triploid genome of M. incognita for TE content analysis using stringent filters (Kozlowski 2020a), and estimated activity by their distribution using a population genomics approach including isolates from different crops and locations. Canonical TE represented around 4.7% of the M. incognita genome of which mostly correspond to TIR (Terminal Inverted Repeats) and MITEs (Miniature Inverted repeat Transposable Elements) followed by Maverick DNA transposons and LTR (Long Terminal Repeats) retrotransposons. The result that most TE found were represented by DNA transposons is similar to the previous studies with the nematode species model Caenorhabditis elegans (Bessereau 2006; Kozlowski 2020b) and other nematodes as well. Canonical TE annotations were highly similar to their consensus sequences containing transposition machinery when TE are autonomous, whereas no genes involved in transposition were found in non-autonomous ones. These findings suggest recent activity of TE in the M. incognita genome. Other relevant result was the significant variation in TE presence frequencies found in more than 3,500 loci across isolates, following a bimodal distribution within isolates. However, variation in TE frequencies was low to moderate between isolates recapitulating the phylogenetic signal of isolates DNA sequences polymorphisms. A detailed analysis of TE frequencies across isolates allowed identifying polymorphic TE loci, some of which might be neo-insertions mostly of TIRs and MITEs (Kozlowski 2020c). Interestingly, the two thirds of the fixed neo-insertions were located in coding regions or in regulatory regions impacting expression of specific genes in M. incognita. Future research on proteomics is needed to evaluate the functional impact that these insertions have on adaptive evolution in M. incognita. In this line, this pioneer research of Kozlowski et al. (2020) is a first step that is also relevant to remark the role that allopolyploidy and reproduction have had on shaping nematode genomes.

References

[1] Bessereau J-L. 2006. Transposons in C. elegans. WormBook. 10.1895/wormbook.1.70.1
[2] Castagnone-Sereno P, Danchin EGJ. 2014. Parasitic success without sex - the nematode experience. J. Evol. Biol. 27:1323-1333. 10.1111/jeb.12337
[3] Faino L, Seidl MF, Shi-Kunne X, Pauper M, Berg GCM van den, Wittenberg AHJ, Thomma BPHJ. 2016. Transposons passively and actively contribute to evolution of the two-speed genome of a fungal pathogen. Genome Res. 26:1091-1100. 10.1101/gr.204974.116
[4] Glémin S, François CM, Galtier N. 2019. Genome Evolution in Outcrossing vs. Selfing vs. Asexual Species. In: Anisimova M, editor. Evolutionary Genomics: Statistical and Computational Methods. Methods in Molecular Biology. New York, NY: Springer. p. 331-369. 10.1007/978-1-4939-9074-0_11
[5] Koutsovoulos GD, Marques E, Arguel M-J, Duret L, Machado ACZ, Carneiro RMDG, Kozlowski DK, Bailly-Bechet M, Castagnone-Sereno P, Albuquerque EVS, et al. 2020. Population genomics supports clonal reproduction and multiple independent gains and losses of parasitic abilities in the most devastating nematode pest. Evol. Appl. 13:442-457. 10.1111/eva.12881
[6] Kozlowski D. 2020a. Transposable Elements prediction and annotation in the M. incognita genome. Portail Data INRAE. 10.15454/EPTDOS
[7] Kozlowski D. 2020b. Transposable Elements prediction and annotation in the C. elegans genome. Portail Data INRAE. 10.15454/LQCIW0
[8] Kozlowski D. 2020c. TE polymorphisms detection and analysis with PopoolationTE2. Portail Data INRAE. 10.15454/EWJCT8
[9] Kozlowski DK, Hassanaly-Goulamhoussen R, Da Rocha M, Koutsovoulos GD, Bailly-Bechet M, Danchin EG (2020) Transposable Elements are an evolutionary force shaping genomic plasticity in the parthenogenetic root-knot nematode Meloidogyne incognita. bioRxiv, 2020.04.30.069948, ver. 4 peer-reviewed and recommended by PCI Evolutionary Biology. 10.1101/2020.04.30.069948
[10] Savary S, Willocquet L, Pethybridge SJ, Esker P, McRoberts N, Nelson A. 2019. The global burden of pathogens and pests on major food crops. Nat. Ecol. Evol. 3:430-439. 10.1038/s41559-018-0793-y
[11] Trudgill DL, Blok VC. 2001. Apomictic, polyphagous root-knot nematodes: exceptionally successful and damaging biotrophic root pathogens. Annu Rev Phytopathol 39:53-77. 10.1146/annurev.phyto.39.1.53

The increasing number of reported cases of antibiotic resistance is one of today’s major public health concerns. Dealing with this threat involves understanding what drives the evolution of antibiotic resistance and investigating whether we can predict (and subsequently avoid or circumvent) it [1].
One of the most illustrative and common models of adaptation (and, hence, resistance evolution) is Fisher’s Geometric Model (FGM). The original model maps phenotypes to fitness, meaning that each point in the fitness landscape corresponds to a phenotype rather than a genotype. However, it has been shown that when mutations are numerous enough, FGM can also describe adaptive walks in genotype space [2]. Nevertheless, limitations have been highlighted, particularly when trying to study complex scenarios such as antibiotic resistance evolution [3].
Harmand et al. [4] incorporated three extensions to the FGM, which allowed them to match the mutational patterns of antibiotic resistance that they obtained from a screen across a gradient of drug concentrations. The implemented extensions took into account that: 1) only a subset of mutations may contribute to traits under selection, reflecting that not all regions in the genome affect the ability to resist antibiotics; 2) mutations that confer a fitness increase in one environment may not reflect a similar increase in others, if the selective constraints are different; and 3) different antibiotic concentrations may either constrain the maximum fitness that populations can reach (changing the height of the fitness peak) or change the rate of fitness increase with each mutation (changing the width/slope of the peak).
Traditionally, most empirical fitness landscape studies have focused on a subset of mutations obtained after laboratory evolution in specific conditions [5, 6]. The results obtained in Harmand et al. [4] indicate a potential shortcoming of studying these small fitness landscapes: rather than having a constrained evolutionary path to a resistant phenotype, as previously observed, their results suggest that antibiotic resistance can be the product of mutations in different regions of the genome. Returning to the fitness landscape perspective, this indicates that there are many alternative paths that can lead to the evolution of antibiotic resistance. This comparison points at a difficult challenge when aiming at developing a predictive framework for evolution: real-time experiments may indicate that evolution is likely to take similar and predictable paths because the strongest and most frequent mutations dictate the outcome, whereas systematic screens of mutants potentially indicate several paths, that may, however, not be relevant in nature. Only a combination of different experimental approaches with motivated theory as presented in Harmand et al. [4] will allow for a better understanding of where in this continuum evolution is taking place in nature, and to which degree we are able to interfere with it in order to slow down adaptation.

References

[1] Palmer AC, and Kishony R. 2013. Understanding, predicting and manipulating the genotypic evolution of antibiotic resistance. Nature Review Genetics 14: 243—248. doi: 10.1038/nrg3351

[2] Tenaillon O. 2014. The utility of Fisher’s geometric model in evolutionary genetics. Annual Review of Ecology, Evolution and Systematics 45: 179—201. doi: 10.1146/annurev-ecolsys-120213-091846

[3] Blanquart F and Bataillon T. 2016. Epistasis and the structure of fitness landscapes: are experimental fitness landscapes compatible with Fisher’s geometric model? Genetics 203: 847—862. doi: 10.1534/genetics.115.182691

[4] Harmand N, Gallet R, Jabbour-Zahab R, Martin G and Lenormand T. 2017. Fisher’s geometrical model and the mutational patterns of antibiotic resistance across dose gradients. Evolution 71: 23—37. doi: 10.1111/evo.13111

[5] de Visser, JAGM, and Krug J. 2014. Empirical fitness landscapes and the predictability of evolution. Nature 15: 480—490. doi: 10.1038/nrg3744

[6] Palmer AC, Toprak E, Baym M, Kim S, Veres A, Bershtein S and Kishony R. 2015. Delayed commitment to evolutionary fate in antibiotic resistance fitness landscapes. Nature Communications 6: 1—8. doi: 10.1038/ncomms8385

Biological invasions are natural experiments, and often show that evolution can affect dynamics in important ways [1-3]. While we often think of invasions as a conservation problem stemming from anthropogenic introductions [4,5], biological invasions are much more commonplace than this, including phenomena as diverse as natural range shifts, the spread of novel pathogens, and the growth of tumors. A major question across all these settings is which set of traits determine the ability of a population to invade new space [6,7]. Traits such as: increased growth or reproductive rate, dispersal ability and ability to defend from predation often show large evolutionary shifts across invasion history [1,6,8]. Are such multi-trait shifts driven by selection on multiple traits, or a correlated response by multiple traits to selection on one? Resolving this question is important for both theoretical and practical reasons [9,10]. But despite the importance of this issue, it is not easy to perform the necessary manipulative experiments [9].
Foucaud et al. [11] tackled this issue by performing experimental evolution on source populations of the invasive ladybug Harmonia axyridis. The authors tested if selection on a single trait could generate correlated responses in other life history traits. Specifically, they used experimental evolution to impose divergent selection on female mass, and reproductive timing. After ten generations, they found that selection for weight did not affect almost any other life history trait. However, nine generations of selection for faster reproduction led to correlated phenotypic changes in developmental, reproduction and survival rate of populations, although not always in the direction we might have expected. Despite this correlated response, none of their selected lines were able to fully recapitulate the trait shifts seen in natural invasions of this species. This implies that selection during natural invasions is operating on multiple traits; a finding in agreement with our growing understanding of how selection acts during introduction and invasion [12,13].
Populations undergoing a colonization process may also be subject to a multitude of different selective pressures [14,15]. The authors expanded their work in this direction by testing whether food availability alters the observed correlations between life history traits. The pervasiveness of genotype by environment interactions observed also points to a role for multiple selective pressures in shaping the suite of life-history shifts observed in wild ladybug populations. The work from Foucaud and colleagues [11] adds to a small but growing list of important studies that use experimental evolution to investigate how life-history traits evolve, and how they evolve during invasions in particular.

References

[1] Sakai, A.K., Allendorf, F.W., Holt, J.S. et al. (2001). The population biology of invasive species. Annual review of ecology and systematics, 32(1), 305-332. doi: 10.1146/annurev.ecolsys.32.081501.114037
[2] Hairston Jr, N. G., Ellner, S. P., Geber, M. A., Yoshida, T. and Fox, J. A. (2005). Rapid evolution and the convergence of ecological and evolutionary time. Ecology letters, 8(10), 1114-1127. doi: 10.1111/j.1461-0248.2005.00812.x
[3] Chuang, A. and Peterson, C. R. (2016). Expanding population edges: theories, traits, and trade‐offs. Global change biology, 22(2), 494-512. doi: 10.1111/gcb.13107
[4] Whitney, K. D. and Gabler, C. A. (2008). Rapid evolution in introduced species,‘invasive traits’ and recipient communities: challenges for predicting invasive potential. Diversity and Distributions, 14(4), 569-580. doi: 10.1111/j.1472-4642.2008.00473.x
[5] Catullo, R. A., Llewelyn, J., Phillips, B. L. and Moritz, C. C. (2019). The Potential for Rapid Evolution under Anthropogenic Climate Change. Current Biology, 29(19), R996-R1007. doi: 10.1016/j.cub.2019.08.028
[6] Suarez, A. V. and Tsutsui, N. D. (2008). The evolutionary consequences of biological invasions. Molecular Ecology, 17(1), 351-360. doi: 10.1111/j.1365-294X.2007.03456.x
[7] Deforet, M., Carmona-Fontaine, C., Korolev, K. S. and Xavier, J. B. (2019). Evolution at the edge of expanding populations. The American Naturalist, 194(3), 291-305. doi: 10.1086/704594
[8] Phillips, B. L., Brown, G. P., and Shine, R. (2010). Life‐history evolution in range‐shifting populations. Ecology, 91(6), 1617-1627. doi: 10.1890/09-0910.1
[9] Colautti, R. I. and Lau, J. A. (2015). Contemporary evolution during invasion: evidence for differentiation, natural selection, and local adaptation. Molecular ecology, 24(9), 1999-2017. doi: 10.1111/mec.13162
[10] Szűcs, M., Melbourne, B. A., Tuff, T., Weiss‐Lehman, C. and Hufbauer, R. A. (2017). Genetic and demographic founder effects have long‐term fitness consequences for colonising populations. Ecology Letters, 20(4), 436-444. doi: 10.1111/ele.12743
[11] Foucaud, J., Hufbauer, R. A., Ravigné, V., Olazcuaga, L., Loiseau, A., Ausset, A., Wang, S., Zang, L.-S., Lemenager, N., Tayeh, A., Weyna, A., Gneux, P., Bonnet, E., Dreuilhe, V., Poutout, B., Estoup, A. and Facon, B. (2020). How do invasion syndromes evolve? An experimental evolution approach using the ladybird Harmonia axyridis. bioRxiv, 849968 ver. 4 peer-reviewed and recommended by PCI Evolutionary Biology. doi: 10.1101/849968
[12] Simons, A. M. (2003). Invasive aliens and sampling bias. Ecology Letters, 6(4), 278-280. doi: 10.1046/j.1461-0248.2003.00430.x
[13] Phillips, B. L. and Perkins, T. A. (2019). Spatial sorting as the spatial analogue of natural selection. Theoretical Ecology, 12(2), 155-163. doi: 10.1007/s12080-019-0412-9
[14] Lavergne, S. and Molofsky, J. (2007). Increased genetic variation and evolutionary potential drive the success of an invasive grass. Proceedings of the National Academy of Sciences, 104(10), 3883-3888. doi: 10.1073/pnas.0607324104
[15] Moran, E. V. and Alexander, J. M. (2014). Evolutionary responses to global change: lessons from invasive species. Ecology Letters, 17(5), 637-649. doi: 10.1111/ele.12262