Skip to main content
  • Research article
  • Open access
  • Published:

Patterns of genetic variation among geographic and host-plant associated populations of the peach fruit moth Carposina sasakii (Lepidoptera: Carposinidae)



Populations of herbivorous insects may become genetically differentiated because of local adaptation to different hosts and climates as well as historical processes, and further genetic divergence may occur following the development of reproductive isolation among populations. Here we investigate the population genetic structure of the orchard pest peach fruit moth (PFM) Carposina sasakii (Lepidoptera: Carposinidae) in China, which shows distinct biological differences when characterized from different host plants. Genetic diversity and genetic structure were assessed among populations from seven plant hosts and nine regions using 19 microsatellite loci and a mitochondrial sequence.


Strong genetic differentiation was found among geographical populations representing distinct geographical regions, but not in host-associated populations collected from the same area. Mantel tests based on microsatellite loci indicated an association between genetic differentiation and geographical distance, and to a lesser extent environmental differentiation. Approximate Bayesian Computation analyses supported the scenario that PFM likely originated from a southern area and dispersed northwards before the last glacial maximum during the Quaternary.


Our analyses suggested a strong impact of geographical barriers and historical events rather than host plants on the genetic structure of the PFM; however, uncharacterized environmental factors and host plants may also play a role. Studies on adaptive shifts in this moth should take into account geographical and historical factors.


Genetic variation among natural populations can develop due to a number of factors that include geographical isolation, ecological isolation and historical processes. Geographical barriers that limit dispersal and consequently lead to isolation by distance (IBD) appear to be particularly important in population divergence in diverse taxa [1]. However, extant patterns of genetic differentiation may also be impacted by historical processes, such as those associated with climate oscillations of the Quaternary, when many species became restricted to refugia in glacial periods, interspersed by range expansions in interglacial periods [2]. And in the last decade, an increasing number of studies have shown that ecological factors also play an important role in shaping genetic differentiation (isolation by environment, IBE) [3,4,5,6,7].

The relative contribution of these factors on population differentiation can be difficult to determine [3]. In empirical and simulation studies, false positives or underestimated correlations between genetic and environmental variations can be generated through the influence of IBD and spatial autocorrelation of ecological variables [5]. In Mantel tests of IBD, hierarchical population structure, which is mostly caused by postglacial recolonisation from multiple refugia, can be confounded with IBE [8]. A better understanding of the complex factors influencing population differentiation needs well designed sampling srategies, and a combined consideration of geography, history and ecology [9].

Herbivorous insects represent a diverse group of species with a wide range of distributions and adaptive potential [10,11,12,13]. Population genetic differentiation of these insects may be influenced by geographical, historical and ecological factors [9, 14]. Host plants represent one obvious form of ecological variation that can play a crucial role in the diversification of herbivorous insect populations [15]. Alternative host-plant species can generate different selection pressures that create ecological barriers to gene flow between insect populations [16,17,18]. Because hosts often differ in traits that are linked ecologically and physiologically to performance (e.g. nutritional quality, recognition cues), fitness trade-offs and divergent selection between plants can occur and contribute to ecological isolation and speciation [19,20,21]. An increasing number of cases of host-associated differentiation have been documented in insects [11, 22,23,24,25,26,27].

The peach fruit moth (PFM), Carposina sasakii Matsumura (Lepidoptera: Carposinidae), is a major phytophagous orchard pest widely distributed in Northeast Asia [28]. Larvae of PFM bore into the fruits of multiple hosts in the Rosaceae and Rhamnaceae, mainly apple, pear, peach, apricot, hawthorn, Chinese quince, jujube, and wild jujube. On different host plants, PFMs vary in performance both under field and laboratory conditions. In the field, peak adult emergence time, oviposition habitat and generation number can vary (Table S1), likely to synchronize developmental stages with hosts. Under laboratory conditions, adult PFMs live significantly longer on jujube than on other hosts, adult females reared from jujube and peach tend to lay more eggs [29], and larval survival also varies with host plant [30]. Phenological isolation associated with host usage may facilitate host-associated adaptation and reduce flow among host-associated populations. Based on biological observations, esterase isozyme patterns [31] and random amplified polymorphic DNA (RAPD) [32], PFMs on different hosts have been proposed as representing host biotypes. Based on mtDNA, two sympatric and cryptic lineages of the PFM were identified in populations from China; however, no association between population variation and host plants was found [33].

Geographical isolation may also contribute to genetic differentiation of PFM. Based on variation in the mtDNA cox1 gene, there is a correlation between genetic differentiation and geographical distance [33]. Using 35 microsatellite loci, genetic differentiation was detected between two geographically distant populations collected from two host plants, Chinese quince and apple [34]. Both host usage and geographic isolation might therefore contribute to genetic differentiation in PFM.

In this study, we simultaneously characterized genetic variation of PFM from both host-associated and geographical populations across China, using microsatellite markers and mtDNA. We hypothesized that populations from different host plants would differ genetically when the influence of geographical influence was removed, and also that populations from different geographical locations would show IBD given the wide distribution range of this species. We therefore estimated the degree of genetic differentiation of PFM associated with different hosts versus geographic distance and also considered historical factors. Our study sheds light on understanding ecological and evolutionary processes that drive divergence of PFM and the possibility of host-associated reproductive isolation in this species.


Specimen collection and DNA extraction

In total 410 PFM larvae were sampled from damaged fruits of host plants in 16 populations with permissions from the orchard owners (Table 1 and Fig. 1). The 10 host-associated populations were collected from seven hosts of apple, pear, hawthorn, apricot, crabapple, Chinese quince (Rosaceae) and jujube (Rhamnaceae). The nine geographical populations cover most of the distribution of PFM in China. To separate geographical distance from host plant effects, nine host-associated populations were collected from Beijing in northern China, as well as a population from Chinese quince collected from Hubei province in southern China, 865 km from Beijing. We included two populations from jujube, apple and apricot in Beijing to evaluate genetic differentiation between populations from the same host plant. The distance among host-associated populations in Beijing ranged from adjacent orchards (BJYQ02X and BJYQ02P) to a distance of 151 km. Samples were obtained from multiple trees at each location, stored in absolute ethanol and frozen at −80 °C prior to DNA extraction. Genomic DNA was extracted from a segment of individual larva using DNeasy Blood & Tissue Kit (QIAGEN, Hilden, Germany).

