Unravelling the invasion history of the Asian tiger mosquito in Europe

Multiple introductions are key features for the establishment and persistence of introduced species. However, little is known about the contribution of genetic admixture to the invasive potential of populations. To address this issue, we studied the recent invasion of the Asian tiger mosquito (Aedes albopictus) in Europe. Combining genome‐wide single nucleotide polymorphisms and historical knowledge using an approximate Bayesian computation framework, we reconstruct the colonization routes and establish the demographic dynamics of invasion. The colonization of Europe involved at least three independent introductions in Albania, North Italy and Central Italy that subsequently acted as dispersal centres throughout Europe. We show that the topology of human transportation networks shaped demographic histories with North Italy and Central Italy being the main dispersal centres in Europe. Introduction modalities conditioned the levels of genetic diversity in invading populations, and genetically diverse and admixed populations promoted more secondary introductions and have spread farther than single‐source invasions. This genomic study provides further crucial insights into a general understanding of the role of genetic diversity promoted by modern trade in driving biological invasions.     

Multiple introductions,admixture and bridgehead invasion characterize the introduction history of Ambrosia artemisiifolia in Europe and Australia

Admixture between differentiated populations is considered to be a powerful mechanism stimulating the invasive success of some introduced species. It is generally facilitated through multiple introductions; however, the importance of admixture prior to introduction has rarely been considered. We assess the likelihood that the invasive Ambrosia artemisiifolia populations of Europe and Australia developed through multiple introductions or were sourced from a historical admixture zone within native North America. To do this, we combine large genomic and sampling data sets analysed with approximate Bayesian computation and random forest scenario evaluation to compare single and multiple invasion scenarios with pre‐ and postintroduction admixture simultaneously. We show the historical admixture zone within native North America originated before global invasion of this weed and could act as a potential source of introduced populations. We provide evidence supporting the hypothesis that the invasive populations established through multiple introductions from the native range into Europe and subsequent bridgehead invasion into Australia. We discuss the evolutionary mechanisms that could promote invasiveness and evolutionary potential of alien species from bridgehead invasions and admixed source populations.     

Genomic admixture increases fitness during a biological invasion

Abstract During biological invasions, multiple introductions can provide opportunities for admixture among genetically distinct lineages. Admixture is predicted to contribute to invasion success by directly increasing fitness through hybrid vigour or by enhancing evolutionary potential within populations . Here, we demonstrate genome‐wide admixture during an invasion that substantially boosted fitness in the cosmopolitan weed, Silene vulgaris. We identified three divergent demes in the native European range that expanded from glacial refugia and experienced historical admixture in a well‐known suture zone. During recent invasion of North America, multiple introductions created additional opportunities for admixture. In common garden experiments, recombinant genotypes from North America experienced a two‐fold increase in fitness relative to nonrecombinants, whereas recombinant genotypes from Europe showed no lasting fitness benefits. This contrast implicates hybrid vigour behind the boost in fitness and supports the hypothesis that admixture can lead to fitness increases that may catapult invasion into a new range.     

Test of colonisation scenarios reveals complex invasion history of the red tomato spider mite Tetranychus evansi

The spider mite Tetranychus evansi is an emerging pest of solanaceous crops worldwide. Like many other emerging pests, its small size, confusing taxonomy, complex history of associations with humans, and propensity to start new populations from small inocula, make the study of its invasion biology difficult. Here, we use recent developments in Approximate Bayesian Computation (ABC) and variation in multi-locus genetic markers to reconstruct the complex historical demography of this cryptic invasive pest. By distinguishing among multiple pathways and timing of introductions, we find evidence for the "bridgehead effect", in which one invasion serves as source for subsequent invasions. Tetranychus evansi populations in Europe and Africa resulted from at least three independent introductions from South America and involved mites from two distinct sources in Brazil, corresponding to highly divergent mitochondrial DNA lineages. Mites from southwest Brazil (BR-SW) colonized the African continent, and from there Europe through two pathways in a "bridgehead" type pattern. One pathway resulted in a widespread invasion, not only to Europe, but also to other regions in Africa, southern Europe and eastern Asia. The second pathway involved the mixture with a second introduction from BR-SW leading to an admixed population in southern Spain. Admixture was also detected between invasive populations in Portugal. A third introduction from the Brazilian Atlantic region resulted in only a limited invasion in Europe. This study illustrates that ABC methods can provide insights into, and distinguish among, complex invasion scenarios. These processes are critical not only in understanding the biology of invasions, but also in refining management strategies for invasive species. For example, while reported observations of the mite and outbreaks in the invaded areas were largely consistent with estimates of geographical expansion from the ABC approach, historical observations failed to recognize the complex pathways involved and the corresponding effects on genetic diversity.     

