Worldwide Alien Invasion: A Methodological Approach to Forecast the Potential Spread of a Highly Invasive Pollinator

The ecological impacts of alien species invasion are a major threat to global biodiversity. The increasing number of invasion events by alien species and the high cost and difficulty of eradicating invasive species once established require the development of new methods and tools for predicting the most susceptible areas to invasion. Invasive pollinators pose serious threats to biodiversity and human activity due to their close relationship with many plants (including crop species) and high potential competitiveness for resources with native pollinators. Although at an early stage of expansion, the bumblebee species Bombus terrestris is becoming a representative case of pollinator invasion at a global scale, particularly given its high velocity of invasive spread and the increasing number of reports of its impacts on native bees and crops in many countries. We present here a methodological framework of habitat suitability modeling that integrates new approaches for detecting habitats that are susceptible to Bombus terrestris invasion at a global scale. Our approach did not include reported invaded locations in the modeling procedure; instead, those locations were used exclusively to evaluate the accuracy of the models in predicting suitability over regions already invaded. Moreover, a new and more intuitive approach was developed to select the models and evaluate different algorithms based on their performance and predictive convergence. Finally, we present a comprehensive global map of susceptibility to Bombus terrestris invasion that highlights priority areas for monitoring.     

Substrate preferences of coexisting invasive amphipods, <Emphasis Type="Italic">Dikerogammarus villosus</Emphasis> and <Emphasis Type="Italic">Dikerogammarus haemobaphes</Emphasis>, under field and laboratory conditions

Two Ponto-Caspian amphipods, Dikerogammarus villosus and Dikerogammarus haemobaphes, have expanded their geographical ranges from eastern Europe into Great Britain in recent years. This study represents one of the first examining the distribution and habitat preferences of coexisting populations of D. haemobaphes and D. villosus via field and laboratory experiments in the UK. Field surveys of a recently invaded lowland reservoir in the UK are complimented with ex situ laboratory mesocosm experiments examining the substrate preferences of coexisting populations of D. villosus and D. haemobaphes. Results from the field study indicated that D. haemobaphes dominated the macroinvertebrate community within the reservoir and demonstrated a strong affinity for large cobble and artificial substrates. D. villosus occurred at lower abundances but displayed a strong preference for coarse cobble substrates. A third invasive amphipod, Crangonyx pseudogracilis, was largely confined to sand/silt habitats. Laboratory mesocosm experiments clearly supported the field observations of D. villosus and D. haemobaphes with both species demonstrating a preference for cobble substrates. Results from the study highlight the importance of characterising physical habitat when investigating biological invasions and suggest that habitat availability may influence the extent and speed at which range expansion of new amphipod invaders occurs.     

Effectiveness of zebra mussels to act as shelters from fish predators differs between native and invasive amphipod prey

Biological invasions cause organisms to face new predators, but also supply new anti-predator shelters provided by alien ecosystem engineers. We checked the level of anti-predator protection provided to three gammarid species by an invasive Ponto-Caspian zebra mussel Dreissena polymorpha, known for its habitat modification abilities. We used gammarids differing in their origin and level of association with mussels: Ponto-Caspian aliens Dikerogammarus villosus (commonly occurring in mussel beds) and Pontogammarus robustoides (not associated with mussels), as well as native European Gammarus fossarum (not co-occurring with dreissenids). The gammarids were exposed to predation of two fish species: the racer goby Babka gymnotrachelus (Ponto-Caspian) and Amur sleeper Perccottus glenii (Eastern Asian). This set of organisms allowed us to check whether the origin and level of association with mussels of both prey and predators affect the ability of gammarids to utilize zebra mussel beds as shelters. We tested gammarid survival in the presence of fish and one of five substrata: sand, macrophytes, stones, living mussels and empty mussel valves. D. villosus survived better than its congeners on all substrata, and its survival was highest in living dreissenids. The survival of the other gammarids was similar on all substrata. Both fish species exhibited similar predation efficiency. Thus, D. villosus, whose affinity to dreissenids has already been established, utilizes them as protection from fish predators, including allopatric predators, more efficiently than other amphipods. Therefore, the presence of dreissenids in areas invaded by D. villosus is likely to help the invader establish itself in a new place.     

Different ammonia tolerances may facilitate spatial coexistence of Gammarus roeselii and the strong invader Dikerogammarus villosus

