
==== Front
Evol Appl
Evol Appl
10.1111/(ISSN)1752-4571
EVA
Evolutionary Applications
1752-4571
John Wiley and Sons Inc. Hoboken

10.1111/eva.70016
EVA70016
EVA-2024-028-OA.R1
Original Article
Original Article
Spatially Varying Wolbachia Frequencies Reveal the Invasion Origin of an Agricultural Pest Recently Introduced From Europe to North America
Lečić Sonja 1
Wolfe Thomas M. 1
Ghosh Animesh 1
Satar Serdar 2
Souza Beraldo Camilla https://orcid.org/0000-0003-3274-375X
3
Smith Emily 4
Dombroskie Jason J. https://orcid.org/0000-0002-7232-6507
5
Jernigan Emily 5
Hood Glen Ray https://orcid.org/0000-0002-9301-6153
4 glenrayhood@wayne.edu

Schuler Hannes https://orcid.org/0000-0001-8307-9831
6 7
Stauffer Christian 1
1 Department of Forest and Soil Sciences Boku University Vienna Austria
2 Department of Plant Protection, Faculty of Agriculture Çukurova University Adana Turkey
3 Organismal and Evolutionary Biology Research Programme, Faculty of Biological and Environmental Sciences The University of Helsinki Helsinki Finland
4 Department of Biological Sciences Wayne State University Detroit Michigan USA
5 Department of Entomology Cornell University Ithaca New York USA
6 Competence Centre for Plant Health Free University of Bozen‐Bolzano Bozen‐Bolzano Italy
7 Faculty of Agricultural, Environmental and Food Sciences Free University of Bozen‐Bolzano Bozen‐Bolzano Italy
* Correspondence:
Glen Ray Hood (glenrayhood@wayne.edu)

20 9 2024
9 2024
17 9 10.1111/eva.v17.9 e7001614 6 2024
12 2 2024
19 8 2024
© 2024 The Author(s). Evolutionary Applications published by John Wiley & Sons Ltd.
https://creativecommons.org/licenses/by/4.0/ This is an open access article under the terms of the http://creativecommons.org/licenses/by/4.0/ License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

ABSTRACT

The introduction of non‐native species across the world represents a major global challenge. Retracing invasion origin is an important first step in understanding the invasion process, often requiring detailed sampling within the native range. Insect species frequently host Wolbachia, a widespread endosymbiotic bacterium that manipulates host reproduction to increase infected female fitness. Here, we draw on the spatial variation in infection frequencies of an actively spreading Wolbachia strain wCer2 to investigate the invasion origin of the European cherry fruit fly, Rhagoletis cerasi. This pest of cherries was introduced from Europe to North America within the last decade. First, we screen the introduced fly population for the presence of Wolbachia. The introduced populations lack the wCer2 strain and the strongly associated mitochondrial haplotype, suggesting strain absence due to founder effects with invading individuals originating from wCer2‐uninfected native population(s). To narrow down geographic regions of invasion origin, we perform spatial interpolation of the wCer2 infection frequency across the native range and predict the infection frequency in unsampled regions. For this, we use an extensive dataset of R. cerasi infection covering 238 populations across Europe over 25 years, complemented with 14 additional populations analyzed for this study. We find that R. cerasi was unlikely introduced from wCer2‐infected populations in Central and Western Europe. We propose wCer2‐uninfected populations from Eastern Europe and the Mediterranean region as the most likely candidates for the invasion origin. This work utilizes Wolbachia as an indirect instrument to provide insights into the invasion source of R. cerasi in North America, revealing yet another application for this multifaceted heritable endosymbiont. Given the prevalence of biological invasions, rapidly uncovering invasion origins gives fundamental insights into how invasive species adapt to new environments.

endosymbiont
European cherry fruit fly
invasive species
Rhagoletis cerasi
spatial interpolation
Animal and Plant Health Inspection Service 10.13039/100009168 AP22PPQS&T00C053 Austrian Science Fund 10.13039/501100002428 31441‐B29 P35327‐B Wayne State University 10.13039/100006710 Department of Innovation, Research and University of the Autonomous Province of Bozen/Bolzano source-schema-version-number2.0
cover-dateSeptember 2024
details-of-publishers-convertorConverter:WILEY_ML3GV2_TO_JATSPMC version:6.4.8 mode:remove_FC converted:20.09.2024
Funding: This work was supported by Wayne State University and USDA‐APHIS‐PPQ (AP22PPQS&T00C053). Additional funding was provided by the FWF projects P35327‐B to CS and 31441‐B29 to HS. Funding for the open access publication charges was provided by the Department of Innovation and Research, and the University of the Autonomous Province of Bozen/Bolzano.

The first two authors contributed equally to this article.

The last three authors share senior authorship.
==== Body
pmc1 Introduction

Facilitated by increasing global connectivity, biological invasions are integral to global environmental change (Early et al. 2016; Pyšek et al. 2020). The growing concerns have steered researchers toward a better understanding of how species adapt to new environments (White et al. 2013; Tepolt et al. 2022) and disentangling the factors that enable their range expansions (Deshpande and Fronhofer 2022). An essential first step in understanding biological invasions is to uncover the invasion origin of introduced species (Dlugosch and Parker 2008; Estoup and Guillemaud 2010). Identifying the source populations relies on a careful and detailed sampling within the native range (Dlugosch et al. 2015) because demographic processes such as founder effect and admixture can confound the inference of invasion origins (Bock et al. 2015; Sherpa and Després 2021).

Due to founder effects, introduced populations may carry a fraction of the genetic diversity compared to populations in their native range (Prentis et al. 2008; Szűcs et al. 2017). Similarly, upon introduction, a fraction of invading individuals can bring their obligate residents into a new range (Hurst and Jiggins 2005; Lu, Hulcr, and Sun 2016) or experience a reduction in endosymbiont diversity (Reuter, Pedersen, and Keller 2005; Shoemaker et al. 2000). Wolbachia, a widespread bacterial endosymbiont among insects (Werren and Windsor 2000; Jiggins et al. 2001), occupies the germline of its insect hosts and is vertically inherited through the maternal line. A common way in which Wolbachia manipulate host reproduction and ensure a reliable transmission in the new host is through cytoplasmic incompatibility (CI) (Turelli and Hoffmann 1991). Here, uninfected females produce few or no offspring when mating with infected males or males infected with another Wolbachia strain (Hoffmann and Turelli, 1988).