Bridgehead Effect in the Worldwide Invasion of the Biocontrol Harlequin Ladybird

Recent studies of the routes of worldwide introductions of alien organisms suggest that many widespread invasions could have stemmed not from the native range, but from a particularly successful invasive population, which serves as the source of colonists for remote new territories. We call here this phenomenon the invasive bridgehead effect. Evaluating the likelihood of such a scenario is heuristically challenging. We solved this problem by using approximate Bayesian computation methods to quantitatively compare complex invasion scenarios based on the analysis of population genetics (microsatellite variation) and historical (first observation dates) data. We applied this approach to the Harlequin ladybird Harmonia axyridis (HA), a coccinellid native to Asia that was repeatedly introduced as a biocontrol agent without becoming established for decades. We show that the recent burst of worldwide invasions of HA followed a bridgehead scenario, in which an invasive population in eastern North America acted as the source of the colonists that invaded the European, South American and African continents, with some admixture with a biocontrol strain in Europe. This demonstration of a mechanism of invasion via a bridgehead has important implications both for invasion theory (i.e., a single evolutionary shift in the bridgehead population versus multiple changes in case of introduced populations becoming invasive independently) and for ongoing efforts to manage invasions by alien organisms (i.e., heightened vigilance against invasive bridgeheads).     

Effects of admixture in native and invasive populations of <Emphasis Type="Italic">Lythrum salicaria</Emphasis>

Intraspecific hybridization between diverged populations can enhance fitness via various genetic mechanisms. The benefits of such admixture have been proposed to be particularly relevant in biological invasions, when invasive populations originating from different source populations are found sympatrically. However, it remains poorly understood if admixture is an important contributor to plant invasive success and how admixture effects compare between invasive and native ranges. Here, we used experimental crosses in Lythrum salicaria, a species with well-established history of multiple introductions to Eastern North America, to quantify and compare admixture effects in native European and invasive North American populations. We observed heterosis in between-population crosses both in native and invasive ranges. However, invasive-range heterosis was restricted to crosses between two different Eastern and Western invasion fronts, whereas heterosis was absent in geographically distant crosses within a single large invasion front. Our results suggest that multiple introductions have led to already-admixed invasion fronts, such that experimental crosses do not further increase performance, but that contact between different invasion fronts further enhances fitness after admixture. Thus, intra-continental movement of invasive plants in their introduced range has the potential to boost invasiveness even in well-established and successfully spreading invasive species.     

Population genomic analyses reveal a history of range expansion and trait evolution across the native and invaded range of yellow starthistle (Centaurea solstitialis)

Identifying sources of genetic variation and reconstructing invasion routes for non‐native introduced species is central to understanding the circumstances under which they may evolve increased invasiveness. In this study, we used genome‐wide single nucleotide polymorphisms to study the colonization history of Centaurea solstitialis in its native range in Eurasia and invasions into the Americas. We leveraged this information to pinpoint key evolutionary shifts in plant size, a focal trait associated with invasiveness in this species. Our analyses revealed clear population genomic structure of potential source populations in Eurasia, including deep differentiation of a lineage found in the southern Apennine and Balkan Peninsulas and divergence among populations in Asia, eastern Europe and western Europe. We found strongest support for an evolutionary scenario in which western European populations were derived from an ancient admixture event between populations from eastern Europe and Asia, and subsequently served as the main genetic ‘bridgehead’ for introductions to the Americas. Introductions to California appear to be from a single source region, and multiple, independent introductions of divergent genotypes likely occurred into the Pacific Northwest. Plant size has evolved significantly at three points during range expansion, including a large size increase in the lineage responsible for the aggressive invasion of the California interior. These results reveal a long history of colonization, admixture and trait evolution in C. solstitialis, and suggest routes for improving evidence‐based management decisions for one of the most ecologically and economically damaging invasive species in the western United States.     