The introduction of invasive species that can replace native species is one of the most critical threats to the biodiversity of aquatic systems. Here we investigated the potential contribution of one factor to the coexistence of the indigenous amphipod Gammarus roeselii and the invasive amphipod Dikerogammarus villosus in the same ecosystem (Lake Constance) within different microhabitats. We quantitatively studied the influence of ambient ammonia concentrations on the distributions of the two amphipod species. We also assessed the ammonia tolerance ranges of both species in laboratory experiments by measuring mortality rate, precopula disruption, egg mortality, and microhabitat choice. The proportion of G. roeselii among the two amphipod species was significantly positively related to the ammonia concentration in the water, which indicated that the distribution of the invasive D. villosus was limited at high ammonia concentrations. Although the mortality rates of the two species did not significantly differ, G. roeselii was more tolerant to ammonia with regard to precopula disruption, egg mortality, and microhabitat choice. The effective ammonia concentrations that led to a significantly reduced direct reproductive success in D. villosus were within the range of the highest field concentrations measured, where only G. roeselii occurred. D. villosus may have a smaller range than the indigenous G. roeselii partially because of its lower tolerance to higher ammonia concentrations, which lead to reduced reproductive success. Beside other habitat parameters differences on ammonia tolerance between the two amphipods might allow their coexistence along a gradient of microhabitats in Lake Constance.     

Getting What Is Served? Feeding Ecology Influencing Parasite-Host Interactions in Invasive Round Goby Neogobius melanostomus

Freshwater ecosystems are increasingly impacted by alien invasive species which have the potential to alter various ecological interactions like predator-prey and host-parasite relationships. Here, we simultaneously examined predator-prey interactions and parasitization patterns of the highly invasive round goby (Neogobius melanostomus) in the rivers Rhine and Main in Germany. A total of 350 N. melanostomus were sampled between June and October 2011. Gut content analysis revealed a broad prey spectrum, partly reflecting temporal and local differences in prey availability. For the major food type (amphipods), species compositions were determined. Amphipod fauna consisted entirely of non-native species and was dominated by Dikerogammarus villosus in the Main and Echinogammarus trichiatus in the Rhine. However, the availability of amphipod species in the field did not reflect their relative abundance in gut contents of N. melanostomus. Only two metazoan parasites, the nematode Raphidascaris acus and the acanthocephalan Pomphorhynchus sp., were isolated from N. melanostomus in all months, whereas unionid glochidia were only detected in June and October in fish from the Main. To analyse infection pathways, we examined 17,356 amphipods and found Pomphorhynchus sp. larvae only in D. villosus in the river Rhine at a prevalence of 0.15%. Dikerogammarus villosus represented the most important amphipod prey for N. melanostomus in both rivers but parasite intensities differed between rivers, suggesting that final hosts (large predatory fishes) may influence host-parasite dynamics of N. melanostomus in its introduced range.     

Assessing the suitable habitat for reintroduction of brown trout (Salmo trutta forma fario) in a lowland river: A modeling approach

Huge efforts have been made during the past decades to improve the water quality and to restore the physical habitat of rivers and streams in western Europe. This has led to an improvement in biological water quality and an increase in fish stocks in many countries. However, several rheophilic fish species such as brown trout are still categorized as vulnerable in lowland streams in Flanders (Belgium). In order to support cost‐efficient restoration programs, habitat suitability modeling can be used. In this study, we developed an ensemble of habitat suitability models using metaheuristic algorithms to explore the importance of a large number of environmental variables, including chemical, physical, and hydromorphological characteristics to determine the suitable habitat for reintroduction of brown trout in the Zwalm River basin (Flanders, Belgium), which is included in the Habitats Directive. Mean stream velocity, water temperature, hiding opportunities, and presence of pools or riffles were identified as the most important variables determining the habitat suitability. Brown trout mainly preferred streams with a relatively high mean reach stream velocity (0.2–1 m/s), a low water temperature (7–15°C), and the presence of pools. The ensemble of models indicated that most of the tributaries and headwaters were suitable for the species. Synthesis and applications. Our results indicate that this modeling approach can be used to support river management, not only for brown trout but also for other species in similar geographical regions. Specifically for the Zwalm River basin, future restoration of the physical habitat, removal of the remaining migration barriers and the development of suitable spawning grounds could promote the successful restoration of brown trout.     

Habitat properties and plant traits interact as drivers of non‐native plant species’ seed production at the local scale

