Introduced species shed friends as well as enemies

Many studies seeking to understand the success of biological invasions focus on species’ escape from negative interactions, such as damage from herbivores, pathogens, or predators in their introduced range (enemy release). However, much less work has been done to assess the possibility that introduced species might shed mutualists such as pollinators, seed dispersers, and mycorrhizae when they are transported to a new range. We ran a cross-continental field study and found that plants were being visited by 2.6 times more potential pollinators with 1.8 times greater richness in their native range than in their introduced range. Understanding both the positive and negative consequences of introduction to a new range can help us predict, monitor, and manage future invasion events.


Far fewer flower visitors
We began by testing the hypothesis that the abundance and taxonomic richness of flower visitors will be lower in the introduced range than in the native range.Interacting with more, or more types of pollinators improves plants' ability to produce more fruits and set more seeds 15 .Therefore, a decrease in the richness of flower visitors is predicted to have a negative impact on successful reproduction, and ultimately invasion.When a plant species invades, it can be integrated into the interaction-web of generalist pollinators 16,17 .But this uptake of some generalists may not be enough to compensate for the decrease in the richness of flower visitors plants receive in their introduced range.Few studies investigating the impact of range on pollinator richness currently exist 10 and often compare introduced species with native species occupying the same habitats 16,[18][19][20] .Our biogeographic comparison of species across their native and introduced ranges enables us to detect patterns that may not be as visible when using a community comparative approach 6,21 .
In a study spanning four countries, two continents, and ten plant species, we observed 2652 flower visits (Fig. 2).We show that plants interact with 2.6 times more flower visitors in the native range than in the introduced OPEN 1   Evolution and Ecology Research Centre, School of Biological, Earth and Environmental Sciences, UNSW Sydney, Sydney, NSW 2052, Australia. 2Kunming Institute of Botany, Chinese Academy of Sciences, Kunming, China. 3Faculty of Arts, Science and Technology, University of Northampton, Northampton, UK. 4 Terrestrial Ecology Group (TEG), Department of Ecology, Institute for Biodiversity and Global Change, Universidad Autónoma de Madrid, 28049 Madrid, Spain. 5Stats Central, Mark Wainwright Analytical Centre, UNSW Sydney, Sydney, NSW 2052, Australia. 6Institute of Ecology and Earth Sciences, University of Tartu, J. Liivi 2, 50409 Tartu, Estonia. 7CSIRO Health & Biosecurity, Brisbane, QLD, Australia. 8Present address: School of Life Sciences, Faculty of Science, University of Technology Sydney, Sydney, NSW 2007, Australia.* email: zoe.xirocostas@uts.edu.aurange (estimate = 0.948, SE = 0.343, P = 0.006; Fig. 3a).Plants also interacted with a 1.8 times greater richness of flower visitors in the native range than in the introduced range (estimate = 0.569, SE = 0.22, P = 0.01; Fig. 3b).One species, Lotus corniculatus, had no observed visitors in the introduced range despite being visited by six different taxa in the native range.Seven of ten species showed evidence for flower visitor assemblage dissimilarity between their native and introduced ranges (P < 0.05; Fig. 4).In plants' native range, we observed 1592 flower visits which were dominated by bees and wasps (Aculeata; 35.3% of interactions), followed by flies (Diptera; 28.1% of interactions), and beetles (Coleoptera; 27.6% of interactions).In the introduced range we observed 1060 interactions which comprised mostly of bees and wasps (39.1% of interactions), flies (38.5% of interactions), and ants (Formicidae; 9.7% of interactions).
Our study provides the most comprehensive test of the Missed Mutualist Hypothesis to date 9,10,12,[22][23][24][25] .The observed reduction in visitor abundance and richness might be substantial enough to impact plants' reproductive potential in their introduced ranges 26 .However, the actual effect of missed mutualists could surpass what our study implies since we only sampled species that had managed to establish populations in a new range.Species reliant on pollinators that are absent in the introduced range would be incapable of establishing viable populations in their new habitat and, therefore, be excluded from our study.
A loss of flower visitors could help to explain the extremely high proportion of unsuccessful plant invasions 10 , as only ~ 25% of plant species successfully take consecutive steps in the invasion process (i.e., introduction, establishment, and spread) 27 .Introduced plants could overcome this disadvantage through evolution of a greater capacity for selfing [28][29][30] .For example, Arctotheca populifolia, a beach daisy native to South Africa, adapted to reproduce asexually in less than 200 years since its introduction to Australia and interaction with fewer pollinators 31 .A global study by Razanajatovo et al. 32 further supports this idea, finding that selfing plants are more likely to establish in new ranges.Similarly, Pyšek et al. 1 found that in Central Europe, selfing was the best reproductive strategy to facilitate invasive plants.Plants missing their mutualist visitors could evolve traits that increase visitation by new taxa 10 .Consistent with this idea, seven of our ten study species showed evidence for visitor assemblage dissimilarity between their native and introduced ranges (P < 0.05).Understanding visitor assemblages between ranges could give important new insights into the factors shaping the reproductive success and spread of introduced plants.
Visitation of introduced plants may not have decreased enough for the plants to be pollen limited, which could explain their ability to thrive despite missing mutualists.There is mixed evidence in the literature for an effect of visitation frequency and pollinator richness on the amount of seed set by plants 15,[33][34][35][36][37][38][39] .While our plants www.nature.com/scientificreports/may be encountering fewer, less diverse, and different pollinators, they could still be setting similar quantities of seeds if the pollinators that do interact with them are highly efficient and transfer high pollen loads (i.e., they are not pollen limited).There is evidence to suggest that some mutualistic interactions may even be enhanced in the www.nature.com/scientificreports/introduced range and promote invasion, however, these studies only consider belowground mutualisms and it is currently unknown whether this concept translates across pollination 40,41 .Our study can prompt future work directed at examining pollen loads/deposition and fruit/seed production across ranges to disentangle specific drivers underpinning the success of introduced plant species.A loss of flower visitors could have evolutionary consequences, possibly even the evolution of increased competitive ability (EICA) 10 .EICA is most commonly associated with enemy release, whereby introduced species gain a competitive edge over native species by reallocating energy from defence to growth and reproduction (Fig. 1), as a result of reduced herbivore pressure 4 .However, a reduction in mutualists, or certain types of mutualists, may also have a similar effect.For example, energy allocated to attracting specialised pollinators in the native range (i.e., costly nectar production, production of showy flowers, or long flower tubes), may be redirected into setting seeds of higher quality or quantity in the introduced range, an advantage not possessed by native competitors.
Our study focuses on entomophilous species whose flowers are visited by an array of generalist pollinators.However, not all introduced plant species adopt this reproductive strategy 42 .Asexually reproducing species (e.g., through rhizomes or stolons) alongside non-outcrossing or self-pollinating plants, do not rely on mutualistic floral interactions to ensure successful reproduction [43][44][45][46] .Coevolved interactions between flowers and visitors may also be highly specialised with successful pollinators constrained to taxa with certain morphologies (e.g., long proboscis to reach down nectar tubes) or belonging to specific clades 47,48 .More research is needed to understand if and how the missed mutualist hypothesis may apply across these different reproductive modes and their subsequent impact on invasion success.