Complementarity of statistical treatments to reconstruct worldwide routes of invasion: the case of the Asian ladybird Harmonia axyridis

Inferences about introduction histories of invasive species remain challenging because of the stochastic demographic processes involved. Approximate Bayesian computation (ABC) can help to overcome these problems, but such method requires a prior understanding of population structure over the study area, necessitating the use of alternative methods and an intense sampling design. In this study, we made inferences about the worldwide invasion history of the ladybird Harmonia axyridis by various population genetics statistical methods, using a large set of sampling sites distributed over most of the species’ native and invaded areas. We evaluated the complementarity of the statistical methods and the consequences of using different sets of site samples for ABC inferences. We found that the H. axyridis invasion has involved two bridgehead invasive populations in North America, which have served as the source populations for at least six independent introductions into other continents. We also identified several situations of genetic admixture between differentiated sources. Our results highlight the importance of coupling ABC methods with more traditional statistical approaches. We found that the choice of site samples could affect the conclusions of ABC analyses comparing possible scenarios. Approaches involving independent ABC analyses on several sample sets constitute a sensible solution, complementary to standard quality controls based on the analysis of pseudo‐observed data sets, to minimize erroneous conclusions. This study provides biologists without expertise in this area with detailed methodological and conceptual guidelines for making inferences about invasion routes when dealing with a large number of sampling sites and complex population genetic structures.     

Using ABC and microsatellite data to detect multiple introductions of invasive species from a single source

The introduction of invasive species to new locations (that is, biological invasions) can have major impact on biodiversity, agriculture and public health. As such, determining the routes and modality of introductions with genetic data has become a fundamental goal in molecular ecology. To assist with this goal, new statistical methods and frameworks have been developed, such as approximate Bayesian computation (ABC) for inferring invasion history. Here, we present a model of invasion accounting for multiple introductions from a single source (MISS), a heretofore largely unexplored model. We simulate microsatellite data to evaluate the power of ABC to distinguish between single and multiple introductions from the same source, under a range of demographic parameters. We also apply ABC to microsatellite data from three invasions of bumblebee in New Zealand. In addition, we assess the performance of several methods of summary statistics selection. Our simulated results suggested good ability to distinguish between one- and two-wave models over much but not all of the parameter space tested, independent of summary statistics used. Globally, parameter estimation was good except for bottleneck timing. For one of the bumblebee species, we clearly rejected the MISS model, while for the other two we found inconclusive results. Since a second wave may provide genetic reinforcement to initial colonists, help relieve inbreeding among founders, or increase the hazard of the invasion, its detection may be crucial for managing invasions; we suggest that the MISS model could be considered as a potential model in future theoretical and empirical studies of invasions.     

Demographic Stochasticity, Environmental Variability, and Windows of Invasion Risk for Bythotrephes Longimanus in North America

Biological invasions are a leading threat to freshwater biodiversity worldwide. A central unanswered question of invasion ecology is why some introduced populations establish while most fail. Answering this question will allow resource managers to increase the specificity and effectiveness of control efforts and policy. We studied the establishment of spiny water flea (Bythotrephes longimanus) in the United States and Canada by modeling introduction failure caused by demographic stochasticity, environmental variation, and seasonal environmental forcing. We compared predicted establishment rates with observed invasions of inland lakes in Ontario, Canada. Our findings suggest that environmental forcing can cause “windows” of invasion opportunity so that timing of introductions might be a greater determinant of population establishment than demographic stochasticity and random environmental variation. We expect this phenomenon to be exhibited by species representing a wide range of life histories. For spiny water flea in North America, a large window of invasion opportunity opens around the fourth week of May, persists through the summer, and closes with decreasing water temperatures in autumn. These results show how timing of introductions with respect to seasonally forced environmental drivers can be a key determinant of establishment success. By focusing on introductions during windows of invasion opportunity, resource managers can more effectively control invasion rates.     

Genetic Bottlenecks in Time and Space: Reconstructing Invasions from Contemporary and Historical Collections