Table 1 Sample collection information for the Carposina sasakii used in this study
Fig. 1
figure 1

Collection sites of Carposina sasakii and BAPS analysis of geographical and host-associated populations based on microsatellite loci and mtDNA. The different colors in each population correspond to the frequency of cluster membership based on the BAPS analysis. Figs a and c show separation in the geographical populations, where one cluster and two separate populations were identified based on microsatellite loci (Fig. a), while based on mtDNA four clusters were identified (Fig. c). Figs b and d show the host-associated populations, which provided no strong evidence of genetic structure based either on microsatellites (Fig. b) or mtDNA (Fig. d)

Microsatellite genotyping and mtDNA sequencing

For the nuclear markers, we genotyped 19 polymorphic microsatellite loci from each individual, developed with methods used in our previous study [34] (Table S2). This involved using PC tail (Primer tail C) modified forward primers and fluorescence-labeled PC tails (FAM, HEX, and ROX) for amplification [35]. For the mitochondrial marker, a fragment of mitochondrial cox1 gene (507 bp) was amplified using primer pair LCO1490 and HCO2198 [36]. Polymerase chain reaction (PCR) was conducted using the Mastercycler pro system (Eppendorf, Germany) with standard PCR conditions and an annealing temperature of 52 °C. Amplified products were purified and sequenced directly from both strands using an ABI 3730xl DNA Analyzer (Applied Biosystems, USA).

Genetic diversity analyses

Prior to population genetic analysis, microsatellites genotyped by GENEMAPPER version 4.0 (Applied Biosystems, USA) were checked for stuttering, scoring error, large allele dropout and presence of null alleles by MICRO-CHECKER [37]. Allele frequencies, number of alleles, observed (H O ) and expected (H E ) heterozygosity, were estimated by macros in Microsatellite Tools [38]. Null allele frequency was estimated using FREENA [39] with 10,000 bootstraps. In addition, deviations from Hardy-Weinberg Equilibrium (HWE) and tests for linkage disequilibrium (LD) were calculated with an exact probability test [40] implemented in GENEPOP version 4.0 [41].

Sequencing results of mtDNA from both strands were assembled. Amino acid sequences were aligned by codons using CLUSTALW [42] implemented in MEGA version 6 [43] under default parameters. Nucleotide sequence alignment was guided by aligned amino acid sequences. The number of polymorphic sites (S), total number of mutations (η), number of haplotypes (H), haplotype diversity (Hd), nucleotide diversity (Pi), nucleotide diversity with Jukes and Cantor correction Pi (JC), Tajima’s D and average number of nucleotide differences (K) were calculated with DnaSP version 5.0 [44].

Population structure analysis

For microsatellites loci, genetic differentiation among 14 populations of PFM was measured by pairwise F ST calculated in FREENA version 4.0 with ENA [41]. For mtDNA, ARLEQUIN suit version 3.5 was used to conduct an exact test of population differentiation based on default parameters [45].

In order to incorporate spatial information into clustering of individuals, the BAPS (Bayesian analysis of population structure) model implemented in software BAPS version 6.0 [46] was used based on microsatellite loci or mtDNA. For microsatellite data, the number of populations (K) ranged from 1 to 20 with 20 iterations per K value, while for mtDNA, 20 runs (K = 20, 15 and 10) were performed to ensure convergence and consistency of the results.

We performed a Discriminant Analysis of Principal Components (DAPC) analysis using adegenet 1.4–2 implemented in R [47] based on microsatellite loci, which plots individuals in space based on genetic similarity without biological assumption.

Isolation by distance and environment

In order to evaluate the effect of geographic distance on genetic differentiation of host-associated populations in Beijing region and assess the level of isolation by distance (IBD) within geographical populations, a Mantel test correlating genetic distance (F ST /(1-F ST )) and geographic distance was undertaken using ade4 version 1.7–4 implemented in R (Daniel et al. 2004) with 999 replicates. The values of F ST were calculated in FREENA version 4.0 with ENA [41] for microsatellite data and ARLEQUIN suit version 3.5 for mtDNA (Excoffier & Lischer 2010).