The European cherry fruit fly, Rhagoletis cerasi (Linnaeus) (Diptera: Tephritidae), is a severe pest in cherry orchards (genus Prunus), but can also be found infesting honeysuckle berries (genus Lonicera). The fly is widespread across its native range in Eurasia (Boller and Prokopy 1976). All native populations of R. cerasi are fixed for a Wolbachia strain wCer1. In addition, R. cerasi can host coinfection by a second strain, wCer2. However, unlike the fixed wCer1 strain, wCer2 is found in spatially varying frequencies in R. cerasi populations across Europe, ranging from (near) fixed in some locations, but completely absent in others (Riegler and Stauffer 2002; Schuler et al. 2016; Bakovic et al. 2018; Schebeck et al. 2019). The wCer2 strain causes strong CI between wCer2‐infected males and wCer2‐uninfected females, with laboratory estimates of CI as high as 98% (Boller et al. 1976; Riegler and Stauffer 2002). In the native range, R. cerasi harbors two mitochondrial haplotypes, denoted HT1 and HT2 (Schuler et al. 2016). The HT1 variant is associated with wCer1‐infected flies, whereas mitochondrial haplotype HT2 is associated with the spreading and CI‐causing wCer2 strain. In Central Europe, wCer2 has spread in R. cerasi through CI, causing haplotype HT2 to hitch‐hike alongside the Wolbachia strain (Schuler et al. 2016). This sweep presumably replaced HT1 in the wCer2‐infected populations, with deviations in transitional populations detected (Schuler et al. 2016; Schebeck et al. 2019).

Classical models of Wolbachia infection frequency dynamics consider the proportion of uninfected eggs produced by infected females (μ), the relative reproductive success of infected and uninfected females (F), and the strength of CI (H) as the three main factors affecting changes in infection frequency (Hoffmann, Turelli, and Harshman 1990). If initially very rare, Wolbachia frequency will either stabilize at 0% or must satisfy F (1–μ) > 1 to spread deterministically from low frequencies, regardless of whether they cause CI (Kriesner et al. 2013; Turelli, Katznelson, and Ginsberg 2022). Once sufficiently common, the effects of strong CI push Wolbachia to high‐equilibrium frequencies (Hoffmann, Turelli, and Harshman 1990; Barton and Turelli 2011). Reproductive manipulations with strong effects pave the way to the rapid spread of new invading Wolbachia strains, exemplified with wRi in Drosophila simulans populations (Turelli and Hoffmann 1991, 1995; Kriesner et al. 2013). A fast spread can also reduce mitochondrial diversity or promote mitochondrial polymorphisms in host populations (Turelli and Hoffmann 1991; Turelli, Hoffmann, and McKechnie 1992; Jiggins and Tinsley 2005; Xiao et al. 2012; Arif et al. 2021). Such rapid spread creates geographical differences in the Wolbachia infection frequency and, sometimes, mitochondrial diversity through the fly's geographic range (Turelli and Hoffmann 1991; Hoshizaki and Shimada 1995; Schuler et al. 2016). Due to this spatial heterogeneity, we hypothesize that in natural conditions Wolbachia can serve as an indirect tool to infer insect invasion routes.

Recent reports indicate that R. cerasi was introduced into the northeastern United States and southeastern Canada sometime within the last decade (Barringer 2018; Wakie, Yee, and Neven 2018). In 2016, the fly was first found in Mississauga, Ontario, Canada, and again in 2017 in the United States in Niagara County, New York (Barringer 2018; Wakie, Yee, and Neven 2018) near or in stands of introduced European honeysuckle. The United States Department of Agriculture (USDA) now considers the fly a threat for farmers and stakeholders in the region where R. cerasi is currently under regulatory control, restricted to a quarantine area comprising five counties in western New York, covering a total of 5187 km2. To further underscore the broader implications of the threat of this pest species, environmental modeling suggests that the fly could spread to and establish in the major cherry‐growing states of Michigan, Washington, Oregon, and California where production is collectively valued at ~$700 million US annually (Wakie, Yee, and Neven 2018). Therefore, there is a need to control the spread of the fly, prevent future invasions, and elucidate the pathways of invasion from the fly's native range in Europe.

In this study, we take advantage of the differences in infection frequencies of the wCer2 strain in R. cerasi across the native range to infer the region(s) from which the fly was introduced into North America. To corroborate the results gleaned from comparing the frequencies of the wCer2 strain in the native in introduced population, we take advantage of the wCer2 association with the mitochondrial haplotype HT2. The invasive population lacked the wCer2 strain and the associated mitochondrial haplotype, suggesting that the fly likely originating from a wCer2‐uninfected native population. We then combine infection frequency information from 14 populations sampled for this study with the wCer2 infection frequency dataset consisting of an additional 238 populations sampled over the last 25 years across the fly's native range (Riegler and Stauffer 2002; Schuler et al. 2016; Bakovic et al. 2018), and perform spatial modeling to predict the invasion origin. Ultimately, our study highlights the utility of endosymbionts as a prospective tool for invasion genetics research and discuss our findings in light of Wolbachia and host mitochondria interactions.

2 Materials and Methods

2.1 Estimating wCer2 Infection Frequencies and mtDNA Haplotyping

To determine the infection frequency of the wCer2 strain in the introduced range, R. cerasi adult flies were caught on yellow sticky traps from several locations in Niagara County, collectively referred to as a single population within the introduced range in 2018 (Figure 1A; Table S1). Adult flies were removed from the sticky traps, cleaned using histoclear, identified to species using distinguishing characteristics in wing banding patterns, sexed, and stored in 70% ethanol. The DNA of whole‐body individuals was extracted using the Qiagen DNeasy Blood and Tissue extraction kit following the manufacturer's protocol.

FIGURE 1 Geographic locations and wCer2 infection frequencies of Rhagoletis cerasi populations and the distribution of wCer2 in (A) the introduced range in Niagara County, New York state, in northeastern United States and (B) the native range across Europe. Due to a very high sampling density in some regions and large sample sizes for several populations, locations from prior studies are represented as single black dots. The infection frequencies of the 13 native populations including the population from Russia (inset) sampled in this study are represented with pie charts showing wCer2 infection frequency (red: wCer2‐infected; grey: wCer2‐uninfected). The size of each pie chart reflects the number of individuals sampled at each location.

To estimate distance‐dependent population variance of the wCer2 infection frequency in the native range, we sampled additional R. cerasi populations from Prunus trees in 14 locations between 2018 and 2022 (Figure 1A; Table S1). Our motivation for the choice of sampling locations was twofold. First, we resampled several regions of the native range to test for patterns of temporal stability in wCer2 frequencies. In this case, observations of a general lack of temporal variation in the uninfected and fixed wCer2 populations would enable us to pool all field estimates of the wCer2 infection frequencies spanning 25 years and perform a spatial population variance analysis to predict infection frequencies in unobserved locations. Second, to improve the accuracy of distance‐dependent population variance estimates of the wCer2 infection frequency, we sampled populations from Turkey, Finland, Sweden, Latvia, and Russia. These populations represent the first field estimates of wCer2 infection frequencies from these regions of the native fly range.