Herbarium accession data offer a useful historical botanical perspective and have been used to track the spread of plant invasions through time and space. Nevertheless, few studies have utilised this resource for genetic analysis to reconstruct a more complete picture of historical invasion dynamics, including the occurrence of separate introduction events. In this study, we combined nuclear and chloroplast microsatellite analyses of contemporary and historical collections of Senecio madagascariensis, a globally invasive weed first introduced to Australia c. 1918 from its native South Africa. Analysis of nuclear microsatellites, together with temporal spread data and simulations of herbarium voucher sampling, revealed distinct introductions to south-eastern Australia and mid-eastern Australia. Genetic diversity of the south-eastern invasive population was lower than in the native range, but higher than in the mid-eastern invasion. In the invasive range, despite its low resolution, our chloroplast microsatellite data revealed the occurrence of new haplotypes over time, probably as the result of subsequent introduction(s) to Australia from the native range during the latter half of the 20^th century. Our work demonstrates how molecular studies of contemporary and historical field collections can be combined to reconstruct a more complete picture of the invasion history of introduced taxa. Further, our study indicates that a survey of contemporary samples only (as undertaken for the majority of invasive species studies) would be insufficient to identify potential source populations and occurrence of multiple introductions.     

Bayesian inference of a complex invasion history revealed by nuclear and chloroplast genetic diversity in the colonizing plant,Silene latifolia

Species invading new ranges are subject to a series of demographic events that can strongly shape genetic diversity. Describing this demographic history is important for understanding where invasive species come from and how they spread, and is critical to testing hypotheses of postinvasion adaptation. Here, we analyse nuclear and chloroplast genetic diversity to study the invasion history of the widespread colonizing weed, Silene latifolia (Caryophyllaceae). Bayesian clustering and PCA revealed strong population structure in the native range of Europe, and although genotypes from multiple native sources were present in the introduced range of North America, the spatial distribution of genetic variance was dramatically reorganized. Using approximate Bayesian computation (ABC), we compared support for different invasion scenarios, including the number and size of independent introduction events and the amount of admixture occurring between sources of introduced genotypes. Our results supported independent introductions into eastern and western North America, with the latter forming a bridgehead for a secondary invasion into the Great Lakes region of central North America. Despite small estimated founder population sizes, the duration of the demographic bottleneck after the initial introduction appeared extremely short‐lived. This pattern of repeated colonization and rapid expansion has effectively eroded the strong population structure and cytonuclear associations present in Europe, but has retained overall high genetic diversity since invasion. Our results highlight the flexibility of the ABC approach for constructing a narrative of the demographic history of species invasions and provide baseline for future studies of evolutionary changes in introduced S. latifolia populations.     

Phylogeography of the widespread marine invader Microcosmus squamiger (Ascidiacea) reveals high genetic diversity of introduced populations and non-independent colonizations

The spread of non-indigenous species into new marine habitats represents an increasing threat to global diversity. Genetic techniques provide basic understanding of the invasion processes. The ascidian Microcosmus squamiger is considered to be native to Australia, having been spread worldwide via transoceanic vessels. It has successfully invaded artificial and natural habitats where it has become a pest. We studied phylogeography and genetic structure of 12 M. squamiger populations, including samples from its native range (Australia) and introduced populations from the Indian, Pacific, and Atlantic oceans, as well as the Mediterranean Sea. We amplified 574 bp of the mitochondrial COI gene in 258 individuals and found a total of 52 haplotypes. A haplotype tree revealed two main groups of haplotypes. The relative frequency of each group of haplotypes, multidimensional scaling, and analysis of molecular variance showed important differences between the western Australia localities and the remaining ones (eastern Australia and introduced populations). Furthermore, we found that the colonization of the different areas by M. squamiger has not occurred independently, as many introduced populations shared some low frequency alleles. A nested clade analysis showed a global pattern of restricted gene flow with isolation by distance, although we found episodes of long-distance dispersal in some clades. A contiguous range expansion was detected between Australian populations. We conclude that M. squamiger is native to Australia and has most likely expanded its range of distribution sequentially through worldwide shipping, especially from the harbours of the more populated eastern Australia. In introduced populations, we found a high genetic diversity which suggests enhanced invasive potential. Consequently, there is a need to control this species, as it outcompetes local biota and is an economic threat.     