Introduced species lose more foe than friends
A parallel study found that plant species experienced an average of 5.8 times more herbivory in their native range than in their introduced range 49 , while the present work revealed only 2.6 times more visitors in their native range (Paired T test mean difference = 0.81, t = 3.19, df = 9, P = 0.01; Fig. 5).That is, introduced plants seem to shed more enemies than mutualists.This may be because pollination is a mutually beneficial interaction posing little risk to naïve mutualists in the introduced range 36 .Conversely, generalist enemies may be more hesitant to ingest leaf material from unfamiliar non-native plants, as they pose a risk of harbouring deadly defences that can reduce herbivore fitness, or even result in herbivore mortality 50,51 .Another possibility is that the selective pressure on species to find alternative pathways for reproduction is more direct than the selective pressure resulting from herbivory.The fact that the benefits of enemy release outweigh the cost of missed mutualisms (Fig. 5) might help to explain the success of introduced plants.www.nature.com/scientificreports/ Ecologists have long focused their attention on understanding how enemy release can facilitate successful invasion.Here, we show evidence that missed mutualisms are also likely to be important in the invasion process.Our findings present a valuable opportunity for ecologists to further explore whether other types of mutualisms are missed in the introduced range (e.g., plants and beneficial fungi or protective mites), and whether the interplay of both enemy release and missed mutualisms could more accurately predict which species, in what locations, would be more likely to invade natural ecosystems.