To sample flies in the native range, we collected and placed cherry fruits sampled at each location in plastic trays at room temperature and stored emerging larvae upon pupation in 95% ethanol. The total DNA of single pupae was extracted using the Qiagen Gentra Puregene extraction kit following the manufacturer's protocol. We then used PCR to assess the infection status (presence vs. absence) of each individual from each location in the native and introduced ranges following Schuler et al. (2016) with wCer2‐specific primers targeting regions of the surface protein wsp (Arthofer et al. 2009). The proportion of infected flies served as our estimate of infection frequency (p) for each location. To score the presence versus absence of wCer2 strain infection in the native and introduced populations, we visualized PCR products via electrophoretic separation on 2% ethidium bromide‐stained agarose gels. To distinguish between mitochondrial haplotype HT1 and HT2 for an individual sample (Table 1), we amplified, sequenced and scored the 546‐bp fragment of the cytochrome oxidase I (COI) gene using the primers following the methods of Schuler et al. (2016).

TABLE 1 Sample size (N), number of wCer2‐infected individuals, infection frequencies p with exact 95% binomial confidence intervals, and mitochondrial haplotype are shown for each Rhagoletis cerasi population sampled from 2018 to 2020 in the native and introduced range.

Population	N	wCer2‐infected	p (confidence intervals)	Haplotype	
Austria	12	12	1 (0.758, 1)	HT2	
Italy	12	12	1 (0.758, 1)	HT2	
France 1	12	12	1 (0.758, 1)	HT2	
France 2	12	10	0.833 (0.552, 0.953)	HT1	
France 3	12	12	1 (0.758, 1)	HT2	
France 4	12	12	1 (0.758, 1)	HT2	
Turkey 1	9	0	0 (0, 0.299)	HT1	
Turkey 2	10	0	0 (0, 0.278)	HT1	
Russia	12	2	0.167 (0.0470, 0.448)	HT2	
Greece	2	0	0 (0, 0.658)	HT1	
Poland	12	10	0.833 (0.552, 0.953)	HT1, HT2	
Sweden	12	11	0.917 (0.646, 0.996)	HT1, HT2	
Latvia	11	5	0.455 (0.213, 0.720)	HT2	
Finland	5	3	0.6 (0.231, 0.882)	HT1, HT2	
USA	40	0	0 (0, 0.087)	HT1	

All statistical analyses were performed in R v. 4.01 (R Core Team 2020). Assuming a binomial distribution for the infection frequency (p), we estimated 95% confidence intervals for p for each population using the “binconf” function in the package Hmisc (Harrell and Dupont 2018). Finally, we performed Fisher's exact tests to determine if pairwise differences in infection frequency exist among R. cerasi populations spatially and/or temporally across the fly's native range.

2.2 Spatial Interpolation of the wCer2 Infection Frequency in the Native Range

To infer the region(s) in Europe from which R. cerasi may have been introduced, we first captured spatial variation in wCer2 frequencies by inferring distance‐dependent population variance of the wCer2 infection frequency in the native range using field estimates of wCer2 infection frequencies gathered over 25 years (Riegler and Stauffer 2002; Arthofer et al. 2009; Schuler et al. 2016; Bakovic et al. 2018). We complemented this dataset with wCer2 infection frequencies from the 14 additional populations sampled for this study (Table S1). The raw dataset consisted of 270 populations sampled from 24 countries. A total of 2911 individuals were genotyped in the native range with a mean per population in the native range was 12.2 (range = 1–68). In the introduced range, we genotyped 40 individuals (Table S2). The estimates of variance for each population are based on independent samples (Wagner et al. 2005), therefore we randomly filtered locations sampled twice, leaving 238 populations for the downstream analysis (Figure 1B).

Next, to estimate distance‐dependent population variance which accounts for spatial autocorrelation (Wagner 2003), we applied variogram analysis, a method robust to sampling variance (Guillot et al. 2009). We fitted the variogram for a single locus (wsp) and analyzed wCer2 infection frequencies obtained for each population. A single locus variogram can be interpreted as the proportion of unlike links against distance (Wagner et al. 2005), giving a range at which statistical dependence disappears (Guillot et al. 2009). First, we computed an empirical semivariogram by the “variogram” function using the gstat package in R (Pebesma 2004) and modelled autocorrelation structure by fitting a theoretical variogram (Wagner et al. 2005). The variogram map and empirical variograms fitted for different isotropy axes confirmed a lack of anisotropy in our data (see Figure S1 and Supporting Information for more details). Therefore, we fit an exponential theoretical variogram to model autocorrelation due to stationary spatial processes following: γr=C0+C11−exp−3rb

where γr is a measure of dissimilarity between pairs of observations and C0 is the proportion of variance that is not spatially structured and C1 is the spatially structured variance component, providing an estimate of the population variance while accounting for spatial autocorrelation (Wagner et al. 2005).

2.3 Predicting Infection Frequencies in Unsampled Regions Across the Native Range

To predict infection frequency in unsampled regions across the R. cerasi native range, we performed spatial interpolation of the wCer2 infection frequency. We first created a grid of coordinates using the sf package in R (Pebesma 2018). Spatial interpolation was then computed using the “krige” function in the gstat package in R. Here, kriging utilizes the fitted theoretical variogram to interpolate values zxi+ri at any unobserved location based on a distance–variance relationship from the observed location, xi returning pairwise prediction values (Table S3). We then applied ordinary kriging with a constant intercept, z ~ 1, where z is the infection frequency of the wCer2, and validated the model (see Supporting Information for more details). The variogram analysis, aided by our 25 year dataset, coupled with predictions of infection frequency in unsampled regions, allows us to infer the region(s) in Europe from which the introduced population of R. cerasi in the United States may have originated.

3 Results

3.1 Rhagoletis cerasi Lacks the wCer2 Strain in the Introduced Range

The 40 R. cerasi sampled from the introduced population located in Niagara County, New York in the northeastern United States were uninfected with the wCer2 strain (p = 0 [0, 0.08]) (Figure 1A; Table 1). The field estimate of p was produced from sampling several locations within the continuous area of outbreak in Niagara County and are, therefore, representative of the region and considered a single population. In accordance with previous studies, the 40 individuals uninfected with the wCer2 strain were associated with the HT1 mitochondrial haplotype (Table 1) (Schuler et al. 2016; Schebeck et al. 2019).

3.2 wCer2 Infection Frequencies Vary Spatially Across the Native Range