The invasion route for an insect pest species: the tobacco aphid in the New World

Biological invasions are rapid evolutionary events in which populations are usually subject to a founder event during introduction followed by rapid adaptation to the new environment. Molecular tools and Bayesian approaches have shown their utility in exploring different evolutionary scenarios regarding the invasion routes of introduced species. We examined the situation for the tobacco aphid, Myzus persicae nicotianae, a recently introduced aphid species in Chile. Using seven microsatellite loci and approximate Bayesian computation, we studied populations of the tobacco aphid sampled from several American and European countries, identifying the most likely source populations and tracking the route of introduction to Chile. Our population genetic data are consistent with available historical information, pointing to an introduction route of the tobacco aphid from Europe and/or from other putative populations (e.g. Asia) with subsequent introduction through North America to South America. Evidence of multiple introductions to North America from different genetic pools, with successive loss of genetic diversity from Europe towards North America and a strong bottleneck during the southward introduction to South America, was also found. Additionally, we examined the special case of a widespread multilocus genotype that was found in all American countries examined. This case provides further evidence for the existence of highly successful genotypes or 'superclones' in asexually reproducing organisms.     

Deciphering the worldwide invasion of the Asian long‐horned beetle: A recurrent invasion process from the native area together with a bridgehead effect

Retracing introduction routes is crucial for understanding the evolutionary processes involved in an invasion, as well as for highlighting the invasion history of a species at the global scale. The Asian long‐horned beetle (ALB) Anoplophora glabripennis is a xylophagous pest native to Asia and invasive in North America and Europe. It is responsible for severe losses of urban trees, in both its native and invaded ranges. Based on historical and genetic data, several hypotheses have been formulated concerning its invasion history, including the possibility of multiple introductions from the native zone and secondary dispersal within the invaded areas, but none have been formally tested. In this study, we characterized the genetic structure of ALB in both its native and invaded ranges using microsatellites. In order to test different invasion scenarios, we used an approximate Bayesian “random forest” algorithm together with traditional population genetics approaches. The strong population differentiation observed in the native area was not geographically structured, suggesting complex migration events that were probably human‐mediated. Both native and invasive populations had low genetic diversity, but this characteristic did not prevent the success of the ALB invasions. Our results highlight the complexity of invasion pathways for insect pests. Specifically, our findings indicate that invasive species might be repeatedly introduced from their native range, and they emphasize the importance of multiple, human‐mediated introductions in successful invasions. Finally, our results demonstrate that invasive species can spread across continents following a bridgehead path, in which an invasive population may have acted as a source for another invasion.     

Genomics reveals the history of a complex plant invasion and improves the management of a biological invasion from the South African–Australian biotic exchange

Many plants exchanged in the global redistribution of species in the last 200 years, particularly between South Africa and Australia, have become threatening invasive species in their introduced range. Refining our understanding of the genetic diversity and population structure of native and alien populations, introduction pathways, propagule pressure, naturalization, and initial spread, can transform the effectiveness of management and prevention of further introductions. We used 20,221 single nucleotide polymorphisms to reconstruct the invasion of a coastal shrub, Chrysanthemoides monilifera ssp. rotundata (bitou bush) from South Africa, into eastern Australia (EAU), and Western Australia (WAU). We determined genetic diversity and population structure across the native and introduced ranges and compared hypothesized invasion scenarios using Bayesian modeling. We detected considerable genetic structure in the native range, as well as differentiation between populations in the native and introduced range. Phylogenetic analysis showed the introduced samples to be most closely related to the southern‐most native populations, although Bayesian analysis inferred introduction from a ghost population. We detected strong genetic bottlenecks during the founding of both the EAU and WAU populations. It is likely that the WAU population was introduced from EAU, possibly involving an unsampled ghost population. The number of private alleles and polymorphic SNPs successively decreased from South Africa to EAU to WAU, although heterozygosity remained high. That bitou bush remains an invasion threat in EAU, despite reduced genetic diversity, provides a cautionary biosecurity message regarding the risk of introduction of potentially invasive species via shipping routes.     

Potential limits to the benefits of admixture during biological invasion