To answer the long‐standing question if we can predict plant invader success based on characteristics of the environment (invasibility) or the invasive species (invasiveness), or the combination of both, there is a need for detailed observational studies in which habitat properties, non‐native plant traits, and the resulting invader success are locally measured. In this study, we assess the interaction of gradients in the environmental and trait space on non‐native species fitness, expressed as seed production, for a set of 10 invasive and noninvasive non‐native species along a wide range of invaded sites in Flanders. In our multidimensional approach, most of the single environmental gradients (temperature, light availability, native plant species diversity, and soil fertility) and sets of non‐native plant traits (plant size, photosynthesis, and foliar chemical attributes) related positively with invader seed production. Yet correlation with seed production was much stronger when several environmental gradients were assessed in interaction, and even more so when we combined plant traits and habitat properties. The latter increased explanatory power of the models on average by 25% for invasive and by 7% for noninvasive species. Additionally, we report a 70‐fold higher seed production in invasive than in noninvasive species and fundamentally different correlations of seed production with plant traits and habitat properties in noninvasive versus invasive species. We conclude that locally measured traits and properties deserve much more attention than they currently get in invasion literature and thus encourage further studies combining this level of detail with the generality of a multiregion and multispecies approach across different stages of invasion.     

Predicting microhabitat selection in juvenile Atlantic salmon Salmo salar by the use of logistic regression and classification trees

1. We compared the capacity of logistic regression (LR) and classification tree (CT) models to predict microhabitat use and the summer distribution of juvenile Atlantic salmon, Salmo salar, in two reaches of a small stream in eastern Quebec. 2. The models predicted the presence or absence of salmon at a location on the basis of habitat features (depth, current velocity, presence of instream and overhead cover, substratum particle size, and distance to stream bank) measured at that location. Models were validated by means of crossover field tests evaluating the performance of models developed for one reach (calibration trials) when applied to the other reach (validation trials). Model performance was evaluated with regard to accuracy, generality and ease of use and interpretation. Prediction maps based on habitat features were also built to compare the observed position of fish with those predicted by LR and CT models. 3. The spatial distribution of active fish differed markedly from that of resting fish, apparently as a result of the selection for water greater than about 30 cm depth by active fish and for the presence of rocky cover by resting fish. 4. All models made accurate predictions, validated by crossover trials. For both LR and CT models, the prediction maps reflected well the actual fish distributions. However, CT models were easier to build and interpret than LR models. CT models also had less variable performance and a smaller decline in predictive capability in crossover trials (for fish at rest), suggesting that they may be more transferable than LR models.     

Simulating the effects of climate change on the distribution of an invasive plant,using a high resolution,local scale,mechanistic approach: challenges and insights

The growing economic and ecological damage associated with biological invasions, which will likely be exacerbated by climate change, necessitates improved projections of invasive spread. Generally, potential changes in species distribution are investigated using climate envelope models; however, the reliability of such models has been questioned and they are not suitable for use at local scales. At this scale, mechanistic models are more appropriate. This paper discusses some key requirements for mechanistic models and utilises a newly developed model (PSS[gt]) that incorporates the influence of habitat type and related features (e.g., roads and rivers), as well as demographic processes and propagule dispersal dynamics, to model climate induced changes in the distribution of an invasive plant (Gunnera tinctoria) at a local scale. A new methodology is introduced, dynamic baseline benchmarking, which distinguishes climate‐induced alterations in species distributions from other potential drivers of change. Using this approach, it was concluded that climate change, based on IPCC and C4i projections, has the potential to increase the spread‐rate and intensity of G. tinctoria invasions. Increases in the number of individuals were primarily due to intensification of invasion in areas already invaded or in areas projected to be invaded in the dynamic baseline scenario. Temperature had the largest influence on changes in plant distributions. Water availability also had a large influence and introduced the most uncertainty in the projections. Additionally, due to the difficulties of parameterising models such as this, the process has been streamlined by utilising methods for estimating unknown variables and selecting only essential parameters.     

Predicting the potential distribution and forest impact of the invasive species Cydalima perspectalis in Europe