Our field estimates of p showed fixation of wCer2 in areas of southwestern and Central Europe, which were associated with HT2 mitochondrial haplotype (Figure 1B; Table 1). One population is southwestern France located in Bayont (France 2, Table 1) showed a wCer2 infection frequency of 83.3% (p = 0.83 [0.552, 0.953]), but was not significantly different from the other three French populations (Fisher's exact test: P = 0.478). The same French population from Bayont, though wCer2 infected in high frequency, was exclusively associated with HT1 (Table 1). Populations in Turkey and Greece in the eastern portion of the range were wCer2 uninfected and were associated with HT1. We did not detect a single wCer2‐uninfected fly associated with the HT2 haplotype. The estimates of infection for population in Poland, Sweden, and Finland varied between p = 0.6 and p = 0.917, where some wCer2‐infected individuals were associated with HT1 instead of HT2. The easternmost population in Russia showed low wCer2 infection frequency (p = 0.167 [0.0470, 0.448]).

To test for temporal variation in the wCer2 frequencies, we compared infection frequencies of population sampled for this study to prior filed estimates of p across sampling years (Figure 2). Our pooled p estimate of French populations from this study was not significantly different from those of pooled regional estimates from 1999 (p = 1; Fisher's exact test: P = 1), 2000 (p = 1; Fisher's exact test: P = 1) and 2001 (p = 1; Fisher's exact test: P = 1). Resampling of the population in Austria showed no significant differences compared to previous sampling years (p = 1; Fisher's exact test: P = 1). The population in Greece uninfected in 1999 remained uninfected in 2000 and 2019, though with large credible intervals due to low sampling size (p = 0 [0, 0.658]). The populations sampled from Poland increased significantly in frequency from 2000 to 2019 (p = 0.286; Fisher's exact test: P = 0.045), but not from 2000 to 2007 (p = 0.5; Fisher's exact test: P = 0.172). Overall, this analysis revealed a general lack of temporal variation of the uninfected and fixed wCer2 populations (Figure 2). Together, these results enabled us to pool all field estimates of the wCer2 infection frequencies over the last 25 years and perform a spatial population prediction of infection frequencies in unobserved locations.

FIGURE 2 The wCer2 infection frequencies for different regions of the native Rhagoletis cerasi range. Data from the present study (orange point) are compared to prior field estimates. Error bars represent 95% binomial confidence intervals. An asterisk indicates a statistically significant difference between the present study and prior years at P < 0.05.

3.3 Identifying the wCer2‐Uninfected Native Regions as Likely Candidates for Invasion Origin

To assess the spatial structure of wCer2 in the native range, we computed distance‐dependent population variance (Figure 3A). The data points of the variogram were grouped into classes over the whole sampling range (0–6300 km) with a distance interval of 60 km. The semivariance first progressively increased up to the distance of about 200 km, after which it plateaued. For scales below 1000 km, the infection frequency was strongly spatially autocorrelated. After the cut‐off distance of 2000 km, additional information did not improve variance estimates.

FIGURE 3 Spatial population variance and interpolation of the wCer2 infection frequency in Rhagoletis cerasi native range predicting the regions from which North American populations of R. cerasi were introduced. (A) Semivariogram of the wCer2 infection frequency. The solid line indicates the fitted exponential variogram model. The dashed vertical line indicates where the curve reaches 95% of the estimated population variance. Data points are shown with a spatial lag distance of 60 km. The size of each point represents the number of point pairs whose separation in the lag interval was used to estimate the semivariance. (B) Spatial interpolation represented with contour lines showing the predicted wCer2 infection frequencies across the native range from 0% infection (dark blue) to 100% infection (yellow), including Russia (inset).

Results of the spatial interpolation analysis of infection frequencies predicted fixation of the wCer2 in R. cerasi populations from Central and Western Europe (Figure 3B) with high confidence (prediction variance ranged from 0.04 to 0.06) (Figure S2; Table S3) due to high sampling density in these regions. The Mediterranean and the eastern portions of the fly's native range including Eastern Europe and Turkey were predicted to be wCer2 uninfected (Figure 3B) with a prediction variance of below 10% (Figure S2; Table S3).

4 Discussion

We investigated the invasion of the European cherry fruit fly R. cerasi in North America using Wolbachia infection frequencies in the fly's native range as a marker to reconstruct the area of invasion origin. Our estimates of wCer2 infection frequencies (p) showed that the introduced R. cerasi population in the northeastern United States did not carry wCer2 and the associated mitochondrial haplotype suggesting a founder effect from wCer2‐uninfected native populations. Field estimates of the wCer2 infection frequency in the R. cerasi native range, combined with the field estimates of p spanning more than 238 populations sampled across 25 years, allowed us to assess distance‐dependent population variance, predict infection frequency in unsampled regions and rule out Central and Western Europe as likely source populations. Our results implicate a number of narrow regions in the Mediterranean and the eastern Europe as the most likely areas from which the flies were introduced into North America (Figure 3B). Our findings provide invaluable information about endosymbiont diversity in the introduced flies, shed light on potential source populations, and highlight that Wolbachia endosymbionts can be used as a tool to infer insect invasion origins.

The wCer2 infection frequencies varied spatially across the R. cerasi native range (Figure 1A; Table 1). We documented the fixation of the wCer2 in Central and Western Europe regions while observing a lack of the wCer2 in Eastern Europe and the Mediterranean (Figure 1A; Table 1). In addition, we observed a lack of temporal variation in infection frequencies across years (Figure 2) in wCer2 uninfected or fixed regions, which appeared to be stable over time. An exception from temporal stability among the fixed populations in Austria was the pooled p estimate from 2008 showing infection frequency less than 100%, though not significantly different from 1998 to 2000, 2015 and 2021 (Figure 2). We acknowledge temporal variation in this region may be due to low Wolbachia titer at concentrations below the detection threshold. However, low densities of Wolbachia have been associated with alternative reproductive strategies of the endosymbiont (Ijichi et al. 2002), often resulting in a lack of CI induction. We observe strong levels of CI in the wCer2 strain, however, suggesting that low titer is unlikely responsible for the pattern. In Drosophila, both environmental factors and/or a variation in maternal transmission have been hypothesized to play a role in temporal fluctuations in Wolbachia infection frequencies (Turelli and Hoffmann 1995; Hague et al. 2020).

Similarly, the 2019 population from Greece harboring only two wCer2‐uninfected individuals had large binomial confidence intervals, but was not significantly different from the two other wCer2‐uninfected Greek populations. Due to low sampling size, it is possible that we failed to detect infection in the 2019 population. However, this scenario is unlikely considering a vast uninfected area surrounding the Greek populations (Figure 2). Our estimate of p for the population from Poland in 2019 was significantly higher compared to the pooled estimate from 2000 (Figure 2). This pattern supports the hypothesis of ongoing spread of the wCer2 in R. cerasi populations in central Europe previously documented by Riegler and Stauffer (2002), Schuler et al. (2016) and quantified by Bakovic et al. (2018). Including the wCer2 infection frequencies in the areas of active spread across the fly's native range further allowed us to “fine tune” the spatial interpolation analysis. Together, the pattern of spatial variation of the wCer2 infection frequencies across the R. cerasi native range and a general temporal stability of the uninfected and fixed wCer2 populations (Figure 2) enabled us to pool all field estimates of the wCer2 infection frequencies over the last 25 years to predict infection frequencies in unobserved locations.

Based on the spatial interpolation prediction of the wCer2 infection frequency, we argue that R. cerasi unlikely originated from Central and Western Europe where fly populations are either fixed for wCer2 or are infected at very high frequencies (Figure 3). A deviation from this pattern is a narrow belt in central Germany (Figure 3B), accounting for a mosaic pattern of the wCer2 infection in this region (Schuler et al. 2016; Schebeck et al. 2019). We propose regions that harbor wCer2‐uninfected individuals as the most likely candidates for invasion source, originating from Eastern Europe and the Mediterranean region. Concurrently, field estimates and spatial interpolation analysis, which take into account sampled regions and predict infection frequencies at unsampled locations, showed that these regions lack the wCer2 strain and could have served as source population(s) (Figure 1A; Figure 3; Table 1). These populations were also associated with HT1 found in R. cerasi in the introduced range (Table 1), strongly associated with the fixed wCer1 strain (Schuler et al. 2016; Schebeck et al. 2019).

In Bakovic et al. (2018), the authors estimated the speed of the wCer2 spatial spread in the native range of R. cerasi. However, the nature of wCer2's local dynamics and spatial spread remains unknown. The mode of Wolbachia local dynamics in natural populations can take two forms depending on the value of effective fecundity. If F (1–μ) < 1, as in Aedes aegypti mosquitoes, Wolbachia will stabilize at 0, when present at very low frequency, following the bistable local dynamics (Hoffmann et al. 2014). Conversely, when F (1–μ) > 1, exemplified by the wRi strain in D. simulans, even very low numbers of Wolbachia are sufficient to push infection frequency to high stable equilibrium following Fisherian dynamics (Hoffmann, Turelli, and Harshman 1990; Barton and Turelli 2011). If we assume that wCer2 infection follows bistable local dynamics, the fact that we failed to detect the presence of the HT2 mitochondrial haplotype associated with the wCer2 strain in the introduced range strongly suggests that wCer2 was unlikely from an infected native population in low numbers that stabilized at p = 0 at the time of sampling. If, however, wCer2 dynamics behave according to Fisherian dynamics, even invading in low initial numbers in the introduced range would have driven wCer2 to high frequencies, unless impeded by differing selection pressures in the newly invaded habitat similar to Wolbachia in the Argentine ant, Linepithema humile (Reuter, Pedersen, and Keller 2005). In the latter case, given the strong association between wCer2 and HT2, we would expect to detect this mitochondrial haplotype in the invasive range as a remnant of such an event. Taken together, evidence strongly supports a founder effect scenario where native wCer2‐uninfected flies served as the founding individuals.

We did detect deviations from wCer2/HT2 association in the native range, namely in transitional populations from France, Poland, Sweden, and Finland, where field estimates of p varied between 0.60 and 0.917. Here, some wCer2‐infected individuals were associated with HT1 (Table 1). Deviations from this association were highlighted in Schuler et al. (2016) who detected wCer2/HT1 co‐occurrence in transitional German populations with a frequency of 21.9% (Schebeck et al. 2019). The incidence of the association between wCer2 and HT1 could be the result of either paternal (e.g. Hoffmann and Turelli 1988) or intraspecific horizontal transmission (Schuler et al. 2016). Using mathematical modeling of Wolbachia and mitochondrial haplotype dynamics, Schuler et al. (2016) showed that the association of wCer2 with HT1 is transient, likely due to weak maternal inheritance of wCer2 and will be lost after complete invasion of wCer2. Their analysis also showed that this association, if transmitted only somatically, would be lost after fixation of the maternally transmitted wCer2. These low levels of misassociation between the Wolbachia strain infecting the fly and fly's mitochondrial haplotype in the native range suggest that the chances that the invasion originated from wCer2/HT1 individuals in transitional populations, followed by a subsequent loss of wCer2, are low. There is a possibility that wCer2 is present in the introduced range, but remained undetected in our sample. We find this scenario unlikely considering a large sample size (n = 40 individuals) along with the fact that our sampling method included several locations within the continuous area of in Niagara County, therefore providing a solid representation of the invaded region.

Our current data lack the sensitivity to detect multiple invasion origins of R. cerasi into North America. Many successful biological invasions are accomplished through multiple introductions to the extent that single introductions are often considered exceptions (Dlugosch et al. 2015; Caizergues et al. 2024). Recurrent introductions can alleviate the adverse effects of demographic bottlenecks associated with founder effects and facilitate the establishment of invasive species in the new range (Verhoeven et al. 2010). However, negative impacts of multiple sources of introduction have also been observed when gene flow between genetically distinct source populations in the invasive range leads to the breakdown of locally adapted genotypes (Rius and Darling 2014; Barker et al. 2019). Future studies will aim to incorporate population‐level genomic variation of the flies across their natural range to corroborate our Wolbachia results.

In a similar study, Zhang et al. (2014) used the obligatory mutualistic endosymbiont Buchnera aphidicola to trace invasion routes in the clonal Russian wheat aphid. A small number of studies have also used Wolbachia to infer demographic history during colonization such as the Argentine ant, Linepithema humile (Reuter, Pedersen, and Keller 2005), and fire ants (Solenopsis spp.) (Shoemaker et al. 2000). This study also utilizes Wolbachia as a marker to focus on the question of invasion origin of R. cerasi in North America, but achieves a high resolution by leveraging an extensive and thorough range‐wide sampling of the native range. In contrast, classic population genetic markers used in invasion biology studies include microsatellites, and mtDNA that often have limited resolution to detect the origin of the invasion (Zhang et al. 2014). In insects, the interpretation of mitochondria‐based biogeographic and colonization history analysis is confounded due to the linkage with maternally transmitted endosymbionts (Hurst and Jiggins 2005; Arif et al. 2021). Our study confirms that when used together, Wolbachia and mtDNA markers can shed light on demographic histories and invasion processes while representing an inexpensive and rapid alternative if genomic data are not available at hand. Given that increasing number of studies report spatial and temporal patterns of Wolbachia infection frequencies (Hoffmann and Turelli 1988; Shoemaker et al. 2003; Kriesner et al. 2016; Hague et al. 2020) and a high prevalence of Wolbachia in insects (Werren and Windsor 2000; Jiggins et al. 2001), as well as other arthropods and nematodes (Weinert et al. 2015), the use of spatially varying Wolbachia‐based markers open up a promising avenue in elucidating invasion origins and a better understanding of the factors governing introductions. Additionally, markers such as SNPs and structural variants that allow researchers to include within‐strain Wolbachia variation (Wolfe et al. 2021) could heighten the resolution and help flesh out invasion scenarios at an even finer scale.

5 Conclusion

Our sampling of the R. cerasi introduced range revealed that the flies lacked the wCer2 strain and the associated mitochondrial haplotype. Our estimates of wCer2 infection frequency varied spatially across the native range. The extensive sampling on the R. cerasi native region enabled us to quantify population variance in infection frequencies, rule out Central and Western Europe of R. cerasi range as the source population, and propose wCer2‐uninfected R. cerasi populations from Eastern Europe and the Mediterranean region as the most likely candidates for the invasion origin. Our study highlights the utility of spatially varying Wolbachia frequencies to infer species invasion origins and demographic histories. More broadly, it puts forward endosymbionts as powerful tools to reconstruct invasion histories.

Conflicts of Interest

The authors declare no conflicts of interest.

Supporting information

Appendix S1.

Appendix S2.

Appendix S3.

Acknowledgments

We thank Norman Barr from the United States Department of Agriculture (USDA) for help conceptualizing the research, Susanne Krumböck (BOKU, Vienna, Austria) for support with molecular analysis, Chris Logue and Margaret Kelly of the New York State Department of Agriculture and Markets (NYSDAM) for access to samples, and the Wayne State University undergraduates Kelvin Selegean and Connor Socrates for help with DNA extractions. The work is part of a cooperative agreement between the Wayne State University and USDA‐APHIS‐PPQ (AP22PPQS&T00C053). Additional funding was provided by the FWF projects P35327‐B to CS and 31441‐B29 to HS. We also thank the Department of Innovation and Research and the University of the Autonomous Province of Bozen/Bolzano for funding the open access publication charges. Mention of trade names or commercial products in this publication is solely for purpose of providing specific information and does not imply recommendation or endorsement by the USDA.

Data Availability Statement

The raw data are available in the supplemental data files. The R code used for statistical analysis, spatial modeling, and constructing the figures can be found in the following link: https://github.com/slecic/USA_invasion_paper.
==== Refs
References

Arif, S. , M. Gerth , W. G. Hone‐Millard , M. D. S. Nunes , L. Dapporto , and T. G. Shreeve . 2021. “Evidence for Multiple Colonisations and Wolbachia Infections Shaping the Genetic Structure of the Widespread Butterfly Polyommatus icarus in the British Isles.” Molecular Ecology 30 , no. 20 : 5196–5213. 10.1111/mec.16126.34402109
Arthofer, W. , M. Riegler , D. Schneider , M. Krammer , W. J. Miller , and C. Stauffer . 2009. “Hidden Wolbachia Diversity in Field Populations of the European Cherry Fruit fly, Rhagoletis cerasi (Diptera, Tephritidae).” Molecular Ecology 18 , no. 18 : 3816–3830. 10.1111/j.1365-294X.2009.04321.x.19732336
Bakovic, V. , M. Schebeck , A. Telschow , C. Stauffer , and H. Schuler . 2018. “Spatial Spread of Wolbachia in Rhagoletis cerasi Populations.” Biology Letters 14 , no. 5 : 20180161. 10.1098/rsbl.2018.0161.29794009
Barker, B. S. , J. E. Cocio , S. R. Anderson , et al. 2019. “Potential Limits to the Benefits of Admixture During Biological Invasion.” Molecular Ecology 28 , no. 1 : 100–113. 10.1111/mec.14958.30485593
Barringer, L. 2018. “First Record of the European Cherry Fruit Fly, Rhagoletis cerasi (Linnaeus) (Diptera: Tephritidae), in North America.” Insecta Mundi 0622 : 1–4. 10.5281/zenodo.5171480.
Barton, N. H. , and M. Turelli . 2011. “Spatial Waves of Advance With Bistable Dynamics: Cytoplasmic and Genetic Analogues of Allee Effects.” American Naturalist 178 , no. 3 : E48–E75. 10.1086/661246.
Bock, D. G. , C. Caseys , R. D. Cousens , et al. 2015. “What We Still Don't Know About Invasion Genetics.” Molecular Ecology 24 , no. 9 : 2277–2297. 10.1111/mec.13032.25474505
Boller, E. F. , and R. J. Prokopy . 1976. “Bionomics and Management of Rhagoletis .” Annual Review of Entomology 21 , no. 1 : 223–246. 10.1146/annurev.en.21.010176.001255.
Boller, E. F. , K. Russ , V. Vallo , and G. L. Bush . 1976. “Incompatible Races of European Cherry Fruit Fly, Rhagoletis cerasi (Diptera: Tephritidae), Their Origin and Potential Use in Biological Control.” Entomologia Experimentalis et Applicata 20 , no. 3 : 237–247. 10.1111/j.1570-7458.1976.tb02640.x.
Caizergues, A. E. , J. S. Santangelo , R. W. Ness , et al. 2024. “Does Urbanization Lead to Parallel Demographic Shifts Across the World in a Cosmopolitan Plant?” Molecular Ecology 33 , no. 7 : E17311. 10.1111/mec.17311.38468155
Deshpande, J. N. , and E. A. Fronhofer . 2022. “Genetic Architecture of Dispersal and Local Adaptation Drives Accelerating Range Expansions.” Proceedings of the National Academy of Sciences of the United States of America 119 , no. 31 : e2121858119. 10.1073/pnas.2121858119.35895682
Dlugosch, K. M. , S. R. Anderson , J. Braasch , F. A. Cang , and H. D. Gillette . 2015. “The Devil is in the Details: Genetic Variation in Introduced Populations and Its Contributions to Invasion.” Molecular Ecology 24 , no. 9 : 2095–2111. 10.1111/mec.13183.25846825
Dlugosch, K. M. , and I. M. Parker . 2008. “Invading Populations of an Ornamental Shrub Show Rapid Life History Evolution Despite Genetic Bottlenecks.” Ecology Letters 11 , no. 7 : 701–709. 10.1111/j.1461-0248.2008.01181.x.18410377
Early, R. , B. A. Bradley , J. S. Dukes , et al. 2016. “Global Threats From Invasive Alien Species in the Twenty‐First Century and National Response Capacities.” Nature Communications 7 , no. 1 : 12485. 10.1038/ncomms12485.
Estoup, A. , and T. Guillemaud . 2010. “Reconstructing Routes of Invasion Using Genetic Data: Why, How and So What?” Molecular Ecology 19 , no. 19 : 4113–4130. 10.1111/j.1365-294X.2010.04773.x.20723048
Guillot, G. , R. Leblois , A. Coulon , and A. C. Frantz . 2009. “Statistical Methods in Spatial Genetics.” Molecular Ecology 18 , no. 23 : 4734–4756. 10.1111/j.1365-294X.2009.04410.x.19878454
Hague, M. T. J. , H. Mavengere , D. R. Matute , and B. S. Cooper . 2020. “Environmental and Genetic Contributions to Imperfect wMel‐Like Wolbachia Transmission and Frequency Variation.” Genetics 215 , no. 4 : 1117–1132. 10.1534/genetics.120.303330.32546497
Harrell, F. E., Jr. , and C. Dupont . 2018. “Hmisc: Harrell Miscellaneous.”
Hoffmann, A. A. , I. Iturbe‐Ormaetxe , A. G. Callahan , et al. 2014. “Stability of the wMel Wolbachia Infection Following Invasion Into Aedes aegypti Populations.” PLoS Neglected Tropical Diseases 8 , no. 9 : e3115. 10.1371/journal.pntd.0003115.25211492
Hoffmann, A. A. , and M. Turelli . 1988. “Unidirectional Incompatibility in Drosophila simulans: Inheritance, Geographic Variation and Fitness Effects.” Genetics 119 , no. 2 : 435–444. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC1203425/.17246430
Hoffmann, A. A. , M. Turelli , and L. G. Harshman . 1990. “Factors Affecting the Distribution of Cytoplasmic Incompatibility in Drosophila simulans .” Genetics 126 , no. 4 : 933–948. 10.1093/genetics/126.4.933.2076821
Hoshizaki, S. , and T. Shimada . 1995. “PCR Based Detection of Wolbachia, Cytoplasmic Incompatibility Microorganisms, Infected in Natural Populations of Laodelphax striatellus (Homoptera: Delphacidae) in Central Japan: Has the Distribution of Wolbachia Spread Recently?” Insect Molecular Biology 4 , no. 4 : 237–243. 10.1111/j.1365-2583.1995.tb00029.x.8825761
Hurst, G. D. D. , and F. M. Jiggins . 2005. “Problems With Mitochondrial DNA as a Marker in Population, Phylogeographic and Phylogenetic Studies: The Effects of Inherited Symbionts.” Proceedings of the Royal Society B: Biological Sciences 272 , no. 1572 : 1525–1534. 10.1098/rspb.2005.3056.
Ijichi, N. , N. Kondo , R. Matsumoto , M. Shimada , H. Ishikawa , and T. Fukatsu . 2002. “Internal Spatiotemporal Population Dynamics of Infection With Three Wolbachia Strains in the Adzuki Bean Beetle, Callosobruchus chinensis (Coleoptera: Bruchidae).” Applied and Environmental Microbiology 68 , no. 8 : 4074–4080. 10.1128/AEM.68.8.4074-4080.2002.12147509
Jiggins, F. M. , J. K. Bentley , M. E. Majerus , and G. D. Hurst . 2001. “How Many Species Are Infected With Wolbachia? Cryptic Sex Ratio Distorters Revealed to be Common by Intensive Sampling.” Proceedings. Biological sciences 268 , no. 1472 : 1123–1126. 10.1098/rspb.2001.1632.11375098
Jiggins, F. M. , and M. C. Tinsley . 2005. “An Ancient Mitochondrial Polymorphism in Adalis bipunctata Linked to a Sex‐Ratio‐Distorting Bacterium.” Genetics 171 , no. 3 : 1115–1124. 10.1534/genetics.105.046342.16079227
Kriesner, P. , W. R. Conner , A. R. Weeks , M. Turelli , and A. A. Hoffmann . 2016. “Persistence of a Wolbachia Infection Frequency Cline in Drosophila melanogaster and the Possible Role of Reproductive Dormancy.” Evolution 70 , no. 5 : 979–997. 10.1111/evo.12923.27076356
Kriesner, P. , A. A. Hoffmann , S. F. Lee , M. Turelli , and A. R. Weeks . 2013. “Rapid Sequential Spread of Two Wolbachia Variants in Drosophila simulans .” PLoS Pathogens 9 , no. 9 : e1003607. 10.1371/journal.ppat.1003607.24068927
Lu, M. , J. Hulcr , and J. Sun . 2016. “The Role of Symbiotic Microbes in Insect Invasions.” Annual Review of Ecology, Evolution, and Systematics 47 , no. 1 : 487–505. 10.1146/annurev-ecolsys-121415-032050.
Pebesma, E. 2018. “Simple Features for R: Standardized Support for Spatial Vector Data.” R Journal 10 , no. 1 : 439. 10.32614/RJ-2018-009.
Pebesma, E. J. 2004. “Multivariable Geostatistics in S: The Gstat Package.” Computers & Geosciences 30 , no. 7 : 683–691. 10.1016/j.cageo.2004.03.012.
Prentis, P. J. , J. R. U. Wilson , E. E. Dormontt , D. M. Richardson , and A. J. Lowe . 2008. “Adaptive Evolution in Invasive Species.” Trends in Plant Science 13 , no. 6 : 288–294. 10.1016/j.tplants.2008.03.004.18467157
Pyšek, P. , P. E. Hulme , D. Simberloff , et al. 2020. “Scientists' Warning on Invasive Alien Species.” Biological Reviews 95 , no. 6 : 1511–1534. 10.1111/brv.12627.32588508
R Core Team . 2020. R: A Language and Environment for Statistical Computing. Vienna, Austria: R Foundation for Statistical Computing.
Reuter, M. , J. S. Pedersen , and L. Keller . 2005. “Loss of Wolbachia Infection During Colonisation in the Invasive Argentine Ant Linepithema humile .” Heredity 94 , no. 3 : 364–369. 10.1038/sj.hdy.6800601.15674380
Riegler, M. , and C. Stauffer . 2002. “ Wolbachia Infections and Superinfections in Cytoplasmically Incompatible Populations of the European Cherry Fruit Fly Rhagoletis cerasi (Diptera, Tephritidae).” Molecular Ecology 11 , no. 11 : 2425–2434. 10.1046/j.1365-294x.2002.01614.x.12406252
Rius, M. , and J. A. Darling . 2014. “How Important is Intraspecific Genetic Admixture to the Success of Colonising Populations?” Trends in Ecology & Evolution 29 , no. 4 : 233–242. 10.1016/j.tree.2014.02.003.24636862
Schebeck, M. , L. Feldkirchner , C. Stauffer , and H. Schuler . 2019. “Dynamics of an Ongoing Wolbachia Spread in the European Cherry Fruit fly, Rhagoletis cerasi (Diptera: Tephritidae).” Insects 10 , no. 6 : 172. 10.3390/insects10060172.31208002
Schuler, H. , K. Köppler , S. Daxböck‐Horvath , et al. 2016. “The Hitchhiker's Guide to Europe: The Infection Dynamics of an Ongoing Wolbachia Invasion and Mitochondrial Selective Sweep in Rhagoletis cerasi .” Molecular Ecology 25 , no. 7 : 1595–1609. 10.1111/mec.13571.26846713
Sherpa, S. , and L. Després . 2021. “The Evolutionary Dynamics of Biological Invasions: A Multi‐Approach Perspective.” Evolutionary Applications 14 , no. 6 : 1463–1484. 10.1111/eva.13215.34178098
Shoemaker, D. D. , M. Ahrens , L. Sheill , M. Mescher , L. Keller , and K. G. Ross . 2003. “Distribution and Prevalence of Wolbachia Infections in Native Populations of the Fire Ant Solenopsis invicta (Hymenoptera: Formicidae).” Environmental Entomology 32 , no. 6 : 1329–1336. 10.1603/0046-225X-32.6.1329.
Shoemaker, D. D. , K. G. Ross , L. Keller , E. L. Vargo , and J. H. Werren . 2000. “ Wolbachia Infections in Native and Introduced Populations of Fire Ants (Solenopsis spp.).” Insect Molecular Biology 9 , no. 6 : 661–673. 10.1046/j.1365-2583.2000.00233.x.11122476
Szűcs, M. , B. A. Melbourne , T. Tuff , C. Weiss‐Lehman , and R. A. Hufbauer . 2017. “Genetic and Demographic Founder Effects Have Long‐Term Fitness Consequences for Colonising Populations.” Ecology Letters 20 , no. 4 : 436–444. 10.1111/ele.12743.28145080
Tepolt, C. K. , E. D. Grosholz , C. E. de Rivera , and G. M. Ruiz . 2022. “Balanced Polymorphism Fuels Rapid Selection in an Invasive Crab Despite High Gene Flow and Low Genetic Diversity.” Molecular Ecology 31 , no. 1 : 55–69. 10.1111/mec.16143.34431151
Turelli, M. , and A. A. Hoffmann . 1991. “Rapid Spread of an Inherited Incompatibility Factor in California Drosophila.” Nature 353 , no. 6343 : 440–442. 10.1038/353440a0.1896086
Turelli, M. , and A. A. Hoffmann . 1995. “Cytoplasmic Incompatibility in Drosophila simulans: Dynamics and Parameter Estimates From Natural Populations.” Genetics 140 , no. 4 : 1319–1338. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC1206697/.7498773
Turelli, M. , A. A. Hoffmann , and S. W. McKechnie . 1992. “Dynamics of Cytoplasmic Incompatibility and mtDNA Variation in Natural Drosophila simulans Populations.” Genetics 132 , no. 3 : 713–723. 10.1093/genetics/132.3.713.1468627
Turelli, M. , A. Katznelson , and P. S. Ginsberg . 2022. “Why Wolbachia‐Induced Cytoplasmic Incompatibility is so Common.” Proceedings of the National Academy of Sciences of the United States of America 119 , no. 47 : e2211637119. 10.1073/pnas.2211637119.36343219
Verhoeven, K. J. F. , M. Macel , L. M. Wolfe , and A. Biere . 2010. “Population Admixture, Biological Invasions and the Balance Between Local Adaptation and Inbreeding Depression.” Proceedings of the Royal Society B: Biological Sciences 278 , no. 1702 : 2–8. 10.1098/rspb.2010.1272.
Wagner, H. H. 2003. “Spatial Covariance in Plant Communities: Integrating Ordination, Geostatistics, and Variance Testing.” Ecology 84 , no. 4 : 1045–1057. 10.1890/0012-9658(2003)084[1045:SCIPCI]2.0.CO;2.
Wagner, H. H. , R. Holderegger , S. Werth , F. Gugerli , S. E. Hoebee , and C. Scheidegger . 2005. “Variogram Analysis of the Spatial Genetic Structure of Continuous Populations Using Multilocus Microsatellite Data.” Genetics 169 , no. 3 : 1739–1752. 10.1534/genetics.104.036038.15654102
Wakie, T. T. , W. L. Yee , and L. G. Neven . 2018. “Assessing the Risk of Establishment of Rhagoletis cerasi (Diptera: Tephritidae) in the United States and Globally.” Journal of Economic Entomology 111 , no. 3 : 1275–1284. 10.1093/jee/toy054.29528438
Weinert, L. A. , E. V. Araujo‐Jnr , M. Z. Ahmed , and J. J. Welch . 2015. “The Incidence of Bacterial Endosymbionts in Terrestrial Arthropods.” Proceedings of the Royal Society B: Biological Sciences 282 , no. 1807 : 20150249. 10.1098/rspb.2015.0249.
Werren, J. H. , and D. M. Windsor . 2000. “ Wolbachia Infection Frequencies in Insects: Evidence of a Global Equilibrium?” Proceedings of the Royal Society B: Biological Sciences 267 , no. 1450 : 1277–1285. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC1690679/.
White, T. A. , S. E. Perkins , G. Heckel , and J. B. Searle . 2013. “Adaptive Evolution During an Ongoing Range Expansion: The Invasive Bank Vole (Myodes glareolus) in Ireland.” Molecular Ecology 22 , no. 11 : 2971–2985. 10.1111/mec.12343.23701376
Wolfe, T. M. , D. J. Bruzzese , L. Klasson , et al. 2021. “Comparative Genome Sequencing Reveals Insights Into the Dynamics of Wolbachia in Native and Invasive Cherry Fruit Flies.” Molecular Ecology 30 , no. 23 : 6259–6272. 10.1111/mec.15923.33882628
Xiao, J.‐H. , N.‐X. Wang , R. W. Murphy , J. Cook , L.‐Y. Jia , and D.‐W. Huang . 2012. “ Wolbachia Infection and Dramatic Intraspecific Mitochondrial DNA Divergence in a Fig Wasp.” Evolution 66 , no. 6 : 1907–1916. 10.1111/j.1558-5646.2011.01561.x.22671555
Zhang, B. , O. Edwards , L. Kang , and S. Fuller . 2014. “A Multi‐Genome Analysis Approach Enables Tracking of the Invasion of a Single Russian Wheat Aphid (Diuraphis noxia) Clone Throughout the New World.” Molecular Ecology 23 , no. 8 : 1940–1951. 10.1111/mec.12714.24617944