Species introductions often bring together genetically divergent source populations, resulting in genetic admixture. This geographic reshuffling of diversity has the potential to generate favourable new genetic combinations, facilitating the establishment and invasive spread of introduced populations. Observational support for the superior performance of admixed introductions has been mixed, however, and the broad importance of admixture to invasion questioned. Under most underlying mechanisms, admixture's benefits should be expected to increase with greater divergence among and lower genetic diversity within source populations, though these effects have not been quantified in invaders. We experimentally crossed source populations differing in divergence in the invasive plant Centaurea solstitialis. Crosses resulted in many positive (heterotic) interactions, but fitness benefits declined and were ultimately negative at high source divergence, with patterns suggesting cytonuclear epistasis. We explored the literature to assess whether such negative epistatic interactions might be impeding admixture at high source population divergence. Admixed introductions reported for plants came from sources with a wide range of genetic variation, but were disproportionately absent where there was high genetic divergence among native populations. We conclude that while admixture is common in species introductions and often happens under conditions expected to be beneficial to invaders, these conditions may be constrained by predictable negative genetic interactions, potentially explaining conflicting evidence for admixture's benefits to invasion.     

Historical environmental change in Africa drives divergence and admixture of Aedes aegypti mosquitoes: a precursor to successful worldwide colonization?

Increasing globalization has promoted the spread of exotic species, including disease vectors. Understanding the evolutionary processes involved in such colonizations is both of intrinsic biological interest and important to predict and mitigate future disease risks. The Aedes aegypti mosquito is a major vector of dengue, chikungunya and Zika, the worldwide spread of which has been facilitated by Ae. aegypti's adaption to human‐modified environments. Understanding the evolutionary processes involved in this invasion requires characterization of the genetic make‐up of the source population(s). The application of approximate Bayesian computation (ABC) to sequence data from four nuclear and one mitochondrial marker revealed that African populations of Ae. aegypti best fit a demographic model of lineage diversification, historical admixture and recent population structuring. As ancestral Ae. aegypti were dependent on forests, this population history is consistent with the effects of forest fragmentation and expansion driven by Pleistocene climatic change. Alternatively, or additionally, historical human movement across the continent may have facilitated their recent spread and mixing. ABC analysis and haplotype networks support earlier inferences of a single out‐of‐Africa colonization event, while a cline of decreasing genetic diversity indicates that Ae. aegypti moved first from Africa to the Americas and then to Asia. ABC analysis was unable to verify this colonization route, possibly because the genetic signal of admixture obscures the true colonization pathway. By increasing genetic diversity and forming novel allelic combinations, divergence and historical admixture within Africa could have provided the adaptive potential needed for the successful worldwide spread of Ae. aegypti.     

Historical range expansion determines the phylogenetic diversity introduced during contemporary species invasion

For a species rapidly expanding its geographic range, such as during biological invasion, most alleles in the introduced range will have their evolutionary origins in the native range. Yet, the way in which historical processes occurring over evolutionary time in the native range contribute to the diversity sampled during contemporary invasion is largely unknown. We used chloroplast DNA (cpDNA) gene genealogies and coalescent methods to study two congeneric plants, Silene latifolia and S. vulgaris. We examined how phylogenetic diversity was shaped by demographic growth and historical range expansions in the native European range, and how this history affected the diversity sampled during their recent invasion of North America. Genealogies from both species depart from neutrality, likely as a result of demographic expansion in the ancestral range, the timing of which corresponds to shortly after each species originated. However, the species differ in the spatial distribution of cpDNA lineages across the native range. Silene latifolia shows a highly significant phylogeographic structure that most likely reflects different avenues of the post-glacial expansion into northern Europe from Mediterranean refugia. By contrast, cpDNA lineages in S. vulgaris have been widely scattered across Europe during, or since, the most recent post-glacial expansion. These different evolutionary histories resulted in dramatic differences in how phylogenetic diversity was sampled during invasion of North America. In S. latifolia, relatively few, discrete invasion events from a structured native range resulted in a rather severe genetic bottleneck, but also opportunities for admixture among previously isolated lineages. In S. vulgaris, lack of genetic structure was accompanied by more representative sampling of phylogenetic diversity during invasion, and reduced potential for admixture. Our results provide clear insights into how historical processes may feed forward to influence the phylogenetic diversity of species invading new geographic ranges.