To check the influence of environmental factors on population genetic differentiation, the presence of isolation by environment (IBE) was tested. Firstly, 19 bioclimatic variables were downloaded from WorldClim database ( using the getData function implemented in R package RASTER. Subsequently, we extracted corresponding bioclimatic values of each location using the getData function. Three vegetation variables (NDVI: normalized difference vegetation index, LAI: leaf area index, and percent tree cover) were downloaded from MODIS landcover database ( and then extracted with ArcGIS multiple version 10.2 (ESRI Inc., Redlands, CA). Then we extracted bioclimatic values and vegetation values of each location using cbind function in R package and ARCGISmultiple. Finally, a principal component analysis was performed to analyze the 22 environmental variables for each locality using prcomp function in R. The first two principal components were used to estimate environmental distances between locations. Environmental distances were compared with genetic distances (F ST /(1-F ST )) based on microsatellite or mtDNA by Mantel tests in R package ade4 version 1.7–4 with 999 replicates.

To investigate the extent of eco-spatial autocorrelation in our data, we performed a Mantel test between the ecological and geographical distance matrices. To further assess the relative contribution of environmental variables and geographical distance, matrix regression with a randomization (MMRR) method implemented in R with 10,000 permutations was used [48].

Haplotype relationship analysis and molecular dating

Haplotype relationships were constructed through the software SPLITSTREE version 4.13.1 [49], while the divergence times for haplotype lineages were estimated using the software BEAST version 1.8.1 [50], as described in [51]. In molecular dating analysis, Carposina fernadana and Carposina hyperlopha were used as outgroups.

Test on scenarios of PFM dispersal

The approximate Bayesian computation (ABC) method implemented in DIYABC version 2.1.0 [52] was followed to compare different dispersal scenarios and infer the ancestral populations in PFM based on microsatellite loci (Fig. 2). Datasets were generated by selecting different populations representing the identified groups of PFM, in order to avoid misleading results and false signals of bottlenecks caused by pooling different samples to identify a group, and simplifying complexity of scenarios to be compared [53, 54]. In total, two datasets were provided in the analysis. Moreover, we assumed two unknown populations as ghost populations divided into two branches. In total, six biologically plausible dispersal scenarios representing the relationships of the three groups were conducted and compared, considering the variation of population size and the split and admixture events. The six scenarios could be split into two categories with or without admixture events. Details of pre-evaluation scenario-prior combinations, estimation of posterior distributions of parameters, model checking, and evaluations of confidence in scenario choice are described in supporting information (Additional file 1: Appendix S1).

Fig. 2
figure 2

Graphical representation of the six scenarios for the three population groups. NE, northeast populations; NO, north populations; SO, south populations. A1 and A2 are two unknown (ghost) populations divided into two branches. Scenario 1, 2 and 3 correspond to possible evolutionary relationships among the three populations without admixture between any two of them. Scenario 4, 5 and 6 assume that one of the three populations is an admixture of two other populations


Genetic diversity and pairwise population differentiation

All microsatellite loci used in the study proved to be polymorphic. The mean number of alleles for each population was high, and the H O was similar to H E in each population. The host-associated and geographically separated populations showed similar values for genetic diversity parameters (Table S3). There was no obvious LD among the 19 microsatellite loci; no loci were significantly linked or departed from HWE across all populations, and no population departed from HWE across all loci.

In total 35 haplotypes were observed (GenBank accession numbers: KY492475-KY492509), among which 14 haplotypes were uniquely represented by one individual, 10 shared among individuals but not across populations, and 11 shared among different populations. Tajima’s D was not significantly different from 0 in all populations (Table S4) after Holm’s correction [55].

Null alleles did not generate bias in estimates of population differentiation (Additional file 1: Appendix S2). Fewer populations were genetically differentiated between pairs of host-associated populations collected in Beijing than between pairs of geographically isolated populations (Table 2). F ST values between the geographically separated Chinese quince population and other populations were mostly higher than those between the other population (Table 2).

Table 2 Pairwise F ST values of comparisons among 14 Carposina sasakii populations based on microsatellite loci (lower triangle) and mtDNA (upper triangle)

Population genetic structure

For the nine geographical populations, BAPS analysis based on microsatellite loci revealed that seven populations clustered into one large group, while one northern population and one southern population were separated from this cluster with minor admixture (Fig. 1a). The analysis based on mtDNA identified four groups, which did not entirely coincide with the microsatellite groups (Fig. 1c). Most individuals in the nine Beijing populations collected from different hosts fell into one major cluster for both types of markers (Fig. 1b and d).

DAPC analyses indicated genetic differentiation between the southern population collected on Chinese quince and other populations within the nine geographical populations (Fig. 3a). No differentiation was found among nine host-associated populations collected from two areas of Beijing (Fig. 3b) or six host-associated populations collected from the Yanqing area of Beijing (Fig. 3c).

Fig. 3
figure 3

Discriminant Analysis of Principal Components (DAPC) in populations of Carposina sasakii. a Nine geographic populations from different regions in China (b) Nine host-associated populations collected from Beijing. c Six host-associated populations collected from Yanqing of Beijing

Isolation by geographical and environmental distances

Both microsatellite and mitochondrial data showed a lack of any association between host-associated population differentiation and distance in the Beijing region (r = −0.056, P = 0.677 for microsatellite loci, r = −0.014, P = 0.502 for mtDNA). For the microsatellite data, Mantel tests indicated the presence of both IBD and IBE when considering the geographically separated populations of PFM. However, a significant correlation between ecological and geographical distance was found (r = 0.767, P = 0.002). The standardized regression coefficient for geographic distance onto genetic distance based on all populations (β D  = 0.476, P = 0.0080) was similar to the equivalent regression coefficient for environmental distance (β D  = 0.443, P = 0.0115), suggesting that IBD was stronger than IBE. For the mtDNA, there was no evidence of either IBD or IBE.

Haplotype network, divergence time and demographic history

SPLITSTREE analysis divided the mitochondrial haplotypes into four major lineages (Fig. 4), mostly corresponding to the four geographical groups identified in the population genetic structure analyses but with some admixture. The southern and western lineages were more closely related to each other than to the other lineages. One haplotype from an eastern population (BJPGL) fell into the western lineage.

Fig. 4
figure 4

The SPLITSTREE network from 14 Carposina sasakii collections based on mtDNA. Four major lineages were found. The largest one included haplotypes mainly from northeastern populations (blue). The second lineage was composed of haplotypes from southern China (green). The remaining two lineages are mainly composed of haplotypes from the northeast population (grey) and the western population (pink), with minor contributions from the other populations. Points in the same color (except for grey) indicate haplotypes from the same population. Points in grey indicate haplotypes shared by populations. The points labeled by hap_1, hap_5, hap_6, hap_7 and hap_16 were haplotypes shared by individuals from northern populations

Molecular clock analysis of the mtDNA indicated that an ancient haplotype diverged from others 1.02 Ma (million years ago) with a 95% highest posterior density (HPD) of 0.43–1.83 (Additional file 1: Figure S1). Two major lineages, corresponding to the southern and western lineages versus northern and northeastern lineages, diverged 0.70 Ma (95% HPD = 0.35–1.15), while the divergence times within the two lineages were 0.52 (95% HPD = 0.24–0.87) and 0.39 (95% HPD = 0.15–0.70) Ma, respectively (Additional file 1: Figure S1).

Dispersal routes

The ABC analyses supported scenario 2 (posterior probability of 0.3768 on average) as the most likely based on microsatellite data (Fig. 2). In this scenario, the southern population and A2 (a ghost population) are from A1 (another ghost population), and A2 is a population established later and linked to the northeastern and northern populations. The choice of scenarios was reliable based on an evaluation of confidence and model checking (Additional file 1: Appendix S1).


Geographical but no host-associated differentiation

Our analyses revealed strong genetic differentiation among geographical populations of PFM, but no host-associated differentiation. The genetic clusters correspond to different geographical regions, indicating a strong effect of geographical barriers on population divergence in PFM. Although some clusters were identified by only a single population, the population genetic structure analysis is congruent with the phylogenetic network pattern. Geographical isolation plays an important role in population divergence in nature [1, 4] and geographically structured populations have been documented in two other orchard insects which also cause heavy damage to fruit [14, 56]. Geographically differentiated pest populations might be more likely to occur in pests of orchards compared to those in ephemeral crop and vegetable fields [51, 57, 58] due to the relatively stable ecosystem provided by orchards [59].

There was a clear lack of host-associated differentiation in PFM based on collections from nine populations in the Beijing area. These populations included two collected from apple with a genetic distance higher than some distances obtained from population pairs from different hosts (Table 2). The lowest F ST value (−0.0005) came from the pair of populations collected from adjacent apple and apricot orchards (BJYQ02P and BJYQ02X, Table 2). The Chinese quince population showed strong genetic differentiation between other host-associated populations. However, this population was geographically separate from the other populations. Two other populations from jujube (NXWZZ and SDTAZ) and one from apple (SXJZP) were also geographically separated (Fig. 1c), and the population from Chinese quince grouped with a population from apple based on mtDNA. Each of the two populations from jujube, apple and apricot separated into two different clusters, further highlighting the lack of plant host effect on genetic differentiation. However, we could not exclude the occasions of locally formed host-associated differentiation out of the study area.

Previous studies have investigated host-associated differentiation of PFM [31,32,33, 60]. Based on esterase isozymes of three host populations, Hua and Hua argued for differentiation between populations collected from apple and populations from jujube and wild jujube; however, the distance between collection locations was about 500 km [31]. RAPD marker data also suggested genetic differentiation linked to hosts and particularly apricot in comparison with apple, hawthorn, peach, cornel, jujube, and wild jujube [32]. This pattern contrasts sharply with the lack of differentiation found here, which might reflect the markers used [61] or the nature of the populations tested. Perhaps esterases detect adaptive differences associated with selection. Nevertheless, the very low genetic differentiation between hosts in adjacent orchards suggests high gene flow between hosts, at least in northern China where this sampling took place. Sympatric host-associated populations from other geographical areas are needed to validate the absence of host-associated differentiation in PFM.

A high level of gene flow in PFM populations may prevent host-related differentiation even if there is host-associated selection. PFM adults tend to be highly mobile [62] and therefore move between adjacent orchards. It is also possible that there is some genetic differentiation among host types, but this was not detected because the markers we used are not involved in host plant adaptation. Biological studies have suggested host differences in the induction of diapause and temperature-dependent development in PFM [62] and these may reflect differentiation at loci under selection. Genes related to biological characteristics, such as circadian clock genes, metabolic arrest, adult eclosion, host selection and oviposition behavior may be differentiated [63]. These might not result in genetic differentiation at neutral markers, particularly if adaptive differentiation is very recent [64, 65]. New technologies targeting genome-wide differentiation may be needed for detecting host-associated genetic differentiation of PFM involving adaptive loci [66].

Historical events and population differentiation

Molecular dating revealed an early divergence of the mtDNA about 1 Ma (within Pleistocene, 2.58–0.0117 Ma), pointing to an influence of climatic vacillation during the Quaternary on PFM. This suggests that PFM may be useful for testing hypotheses about the historical effects of the Quaternary on phylogeographical patterns in China, with patterns found so far contrasting with those of well-studied areas of Europe and North America [2, 67]. During the glacial periods of Quaternary, no unified ice sheet had developed in China [68], providing opportunities for divergence and even regional expansion of organisms before the last glacial maximum (LGM, 0.018–0.025 Ma) [14, 69]. Few studies on insects have traced these patterns of divergence [14], unlike plant studies that have shown evidence for multiple refugees in China (east Asia) during the Quaternary [67], mostly located in the southern region [70], but also in the northern region [71].

We explored all possible hypotheses on the origin and dispersal of PFM based on the identified genetic groups using the ABC method. This method is suitable to test complex scenarios in population genetics [72], and has been used in recent work [14, 53, 73, 74]. Our analyses support the notion that the PFM originated from southern China followed by dispersal from south to north. In terms of pest management, PFM was considered as a major pest of deciduous fruit trees in northern China, although damage was occasionally found in southern China. Southern China was warmer than northern regions during the Quaternary, likely allowing species to persist there. A similar pattern of origin and dispersal has been reported in another orchard pest, Grapholita molesta [14].

Molecular dating showed the divergence time of major lineages before 0.39 Ma, indicating colonization of northern China by PFM before LGM. This is congruent with the pattern for G. molesta in China [14].

Isolation by distance and environment

Apart from IBD, habitats can contribute to genetic divergence by creating barriers to gene flow [75], resulting in IBE [3]. Mantel tests based on microsatellite data showed the presence of both IBD and IBE in populations of PFM. High false positive rates for Mantel tests of IBE can arise when high levels of IBD and eco-spatial autocorrelation occurs [76], and Mantel tests showed strong correlations in our data. However, while the MMRR analysis suggested that the effect of geographical isolation on genetic differentiation was stronger than environmental factors, there was also an effect of IBE in PFM populations, in support of patterns in the literature that suggest IBE is common [4, 5]. In PFM, the emergence of adults, development rate and voltinism depend on temperature [77, 78], and the sampled populations cover a wide geographical range along a temperature gradient. With temperature affecting life history traits, relatively higher rates of gene flow might be expected across populations sharing a similar thermal environment. This might be tested further by comparing patterns of gene flow across topographically complex areas where a high degree of local temperature variation might be present.

It is unclear why there was an apparent hierarchical structure of mtDNA variation in PFM, which did not appear connected to IBD or IBE. Incongruent population structure between mitochondrial and nuclear genes has been noted in many studies [79, 80], and could be due to an incomplete natural history of the mitochondrial genome due to a range of factors such as a small effective population size, high mutation rate and patterns of introgression [81]. High differentiation based on mitochondrial genes but low based on microsatellite loci indicated complicated population history, such as the existence of multiple refugia populations during glacial periods followed by admixture in the interglacial periods, as reported in other species [14, 54]. Admixture of clusters was noted in several populations identified by BAPS analysis, suggesting ongoing introgression or incomplete lineage sorting. This was further supported by the SPLITSTREE analysis on mtDNA, in which haplotypes from the southern population (HBYCM) and one of its nearby population (SXJZZ) clustered in the same lineage.


Based on microsatellite loci and mtDNA, we found strong genetic differentiation in populations of PFM, but no obvious evidence for host-associated differentiation in PFM involving its common plant hosts, even though these plants alter PFM phenology and life history. Our study suggests that the geographical isolation and historical events in the Quaternary had a strong impact on current genetic differentiation of PFM in China. These strong effects may conceal other factors such as host-associated adaption and the impact of local environmental conditions. While host-associated adaptation of PFM might be present, it is not sufficient to generate separate gene pools of PFM that might reflect incipient speciation. Our study also suggests that geographical and historical factors need to be considered in experimental designs when attempting to assess adaptive divergence in PFM. Taking advantage of genomic tools, there are opportunities to investigate these processes further by incorporating a high density of markers across the genome that might include markers linked to loci under selection [82, 83].



approximate Bayesian computation


Bayesian analysis of population structure

cox1 :

cytochrome c oxidase subunit I


discriminant analysis of principal components


excluding null alleles


Hardy-Weinberg Equilibrium




isolation by distance


isolation by environment




linkage disequilibrium


last glacial maximum


multiple matrix regression with randomization


mitochondrial DNA


peach fruit moth


  1. Jenkins DG, Carey M, Czerniewska J, Fletcher J, Hether T, Jones A, Knight S, Knox J, Long T, Mannino M. A meta-analysis of isolation by distance: relic or reference standard for landscape genetics? Ecography. 2010;33(2):315–20.

    Google Scholar 

  2. Hewitt G. The genetic legacy of the quaternary ice ages. Nature. 2000;405(6789):907–13.

    Article  CAS  PubMed  Google Scholar 

  3. Orsini L, Vanoverbeke J, Swillen I, Mergeay J, Meester L. Drivers of population genetic differentiation in the wild: isolation by dispersal limitation, isolation by adaptation and isolation by colonization. Mol Ecol. 2013;22(24):5983–99.

    Article  PubMed  Google Scholar 

  4. Sexton JP, Hangartner SB, Hoffmann AA. Genetic isolation by environment or distance: which pattern of gene flow is most common? Evolution. 2014;68(1):1–15.

    Article  CAS  PubMed  Google Scholar 

  5. Shafer A, Wolf JB. Widespread evidence for incipient ecological speciation: a meta-analysis of isolation-by-ecology. Ecol Lett. 2013;16(7):940–50.

    Article  PubMed  Google Scholar 

  6. Yassin A, Debat V, Bastide H, Gidaszewski N, David JR, Pool JE. Recurrent specialization on a toxic fruit in an island Drosophila population. Proc Natl Acad Sci U S A. 2016;113(17):4771–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Wang IJ, Bradburd GS. Isolation by environment. Mol Ecol. 2014;23(23):5649–62.

    Article  PubMed  Google Scholar 

  8. Meirmans PG. The trouble with isolation by distance. Mol Ecol. 2012;21(12):2839–46.

    Article  PubMed  Google Scholar 

  9. Bagley RK, Sousa VC, Niemiller ML, Linnen CR. History, geography, and host use shape genome-wide patterns of genetic variation in the redheaded pine sawfly (Neodiprion lecontei). Mol Ecol. 2016;26(4):1022–44.

    Article  Google Scholar 

  10. Darwell CT, Fox KA, Althoff DM. The roles of geography and founder effects in promoting host-associated differentiation in the generalist bogus yucca moth Prodoxus decipiens. J Evol Biol. 2014;27(12):2706–18.

    Article  CAS  PubMed  Google Scholar 

  11. Hernández-Vera G, Mitrovic M, Jovic J, Tosevski I, Caldara R, Gassmann A, Emerson BC. Host-associated genetic differentiation in a seed parasitic weevil Rhinusa antirrhini (Coleptera: Curculionidae) revealed by mitochondrial and nuclear sequence data. Mol Ecol. 2010;19(11):2286–300.

    Article  PubMed  Google Scholar 

  12. Agunbiade TA, Coates BS, Kim KS, Forgacs D, Margam VM, Murdock LL, Ba MN, Binso-Dabire CL, Baoua I, Ishiyaku MF, et al. The spatial genetic differentiation of the legume pod borer, Maruca vitrata F. (Lepidoptera: Crambidae) populations in West Africa. Bull Entomol Res. 2012;102(5):589–99.

    Article  CAS  PubMed  Google Scholar 

  13. Mari-Mena N, Lopez-Vaamonde C, Naveira H, Auger-Rozenberg M-A, Vila M. Phylogeography of the Spanish moon moth Graellsia isabellae (Lepidoptera, Saturniidae). BMC Evol Biol. 2016;16:139.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Wei SJ, Cao LJ, Gong YJ, Shi BC, Wang S, Zhang F, Guo XJ, Wang YM, Chen XX. Population genetic structure and approximate Bayesian computation analyses reveal the southern origin and northward dispersal of the oriental fruit moth Grapholita molesta (Lepidoptera: Tortricidae) in its native range. Mol Ecol. 2015;24(16):4094–111.

    Article  PubMed  Google Scholar 

  15. Ehrlich PR, Raven PH. Butterflies and plants: a study in coevolution. Evolution. 1964;18(4):586–608.

    Article  Google Scholar 

  16. Pashley DP. Host-associated genetic differentiation in fall armyworm (Lepidoptera: Noctuidae): a sibling species complex? Ann Entomol Soc Am. 1986;79(6):898–904.

    Article  Google Scholar 

  17. Feder JL, Hunt TA, Bush L. The effects of climate, host plant phenology and host fidelity on the genetics of apple and hawthorn infesting races of Rhagoletis pomonella. Entomologia Experimentalis Et Applicata. 1993;69(2):117–35.

    Article  Google Scholar 

  18. Nosil P, Crespi BJ. Ecological divergence promotes the evolution of cryptic reproductive isolation. Proc R Soc B. 2006;273(1589):991–7.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Berlocher SH, Feder JL. Sympatric speciation in phytophagous insects: moving beyond controversy? Annu Rev Entomol. 2003;47(1):773–815.

    Article  Google Scholar 

  20. Matsubayashi KW, Issei O, Patrik N. Ecological speciation in phytophagous insects. Entomologia Experimentalis Et Applicata. 2010;134(1):1–27.

    Article  Google Scholar 

  21. Stireman JO, Devlin H, Abbot P. Rampant host- and defensive phenotype-associated diversification in a goldenrod gall midge. J Evol Biol. 2012;25(10):1991–2004.

    Article  PubMed  Google Scholar 

  22. Althoff DM, Svensson GP, Pellmyr O. The influence of interaction type and feeding location on the phylogeographic structure of the yucca moth community associated with Hesperoyucca whipplei. Mol Phylogenet Evol. 2007;43(2):398–406.

    Article  CAS  PubMed  Google Scholar 

  23. Funk DJ, Filchak KE, Feder JL. Herbivorous insects: model systems for the comparative study of speciation ecology. Genetica. 2002;116:251–67.

    Article  PubMed  Google Scholar 

  24. Heard SB. Use of host-plant trait space by phytophagous insects during host-associated differentiation: the gape-and-pinch model. Int J Ecol. 2012;2012:15.

    Article  Google Scholar 

  25. Antwi JB, Sword GA, Medina RF. Host-associated differentiation in a highly polyphagous, sexually reproducing insect herbivore. Ecol Evol. 2015;5(13):2533–43.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Simon JC, Carre S, Boutin M, Prunier-Leterme N, Sabater-Munoz B, Latorre A, Bournoville R. Host-based divergence in populations of the pea aphid: insights from nuclear markers and the prevalence of facultative symbionts. Proc R Soc Lond [Biol]. 2003;270(1525):1703–12.

    Article  Google Scholar 

  27. Dopman EB, Sword GA, Hillis DM. The importance of the ontogenetic niche in resource-associated divergence: evidence from a generalist grasshopper. Evolution. 2002;56(4):731–40.

    Article  PubMed  Google Scholar 

  28. Kim DS, Lee JH, Yiem MS. Spring emergence pattern of Carposina sasakii (Lepidoptera: Carposinidae) in apple orchards in Korea and its forecasting models based on degree-days. Environ Entomol. 2000;29(6):1188–98.

    Article  Google Scholar 

  29. Lei X, Li D, Li Z, Zalom FG, Gao L, Shen Z. Effect of host plants on developmental time and life table parameters of Carposina sasakii (Lepidoptera: Carposinidae) under laboratory conditions. Environ Entomol. 2012;41(2):349–54.

    Article  PubMed  Google Scholar 

  30. Kim DS, Lee JH. Egg and larval survivorship of Carposina sasakii (Lepidoptera: Carposinidae) in apple and peach and their effects on adult population dynamics in orchards. Environ Entomol. 2002;31(4):686–92.

    Article  Google Scholar 

  31. Hua L, Hua BZ. Preliminary study on the host-biotypes of peach fruit borer. Acta Phytophylacica Sinica. 1995;22(2):165–70.

    Google Scholar 

  32. Xu QG, Hua BZ. RAPD analysis on the speciation in host races of Carposina sasakii Matsumura (Lepidoptera: Carposinidae). Acta Entomol Sin. 2004;47(3):379–83.

    CAS  Google Scholar 

  33. Wang J, Yu Y, Li LL, Guo D, Tao YL, Chu D. Carposina sasakii (Lepidoptera: Carposinidae) in its native range consists of two sympatric cryptic lineages as revealed by mitochondrial COI gene sequences. J Insect Sci. 2015;15(1):1–6.

    Article  Google Scholar 

  34. Wang YZ, Cao LJ, Zhu JY, Wei SJ. Development and characterization of novel microsatellite markers for the peach fruit moth Carposina sasakii (Lepidoptera: Carposinidae) using next-generation sequencing. Int J Mol Sci. 2016;17(3):362.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Blacket MJ, Robin C, Good RT, Lee SF, Miller AD. Universal primers for fluorescent labelling of PCR fragments — an efficient and cost-effective approach to genotyping by fluorescence. Mol Ecol Resour. 2012;12(3):456–63.

    Article  CAS  PubMed  Google Scholar 

  36. Folmer O, Black M, Hoeh W, Lutz R, Vrijenhoek R. DNA primers for amplification of mitochondrial cytochrome coxidase subunit I from diverse metazoan invertebrates. Mol Mar Biol Biotechnol. 1994;3(5):294–9.

    CAS  PubMed  Google Scholar 

  37. Van Oosterhout C, Hutchinson WF, Wills DPM, Shipley P. Micro-checker: software for identifying and correcting genotyping errors in microsatellite data. Mol Ecol Notes. 2004;4(3):535–8.

    Article  CAS  Google Scholar 

  38. Park SDE. Trypanotolerance in west African cattle and the population genetic effects of selection. Ph D thesis: University of Dublin; 2001.

    Google Scholar 

  39. Chapuis MP, Estoup A. Microsatellite null alleles and estimation of population differentiation. Mol Biol Evol. 2007;24(3):621–31.

    Article  CAS  PubMed  Google Scholar 

  40. Guo S, Thompson EA. Performing the exact test of hardy-Weinberg proportion for multiple alleles. Biometrics. 1992;48(2):361–72.

    Article  CAS  PubMed  Google Scholar 

  41. Rousset F. Genepop ‘007: a complete re-implementation of the Genepop software for windows and Linux. Mol Ecol Resour. 2008;8(1):103–6.

    Article  PubMed  Google Scholar 

  42. Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, Valentin F, Wallace IM, Wilm A, Lopez R, et al. Clustal W and Clustal X version 2.0. Bioinformatics. 2007;23(21):2947–8.

    Article  CAS  PubMed  Google Scholar 

  43. Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30(12):2725–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Librado P, Rozas J. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics. 2009;25(11):1451–2.

    Article  CAS  PubMed  Google Scholar 

  45. Excoffier L, Lischer HE. Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and windows. Mol Ecol Resour. 2010;10(3):564–7.

    Article  PubMed  Google Scholar 

  46. Cheng L, Connor TR, Sirén J, Aanensen DM, Corander J. Hierarchical and spatially explicit clustering of DNA sequences with BAPS software. Mol Biol Evol. 2013;30(5):1224–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Jombart T, Devillard S, Dufour AB, Pontier D. Revealing cryptic spatial patterns in genetic variability by a new multivariate method. Heredity. 2008;101(1):92–103.

    Article  CAS  PubMed  Google Scholar 

  48. Wang IJ. Examining the full effects of landscape heterogeneity on spatial genetic variation: a multiple matrix regression approach for quantifying geographic and ecological isolation. Evolution. 2013;67(12):3403–11.

    Article  PubMed  Google Scholar 

  49. Huson DH, Bryant D. Application of phylogenetic networks in evolutionary studies. Mol Biol Evol. 2006;23(2):254–67.

    Article  CAS  PubMed  Google Scholar 

  50. Drummond AJ, Rambaut A. BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol Biol. 2007;7:214.

    Article  PubMed  PubMed Central  Google Scholar 

  51. Wei SJ, Shi BC, Gong YJ, Jin GH, Chen XX, Meng XF. Genetic structure and demographic history reveal migration of the diamondback moth Plutella xylostella (Lepidoptera: Plutellidae) from the southern to northern regions of China. PLoS One. 2013;8(4):e59654.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Cornuet JM, Pudlo P, Veyssier J, Dehne-Garcia A, Gautier M, Leblois R, Marin JM, Estoup A. DIYABC v2.0: a software to make approximate Bayesian computation inferences about population history using single nucleotide polymorphism, DNA sequence and microsatellite data. Bioinformatics. 2014;30(8):1187–9.

    Article  CAS  PubMed  Google Scholar 

  53. Lombaert E, Guillemaud T, Lundgren J, Koch R, Facon B, Grez A, Loomans A, Malausa T, Nedved O, Rhule E, et al. Complementarity of statistical treatments to reconstruct worldwide routes of invasion: the case of the Asian ladybird Harmonia axyridis. Mol Ecol. 2014;23(24):5979–97.

    Article  PubMed  Google Scholar 

  54. Cao LJ, Wei SJ, Hoffmann AA, Wen JB, Chen M, Vaclavik T. Rapid genetic structuring of populations of the invasive fall webworm in relation to spatial expansion and control campaigns. Divers Distrib. 2016;22(12):1276–87.

    Article  Google Scholar 

  55. Holm S. A simple sequential rejective method procedure. Scand Stat Theory Appl. 1979;6:65–70.

    Google Scholar 

  56. Duan X, Li Y, Men Q, Zhang M, Qiao X, Harari A, Chen M. Limited gene flow among Cydia pomonella (Lepidoptera: Tortricidae) populations in two isolated regions in China: implications for utilization of the SIT. Fla Entomol. 2016;99(sp1):23–9.

    Article  Google Scholar 

  57. Sun JT, Wang MM, Zhang YK, Chapuis MP, Jiang XY, Hu G, Yang XM, Ge C, Xue XF, Hong XY. Evidence for high dispersal ability and mito-nuclear discordance in the small brown planthopper, Laodelphax striatellus. Sci Rep. 2015;5:8045.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  58. Xun H, Li H, Li S, Wei S, Zhang L, Song F, Jiang P, Yang H, Han F, Cai W. Population genetic structure and post-LGM expansion of the plant bug Nesidiocoris tenuis (Hemiptera: Miridae) in China. Sci Rep. 2016;6:26755.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  59. Elameen A, Eiken HG, Knudsen GK. Genetic diversity in apple fruit moth indicate different clusters in the two most important apple growing regions of Norway. Diversity. 2016;8(2):10.

    Article  Google Scholar 

  60. Hua L, Wang JF, Hui XX. The esterase isozymes of the Carposina sasakii from different hosts. Acta Agriculturae Boreali-Occidentalis Sinica. 1995;4(3):95–6.

    Google Scholar 

  61. Grover A, Sharma PC. Development and use of molecular markers: past and present. Crit Rev Biotechnol. 2016;36(2):290–302.

    Article  CAS  PubMed  Google Scholar 

  62. Ishiguri Y, Shirai Y. Flight activity of the peach fruit moth, Carposina sasakii (Lepidoptera: Carposinidae), measured by a flight mill. Appl Entomol Zool. 2004;39(1):127–31.

    Article  Google Scholar 

  63. Hand SC, Denlinger DL, Podrabsky JE, Roy R. Mechanisms of animal diapause: recent developments from nematodes, crustaceans, insects, and fish. Am J Phys Regul Integr Comp Phys. 2016;310(11):R1193–R211.

    Google Scholar 

  64. Hoban S, Kelley JL, Lotterhos KE, Antolin MF, Bradburd G, Lowry DB, Poss ML, Reed LK, Storfer A, Whitlock MC. Finding the genomic basis of local adaptation: pitfalls, practical solutions, and future directions. Am Nat. 2016;188(4):379–97.

    Article  PubMed  PubMed Central  Google Scholar 

  65. Bernardi G, Azzurro E, Golani D, Miller MR. Genomic signatures of rapid adaptive evolution in the bluespotted cornetfish, a Mediterranean Lessepsian invader. Mol Ecol. 2016;25(14):3384–96.

    Article  CAS  PubMed  Google Scholar 

  66. Narum SR, Buerkle CA, Davey JW, Miller MR, Hohenlohe PA. Genotyping-by-sequencing in ecological and conservation genomics. Mol Ecol. 2013;22(11):2841–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  67. Liu JQ, Sun YS, Ge XJ, Gao LM, Qiu YX. Phylogeographic studies of plants in China: advances in the past and directions in the future. J Syst Evol. 2012;50(4):267–75.

    Article  Google Scholar 

  68. Shi YF, Zhao JD, Wang J. New understanding of quaternary glaciations in China. Shanghai: Shanghai Popular Science Press; 2011.

    Google Scholar 

  69. Yang F, Li YF, Ding X, Wang XQ. Extensive population expansion of Pedicularis longiflora (Orobanchaceae) on the Qinghai-Tibetan plateau and its correlation with the quaternary climate change. Mol Ecol. 2008;17(23):5135–45.

    Article  PubMed  Google Scholar 

  70. Wang J, Gao P, Kang M, Lowe AJ, Huang H. Refugia within refugia: the case study of a canopy tree (Eurycorymbus cavaleriei) in subtropical China. J Biogeogr. 2009;36(11):2156–64.

    Article  Google Scholar 

  71. Zhang H, Yan J, Zhang G, Zhou K. Phylogeography and demographic history of Chinese black-spotted frog populations (Pelophylax nigromaculata): evidence for independent refugia expansion and secondary contact. BMC Evol Biol. 2008;8(1):1.

    Article  CAS  Google Scholar 

  72. Beaumont MA, Zhang WY, Balding DJ. Approximate Bayesian computation in population genetics. Genetics. 2002;162(4):2025–35.

    PubMed  PubMed Central  Google Scholar 

  73. Lye G, Lepais O, Goulson D. Reconstructing demographic events from population genetic data: the introduction of bumblebees to New Zealand. Mol Ecol. 2011;20(14):2888–900.

    Article  CAS  PubMed  Google Scholar 

  74. Havill NP, Shiyake S, Lamb Galloway A, Foottit RG, Yu G, Paradis A, Elkinton J, Montgomery ME, Sano M, Caccone A. Ancient and modern colonization of North America by hemlock woolly adelgid, Adelges tsugae (Hemiptera: Adelgidae), an invasive insect from East Asia. Mol Ecol. 2016;25(9):2065–80.

    Article  PubMed  Google Scholar 

  75. Nosil P, Funk DJ, Ortiz-Barrientos D. Divergent selection and heterogeneous genomic divergence. Mol Ecol. 2009;18(3):375–402.

    Article  PubMed  Google Scholar 

  76. Guillot G, Rousset F. Dismantling the mantel tests. Methods Ecol Evol. 2013;4(4):336–44.

    Article  Google Scholar 

  77. Kim DS, Lee JH, Yiem MS. Temperature-dependent development of Carposina sasakii (Lepidoptera: Carposinidae) and its stage emergence models. Environ Entomol. 2001;30(2):298–305.

    Article  Google Scholar 

  78. Zhang B, Peng Y, Zhao XJ, Hoffmann AA, Li R, Ma CS. Emergence of the overwintering generation of peach fruit moth (Carposina sasakii) depends on diapause and spring soil temperatures. J Insect Physiol. 2016;86:32–9.

    Article  PubMed  Google Scholar 

  79. Bergamo LW, Fresia P, Azeredo-Espin AML. Incongruent nuclear and mitochondrial genetic structure of new world screwworm fly populations due to positive selection of mutations associated with dimethyl-and diethyl-organophosphates resistance. PLoS One. 2015;10(6):e0128441.

    Article  PubMed  PubMed Central  Google Scholar 

  80. Egger B, Koblmüller S, Sturmbauer C, Sefc KM. Nuclear and mitochondrial data reveal different evolutionary processes in the Lake Tanganyika cichlid genus Tropheus. BMC Evol Biol. 2007;7(1):137.

    Article  PubMed  PubMed Central  Google Scholar 

  81. Ballard JWO, Whitlock MC. The incomplete natural history of mitochondria. Mol Ecol. 2004;13(4):729–44.

    Article  PubMed  Google Scholar 

  82. Savolainen O, Lascoux M, Merilä J. Ecological genomics of local adaptation. Nat Rev Genet. 2013;14(11):807–20.

    Article  CAS  PubMed  Google Scholar 

  83. Storz JF. Using genome scans of DNA polymorphism to infer adaptive population divergence. Mol Ecol. 2005;14(3):671–88.

    Article  CAS  PubMed  Google Scholar 

Download references


We thank Li-Na Sun, Zong-Jiang Kang, Chun-Yan Jiang, Yu-Chun Zhu, Jiu-Yuan Ji, Liang Zhu, Long-Long Zhao, Rui-Qing Ma, Wan-Da Liu, Gui-Sheng Qiu, Li-Li Li, Jiang-Wen Qin and You-Gui Wang for their help with sample collection.


The research was funded by the Natural Science Foundation of Beijing Municipality (6162010), the National Natural Science Foundation of China (31472025), the National Basic Research Program of China (2013CB127600) and the Beijing Key Laboratory of Environmentally Friendly Pest Management on Northern Fruits (BZ0432). The funding bodies had no role in design or conclusion of the study.

Availability of data and materials

The data supporting the conclusions of this article are included within the additional files and is available at Dryad:

Author information

Authors and Affiliations



SJW conceived the study; SJW and YJG organized the collection of specimens; YZW and BYL, performed the molecular analyses; YZW, LJC, WS and SJW analyzed the data; YZW, SJW, AH, JYZ, YJG and LJC discussed the results and wrote the paper. All authors have read and approved the manuscript.

Corresponding author

Correspondence to Shu-Jun Wei.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional file

Additional file 1: Appendix S1.

Analysis of dispersal scenarios of Carposina sasakii using approximate Bayesian computation (DIYABC). Appendix S2. Power analysis of the microsatellite markers used in the study and effect of null alleles on population differentiation estimation. Figure S1. Phylogenetic trees of the haplotypes of the peach fruit moth Carposina sasakii based on mtDNA. The population name(s) followed by the haplotype indicates the haplotype found in the population(s). The value near the node indicates the corresponding divergence time (million years ago) of the two branches. In addition, the values of divergence times in red are given with a 95% highest posterior density (HPD). Table S1. Summary information on the biology of Carposina sasakii on different host plants. Table S2. Microsatellite loci used in this study. Table S3. Summary statistics for diversity of the 19 microsatellite loci examined in 14 populations of Carposina sasakii. Table S4. Genetic diversity of Carposina sasakii populations based on mitochondrial cox1 gene. (DOCX 741 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Wang, YZ., Li, BY., Hoffmann, A.A. et al. Patterns of genetic variation among geographic and host-plant associated populations of the peach fruit moth Carposina sasakii (Lepidoptera: Carposinidae). BMC Evol Biol 17, 265 (2017).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: