Skip to main content

Worldwide population genetic structure of the oriental fruit moth (Grapholita molesta), a globally invasive pest



Invasive pest species have large impacts on agricultural crop yields, and understanding their population dynamics is important for ensuring food security. The oriental fruit moth Grapholita molesta is a cosmopolitan pest of stone and pome fruit species including peach and apple, and historical records indicate that it has invaded North and South America, Europe, Australia and Africa from its putative native range in Asia over the past century.


We used 13 microsatellite loci, including nine newly developed markers, to characterize global population structure of G. molesta. Approximately 15 individuals from each of 26 globally distributed populations were genotyped. A weak but significant global pattern of isolation-by-distance was found, and G. molesta populations were geographically structured on a continental scale. Evidence does not support that G. molesta was introduced to North America from Japan as previously proposed. However, G. molesta was probably introduced from North America to The Azores, South Africa, and Brazil, and from East Asia to Australia. Shared ancestry was inferred between populations from Western Europe and from Brazil, although it remains unresolved whether an introduction occurred from Europe to Brazil, or vice versa. Both genetic diversity and levels of inbreeding were surprisingly high across the range of G. molesta and were not higher or lower overall in introduced areas compared to native areas. There is little evidence for multiple introductions to each continent (except in the case of South America), or for admixture between populations from different origins.


Cross-continental introductions of G. molesta appear to be infrequent, which is surprising given its rapid worldwide expansion over the past century. We suggest that area-wide spread via transport of fruits and other plant materials is a major mechanism of ongoing invasion, and management efforts should therefore target local and regional farming communities and distribution networks.


Global interest in the ecology of agricultural systems has driven a substantial amount of fundamental research into plant-herbivore interactions, modes of plant (i.e. weed) and invertebrate dispersal and adaptation, and the role of multi-trophic interactions in community dynamics[1]. A growing need to understand the potential effects of climate change on cash crop production systems[2], combined with the necessity of increasing worldwide crop yields in response to projected global population growth[3], continue to motivate research into fundamental and applied questions in ecology and evolutionary biology.

The mechanisms that facilitate the establishment and spread of invasive pests are of particular interest to ecologists, evolutionary biologists and policy makers for a number of reasons; agricultural pests are responsible for substantial yield losses in agricultural systems, and there are considerable economic incentives to understand mechanisms of pest invasions. Also, because invasions often occur relatively rapidly, they provide useful case studies to understand the ecological factors (both biotic and abiotic) that mediate population dynamics and the establishment of range limits[4, 5]. Finally, pest species often adapt to new selection pressures in their invaded ranges, and therefore represent contemporary examples of rapid evolutionary change[5, 6].

The oriental fruit moth Grapholita (=Cydia) molesta Busck (Lepidoptera: Tortricidae) is one example of an insect pest species that continues to cause significant economic damage to crops on a global scale[7]. Grapholita molesta is a major pest of stone and pome fruit species mainly belonging to the Rosaceae family, including peach, apple, pear, nectarines, cherries, quince, and persimmons (Ebenaceae)[8, 9]. Assumed to be native to China[7], G. molesta is now distributed throughout temperate regions of Asia, Europe, The Americas, Africa, and Australia.

The global demographic history of G. molesta has not been accurately traced, although historical records permit the reconstruction of a rough timeline of continental invasions. In Europe, G. molesta was first recorded in Slovenia in 1920[10], in southeastern France and much of north-central Italy in the early 1920s[11], and has since dispersed throughout eastern, southern and western Europe, where stone fruit trees are grown. It is commonly assumed that the species was introduced to North America via a fruit shipment to Washington D.C. from Japan around 1913[12], however this assumption is based on the anecdotal report of a single G. molesta specimen recovered from a shipment of Japanese pears[12]. After its introduction to Washington D.C., G. molesta is presumed to have dispersed to the neighboring states of Virginia and Maryland, and northward to Ontario, Canada. By the mid-1940s, it had spread from the mid-West to California, Washington state, and Oregon[13] and references therein]. In Australia, G. molesta was introduced around 1910, and has since spread to most of the stone and pome fruit growing regions on that continent[14] and references therein]. The species is assumed to have been present in southern Brazil since at least the 1940s[15]. In other parts of its range, introductions are thought to be more recent; G. molesta was first identified in New Zealand in 1976[16], and in South Africa in 1980[17].

A number of regional studies have yielded some information about local factors that contribute to the dispersal and population structuring of G. molesta over regional spatial scales. Population structure of G. molesta in a major fruit-growing region of Northern Italy was relatively low[18], and could be explained by both natural dispersal[19] and anthropogenic displacement of individuals between orchards[20]. Based on molecular genetic analysis of six populations in South Africa, it was also concluded that anthropogenic movement of fruit, bins, and nursery material likely drives range expansion at a regional level[21]. Another study showed that landscape features such as rivers might act as ecological barriers to regional dispersal[22]. In spite of substantial effort to understand the lifecycle of G. molesta and to develop strategies for its control, little is known about its ongoing global dispersal patterns.

In this study, we used a set of 13 microsatellite loci (including nine newly developed ones) to characterize global population structure of G. molesta, in order to retrace invasion routes, and to investigate global population dynamics. We hypothesized that i) patterns of population structure will reflect evidence of multiple introductions via human-mediated dispersal in the invaded range of the species, and ii) patterns of isolation-by-distance will be absent or weak, since international trade routes likely determine patterns of global population structure, and iii) genetic diversity will be reduced in introduced compared to native populations. We also tested whether populations are structured according to host species (predominantly peach and apple).