Data collection
To determine whether plants are interacting with fewer, less diverse pollinators in their introduced range we chose target species that were: 1. Biotically pollinated.2. Present across a wide area in their native (Europe) and introduced ranges (Australia).
Using these criteria, we identified 15 plant species.Some species, however, were unable to be located or observed at least once in both ranges despite our best efforts in the field.Our final dataset is comprised of observations from the ten herbaceous plant species belonging to nine families and eight orders for which we were able to take observations in at least one site across both native and introduced ranges (Supplementary information).
We conducted floral observations at nine sites within the native and introduced ranges of ten plant species.For most of our study species in Australia, the exact source population, or populations for multiple introductions (as in Hypericum perforatum 52 and Trifolium repens 53 ), are not known.Further, most of the introduced species present in Australia have very wide home ranges [54][55][56][57][58][59][60][61][62][63] .Thus, studying a single native population and comparing it with a single introduced population could give misleading results.We therefore selected a broad range of sites and climatic conditions to get a general understanding of plant-pollinator interactions within and across native and introduced ranges (Fig. 6).We prioritised maximising the latitudinal range and landscape diversity in each range when selecting our study sites.We also considered the presence of our target species when choosing sites as we favoured places that would increase our sampling potential.Prior to choosing site locations we used online resources such as the Global Biodiversity Information Facility (gbif.org)and the Atlas of Living Australia (ala.org.au) to check the presence of our target species.However, not all study species were present at each site (i.e., city or region where sampling took place) (see Supplementary information).
Floral observations were made in early to mid-summer of 2019 from May to August in Europe and from September to December in Australia.When choosing individuals, we looked for those that were actively in flower www.nature.com/scientificreports/and then randomly selected individual plants from that subset population.This was done using a compass and random number generator to determine the observer's direction of movement, the first individual encountered (or nearby) when walking in this direction was sampled.If ten or fewer individuals were present, all were sampled.All observations were conducted in daylight hours between 10am and 3 pm on days with no rain and minimal wind.Each observation was timed at 15 min and conducted with the observer placed 1-2 m from the target plant.We aimed to repeat these timed observations for at least 10 individuals of each species at every site.
In total, we conducted 250 timed observations in the native range and 236 in the introduced range.
To be considered a pollinator, an animal must act as a vector for transferring pollen from one flower to another 66 .As we could not visually ensure successful pollination occurring during our observations, we instead quantified the number of times a flower was visited by invertebrates (potential pollinators).We defined a 'visit' as any time an invertebrate came into contact with floral reproductive organs (e.g., anthers or stigma), as implemented across similar field studies 67,68 .The duration of a visit was not recorded, and longer visits were considered the same as shorter ones (e.g., if a bee interacted with a flower for 5 min or 5 s it was still considered one visit).Visitors were categorised according to their taxonomic group as bees/wasps (Aculeata), flies (Diptera), ants (Formicidae), beetles (Coleoptera), butterflies/moths (Lepidoptera), spiders (Araneae), true bugs (Hemiptera), thrips (Thysanoptera), dragonflies (Odonata), and grasshoppers (Orthoptera).The total number of "floral units" (defined for each species in the Supplementary information) where visitors were observed on the target individual were also counted.If an invertebrate visited the same floral unit more than once (i.e., the animal left and then returned) or visited multiple floral units on the same individual, they were counted as multiple visits.
Visitor abundance was calculated as the total number of visits from each taxonomic group for each observation.Visitor richness was defined as the number of taxonomic groups that were recorded for each observation.
Estimates of enemy release for our comparative analysis were sourced from Xirocostas et al. 49 .

