Codling moth, Cydia pomonella (L.) (Lepidoptera: Tortricidae), a highly invasive species, recently became established in the Hexi Corridor, which is a long narrow passage area with many oases surrounded by deserts and tall mountains in Gansu province, China. The corridor is an important temperate fruit growing region in northwestern China as well as a natural barrier to prevent C. pomonella from invading other fruit growing areas of the country. Since the codling moth was firstly reported, pome fruit damage in this corridor has been severe. The sterile insect technique (SIT) is considered a possible effective control tactic for integration in a future area-wide integrated pest management (AW-IPM) program against C. pomonella in the corridor. Knowledge of population genetics and more specifically of genetic differentiation and gene flow patterns may be important for developing AW-IPM strategies that include the SIT. In the current study, we collected C. pomonella samples from 8 populations distributed across 2 adjacent regions in the Hexi corridor that are geographically separated by stone deserts and high mountains. Eight microsatellite loci were used to investigate the genetic diversity, structure and differentiation of these 8 populations. Significant genetic differentiation was found between populations of each of the 2 regions, whereas populations within each region showed a similar genetic structure, demonstrated by higher Nm and lower FST values for population pairs within the same region than in pairs between the regions. Our findings indicate limited gene flow of C. pomonella between the 2 regions, which suggests that SIT can be implemented to control the pest in the Hexi Corridor of China.
The codling moth, Cydia pomonella (L.) (Lepidoptera: Tortricidae), is one of the most destructive pests of pome fruits (apple, pear and quince) and walnut orchards in most temperate regions of the world (Barnes 1991; Willett et al. 2009; Vreysen et al. 2010). It is generally assumed that the codling moth originated from southeastern Europe from where it spread around the world as a result of expanding apple and pear culture (Shel′Deshova 1967; Boivin et al. 2004; Franck et al. 2007; Meraner et al. 2008; Thaler et al. 2008). Codling moth is a serious invasive pest in China causing 20%-30% damage to pome fruit in those areas where it is distributed (Zhang et al. 2012). In 1957, codling moth was first detected in Xinjiang Province in Western China, which covers an area of more than 1.6 million km2 consisting mostly of mountains and deserts (Zhang 1957). Although C. pomonella has been the object of both internal and external quarantines in China, it still managed to spread from Xinjiang Province to the neighboring Hexi Corridor of Gansu Province in 1989 as a direct result of dense traffic and commercial transport (Shi et al. 2008). Thereafter, it took 20 years for codling moth to spread along the Hexi Corridor, which is composed of a1,000 km long string of oases along the northern edge of the high Tibetan Plateau, and is bordered on the north by the Gobi Desert and the grasslands of Outer Mongolia (Fig. 1). Irrigated by water from the snow-capped mountains, the oases of the Hexi Corridor are the main farming regions that supply most of the food including fruit to the arid Gansu Province. Apple and pear trees are planted widely in the oases and in the last several decades approximately 40% apple losses due to the codling moth have been reported in the corridor (Zhai et al. 2010).
Intensive spraying of chemical insecticides during the growing season is the main measure for controlling C. pomonella in China. However, the continuous use of broad-spectrum insecticides has not only become a serious concern for the environment and human health, but has also resulted in increased insecticide resistance of codling moth and in the loss of natural pest control (Bloem & Carpenter 2001; Reyes et al. 2009; Vreysen et al. 2010). The use of more environment-friendly pest management strategies such as the sterile insect technique (SIT) would have several advantages for the control of C. pomonella in China (Liu et al. 2012). The SIT is highly species-specific and friendly to the environment (Klassen 2005; Vreysen et al. 2010), and has been used successfully against a number of pest species, such as the Mediterranean fruit fly, Ceratitis capitata (Wiedemann) (Diptera: Tephritidae), the melon fly, Bactrocera cucurbitae (Coquilett) (Diptera: Tephritidae), the pink bollworm, Pectinophora gossypiella (Saunders) (Lepidoptera: Gelechiidae), the false codling moth, Thaumatotibia leucotreta (Meyrick) (Lepidoptera: Tortricidae), the New World screwworm, Cochliomyia hominivorax (Coquerel) (Diptera: Calliphoridae) and the tsetse fly Glossina austeni Newstead (Diptera: Glossinidae) (Tan 2000; Vreysen et al. 2000; Wyss 2000; Bloem & Carpenter 2001; Hendrichs et al. 2005; Klassen & Curtis 2005; Parker & Mehta 2007). The successful suppression of the codling moth in the Okanagan Valley of British Columbia, Canada for the last 20 years is arguably one of the most successful AW-IPM programs that integrated the SIT with other control tactics (Klassen 2005; Bloem et al. 2007 b; Vreysen et al. 2010). Given that the invasion of the Hexi Corridor by the codling moth is fairly recent, and that the oases of the corridor are isolated, the SIT could be considered as a potential part of an effective strategy of managing codling moth in the Hexi Corridor (Li 2013).
Population genetics data and more specifically genetic differentiation and gene flow patterns of target insect populations can provide sound baseline information to develop efficient insect pest-management strategies including the SIT (Klassen 2005; Krafsur 2005; Timm et al. 2006; Chen & Dorn 2010; Vreysen et al. 2010; Zheng et al. 2013). These data can reveal the distance over which members of a given species typically disperse, estimate the degree of genetic isolation of target populations from each other and from untargeted populations, determine whether a sub-population or sibling species exists, and establish the origin of outbreaks or reintroductions (Krafsur 2005; Chen & Dorn 2010). In practice, if immigration of gravid female moths into a given agro-ecosystem can be kept to a minimum, the SIT would be more effective. As moths released in AW-IPM programs with a SIT component should be capable of mating with females from different geographic regions, it is important to know to which degree an insect species migrates between different field populations (Krafsur 2005; Chen & Dorn 2010). Moreover, knowledge on the genetic differentiation and rate of gene flow among populations can provide an estimate of the rate of invasion or reinvasion of eradicated areas by unchallenged populations (Krafsur 2005).
Although the Hexi Corridor has always been considered an important natural barrier that prevents C. pomonella from invading the more southerly major fruit growing areas, C. pomonella has been established in several oases of the Hexi Corridor. However, the genetic differentiation and gene flow of the different codling moth populations is still unknown. In this study, we collected C. pomonella samples from 2 geographically adjacent regions (Jiuquan and Zhangye) in the Hexi Corridor separated by unpopulated tall mountains and stone deserts serving as a geographic barrier. Using microsatellites and the method developed by Men et al. (2013) we aimed to characterize the genetic structure and gene flow pattern of the C. pomonella populations in the 2 regions. The findings are of interest for the development and implementation of a rational AW-IPM strategy that includes a SIT component.
Materials and Methods
Cydia pomonella was sampled in 2 isolated adjacent regions (Jiuquan and Zhangye) of the Hexi Corridor (Fig. 1). Adult moths were sampled with 6 pheromone traps deployed for 5 consecutive days in apple orchards. Traps were replaced every 24 h and relocated to a new site within each sublocation. This resulted in a total of 30 traps per sublocation. To avoid sibling sampling, only one individual was taken from each trap and preserved in a 10 mL Falcon tube filled with ethanol. The same sampling procedure was repeated 3 times in each sublocation during the apple growing season giving a total of 45-65 moths collected in each sublocation. Around half of the samples of each region were randomly selected for genotyping in the current study. We use the term “population” for C. pomonella moths sampled from a single orchard (sampling unit). Samples from 4 locations in each of the 2 regions were used (Fig. 1).
DNA EXTRACTION AND MICROSATELLITE GENOTYPING
Genomic DNA was extracted from 8-10 mg of insect material using the Easy Pure TM Genomic DNA Kit (TransGen Biotech Co., Ltd. Beijing, China) according to the manufacturer's protocol for animal tissues. DNA was eluted in 200 μL of ultra pure water and stored at −20 °C. All samples were genotyped at 8 polymorphic microsatellite loci (Cp1.62, Cp2.39, Cp2.P, Cp3.56, Cp3.169, Cp3.K, Cp4.56 and Cp4.129) isolated from Cydia pomonella (Franck et al. 2005). Three primers, which included the forward primer synthesized with an M13 (−21) at the 5′ end, the reverse primer and the FAM fluorescent dye labeled M13 (−21) primer, were used for amplification each microsatellite locus in a PCR reaction (Schuelke 2000). PCR reactions for the 8 loci were carried out in a mix of 12.5 μL 2 × Taq Master Mix (containing 0.05 U/μL Taq DNA Polymerase, 2 × Taq PCR Buffer, 3 mM MgCl2 and 400 μM dNTP mix) (Beijing CoWin Biotech Co, Ltd. Beijing, China), 0.2 μM of each forward primer, 0.8 μM of each reverse primer, 0.8 μM M13 primer, and 1 μL genomic DNA (10-30 ng/μL). The PCR involved denaturation at 95 °C for 10 min, followed by 30 amplification cycles consisting of 95 °C for 30 s, 45 s at the primer-specific annealing temperature (57 °C for Cp1.62, 60 °C for Cp2.39, 58 °C for Cp2.P, 60 °C for Cp3.56, 54 °C for Cp3.169, 48 °C for Cp3.K, 60 °C for Cp4.129 and 56 °C for Cp4.56), 72 °C for 45 s, then 8 cycles consisting of 95 °C for 30 s, 53 °C for 45 s and 72 °C for 45 s, and a final step at 72 °C for 10 min. Genotyping was carried out on an ABI3730XL automated DNA sequencer (Applied Biosystems, Foster City, California, USA) and GENEMAPPER version 4.0 (Applied Biosystems, Foster City, California, USA).
Null allele and frequency were checked with Micro-Checker version 2.2.3 (Van Oosterhout et al. 2004). FSTAT Version 188.8.131.52 (Goudet 2001) was utilized for calculating the number of alleles (Na), and for calculating the mean allelic richness per locus (R) for all the sampling populations with a rarefaction method that adjusted the observed number of alleles to a common size for each population sample. Observed heterozygosity (H o) and expected heterozygosity (H e) were calculated using microsatellite analyzer (MSA) version 3.15 (Dieringer & Schlötterer 2003). The Hardy-Weinberg equilibrium (HWE) and P value for Hardy-Weinberg equilibrium (HWE-P) were calculated using GENEPOP version 4.0.1 (Rousset 2008).
Analysis of molecular variance (AMOVA), inbreeding coefficients (F IS), fixation index among groups (F CT), fixation index among populations within groups (F SC) and fixation index within populations (FST) were performed using ARLEQUIN version 184.108.40.206 (Weir & Cockerham 1984; Excoffier &Lischer 2010). We calculated the number of migrants moving between each pair of populations (N m) per generation as (1/F ST − 1)/4 (Excoffier &Lischer 2010; Munshi-South 2012). For AMOVA, samples were arbitrarily grouped according to the collection regions: (1) Jiuquan and (2) Zhangye. The matrices of genetic distance F ST/(1 − F ST) and the geographic distance (ln) between 8 C. pomonella populations were used to test isolation by distance (IBD) with 10,000 permutations (Mantel 1967). The IBD analysis was performed using the ZT software package (Bonnet &Van der Peer 2002).
The Bayesian clustering approach was used to estimate the population structure in STRUCTURE version 2.3.3 (Pritchard et al. 2000) with a burn-in period of 50,000 iterations and 1 million Markov chain Monte Carlo (MCMC) repetitions. We set the number of the cluster (K) from 1 to 10 and repeated 10 times, and the most likely number of clusters was calculated by delta K methods (Evanno et al. 2005). The graphical display of genetic structure was produced by DISTRUCT (Rosenberg 2004). A neighbor-joining (NJ) tree based on Nei's genetic distances among populations was constructed by MEGA version 5 (Tamura et al. 2011).
A total of 214 C. pomonella individuals were genotyped across all 8 microsatellite loci. All 8 loci were polymorphic and the number of alleles at each locus ranged from 7 to 25. The frequency of null alleles ranged from 0.049 to 0.241, which was similar to previous reports of C. pomonella (Franck et al. 2007; Men et al. 2013) and is typical for Lepidoptera (Dakin & Avise 2004; Meglécz et al. 2004). The average number of alleles per population ranged from 4.3 to 6.9 (Table 1). The allele richness value per locus per population was between 4.07 and 6.64. The value of observed heterozygosity (Ho) was between 0.190 and 0.534, whilst the expected heterozygosity (He) values ranged from 0.311 to 0.569. Three populations revealed significant departures from HWE, together with large mean values of the inbreeding coefficient, FIS, ranging from 0.422 to 0.543. This indicates the existence of heterozygote deficiencies in these 3 populations. The genetic indices indicated that the 8 populations showed similar genetic diversity over the 8 microsatellite loci.
AMOVA results (Table 2) indicated significant genetic differentiation at the 3 hierarchical levels (F CT = 0.097, P < 0.001; F SC = 0.085, P < 0.001; F ST = 0.173, P < 0.001) with populations grouped according to the 2 sampling regions. Although the proportion of explained variance was low (9.7%), significant genetic differentiation existed among populations from the 2 regions.
Gene flow over generations as revealed by N m is shown in Table 3. The results indicate that N m values within populations from Jiu quan region were between 2.103 and 133.276, with an average value of 25.561, whilst N m values within populations from Zhangye region ranged from 1.108 to 4.765, with an average value of 2.495. The N m values of the population pairs from the 2 different regions (Jiuquan and Zhangye) averaged 1.386, with the values ranging between 0.649 and 1.915.
Sample information and population statistics for Cydia pomonella screened with 8 microsatellites.
Pairwise F ST values for each pair of populations over the 8 loci ranged from 0.002 to 0.259. The results show significant genetic variance between all pairs of populations except one pair (JJT and JYD) that was sampled from 2 neighboring sublocations in Jiuquan region (Table 3). The average F ST value of population pairs from the Jiuquan region was 0.059, versus 0.110 for population pairs from the Zhangye region. The average F ST value of the population pairs from the 2 different regions (Jiuquan and Zhangye) was 0.178. The Mantel test for the 8 C. pomonella populations from the 2 regions revealed no correlation between genetic distances and geographic distances (r = 0.315, P = 0.133).
In the Bayesian analysis performed with the software STRUCTURE 2.3.3, the most likely value of K was 2 (Pritchard et al. 2000; Evanno et al. 2005), indicating that the 8 populations included in this study can be assigned to 2 distinct clusters, hereafter referred to as cluster 1 and cluster 2. The proportions of each population that contributed to each of the 2 clusters are shown in Fig. 2. Cluster 1 was mainly made up of individuals sampled from 4 sublocations of the Jiuquan region (JJT, 0.981 implying that 98.1% of the individuals from the JJT sublocation contributed to this cluster, and that corresponding contributions from the remaining 3 sublocations in this region were: JYD, 0.966; JXD, 0.944 and JZZ, 0.774). Likewise, cluster 2 consisted largely of individuals of the 4 populations collected in the 4 sublocations of Zhangye region and with corresponding contributions as follows: ZLT, 0.885; ZNJ, 0.936; ZMY, 0.906; and ZXM, 0.938. Populations within the same region showed similar genetic structures, but were different from the genetic structures of populations in the other region.
The NJ tree based on the Nei's genetic distances is shown in Fig. 3. The 8 populations from the 2 separate regions were divided into 2 clades in the NJ tree. All 4 populations (JJT, JYD, JXD and JZZ) from the Jiuquan region made up one clade (clade 1), and the 4 populations (ZLT, ZNJ, ZMY and ZXM) from the Zhangye region constituted the other clade (clade 2). Both the results of the NJ tree and the STRUCTURE analysis consistently and clearly indicated a distinct clustering according to region.
Analysis of molecular variance of populations from Jiuquan and Zhangye regions of the Hexi Corridor.
We analyzed the genetic variation at 8 microsatellite loci in 214 individuals from the Jiuquan and Zhangye regions of the Hexi Corridor that were separated by mountains and stone deserts. There was significant genetic differences between the various populations of the 2 regions. Populations within the same region showed similar genetic structures, but their genetic structures differed from those of the populations from the other region. When compared with population pairs from the other region, population pairs within a region showed greater gene flow (N m) and smaller F ST values. Our results suggest only limited gene flow between populations of the Jiuquan region and the Zhangye region.
The flight capacity of a large proportion of a C. pomonella population is low, and only a small proportion of a population is capable of long distance flights (Schumacher et al. 1997; Dorn et al. 1999). In the current study, AMOVA analyses among the 8 populations from the 2 sampling regions showed significant genetic differentiation. Based on cluster analysis, populations sampled in the 2 regions differed in genetic structure, and they formed 2 distinct clades in the NJ tree. The codling moth has traditionally been regarded as a sedentary insect (Keil et al. 2001; Timm et al. 2006; Thaler et al. 2008; Chen & Dorn 2010), and the limited gene flow implied by the significantly different population genetic characteristics of the 2 separate regions is likely the result of the limited dispersal capacity of C. pomonella. This was verified by mark-release-recapture studies in the field and flight mill experiments in the laboratory (Schumacher et al. 1997; Dorn et al. 1999; Keil et al. 2001). Our results are consistent with findings in South Africa, northern Italy and Switzerland, where limited gene flow and significant genetic differentiation between different local geographic regions were found (Timm et al. 2006; Thaler et al. 2008; Chen & Dorn 2010). Similar results were also reported in the closely related species with small flight capacity, Grapholita molesta (Timm et al. 2008; Torriani et al. 2010; Kirk et al. 2013 a, b; Zheng et al. 2013 ), as well as in aphids (Loxdale et al. 1993; Miller et al. 2003).
Pairwise population differentiation estimates (FST) averaging 8 loci between 8 populations of Cydia pomonella with P value (below the diagonal) and gene flow (Nm) values (above the diagonal).
We found similar genetic diversities and similar genetic structures within each region in this study. In the current analysis, Nm values were greater and FST values were smaller within population pairs from the same region in comparison to those found in population pairs from the 2 regions, indicating limited gene flow between the regions. Dispersal caused by humans might have contributed to these results (Mazzi & Dorn 2012; Kirk et al. 2013 a, b). Population genetics of local populations of sedentary species such as C. pomonella and G. molesta can obviously be affected by anthropogenic influences including neurotoxic pesticide spraying, regional fruit transportation, the recent spread of modern fruit culture and regional integrated production and marketing programs (Chen & Dorn 2008; Torriani et al. 2010; Mazzi & Dorn 2012). Anthropogenic influences are considered important reasons for the small genetic differentiation of C. pomonella as found in populations in France, Chile and Croatia (Franck et al. 2007; Fuentes-Contreras et al. 2008; Voudouris et al. 2012). Actually, local Chinese governments promote centralized purchasing of fruit and distribution of apple and pear cultivars that might promote gene flow and facilitate genetic exchange of G. molesta populations within regions, as is also the case for cooperatives in Italy (Torriani et al. 2010). Two populations (JJT and JYD) from the Jiuquan region were very similar in genetic diversity and showed very large Nm values, which can be explained by the frequent fruit transportation between the 2 regions offering ample opportunity for interbreeding between the populations of the 2 sublocations.
The SIT has been successfully used to suppress C. pomonella in the fruit growing Okanagan Valley of British Columbia in Canada (Klassen 2005; Bloem et al. 2007b; Vreysen et al. 2010). The Hexi Corridor is considered an important natural barrier that prevents the codling moth from invading into the more southerly major fruit growing areas of China. Our results show that gene flow between the 2 investigated regions is limited, which creates the possibility of developing and implementing a sequential eradication strategy for the sustainable removal of the codling moth populations in the 2 regions. The use of the SIT is considered a possibility for use in an AW-IPM strategy in view of the absence of mating barriers between codling moth populations worldwide (Taret et al. 2010; Liu et al. 2012; Li 2013).
We thank Yong Zhang for help with the insect sampling. This work was supported by grants of the National Natural Science Foundation of China (No. 31071687 and No. 31272036), the National Key Technology R&D Program (No. 2012BAK11B03) from the Ministry of Science and Technology of China, and the Open Project of Fund of Hubei Key Laboratory of Animal Embryo Engineering and Molecular Breeding (No. 2012-2D-128). This work was part of the FAO/IAEA Coordinated Research Project on Increasing the Efficiency of Lepidoptera SIT by Enhanced Quality Control.