Invasive species have considerably increased in recent decades due to direct and indirect effects of ever‐increasing international trade rates and new climate conditions derived from global change. We need to better understand how the dynamics of early species invasions develop and how these result in impacts on the invaded ecosystems. Here we studied the distribution and severe defoliation processes of the box tree moth (Cydalima perspectalis W.), a tree defoliator insect native to Asia and invasive in Europe since 2007, through the combination of species distribution models based on climate and landscape composition information. The results showed that the combination of data from the native and the invaded areas was the most effective methodology for the appropriate invasive species modeling. The species was not influenced by overall landscape factors, but only by the presence of its host plant, dispersal capacity, and climate suitability. Such climate suitability was described by low precipitation seasonality and minimum annual temperatures around 0°C, defining a continentality effect throughout the territory. We emphasize the need of studying distribution and severe defoliation processes separately because we identified that climate suitability was slightly involved in limiting species spread processes but strongly constrained ecosystem impact in terms of defoliation before the species reaches equilibrium with the new environment. New studies on habitat recovery after disturbance, ecological consequences of such impact, and community dynamics in a context of climate change are required for a better understanding of this invasive species.     

Invader–invader interactions in relation to environmental heterogeneity leads to zonation of two invasive amphipods, Dikerogammarus villosus (Sowinsky) and Gammarus tigrinus Sexton: amphipod pilot species project (AMPIS) report 6

As biological invasions continue, interactions occur not only between invaders and natives, but increasingly new invaders come into contact with previous invaders. Whilst this can lead to species replacements, co-existence may occur, but we lack knowledge of processes driving such patterns. Since environmental heterogeneity can determine species richness and co-existence, the present study examines habitat use and its mediation of the predatory interaction between invasive aquatic amphipods, the Ponto-Caspian Dikerogammarus villosus and the N. American Gammarus tigrinus. In the Dutch Lake IJsselmeer, we found broad segregation of D. villosus and G. tigrinus by habitat type, the former predominating in the boulder zone and the latter in the soft sediment. However, the two species co-exist in the boulder zone, both on the short and longer terms. We used an experimental simulation of habitat heterogeneity and show that both species utilize crevices, different sized holes in a plastic grid, non-randomly. These amphipods appear to optimise the use of holes with respect to their ‘C-shape’ body size. When placed together, D. villosus adults preyed on G. tigrinus adults and juveniles, while G. tigrinus adults preyed on D. villosus juveniles. Juveniles were also predators and both species were cannibalistic. However, the impact on G. tigrinus of the superior intraguild predator, D. villosus, was significantly reduced where experimental grids were present as compared to absent. This mitigation of intraguild predation between the two species in complex habitats may explain the co-existence of these two invasive species.     

Evidence for pollen limitation of a native plant in invaded communities

Animal-pollinated invasive species have frequently been demonstrated to outcompete native species for pollinator attention, which can have detrimental effects on the reproductive success and population dynamics of native species. Many animal-pollinated invasive species exhibit showy flowers and provide substantial rewards, allowing them to act as pollinator ‘magnets’, which, at a large scale, can attract more pollinators to an area, but, at a smaller scale, may reduce compatible pollen flow to local native species, possibly explaining why most studies detect competition. By performing pollen limitation experiments of populations in both invaded and uninvaded sites, we demonstrate that the invasive plant Lythrum salicaria appears to facilitate, rather than hinder, the reproductive success of native confamilial Decodon verticillatus, even at a small scale, in a wetland habitat in southeastern Ontario. We found no evidence for a magnet species effect on pollinator attraction to invaded sites. Germination experiments confirmed that seeds from invaded sites had similar germination rates to those from uninvaded sites, making it unlikely that a difference in inbreeding was masking competitive effects. We describe several explanations for our findings. Notably, there were no differences in seed set among populations at invaded and uninvaded sites. Our results underscore the inherent complexity of studying the ecological impacts of invasive species on natives.     

Roe deer population structure in a highly fragmented landscape

Northern Belgium (Flanders) is one of the most densely populated and urbanized regions in Europe. Many species are therefore likely to suffer from anthropogenic pressure and habitat destruction and fragmentation. Although many large mammals are recolonizing in parts of Europe, including Belgium, due to adaptation, a relaxation of persecution and habitat restoration, we have little actual data concerning the effects of landscape features on their population structure. We analysed the genetic structure of discrete roe deer (Capreolus capreolus) populations in the Eastern part of Flanders, with special emphasis on the impact of habitat fragmentation and anthropogenic barriers. The sampled populations were clearly genetically differentiated. Genetic structure could be explained by purely distance-based landscape modelling, but a simpler model focusing solely on barrier effects of large transportation infrastructure explained nearly as much genetic variance. In contrast, analyses based on least-cost landscape modelling failed to yield a significant effect. Overall, the results suggest considerable landscape-level effects of transportation infrastructure.     

Predicting the current distribution and potential spread of the exotic grass Eragrostis plana Nees in South America and identifying a bioclimatic niche shift during invasion