Data analysis
All statistical analyses were performed in RStudio version 4.2.0 69 .We included all groups of arthropods observed on floral units in our analyses, including those not (currently) known to be pollinators such as Odonata and Araneae 66 , because although they may not act as pollinators in Europe, it is not known if they play this role in Australia.
To test the hypothesis that plants in their introduced range will receive fewer visits to their flowers than conspecifics in the native range, we performed an overall comparison of visitor abundance between ranges for each species with Generalized Linear Mixed Models using Template Model Builder 70 .Our response variable was visitor abundance, and the predictor variable was range.We also included random effects terms for site and species.We offset our model by the log-transformed number of floral units per observation to account for this variance in affecting our abundance counts and used a negative binomial family as our data were over-dispersed.
To determine whether visitor richness was higher in plants' native range than in the introduced range, we also fit a Generalized Linear Mixed Model.Our response variable was visitor richness and predictor variable was range; we also included random effects terms for site, species, and the polynomial-transformed number of floral units.We offset our model by the log-transformed number of floral units to account for this variance in affecting our richness counts and used a Poisson distribution as our data were not as over dispersed as the abundance data.
To determine whether the visitor assemblage of flower visitors differs between ranges we performed separate comparisons for each species across their native and introduced range.For each target species we created a community matrix of all observed visitor taxa using the mvabund function in the mvabund package 71 .Then we ran alternate and null hypothesis generalised linear models using the manyglm function.Our alternate models used the community matrix as the response variable, range as the predictor variable, and were offset by the log-transformed number of floral units to account for this variance across observations.Our null models were similar, except that our predictor variable (range) was removed and replaced with an intercept term (1).We then ran an ANOVA to calculate the distribution of our test statistic under the null hypothesis of no effect of range.Adjusted P values were calculated using 999 iterations via case block resampling with a stratified cluster bootstrap to account for correlation due to site.
To test whether plants experience differing degrees of enemy release in comparison to pollinator loss, we compared our floral visitation data with enemy release data from a parallel study 49 conducted on the same plant species.For each species, herbivory data were collected by visually estimating leaf damage 72 on a percentage scale of zero (no damage) to one hundred (complete consumption) on 10 leaves for at least 12 randomly selected individuals per site [full protocol in Xirocostas et al. 49 ].We then performed a paired samples t-test on the coefficients generated from the generalised linear mixed models performed in Xirocostas et al. 49 and our study.These model coefficients represent the magnitude of the effect (log response ratio), for enemy release (herbivore damage) and missed mutualisms (flower visitor abundance) for each of our target species.

Figure 1 .
Figure 1.Conceptual illustration of the differences in key interactions (non-exhaustive) explained by the Enemy Release Hypothesis and Missed Mutualist Hypothesis in plants' native and introduced ranges.In this paper we consider interactions between plants and herbivores/flower visitors only.

Figure 3 .
Figure 3.Comparison of (a) the mean number of flower visits and (b) mean number of visiting taxa per 15-min observation for each species in the native (yellow) and introduced (pink) ranges.To calculate this metric, visits were first divided by the number of floral units per observation before taking means (per species per site) [see Supplementary information for variance across sites and model coefficients/confidence intervals per species].Images of plant species observed in this study are displayed to the right of both graphs accompanied by letters A-J which correspond to letters displayed on the graphs (all images credited to Z Xirocostas).

Figure 4 .
Figure 4. Comparison between ranges of the community composition of visitors per floral unit for each plant species.Significant differences (P < 0.05) are determined by multivariate abundance analyses and are denoted by asterisks (*) (values listed in Supplementary information).

Figure 6 .
Figure 6.Maps of Europe (native range) and Australia (introduced range) where floral observations took place.[Bottom left] Sites in Europe include Madrid (Spain), Northampton (United Kingdom) and Tartu (Estonia).[Bottom right] Sites in Australia include Hobart (Tasmania), Melbourne (Victoria), Cooma (New South Wales), Canberra (Australian Capital Territory), Robertson (New South Wales), and Sydney (New South Wales).Maps are shaded according to mean temperature of the warmest quarter from WorldClim version 2.1 climate data for 1970-2000 64 , as it is the time of year when most pollination occurs 65 .