Insect sampling

A total of 376 G. molesta larvae (sampled from fruits or shoots) or adult males (captured from pheromone traps) were collected from fruit trees at 26 worldwide locations between April and October 2011 (see Table 1; Figure 1). DNA has been stored, and voucher samples are deposited at the Entomological Collection of ETH Zurich. We sampled 15 individuals per site in most cases (Table 1). Upon sampling, geographical coordinates and host plant species were recorded. To maximize the genetic variability of larvae sampled within sites, we collected larvae from as many different trees as possible. Samples were stored in 70% ethanol.

Table 1 Sampling information
Figure 1

Locations of 26 globally sampled Grapholita molesta populations. Approximately 15 individuals were sampled from each population, and were genotyped at 13 microsatellite loci.

DNA Extraction and Simple Sequence Repeat (SSR) genotyping

Genomic DNA was extracted from 8–10 mg larvae or adult male moths (heads and thoraxes) (see Table 1) using the Nucleospin Genomic DNA from Tissue Kit (Macherey-Nagel AG, Oensingen, Switzerland). Extraction was performed according to the bench protocol for animal tissues. In the final step, DNA was eluted in 200 μL TE buffer and stored at −20°C.

A total of 13 SSR loci were used for genotyping. Primer sequences and reaction conditions for the amplification of four of these loci were previously described[18], and are given in Table 2. A new set of nine SSR markers was developed by ecogenics GmbH (Zurich-Schlieren, Switzerland) using the high-throughput genomic sequencing approach described by[23]. 1.25 μL of genomic DNA were analyzed on a Roche 454 GS-FLX platform (Roche, Basel, Switzerland) using a 1/16th run and the GS-FLX titanium reagents. Of these, 2 183 contained a microsatellite insert with a tetra- or a trinucleotide of at least six repeat units or a dinucleotide of at least 10 repeat units. Suitable primer design was possible in 564 reads, of which 36 were tested for polymorphism in fifteen unrelated individuals derived from geographically disparate populations. Of these, 27 were eliminated because they showed failures or high null allele frequencies, were not polymorphic, or had complex peak patterns that could not be reliably interpreted. Characteristics of the newly developed microsatellite DNA markers are given in Table 2.

Table 2 Characteristics of 13 microsatellite loci in Grapholita molesta

All polymerase chain reactions (PCRs) were performed as singleplex reactions. Reactions were performed in a Labcycler thermocycler (Sensoquest GmbH, Göttingen, Germany) in a total reaction volume of 10 μL. The forward primers of the loci GM02, GM05, GM07, and GM10 were fluorescently labeled with joe, fam, fam, and hex respectively. Reactions for these four loci were carried out in a mix of 5.4 μL ddH2O, 1.0 μL PCR buffer (Qiagen, Hombrechtikon, Switzerland), 2.0 μL dNTPs (2.5 mM, Qiagen), 0.5 μL of each forward and reverse primer (2 μM), 0.1 μL HotstarTaq (5U, Qiagen) and 0.5 μL DNA (10–35 ng⁄ μL). Reactions for the remaining nine newly developed loci were carried out in a mix of 5.7 μL ddH2O, 1.0 μL PCR buffer (Qiagen), 1.0 μL dNTPs (2.5 mM, Qiagen), 0.2 μL of each forward primer (2 μM), 0.5 μL of each reverse primer (2 μM), 0.5 μL M13 primer (2 μM, fluorescently labeled with fam, Microsynth, Balgach, Switzerland), 0.1 μL HotstarTaq (5U, Qiagen) and 0.5 μL DNA (10–35 ng⁄ μL). Locus-specific PCR conditions are reported in Additional file1: Appendix 1.