Eragrostis plana (Poaceae) is a perennial grass introduced from South Africa to the state of Rio Grande do Sul in southern Brazil. Currently, it is considered an invasive grass in several regions of the world, including South America, where it has caused negative ecological and socio‐economic impacts. Ecological niche models, using bioclimatic variables, are often used to predict the potential distribution of invasive species. In this study we prepared two bioclimatic models for E. plana using the Genetic Algorithm for Rule‐set Production, the first based on data from its native region (South Africa) and the second on data from both the native and invaded (South America) regions. We then projected each model onto South America to identify regions vulnerable to invasion by the species, and compared our results with available records of the species in South America. Finally, we explored the model's predictions for the existence of a bioclimatic niche shift during the invasion process of E. plana in South America, using multivariate statistical analysis. The model created with native distribution data was only able to predict (with highly suitable habitat) the region of introduction of E. plana in South America. However, the current distribution, as well as the region of introduction of the species, was reliably predicted by the model created with data from both native and invaded regions. Our multivariate analysis supports a hypothesis of bioclimatic niche shift during the invasion process of E. plana in South America.     

Molecular markers for systematic identification and population genetics of the invasive Ponto‐Caspian freshwater gammarid Dikerogammarus villosus (Crustacea,Amphipoda)

The Ponto‐Caspian amphipod, Dikerogammarus villosus, is an invasive species of many European rivers. First, we show that size difference of nrDNA ITS1 allows discriminating D. villosus from Dikerogammarus bispinosus, a closely related but morphologically hardly distinguishable species. Second, we present two types of polymorphic markers for D. villosus, three microsatellites and two single nucleotide polymorphisms (SNPs) of mtDNA COI gene, which were scored by polymerase chain reaction‐single strand conformational polymorphism (PCR‐SSCP). These markers will be very useful in studying population genetics of D. villosus.     

Modelling the spread of European buckthorn in the Region of Waterloo

Established invasive species, such as European buckthorn (Rhamnus cathartica), pose a challenging problem for land managers who must decide when and how to control them. In order to make an informed decision land managers need to be able to predict the spread of these invasive species at local scales and without the need for excessive sampling. Our approach uses a hybrid model, combining habitat suitability and the presence of the invasive in neighbouring cells to predict the probability of a cell being invaded over time. A generalized linear mixed-effects model was used to create a habitat suitability model and a spread model. The habitat suitability model predicts the presence of buckthorn based on environmental characteristics and the results are used in the spread model. The spread model indicates that the invasion of buckthorn is influenced by the suitability of habitat and the presence of buckthorn in neighbouring cells. The success of the spread model suggests that this approach can be used to create a spatiotemporally explicit model with limited sampling effort.     

Niche differentiation between sympatric alien aquatic crustaceans: An isotopic evidence

Among the mechanisms that allow competing species to coexist are resource partitioning and dietary segregation. The current study uses multiple stable isotopes, carbon (δ¹³C), nitrogen (δ¹⁵N) and sulphur (δ³⁴S), to test the hypothesis that dietary segregation in cohabiting invasive mysids (Limnomysis benedeni and Katamysis warpachowskyi) and gammarids (Dikerogammarus villosus and Gammarus roeselii) will be reflected by differences in isotope values. Furthermore, IsoError mixing models were used to estimate the relative contributions of periphyton and seston to the invaders’ diets. Whole tissue δ¹³C, δ¹⁵N and δ³⁴S analysis in L. benedeni and K. warpachowskyi imply that these sympatric, non-native mysids maintain differentiated feeding niches or resource partitioning by feeding on distinct components of the available food resources (predominantly seston by L. benedeni and periphyton by K. warpachowskyi). By contrast, the gammarids D. villosus (‘killer shrimp’) and G. roeselii exhibited no significant difference in δ¹³C and δ¹⁵N, indicating a considerable overlap between the dietary sources of these sympatric invaders. Feeding niche differentiation, irrespective of season or the nature of habitat invaded (lake or river), might facilitate the coexistence of invasive mysids in their ‘new’ environment by minimizing direct resource competition. The mutual interaction by the invasive gammarids, coupled with voracious behavior, could assist their success at co-invasion with serious implications for local biodiversity including the potential extinction of native species.     

Current distribution and potential extent of the most invasive alien plant species on La Réunion (Indian Ocean,Mascarene islands)