Genotyping was carried out by denaturing 1 μL PCR product in 9 μL formamide (deionized, 99.5% minimum, Sigma-Aldrich, Buchs, Switzerland) with 0.09 μL Gene Scan 500 LIZ size standard (Applied Biosystems, Foster City, CA, USA), and running the products on a laser detection system (3730xl DNA Analyzer, Applied Biosystems). DNA sizing and allele definitions were performed using Gene Mapper 4.0 software (Applied Biosystems). SSR data are archived at the Dryad repository (

Statistical Analyses

We identified departures from Hardy–Weinberg equilibrium (HWE) using Gen Al Ex v. 6.41[24]. We calculated null allele frequencies using the R package Geneland 3.1.4[25, 26]. We tested for linkage disequilibrium between all pairs of loci using Arlequin 3.01[27], according to[28].

Population genetic structure

Genetic and geographic distances were calculated in Gen Al Ex using a Euclidean distance metric according to[29], and missing data were interpolated by the software. A test of isolation-by-distance was carried out by applying a Mantel’s test to matrices of genetic and geographic distances in Gen Al Ex using 999 permutations.

To better understand the partitioning of regional population genetic structure, we carried out an AMOVA in Gen Al Ex, using 999 permutations, and also calculated global (mean) FST value between populations[30]. We defined populations according to the sample sites listed in Table 1, and we defined eight separate regions: North America, Europe (excluding The Azores), Asia (excluding Japan), South America, Australia, South Africa, Japan, and the Azores. Japan and the Azores were included as separate regions because we hypothesized that island populations may be significantly differentiated from mainland populations. In order to determine whether host plant species accounted for a portion of the genetic variation, we carried out an additional AMOVA for which we defined “region” as host plant (apple or peach). We removed three samples that were collected from persimmon (South Africa), and 45 samples that were collected from pear. The latter samples were collected from three populations in China, and are therefore not representative of a broad geographic sampling range. The final analysis contained 328 samples, 256 of which were collected from peach, and 72 of which were collected from apple.

A Bayesian model-based analysis was performed to infer population structure using Structure version 2.3[31, 32]. Although Structure assumes HWE and no linkage disequilibrium between loci, the software is robust to deviations from these assumptions, particularly if no spurious populations are observed across multiple runs[33]. Data were analyzed using an admixture ancestry model with correlated allele frequencies, to estimate the posterior probabilities L(K) of K groups and the individual percentages of membership assigned to them according to their molecular multilocus profiles[32]. We examined the probabilities for a range of K (K = 1–10), using a burn-in period and a run length of the Markov Chain Monte Carlo (MCMC) of 30 000 and 100 000 iterations, respectively. Based on a trial run, we found that incorporating higher values of K, and longer burn-ins or MCMC did not appreciably change the results. Five runs were carried out for each value of K. We calculated ΔK according to[34]. Figure 2 was drawn using the program Distruct[35].

Figure 2

Bayesian clustering according to the software program S tructure ; results for K =4. Each individual is represented by a line, which is partitioned into K colored segments according to the individuals’ estimated membership fractions in each of the K clusters. Graphs show mean L(K) (±SD) over 5 runs for each value of K between 1 and 10, and ΔK calculated according to [34].

We also carried out an independent analysis of spatial structure using the R package Geneland 3.1.4[25, 26], which has been explicitly tested for robustness to the presence of null alleles. Like Structure, the software uses a MCMC strategy[36] to determine the most likely number of populations (K), and assigns individuals to the most appropriate population basing on individual multilocus genotypes. We carried out ten independent MCMC simulations (100 000 iterations, thinning of 100 iterations) using the non-spatial approach, during which we allowed K to vary between 1 and 20. Afterwards, the model with the highest mean logarithm of posterior probability was post-processed. Geneland was used to estimate inbreeding coefficients (FIS) of inferred clusters.

We also attempted to estimate the magnitude and direction of historical gene flow between our 26 sample populations using a Bayesian coalescent approach implemented in the software program Migrate-N[37]. However the lower 0.025 posterior distribution value of the pairwise migration estimates (M) never differed from 0. In other words, migration rates between populations never differed significantly from zero, and the direction of gene flow between sampling locations could not be inferred using this approach (results not shown).

Genetic diversity

We calculated expected heterozygosity (HE), observed heterozygosity (HO), and allelic richness (NA) using Gen Al Ex. Gen Al Ex was used to identify private alleles from each continent, and from Japan and the Azores.


All of the nine newly developed markers (Table 2) proved to be polymorphic and informative. Among the 13 microsatellite loci used for genotyping, the number of alleles per locus ranged from 7 to 28 (Table 2). There were significant deviations from HWE at multiple loci from all sampling locations (data not shown), and in all cases these were the result of heterozygote deficiencies (Table 3). Null allele frequencies varied from 0.05 to 0.33 among the 13 microsatellite loci (Table 2), which is typical for lepidopteran DNA loci[38], and lower than or comparable to other studies that have used microsatellite markers to determine the genetic structure of populations[39, 40]. Significant linkage disequilibrium was present in 15 out of 78 pairs of loci overall, but was unlikely due to physical linkage; patterns of significant allelic associations were not consistently restricted to certain pairs of loci within the 26 populations (data not shown).

Table 3 Number of sampled individuals ( N ), inbreeding coefficient (F IS ), expected heterozygosity ( H E ), observed heterozygosity ( H O ), and allelic diversity ( N A ; corrected for sample size) of Grapholita molesta individuals from 26 sampling sites

Population genetic structure

Global FST across the data set was 0.219. AMOVA results showed that variation within populations, between populations, and among regions accounted for 72%, 18%, and 10% of the total variation, respectively (P < 0.001). Incorporating host species (peach or apple) rather than geographic region as the “regional” factor revealed that host species did not account for any of the variation in the data set, while variation within and between populations accounted for 70% and 30% of the variation, respectively (P < 0.001). A significant pattern of isolation-by-distance was detected (P < 0.001), although this pattern explained a low proportion of the variation in the data set (R2 = 0.111).

Bayesian analysis of population genetic structure conducted with the program Structure yielded a modal value of ΔK at K = 4 (Figure 2), and different runs at K = 4 produced consistent clustering solutions (Figure 2). Populations from the Azores, South Africa, Argentina, and Chile clustered with populations from North America. Australian and Japanese populations clustered with populations from mainland Asia. Most European populations clustered together, except for single populations from France (Lafitte-sur-Lot) and Spain, which clustered with populations from Brazil. There was little evidence for substantial admixture between populations that clustered separately (Figure 2).

Based on analysis in Geneland, the most likely number of inferred populations was K = 14 consistently across all ten independent runs (Table 3). Although the results from Geneland indicated a greater degree of population structure compared to the Structure results, there were no qualitative disagreements between the two analyses. In other words, populations that clustered together based on the Geneland analysis (Table 3) always clustered together in the Structure analysis (Figure 2). Pairwise FST values between clusters derived from the Geneland analysis are reported in Additional file2: Appendix 2.

In order to further investigate the relationship between the two Brazilian and the two European populations that clustered together based on the Structure analysis, we used the private allele function in Gen Al Ex to check whether the sets of alleles from the two Brazilian populations were a subset of the alleles found in France (Lafitte-sur-Lot) and Spain, or vice versa, which would allow us to infer the direction of gene flow between continents. We found a similar number of private alleles among the Brazilian versus the French/Spanish populations; compared to Brazilian populations, there were 11 alleles that were specific to France (Lafitte-sur-Lot)/Spain, with allele frequencies that varied between 0.033, and 0.310. Similarly, there were 13 alleles that were specific to Brazil compared to France (Lafitte-sur-Lot)/Spain, with allele frequencies that varied between 0.017 and 0.433. The direction of gene flow between Brazil and Europe could therefore not be unambiguously traced using molecular genetic data.

Genetic diversity and inbreeding

Measures of genetic diversity were high in native as well as in many introduced populations (Table 3). Measures of genetic diversity were substantially lower in one population from Brazil (Geneland Cluster 6), one population from Spain (Geneland Cluster 4), and the South African population (Geneland Cluster 13).

The number of private alleles was more than six times greater among Asian populations compared to European and North American populations (Table 4), even though fewer individuals and populations were sampled from Asia compared to either of the latter two continents. FIS estimates corrected for null alleles (computed using Geneland), were high across the range of sampled populations, and did not differ appreciably among populations from the introduced versus invaded ranges (Table 3).

Table 4 Private alleles from each geographic region


Grapholita molesta populations are structured on a global scale, and geographic structure generally reflects continental divisions; these data, combined with knowledge gleaned from historical records, provide important information regarding dispersal and population genetic structure of this invasive species throughout its range.

Population genetic structure and migration

Populations were highly differentiated, with a global FST value of 0.219. In comparison, population genetic structure of the codling moth (Cydia pomonella), which belongs to the same tribe (Grapholitini) as G. molesta[41] and is also a major pest of fruit crops, varies from absent/low (France[42]; Chile[43]) to moderate (Switzerland[44]) to high (South Africa[45]; Italy[46]) among different geographic areas. However, the present study was carried out at a global scale, while individual studies of C. pomonella were carried out at scales varying from a few to several hundred kilometers, impairing comparisons between the two species. Nonetheless, low levels of population structure among C. pomonella populations from some regions relative to G. molesta could be the combined result of higher dispersal abilities of the former species compared to the latter[19, 4749] and the frequent use of pesticides against C. pomonella in some countries, which can stimulate flight in this species (reviewed by[44]).

In the present study, clustering outcomes from two distinct Bayesian approaches were largely congruent, except that the analysis from Geneland suggested greater population substructure compared to the analysis from Structure. Evidence indicates that G. molesta was introduced to Australia from mainland Asia. The data also support a North American source of G. molesta in Argentina and Chile, The Azores and South Africa, since historical records indicate that G. molesta was present in North America from the beginning of the 20th century, but arrived much later in the latter three regions. However, a previous study[21] suggested that South African populations were unlikely to be derived from Canadian populations, since the genetic distance estimates from these two regions were relatively high.

Unexpectedly, Brazilian populations clustered together with two European populations (Structure), and shared ancestry between Western European and Brazilian populations was also supported by results from Geneland. The two European populations that shared ancestry with populations from Brazil showed little evidence of admixture with other European populations. It is possible that these two populations from Europe were derived from the same source as other European populations, but were isolated from other populations on the same continent, and therefore became differentiated over time. If this is the case, Brazilian populations are likely derived from European origins. However, it is also possible that G. molesta was recently re-introduced to Europe from the South American part of its range. Chile is the primary source of off-season peach and nectarine imports to Europe, and Brazil exports nearly 90% of its apple crop to Europe[50, 51]. In contrast, stone and pome fruit exports from Europe to South America are negligible[50, 51]. This latter interpretation of our data is in line with recent findings[42] of C. pomonella individuals with putative South American origin in France, suggesting that South America may be a common source of introduction for important stone and pome fruit pests to Western Europe. Moreover, other invertebrate pests of distantly related plant host species have recently been introduced from South America to Europe, including the pine-infesting woodwasp Sirex noctilio, which was introduced from Chile to Switzerland[52]. Although most invasive species in Europe are assumed to originate from North America or Asia, there is growing evidence that a number of invasive invertebrates in Europe come from South America (e.g.[42],[52], this manuscript).

Although it has been suggested that G. molesta was introduced to North America from Japan[7, 12, 13, 53], we did not find evidence to support this hypothesis. The single Japanese population incorporated in this study clustered with other Asian populations, and did not show evidence of shared ancestry with North American populations. However, it is possible that limited sampling of populations in Japan prevented us from capturing a likely source population from that region, or that a historical bottleneck followed by high rates of mutation since the North American introduction have generated divergence levels that prevent the detection of shared ancestry. Our data did not provide an alternative hypothesis regarding the source of North American G. molesta populations.

We found an overall (albeit weak) pattern of isolation-by-distance, suggesting that regional introductions occur in a relatively stepwise manner. Also, we found little evidence for substantial admixture between populations sampled from different continents, which, combined with the clustering outcomes, implies that multiple cross-continental introductions are unlikely. We note however, that sampling from Australia and Africa was limited to one population each, and we cannot make definitive continent-wide conclusions for these regions based on these sample sizes.

It has earlier been suggested that G. molesta disperses primarily through the movement of fruit, bins and plant material between orchards[18, 21]. The natural dispersal ability of both males and females is limited, even though certain environmental factors can increase flight capacity[54]. Flights between non-contiguous orchards are possible, however they are generally short[19, 47], and human-mediated dispersal is implicated as the main mode of dispersal[20]. The poor natural dispersal ability of G. molesta likely accounts for the high levels of inbreeding observed in our study, since close relatives are probably geographically constrained to the same or neighboring orchards. Combined, these findings suggest that improved management of this invasive pest at regional and community levels (e.g. by area-wide pest management) and at the national level (e.g. by quarantine regulations) may be efficient strategies to limit its ongoing spread.

We did not find any evidence that G. molesta populations are structured according to host plant species. Although G. molesta has been shown to perform better on peach, its primary host, compared to apple[55], serial generations of a population sometimes infest different hosts as the growing season progresses[56]. Such seasonal host shifts may preclude adaptation of populations to any specific host, particularly if several suitable hosts are available within close geographic proximity.

Genetic diversity and inbreeding

Estimates of genetic diversity among different geographic regions were high overall. A six times higher number of private alleles among Asian populations compared to populations from other continents provides, for the first time, convincing evidence that Asia is indeed the native range of this species, as is frequently suggested in the literature. We found that genetic diversity levels are surprisingly high across most of the invaded range, except for a few populations (Brazil, The Azores, and South Africa). It is therefore unlikely that founder effects impair the adaptive potential of G. molesta in most non-native areas, as is suggested to occur in other invasive species. Levels of genetic diversity were reported from the non-native range of 14 species of invasive insects[6], and of these only four species (Ceratitis rosa, Drosophila pseudoobscura, Polistes dominulus, and Rhagoletis completa) exhibited HE values above 0.3. In our study, HE exceeded 0.5 in most populations from both the native and invaded range. High levels of genetic diversity observed in this study could be a result of high microsatellite mutation rates in lepidopterans[57], and/or may result from the introduction of many individuals at each founding event.

High levels of genetic diversity may account for high levels of heritable phenotypic variation within G. molesta populations, and may contribute to the ability of G. molesta to adapt to pest management regimes. For example, there is considerable genetic variation within G. molesta populations with regard to olfactory response of female moths to their host plants[58]. Similarly, intra- and inter-population variation in other traits may enable adaptation to new environments and/or pest control regimes, although heritable variation in such traits has not yet been quantified. Nonetheless, a number of authors have demonstrated that G. molesta has developed insecticide resistance in at least some parts of its invaded range[59, 60].


Unlike many other invasive pests (reviewed by[61]), which have high dispersal abilities and/or are otherwise predisposed to multiple introductions, we found little evidence for frequent introductions between different continents in the case of G. molesta, with the possible exception of a reintroduction to Europe from South America, or two independent introductions to South America. This is likely a result of the poor natural dispersal abilities of G. molesta, combined with stringent international trade regulations that aim to prevent its movement across international borders. Although bottlenecks have occurred in a subset of introduced populations, we argue that, overall, reductions in genetic diversity will not prevent G. molesta from adapting to new control measures or changes in climatic conditions, and we therefore suggest that improved community-based control measures, in addition to more stringent national and intracontinental regulations regarding the regional transport of plant material, will be most effective for containing its ongoing spread.


  1. 1.

    Thrall PH, Oakeshott JG, Fitt G, Southerton S, Burdon JJ, Sheppard A, Russell RJ, Zalucki M, Heino M, Denison RF: Evolution in agriculture: the application of evolutionary approaches to the management of biotic interactions in agro-ecosystems. Evol Appl. 2011, 4: 200-215. 10.1111/j.1752-4571.2010.00179.x.

    PubMed Central  Article  PubMed  Google Scholar 

  2. 2.

    Gornall J, Betts R, Burke E, Clark R, Willett K, Wiltshire A: Implications of climate change for agricultural productivity in the early twenty-first century. Phil Trans Roy Soc B. 2010, 365: 2973-2989. 10.1098/rstb.2010.0158.

    Article  Google Scholar 

  3. 3.

    Godfray HCJ, Crute IR, Haddad L, Lawrence D, Muir JF, Nisbett J, Pretty J, Robinson S, Toulmin C, Whiteley R: The future of the global food system. Phil Trans Roy Soc B. 2010, 365: 2769-2777. 10.1098/rstb.2010.0180.

    Article  Google Scholar 

  4. 4.

    Gilchrist AS, Meats AW: The genetic structure of populations of an invading pest fruit fly, Bactrocera tryoni, at the species climatic range limit. Heredity. 2010, 105: 165-172. 10.1038/hdy.2009.163.

    CAS  Article  PubMed  Google Scholar 

  5. 5.

    Midamegbe A, Vitalis R, Malausa T, Delava E, Cros-Arteil S, Streiff R: Scanning the European corn borer (Ostrinia spp.) genome for adaptive divergence between host-affiliated sibling species. Mol Ecol. 2011, 20: 1414-1430. 10.1111/j.1365-294X.2011.05035.x.

    Article  PubMed  Google Scholar 

  6. 6.

    Dlugosch KM, Parker IM: Founding events in species invasions: genetic variation, adaptive evolution, and the role of multiple introductions. Mol Ecol. 2008, 17: 431-449. 10.1111/j.1365-294X.2007.03538.x.

    CAS  Article  PubMed  Google Scholar 

  7. 7.

    Rothschild GHL, Vickers RA: Biology, ecology and control of the oriental fruit moth. Tortricid pests: their biology, natural enemies and control. Edited by: van der Geest LPS, Evenhuis HH. 1991, Amsterdam: Elsevier, 389-412.

    Google Scholar 

  8. 8.

    Myers CT, Hull LA, Krawczyk G: Effects of orchard host plants (apple and peach) on development of oriental fruit moth (Lepidoptera: Tortricidae). J Econ Entomol. 2007, 100: 421-430. 10.1603/0022-0493(2007)100[421:EOOHPA]2.0.CO;2.

    Article  PubMed  Google Scholar 

  9. 9.

    Piñero JC, Dorn S: Response of female oriental fruit moth to volatiles from apple and peach trees at three phenological stages. Entomol Exp Appl. 2009, 131: 67-74. 10.1111/j.1570-7458.2009.00832.x.

    Article  Google Scholar 

  10. 10.

    Lopez-Vaamonde C, Agassiz D, Augustin S, De Prins J, De Prins W, Gomboc S, Ivinskis P, Karsholt O, Koutroumpas A, Koutroumpa F, Laštůvka Z, Marabuto E, Olivella E, Przybylowicz L, Roques A, Ryrholm N, Sefrova H, Sima P, Sims I, Sinev S, Tomov R, Zilli A, Lees DC: Lepidoptera Chapter 11. BioRisk. 2010, 4: 603-668.

    Article  Google Scholar 

  11. 11.

    Balachowsky AS, Mesnil L: Les insectes nuisibles aux plantes cultivées, leurs moeurs, leur destruction: traité d’entomologie agricole concernant la France, la Corse, l’Afrique du Nord et les régions limitrophes. 1935, Paris: Busson

    Google Scholar 

  12. 12.

    Quaintance AL, Wood WB: Laspeyresia molesta, an important new insect enemy of the peach. J Agr Res. 1916, 7: 373-378.

    Google Scholar 

  13. 13.

    Rings RW: Economic aspects of the biology and control of the oriental fruit moth, Grapholita molesta Busck, in the United States. Ohio J Sci. 1970, 70: 58-61.

    Google Scholar 

  14. 14.

    Bailey P: An attempt to control oriental fruit moth, Cydia molesta Busck by mass releases of Macrocentrus ancylivorus Rohwer (Hymenoptera: Braconidae). J Aust Entomol Soc. 1979, 18: 211-212.

    Article  Google Scholar 

  15. 15.

    Reis W, Nora I, Melzer R: Population dynamics of Grapholita molesta, Busck, 1916, and its adaptation on apple in South Brazil. Acta Hortic. 1988, 232: 204-208.

    Article  Google Scholar 

  16. 16.

    Cox JM, Dale PS: New records of plant pests in New Zealand No. 2. New Zeal J Agr Res. 1977, 20: 109-111. 10.1080/00288233.1977.10427312.

    Article  Google Scholar 

  17. 17.

    Blomefield TL, Barnes BN: Integrated control of oriental fruit moth, Grapholita molesta, on peaches using a spray-date prediction model. Acta Hortic. 2000, 525: 161-168.

    Article  Google Scholar 

  18. 18.

    Torriani MVG, Mazzi D, Hein S, Dorn S: Structured populations of the oriental fruit moth in an agricultural ecosystem. Mol Ecol. 2010, 19: 2651-2660. 10.1111/j.1365-294X.2010.04711.x.

    CAS  Article  PubMed  Google Scholar 

  19. 19.

    Hughes J, Dorn S: Sexual differences in the flight performance of the oriental fruit moth, Cydia molesta. Entomol Exp Appl. 2002, 103: 171-182. 10.1046/j.1570-7458.2002.00967.x.

    Article  Google Scholar 

  20. 20.

    Mazzi D, Dorn S: Movement of insect pests in agricultural landscapes. Ann Appl Biol. 2012, 160: 97-113. 10.1111/j.1744-7348.2012.00533.x.

    Article  Google Scholar 

  21. 21.

    Timm AE, Geertsema H, Warnich L: Population genetic structure of Grapholita molesta (Lepidoptera: Tortricidae) in South Africa. Ann Entomol Soc Am. 2008, 101: 197-203. 10.1603/0013-8746(2008)101[197:PGSOGM]2.0.CO;2.

    Article  Google Scholar 

  22. 22.

    Sciarretta A, Trematerra P: Geostatistical characterization of the spatial distribution of Grapholita molesta and Anarsia lineatella males in an agricultural landscape. J Appl Entomol. 2006, 130: 73-83. 10.1111/j.1439-0418.2006.01034.x.

    Article  Google Scholar 

  23. 23.

    Abdelkrim J, Robertson BC, Stanton J-AL, Gemmell NJ: Fast, cost-effective development of species-specific microsatellite markers by genomic sequencing. Biotechniques. 2009, 46: 185-191. 10.2144/000113084.

    CAS  Article  PubMed  Google Scholar 

  24. 24.

    Peakall R, Smouse PE: GenAlEx 6: genetic analysis in Excel. Population genetic software for teaching and research. Mol Ecol Notes. 2006, 6: 288-295. 10.1111/j.1471-8286.2005.01155.x.

    Article  Google Scholar 

  25. 25.

    Guillot G, Estoup A, Mortier F, Cosson JF: A spatial statistical model for landscape genetics. Genetics. 2005, 170: 1261-1280. 10.1534/genetics.104.033803.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  26. 26.

    Guillot G, Mortier F, Estoup A: Geneland: a computer package for landscape genetics. Mol Ecol Notes. 2005, 5: 712-715. 10.1111/j.1471-8286.2005.01031.x.

    CAS  Article  Google Scholar 

  27. 27.

    Excoffier L, Laval G, Schneider S: Arlequin (version 3.0): an integrated software package for population genetics data analysis. Evol Bioinf Online. 2005, 1: 47-50.

    CAS  Google Scholar 

  28. 28.

    Slatkin M: Linkage disequilibrium in growing and stable populations. Genetics. 1994, 137: 331-336.

    PubMed Central  CAS  PubMed  Google Scholar 

  29. 29.

    Huff DR, Peakall R, Smouse PE: RAPD Variation within and among natural populations of outcrossing buffalo grass (Buchloe dactyloides (Nutt.) Engelm.). Theor Appl Genet. 1993, 86: 927-934.

    CAS  Article  PubMed  Google Scholar 

  30. 30.

    Weir BS, Cockerham CC: Estimating F-statistics for the analysis of population structure. Evolution. 1984, 38: 1358-1370. 10.2307/2408641.

    Article  Google Scholar 

  31. 31.

    Pritchard JK, Stephens M, Donnelly P: Inference of population structure using multilocus genotype data. Genetics. 2000, 155: 945-959.

    PubMed Central  CAS  PubMed  Google Scholar 

  32. 32.

    Falush D, Stephens M, Pritchard JK: Inference of population structure using multilocus genotype data: dominant markers and null alleles. Mol Ecol Notes. 2007, 7: 574-578. 10.1111/j.1471-8286.2007.01758.x.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  33. 33.

    Falush D, Stephens M, Pritchard JK: Inference of population structure using multilocus genotype data: linked loci and correlated allele frequencies. Genetics. 2003, 164: 1567-1587.

    PubMed Central  CAS  PubMed  Google Scholar 

  34. 34.

    Evanno G, Regnaut S, Goudet J: Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol Ecol. 2005, 14: 2611-2620. 10.1111/j.1365-294X.2005.02553.x.

    CAS  Article  PubMed  Google Scholar 

  35. 35.

    Rosenberg NA: DISTRUCT: a program for the graphical display of population structure. Mol Ecol Notes. 2004, 4: 137-138.

    Article  Google Scholar 

  36. 36.

    Guillot G, Santos F, Estoup A: Analysing georeferenced population genetics data with Geneland: a new algorithm to deal with null alleles and a friendly graphical user interface. Bioinformatics. 2008, 24: 1406-1407. 10.1093/bioinformatics/btn136.

    CAS  Article  PubMed  Google Scholar 

  37. 37.

    Beerli P: Comparison of Bayesian and maximum likelihood inference of population genetic parameters. Bioinformatics. 2006, 22: 341-345. 10.1093/bioinformatics/bti803.

    CAS  Article  PubMed  Google Scholar 

  38. 38.

    Zhang DX: Lepidopteran microsatellite DNA: redundant but promising. Trends Ecol Evol. 2004, 19: 507-509. 10.1016/j.tree.2004.07.020.

    Article  PubMed  Google Scholar 

  39. 39.

    Gilabert A, Simon JC, Mieuzet L, Halkett F, Stoekel S, Plantegenest M, Dedryver C-A: Climate and agricultural context shape reproductive mode variation in an aphid crop pest. Mol Ecol. 2009, 18: 3050-3061. 10.1111/j.1365-294X.2009.04250.x.

    CAS  Article  PubMed  Google Scholar 

  40. 40.

    Chapuis MP, Popple JAM, Berthier K, Simpson SJ, Deveson E, Spurgin P, Steinbauer MJ, Sword GA: Challenges to assessing connectivity between massive populations of the Australian plague locust. Proc Roy Soc B. 2011, 278: 3152-3160. 10.1098/rspb.2010.2605.

    Article  Google Scholar 

  41. 41.

    Chen MH, Dorn S: Cross-amplification of microsatellites from the codling moth Cydia pomonella to three other species of the tribe Grapholitini (Lepidoptera: Tortricidae). Mol Ecol Res. 2010, 10: 1034-1037. 10.1111/j.1755-0998.2010.02837.x.

    CAS  Article  Google Scholar 

  42. 42.

    Franck P, Reyes M, Olivares J, Sauphanor B: Genetic architecture in codling moth populations: comparison between microsatellite and insecticide resistance markers. Mol Ecol. 2007, 16: 3554-3564. 10.1111/j.1365-294X.2007.03410.x.

    CAS  Article  PubMed  Google Scholar 

  43. 43.

    Fuentes-Contreras E, Espinoza JL, Lavandero B, Ramírez CC: Population genetic structure of codling moth (Lepidoptera: Tortricidae) from apple orchards in central Chile. J Econ Entomol. 2008, 101: 190-198. 10.1603/0022-0493(2008)101[190:PGSOCM]2.0.CO;2.

    Article  PubMed  Google Scholar 

  44. 44.

    Chen MH, Dorn S: Microsatellites reveal genetic differentiation among populations in an insect species with high genetic variability in dispersal, the codling moth, Cydia pomonella (L.) (Lepidoptera: Tortricidae). B Entomol Res. 2010, 100: 75-85. 10.1017/S0007485309006786.

    CAS  Article  Google Scholar 

  45. 45.

    Timm AE, Geertsema H, Warnich L: Gene flow among Cydia pomonella (Lepidoptera: Tortricidae) geographic and host populations in South Africa. J Econ Entomol. 2006, 99: 341-348. 10.1603/0022-0493-99.2.341.

    CAS  Article  PubMed  Google Scholar 

  46. 46.

    Thaler R, Brandstaetter A, Meraner A, Chabicovski M, Parson W, Zelger R, Dalla Via J, Dallinger R: Molecular phylogeny and population structure of the codling moth (Cydia pomonella) in Central Europe: II. AFLP analysis reflects human-aided local adaptation of a global pest species. Mol Phylogenet Evol. 2008, 48: 838-849. 10.1016/j.ympev.2008.05.027.

    CAS  Article  PubMed  Google Scholar 

  47. 47.

    Yetter WP, Steiner LF: Efficiency of bait traps for the oriental fruit moth as indicated by the release and capture of marked adults. J Econ Entomol. 1932, 25: 106-116.

    CAS  Article  Google Scholar 

  48. 48.

    Mani E, Wildbolz T: Dispersal of male codling moths (Laspeyresia pomonella L) in upper Rhine valley. J Appl Entomol. 1977, 83: 161-168.

    Google Scholar 

  49. 49.

    Schumacher P, Weyeneth A, Weber DC, Dorn S: Long flights in Cydia pomonella L. (Lepidoptera: Tortricidae) measured by a flight mill: influence of sex, mated status and age. Physiol Entomol. 1997, 22: 149-160. 10.1111/j.1365-3032.1997.tb01152.x.

    Article  Google Scholar 

  50. 50.

    USDA Foreign Agricultural Service: Global Agricultural Information Network (Gain) Report: Brazil, Annual Fresh Deciduous Fruit. 2011

    Google Scholar 

  51. 51.

    USDA Foreign Agricultural Service: Global Agricultural Information Network (Gain) Report: EU-27, Stone Fruit Annual. 2011

    Google Scholar 

  52. 52.

    Boissin E, Hurley B, Wingfield MJ, Vasaitis R, Stenlid J, Davis C, De Groot P: Retracing the routes of introduction of invasive species: the case of the Sirex noctilio woodwasp. Mol Ecol. 2012, 21: 5728-5744. 10.1111/mec.12065.

    CAS  Article  PubMed  Google Scholar 

  53. 53.

    Bellerose S, Chouinard G, Roy M: Occurrence of Grapholita molesta (Lepidoptera: Tortricidae) in major apple-growing areas of southern Quebec. Can Entomol. 2011, 139: 292-295.

    Article  Google Scholar 

  54. 54.

    Hughes J, Hern A, Dorn S: Preimaginal environment influences adult flight in Cydia molesta (Lepidoptera: Tortricidae). Environ Entomol. 2004, 33: 1155-1162. 10.1603/0046-225X-33.5.1155.

    Article  Google Scholar 

  55. 55.

    Silva OABNE, Botton M, Garcia MS, Bisognin AZ, Nava DE: Development and reproduction of the oriental fruit moth in apple and peaches. Pesquisa Agropecuária Brasileira. 2010, 45: 1082-1088.

    Google Scholar 

  56. 56.

    Zhao ZR, Wang YG, Yan GY: A preliminary report on the oriental fruit moth in north Jiangsu. Insect Knowledge. 1989, 26: 17-19.

    Google Scholar 

  57. 57.

    Tay WT, Behere GT, Batterham P, Heckel DG: Generation of microsatellite repeat families by RTE retrotransposons in lepidopteran genomes. BMC Evol Biol. 2010, 10: 144-10.1186/1471-2148-10-144.

    PubMed Central  Article  PubMed  Google Scholar 

  58. 58.

    Najar-Rodriguez A, Schneeberger M, Bellutti N, Dorn S: Variation in attraction to host plant odors in an invasive moth has a genetic basis and is genetically negatively correlated with fecundity. Behav Genet. 2012, 42: 687-697. 10.1007/s10519-012-9539-9.

    CAS  Article  PubMed  Google Scholar 

  59. 59.

    Kanga LHB, Pree DJ, van Lier JL, Walker GM: Management of insecticide resistance in oriental fruit moth (Grapholita molesta; Lepidoptera: Tortricidae) populations from Ontario. Pest Manag Sci. 2003, 59: 921-927. 10.1002/ps.702.

    CAS  Article  PubMed  Google Scholar 

  60. 60.

    Siegwart M, Monteiro LB, Maugin S, Olivares J, Carvalho SM, Sauphanor B: Tools for resistance monitoring in oriental fruit moth (Lepidoptera: Tortricidae) and first assessment in Brazilian populations. J Econ Entomol. 2011, 104: 636-645. 10.1603/EC10302.

    CAS  Article  PubMed  Google Scholar 

  61. 61.

    Handley LJL, Estoup A, Evans DM: Ecological genetics of invasive alien species. Biocontrol. 2011, 56: 409-428. 10.1007/s10526-011-9386-2.

    Article  Google Scholar 

Download references


We thank our collaborators, listed in Table 1, who kindly collected and sent G. molesta samples for analysis, and Claudio Sedivy and two reviewers for providing useful comments on the manuscript. Data analyzed in this paper were generated in the Genetic Diversity Centre of ETH Zurich.

Author information



Corresponding author

Correspondence to Dominique Mazzi.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

SD and DM conceived of and designed the study. HK carried out the molecular genetic work, analyzed the data, and drafted the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material


Additional file 1:Appendix 1. Locus specific PCR conditions. Description: Detailed description of PCR conditions. (DOCX 80 KB)


Additional file 2:Appendix 2. Pairwise FST values. Description: Pairwise FST values from 14 clusters derived from analysis in GENELAND. (DOCX 19 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and Permissions

About this article

Cite this article

Kirk, H., Dorn, S. & Mazzi, D. Worldwide population genetic structure of the oriental fruit moth (Grapholita molesta), a globally invasive pest. BMC Ecol 13, 12 (2013).

Download citation


  • Microsatellites
  • Human-mediated dispersal
  • Isolation-by-distance
  • Invasive pest
  • Genetic structure
  • Population differentiation
  • Oriental fruit moth