Abstract La Réunion Island has the largest area of intact vegetation of the islands in the Mascarene archipelago. Biological invasions are the primary threat to biodiversity in the intact habitats of the island (those not already transformed by agriculture and urbanization). Our study aimed to identify areas to prioritize in managing invasive alien plants for biodiversity conservation. We used extensive surveys of 238 distinct untransformed areas on La Réunion to define the current distribution patterns of all invasive species. Using expert knowledge, we compiled maps of the current distribution of the 46 most widespread/important invasive plants at the habitat scale (identified according to vegetation structure). Data from 440 botanical relevés for the 20 most threatening invasive alien plant species across the island and climatic envelope models were used to derive climatic suitability surfaces; these were used to map potential distributions for these species. More than 10 species invade 16.7% of the remaining habitat. Five habitat types are invaded by 25 or more species, and eight have fewer than 10 invasive alien plant species. Cluster analysis based on presence/absence of species in the 18 habitat types produced eight groups of species that invade particular habitats. Potential distribution models show that some species have invaded large parts of their potential range (e.g. Fuchsia magellanica, Furcraea foetida, Hiptage benghalensis), whereas others have the potential to increase their range substantially (e.g. Clidemia hirta, Strobilanthes hamiltonianus, Ulex europaeus). Management implications are identified for both groups. Three broad groups of habitats were identified: (i) intact habitats with a low level of invasion (e.g. subalpine shrubland); (ii) moderately invaded habitats with varying levels of intactness (ranging from windward submountain rainforest to the Acacia heterophylla forest); and (iii) habitats with little remaining intact area and high levels of invasion (e.g. lowland rainforest). Different management interventions are appropriate for these three groups.     

Parasitism may alter functional response comparisons: a case study on the killer shrimp <Emphasis Type="Italic">Dikerogammarus villosus</Emphasis> and two non-invasive gammarids

The Ponto-Caspian freshwater amphipod Dikerogammarus villosus has colonized most of the water bodies of continental Europe where it causes strong structural alterations in recipient communities that can lead to changes in ecosystem-level processes, mainly because of a strong predatory behaviour. Most of the D. villosus populations from the invaded range have been found infected with the co-introduced microsporidian parasite Cucumispora dikerogammari, known to decrease the predation rate of its host. Infection might thus mitigate the ecological impact of D. villosus and we wanted to test this assumption using the comparative functional response approach. We compared the relationship between resource use and resource availability (i.e. the functional response, FR) of D. villosus, either with infected individuals or not, to that of two non-invasive gammarids: Gammarus pulex and Echinogammarus berilloni. With infected individuals included, D. villosus displayed a higher FR than the two non-invasive gammarids. Although this effect was not significant, C. dikerogammari infection tended to alter the FR of D. villosus with a slight decrease in attack rate and handling time, resulting in a less steep initial slope and a higher asymptote, respectively. Removing infected D. villosus from the dataset did not affect the FR comparison with G. pulex but suppressed the difference in FR with E. berilloni. Although we cannot exclude the role of sample size reduction in this effect, this suggests that C. dikerogammari infection might increase the predation pressure on local prey populations in case of species replacement between D. villosus and E. berilloni. From a more general perspective, our study illustrates how parasites may alter our capacity to predict invasive species impacts from FR comparisons.     

Two colonisation stages generate two different patterns of genetic diversity within native and invasive ranges of Ulex europaeus

Genetic diversity and the way a species is introduced influence the capacity of populations of invasive species to persist in, and adapt to, their new environment. The diversity of introduced populations affects their evolutionary potential, which is particularly important for species that have invaded a wide range of habitats and climates, such as European gorse, Ulex europaeus. This species originated in the Iberian peninsula and colonised Europe in the Neolithic; over the course of the past two centuries it was introduced to, and has become invasive in, other continents. We characterised neutral genetic diversity and its structure in the native range and in invaded regions. By coupling these results with historical data, we have identified the way in which gorse populations were introduced and the consequences of introduction history on genetic diversity. Our study is based on the genotyping of individuals from 18 populations at six microsatellite loci. As U. europaeus is an allohexaploid species, we used recently developed tools that take into account genotypic ambiguity. Our results show that genetic diversity in gorse is very high and mainly contained within populations. We confirm that colonisation occurred in two stages. During the first stage, gorse spread out naturally from Spain towards northern Europe, losing some genetic diversity. During the second stage, gorse was introduced by humans into different regions of the world, from northern Europe. These introductions resulted in the loss of rare alleles but did not significantly reduce genetic diversity and thus the evolutionary potential of this invasive species.