Genetic structure and diversity of a newly invasive ... - Springer Link

4 downloads 66364 Views 529KB Size Report
Aug 8, 2012 - the historical events, such as bottleneck or founder effects ... collected from 12 apple orchards in the C. pomonella .... 19.0 (SPSS Inc. Chicago).
Biol Invasions (2013) 15:447–458 DOI 10.1007/s10530-012-0299-5

ORIGINAL PAPER

Genetic structure and diversity of a newly invasive species, the codling moth, Cydia pomonella (L.) (Lepidoptera: Tortricidae) in China Qiu-Lei Men • Mao-Hua Chen Ya-Lin Zhang • Ji-Nian Feng



Received: 14 November 2011 / Accepted: 25 July 2012 / Published online: 8 August 2012 Ó Springer Science+Business Media B.V. 2012

Abstract Cydia pomonella (L.) was firstly reported in China in the 1950s and considered as one of the most serious invasive pest in fruit orchards of China. It spread rapidly from the original site in Xinjiang to other northwestern regions. The pest has further penetrated northeastern China since 2006. With its rapid invasion rate, most pome fruit production areas of China are being threatened. As yet there has been no research into the genetic diversity and structure of the codling moth population in China. We investigated the genetic variations of 12 C. pomonella populations collected from the main distribution regions (Xinjiang, Gansu and Heilongjiang Provinces) in China and compared them with one German and one Swiss population using eight microsatellites loci to infer the characteristics of genetic diversity and genetic structure. We observed sequential loss of genetic diversity and significant structuring associated with distribution but no significant correlation between genetic distance and geographic distance among northwestern populations. There was no genetic evidence for bottleneck effects in any of the populations. The results suggest that the loss of genetic diversity in C. pomonella

Q.-L. Men  M.-H. Chen  Y.-L. Zhang  J.-N. Feng (&) Key Laboratory of Plant Protection Resources and Pest Management of Ministry of Education, Entomological Museum, Department of Entomology, Northwest A & F University, Yangling 712100, Shaanxi, China e-mail: [email protected]

populations resulted from the successive colonization of founder populations. Recent invasion history led to the lack of any bottleneck effect. The high level of population genetic structuring is related to the weak flight capacity of the codling moth and the humanaided dispersal rather than to geographic distance. These genetic data not only provide us with an understanding of the micro-evolutionary processes related to successful biological invasions, but also provide guidance for pest management strategies. Keywords Cydia pomonella  Invasive  Microsatellites  Genetic structure  Genetic diversity  China

Introduction Many factors influence the genetic structure of animal populations. For an invasive species, the factors which influence population genetic structure are often related to its history of invasion (Grapputo et al. 2005; Watts et al. 2010). Populations that have been in stable environments and are connected by dispersal over long periods of time will reach a genetic equilibrium where the loss of alleles as a result of drift is balanced by the introduction of new alleles through migration (Wright 1951). In contrast, in newly established populations, genetic structure is often determined by the historical events, such as bottleneck or founder effects, vicariance events and range of dispersal, rather

123

448

than gene flow and genetic drift (Bohonak 1999). The levels of genetic diversity and genetic structure for the newly invasive populations are affected by the recent expansion of the species, the numbers of founders and passive dispersal by anthropic factors (Ramstad et al. 2004; Dlugosch and Parker 2008; Watts et al. 2010). The codling moth, Cydia pomonella (L.), is a destructive fruit pest in fruit orchards (Barnes 1991; Willett et al. 2009). This moth’s rapid spread to nearly all apple growing areas throughout the world has largely been due to a combination of the growing popularity of pome fruit culturing and accidentally human-aided dispersal (Boivin et al. 2004; Meraner et al. 2008). Host plants for the codling moth include apple, pear, apricot, plum, peach, nectarine, cherry, quince and walnut (Barnes 1991). The first reported sightings of the codling moth in China were in 1957 (Zhang 1957). It has now been observed throughout Xinjiang which covers more than 1.6 million square kilometers of mostly mountains, deserts and unpopulated areas. From this first site, it took more than 20 years to reach Gansu, spreading along the Hexi Corridor. In 2006, new invasion sites were found in Heilongjiang, the first report in northeastern China. As a notoriously invasive pest, the codling moth has caused a great deal of damage to Chinese production of pomes (primarily apple and pear) since it was first reported, giving rise to substantial economic losses for affected farmers every year (Wan et al. 2009). A better understanding of the micro-evolutionary processes and ecological strategies related to successful biological invasions is necessary in order to develop effective pest management practices (Dorn et al. 1999; Miller et al. 2003). Previous researchers focused on the genetic diversity and genetic structure of older C. pomonella populations from Europe and Africa (Timm et al. 2006; Espinoza et al. 2007; Franck et al. 2007; Fuentes-Contreras et al. 2008; Chen and Dorn 2009; Franck and Timm 2010), but there has been no investigation of the population genetic structure and genetic diversity of the codling moth in China. Molecular markers have greatly facilitated the study of population genetic structure and diversity. Microsatellites, AFLP and mitochondrial genes have been successfully used to characterize many insect species (De Barro 2005; Paupy et al. 2005; Mozaffarian et al. 2007; Meng et al. 2008; Orsini et al. 2008; Ahern et al. 2009; Cristescu et al. 2010). In this study, we use eight microsatellite loci to investigate the

123

Q.-L. Men et al.

genetic diversities and genetic structures of 12 populations collected from the main distribution regions of C. pomonella in China. Our main objectives were to analyze the genetic diversity and genetic structure of codling moths in newly-invaded ecosystems in China.

Materials and methods Sampling A total of 339 Cydia pomonella individuals were collected from 12 apple orchards in the C. pomonella distribution regions of China from 2008 to 2010. In Xinjiang, an infested fruit was collected from each fruit tree per orchard and 1 second or third instar larva per fruit was used. In the newly invaded Gansu and Heilongjiang, more than 30 pheromone traps were used in each apple orchard ([5 ha), and only one moth from each trap was used. All the sample locations are illustrated in Fig. 1. For making comparisons with populations in China, C. pomonella samples collected by pheromone traps in Germany (samples from Karlsruhe in July, 2008) and Switzerland (samples from Chur in June, 2009) were included in the analysis. All the samples were stored at -20 °C prior to genetic analyses. DNA extraction and microsatellite genotyping Genomic DNA was extracted from the thorax and leg of each individual sample using DNeasy Tissue Kit (QIAGEN, Hilden, Germany). Extraction was performed according to the bench protocol for animal tissues. Before extraction, xylene was used to dissolve the glue covering specimens taken from pheromone traps. Eight loci (Cp1.62, Cp2.39, Cp2.P, Cp3.56, Cp3.169, Cp3.K, Cp4.56 and Cp4.129) were selected for the present study, based on their characterization in the literature as highly polymorphic (Franck et al. 2005). Each locus was amplified using fluorescently-labeled primers (FAM) (Schuelke 2000). Polymerase chain reactions (PCR) were performed using a S1000 Thermal Cycler (BIO-RAD, Hercules, CA, USA) in a total volume of 25 ll, containing 19 PCR amplification buffer (Takara, Dalian, China), 2 mM MgCl2, 0.2 mM of each dNTP (Takara, Dalian, China), 0.2 lM of each forward primer, 0.8 lM of each reverse primer, 0.8 lM M13 primer, 1.0 U Taq polymerase (Takara, Dalian, China), and 1 ll genomic DNA (10–30 ng/ll).

Genetic structure and diversity of a newly invasive species

449

Fig. 1 a Distribution regions and sample locations for the 12 C. pomonella populations from China, the codes for C. pomonella populations are explained in Table 1; the light gray region represents the geographic distribution of C. pomonella in Xinjiang; the dark gray region represents the geographic distribution in Gansu; the black region represents the

geographic distribution in Heilongjiang; b sample locations of the two populations from Heilongjiang; c sample locations of the 10 populations from Xinjiang and Gansu. Pie charts show the mean membership fractions to each of the four genetic clusters in the populations by cluster analysis using STRUCTURE software

PCR amplification was employed with 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, 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. To examine the length and genotype of the amplified PCR products, an ABI3730XL automated DNA sequencer (Applied Biosystems, Foster City, CA, USA) and GENESCAN version 4.0 (Applied Biosystems, Foster City, CA, USA) were used.

Data analysis Population genetic diversity We used Micro-Checker version 2.2.3 (Van Oosterhout et al. 2004) to check the data for null alleles, stutter-errors and large allele drop-out. Using MICROSATELLITE ANALYSER (MSA) version 3.15 (Dieringer and Schlo¨tterer 2003), we calculated number of alleles, and observed versus expected heterozygosity. Using GENEPOP version 4.0.1

123

450

(Rousset 2008), we estimated the Hardy–Weinberg equilibrium (HWE), HWE-P and genotypic linkage disequilibrium. Allelic Diversity Analyzer (ADZE) version 1.0 (Szpiech et al. 2008) was utilized for calculating the mean allelic richness per locus for all the sampling populations with a rarefaction method that adjusted the observed number of alleles to a common size for each population sample. We used BOTTLENECK version 1.2.02 (Cornuet and Luikart 1996; Piry et al. 1999) to assess deviation from expected heterozygote excess relative to allelic diversity across all eight loci, based on the loss of allelic diversity exceeding heterozygosity during a bottleneck effect (Nei et al. 1975; Hufbauer et al. 2004; Pruett and Winker 2005). Two appropriate models, stepwise mutation model (SMM) and twophase model (TPM) were utilized for this analysis (Di Rienzo et al. 1994; Primmer et al. 1998). Probability values were determined using one-tailed Wilcoxon tests for heterozygote excess (Luikart et al. 1998; Piry et al. 1999). Values of P less than 0.05 were considered significant. The significances were corrected for multiple comparisons by the Bonferroni method (Rice 1989) using SPSS version 19.0 (SPSS Inc. Chicago). Population genetic structure Population structure was analyzed using STRUCTURE version 2.3.3 (Pritchard et al. 2000) based on Bayesian clustering approach with a burn-in period of 50,000 iterations and one million Markov chain Monte Carlo (MCMC) repetitions. We used the admixture ancestry model and the correlated allele frequency model. We performed 20 independent runs for each K to confirm consistency across runs, tested K from 1 to 10 and DK methods (Evanno et al. 2005) was used to calculate the number of genetic clusters (K). The graphical display of genetic structure was produced by DISTRUCT (Rosenberg 2004). Analysis of molecular variance (AMOVA) was performed using ARLEQUIN version 3.0 (Excoffier et al. 2005), along with calculating the inbreeding coefficients (FIS), the pair fixation indices (FST) and their significance with 10,000 permutations (Weir and Cockerham 1984). The significances of fixation indices were corrected for multiple comparisons by the Bonferroni method (Rice 1989) using SPSS version 19.0 (SPSS Inc. Chicago). For AMOVA, samples were

123

Q.-L. Men et al.

arbitrarily grouped according to the collection regions: (a) Xinjiang Province; (b) Gansu Province; (c) Heilongjiang Province; and (d) western Europe. Referring to the criterion for genetic differentiation by Wright (1978), we defined genetic differentiation as low for FST \ 0.05, moderate for 0.05 \ FST \ 0.15, high for 0.15 \ FST \ 0.25, and very high for FST [ 0.25. In order to test for isolation by distance (IBD), the matrices of genetic distance FST/(1-FST) and the geographic distance (ln) between 10 populations from northwestern China as well as all 12 Chinese sampling populations were compared using the Mantel test with 10,000 permutations (Mantel 1967). This analysis was performed using the ZT software package (Bonnet and Van der Peer 2002).

Results Microsatellite markers and genetic diversity We genotyped 339 individuals for eight microsatellite loci. Population statistics for all investigated populations are given in Table 1. All eight selected loci were polymorphic and the number of alleles per locus ranged from four to 16. A 21.9 % disequilibrium was found in 392 pairs of loci, but did not appear at certain pairs of loci, which was not considered to be a result of physical linkage. The frequency of null alleles ranged from 0.010 to 0.203, values typical for lepidopterans (Dakin and Avise 2004; Megle´cz et al. 2004). Twelve samples with two non-amplified loci were omitted. The mean number of alleles per locus ranged from 4.3 to 12.6 (Table 2). Two populations from Heilongjiang Province in northeastern China had the largest number of alleles (Don, 12.6; Mud, 10.6). Of populations from northwestern China, the Ili population showed the highest value of mean number of alleles (9.6), followed by the Urumqi population (8.6) and Kuytun population (8.4). In comparison with the Chinese codling moth populations, the German and Swiss population had a relatively high mean number of alleles: 8.6 and 10.4 respectively. The mean value of observed heterozygosity was between 0.201 and 0.550, as compared to the expected heterozygosity of between 0.395 and 0.821. The two populations from Heilongjiang had the largest observed heterozygosity (Mud, 0.517; Don, 0.482), followed by populations from the north of Xinjiang (Kuy, 0.387; Ili, 0.356; Jin, 0.321).

Genetic structure and diversity of a newly invasive species

451

Table 1 Population details of 12 Cydia pomonella populations collected from the main distribution regions in China and two populations from other countries including collection region, location, location code, longitude, latitude and collection date Collection region

Location

Location code

Longitude (°E)

Latitude (°N)

Collection date

Gansu Province

Zhangye

Zha

100.38

38.85

Mar. 2010

Jiuquan

Jiu

98.42

39.67

Mar. 2010

Dunhuang

Dun

95.50

41.58

Jul. 2009

Ili

Ili

81.32

43.95

Aug. 2008

Jinghe

Jin

83.53

45.12

Aug. 2008

Kuytun

Kuy

84.93

44.45

Aug. 2008

Urumqi

Uru

87.47

43.90

Aug. 2008

Kumul

Kum

93.67

42.83

Aug. 2008

Korla

Kor

86.13

41.60

Aug. 2008

Kashgar

Kas

76.02

39.53

Aug. 2008

Dongning

Don

131.23

44.10

Aug. 2010

Mudanjiang

Mud

129.57

44.58

Aug. 2010

Switzerland (Chur)

Swi

9.52

46.85

Jun. 2009

Germany (Karlsruhe)

Ger

8.40

49.02

Jul. 2008

Xinjiang Province

Heilongjiang Province Other countries

Table 2 Population statistics for Cydia pomonella investigated using eight microsatellite loci Population code

N

NA

r

FIS

Ho

He

HWE-P

TPM

SMM

Zha

25

5.3

5.14

0.484

0.201

0.416

0.006

0.99023

1.00000

Jiu

25

4.3

4.17

0.298

0.272

0.434

0.118

0.90234

0.98633

Dun Kum

33 33

5.6 6.0

5.50 6.13

0.400 0.257

0.220 0.298

0.460 0.495

0.006 0.004

0.97266 0.99609

0.99414 0.99609

Kas

25

5.5

5.43

0.330

0.232

0.395

0.020

0.90234

0.98047

Kor

11

4.9

4.46

0.097

0.290

0.475

0.268

0.96094

0.98828

Uru

26

8.6

8.30

0.486

0.262

0.552

0.125

0.98633

0.99609

Jin

21

6.9

6.82

0.267

0.321

0.519

0.211

0.99414

0.99805

Ili

27

9.6

9.16

0.354

0.356

0.616

0.064

0.98633

0.99805

Kuy

26

8.4

8.08

0.275

0.387

0.537

0.054

0.99805

1.00000

Mud

24

10.6

10.38

0.295

0.517

0.821

0.084

0.62891

0.99609

Don

24

12.6

12.23

0.293

0.482

0.791

0.160

0.98047

0.98047

Ger

15

8.6

7.96

0.259

0.550

0.770

0.192

0.52734

0.98047

Swi

24

10.4

9.63

0.304

0.502

0.801

0.105

0.52734

0.99023

This table includes population code, sample size (N), allelic richness (r), mean number of alleles per locus (NA), multilocus estimate of inbreeding coefficient (FIS), observed heterozygosity (Ho), excepted heterozygosity (He), P value for Hardy–Weinberg equilibrium (HWE) analysis, and P value for heterozygote excess using two phase model (TPM) and stepwise mutation model (SMM). Significant departures from Hardy–Weinberg equilibrium are given in bold, P \ 0.05. NA, FIS, Ho, He and HWE-P are all indicated by mean values over eight loci. Probability values for heterozygote excess were determined using one-tailed Wilcoxon tests, and P values less than 0.05 were considered significant. The significances were tested for multi comparisons by the Boferroni method for k = 14, P \ 0.05

The German and Swiss population also showed relatively high level of observed heterozygosity (Ger, 0.550; Swi, 0.502). Four of the 14 populations revealed

significant departures from HWE, together with high mean values of FIS ranging from 0.097 to 0.486, indicating the existence of heterozygote deficiencies (Table 2).

123

452

Q.-L. Men et al.

the three populations from the Hexi Corridor were at the same level as the population from the adjacent region of Xinjiang (Kum). Probability values testing for a significant heterozygote excess based on SMM and TPM models of microsatellite allele mutation using BOTTLENECK software are listed in Table 2. Results show no significant heterozygosity excess in any population. All P values remained insignificant (P \ 0.05) after Bonferroni correction. Genetic structure Fig. 2 The mean number of distinct alleles per locus for 14 sample populations. The lines represent the trend in the fluctuations of allelic richness accompanied by the change of sample size for each population

Among 12 Chinese populations, two populations from the Heilongjiang had the highest mean allelic richness (Mud, 10.38; Don, 12.23), followed by four populations from northern Xinjiang (Ili, 9.16; Uru, 8.30; Kuy, 8.08; Jin 6.82). Two populations from southern Xinjiang (Kor, 4.46; Kas, 5.43) and three populations from Gansu (Zha, 5.14; Dun, 5.50; Jiu, 4.17) had the relatively low mean allelic richness (Fig. 2). Based on the mean number of alleles per locus, the observed heterozygosity and the fluctuation of allelic richness, the populations from Heilongjiang (Mud and Don) had higher genetic diversity than populations from Xinjiang and Gansu. Populations from northern Xinjiang (Ili, Jin, Kuy and Uru) had similar genetic diversities, and the two populations from southern Xinjiang (Kas and Kor) were similar. Genetic diversities of populations from northern Xinjiang were higher than populations in southern Xinjiang. Meanwhile, population genetic diversities of

Fig. 3 Proportion of membership coefficient with K = 4 for 14 populations of codling moth obtained using STRUCTURE software

123

We performed cluster analysis using STRUCTURE version 2.3.3, which simultaneously identifies clusters and assigns individuals to populations using a Bayesian approach. Following the DK methods (Evanno et al. 2005), we obtained the most likely K value of four. All populations in this study were assigned to four clusters which were indicated by different colors and named as cluster 1, cluster 2, cluster 3 and cluster 4. The proportions of each population that contributed to each of the four clusters are shown in Fig. 3 and pie charts in Fig. 1b and c. Cluster 2 had a large contribution from three populations from northern Xinjiang in the northwest of China (Ili, 0.526; Jin, 0.869 and Kuy, 0.577). Cluster 4 mainly consisted of individuals from two population from northeastern China and the two populations from western Europe (Don, 0.861; Mud, 0.960; Swi, 0.901; Ger, 0.963). Cluster 3 primarily contained individuals from the two populations from Gansu (Dun, 0.674 and Jiu, 0.726) and one population from the boundary region between Gansu and Xinjiang (Kum, 0.540). Cluster 1 was mostly made-up of

Genetic structure and diversity of a newly invasive species

453

Table 3 Analysis of molecular variance of populations from Xinjiang, Gansu, Heilongjiang and other countries Source of variation

Among groups Among populations within groups Within populations

d.f.

Sum of squares

Variance components

Percentage of variation

P value

Fixation indices

3

155.756

0.23546Va

9.64

P \ 0.0001

FCT = 0.09644

10

165.522

0.30390Vb

12.45

P \ 0.0001

FSC = 0.13776

666

1,266.820

1.90213Vc

77.91

P \ 0.001

FST = 0.22091

samples from the two populations of southern Xinjiang (Kor, 0.760 and Kas, 0.655) and one population from Gansu (Zha, 0.805). The Urumqi population contributed almost equally to cluster 1 and cluster 3, with cluster 1 being 0.412 and cluster 3 being 0.394; this presents a unique genetic structure. AMOVA results (Table 3) indicated that the major portion of the molecular genetic variation was found within populations (77.91 %) with 9.64 % among the four groups and 12.45 % among the populations within groups. Exact tests showed significant genetic variance on all of these three levels. The values of pairwise FST ranged from 0.040 to 0.337 (Table 4, ‘‘Appendix’’). Sixty-eight comparisons out of ninety-one showed high or very high genetic differentiation. The pairwise FST values among Ili, Jin and Kuy populations from the north of Xinjiang were less than 0.10, while other pairwise FST values were above 0.10 for populations from Xinjiang. The pairwise FST values between the populations from the northeast (Heilongjiang) and northwest (Xinjiang and Gansu) of China ranged from 0.157 to 0.324, indicating high or very high genetic differentiation (Wright 1978). All P values remained significant (P \ 0.05) after Bonferroni correction. The Mantel test for 10 populations from northwestern China revealed no correlation between genetic distances and geographic distances (r = 0.164, P = 0.176). However, the test showed a strong isolation-by-distance effect present among populations (r = 0.399, P = 0.015) when we added the two populations from Heilongjiang.

Discussion By genotyping 12 populations from the main C. pomonella distribution regions in China, we have

documented successive losses of genetic diversity and a high level of genetic structure in these moth populations. Genetic variation among populations from different regions was significant, and we found no evidence for a bottleneck effect. Genetic diversity in the Cydia pomonella population from Ili was highest among the populations from northwestern China. The more eastern populations of Kum, Dun, Jiu and Zha, and the more southern populations of Kor and Kas showed relatively lower genetic diversity. Gradual loss of genetic diversity was observed in other insect species (Ciosi et al. 2008; Margaritopoulos et al. 2009; Watts et al. 2010). For invasive species, genetic diversity is expected to decrease with range expansion and colonization of new areas via sequential founder populations or bottleneck effects (Ramachandran et al. 2005; Herborg et al. 2007; Dlugosch and Parker 2008). However, the BOTTLENECK analysis revealed no clear genetic evidence for a bottleneck effect in any of the studied populations. Lack of a bottleneck effect may be attributed to the recent invasion history and successful colonisation by a large numbers of migrants (Cham 2002; Parr 2005). Comparison of genetic variations between old and newly-introduced populations can provide a useful clue for deducing the origins of an introduced population or earliest founder population (Dlugosch and Parker 2008). The high level of genetic diversity in the Ili population together with the strategic location of Ili which is situated close to Horgost, the important land port to central Asia, both suggest that the Ili region is most likely the origin of codling moth introduction into northwestern China. The Ili River valley, with its extensive apple orchards provides a suitable habitat for successful invasion of codling moth. However, two populations (Kas and Kor) from southern Xinjiang had low genetic diversity. This

123

454

result can be explained as populations expressing a fraction of the genetic variation that occurred in the source populations (Nei et al. 1975). A similar situation was also observed in populations from Kumul, Dunhuang, Jiuquan and Zhangye. Their genotypes are displaced and mixed through active or passive dispersal of individuals during the course of invasion. Selection of beneficial combinations related to local microclimatic and ecological conditions then improves local fitness and leads to new infestations (Lee 2002). As a result of multiple introduction events from distinct genetic sources, the newly-introduced populations often have higher genetic diversity than native or older populations (Novak and Mack 1993; Roman and Darling 2007; Dlugosch and Parker 2008; Brown and Stepien 2009; Gillis et al. 2009; Ghabooli et al. 2011). Since Cydia pomonella was reported in the northeast of China only 6 years ago, the two populations from Heilongjiang showing the highest genetic diversities among the studied populations could be explained by the new invasion event drawing from different genetic resources. Heilongjiang is around 3,000 km and 5,000 km away from Gansu and Xinjiang, respectively. In spite of the considerable distance from Heilongjiang to Xinjiang and Gansu, no occurrence of this pest has been recorded in the large area between them despite continuous monitoring. Taking together the higher genetic diversity and the distinct genetic structure revealed by the cluster analysis, we deduce that C. pomonella were not established by populations invading from western China. However, the Urumqi moths, as a relatively older population in China, showed a high level of genetic diversity and distinct structure. As the capital of Xinjiang, Urumqi could provide a corridor for accidental human-aided dispersal of C. pomonella, resulting in the intermixing of populations from northern and southern Xinjiang after years of successful invasion. The codling moth populations from Germany and Switzerland displayed a higher level of genetic diversity compared with populations from northwestern China. This aligns with the codling moth having a relatively shorter history in China than in southeastern Europe, the moth’s documented origin (Shel‘Deshova 1967). There was strong genetic structuring between populations based on the Bayesian clustering approach, AMOVA and pairwise FST. Previous studies

123

Q.-L. Men et al.

of genetic variation showed that Cydia pomonella populations exhibited distinct characteristics of genetic structure in different habits, and that this is governed by a number of factors. Fuentes-Contreras et al. (2008) reported the genetic variation in central Chile and concluded that the low level of genetic differentiation observed was due to the adjacent applegrowing areas, the domestic culturing of fruits, the recent spread of modern fruit culture and the recent introduction of C. pomonella. Chen and Dorn (2009) analyzed genetic variation of codling moth populations in Switzerland and found that host specificity, geographic isolation, intrinsic flight capacity and anthropogenic measures could shape the population structure. Franck and Timm (2010) concluded that the population genetic structure of codling moths is primarily affected by hosts, geography and time. In the present study, the results of the Mantel test indicated that the pattern of variation was weakly associated with geographic distance among populations from northwestern China. Considering that the history of the codling moth in China spans only a little more than 50 years, this pattern may be the result of non-equilibrium from migration and drift, and due to the recent range expansion associated with human intervention (Slatkin 1993). The population from Zhangye had a very similar genetic structure to the two populations from southern Xinjiang. As the samples of Zhangye were collected in an apple and pear breeding center, where the famous Korla pear, Pyrus sinkiangensis Yu, has been frequently exported from Korla over the last 10 years, it is very possible that the similarity in genetic structure is due to human activity associated with the movement of nursery material. Human intervention in the form of fruit and seedling transport as well as the movement of bins has played an important role in the mixing of populations from distant geographic regions (Higbee et al. 2001). Similar circumstances were conjectured in northern Italy, Grapholita molesta populations (Torriani et al. 2010) and for codling moth populations from South Africa (Timm et al. 2006), which were genetically little differentiated despite the large geographic distance ([150 km) between some of the sampled orchards. Furthermore, C. pomonella has a limited dispersal ability, which would be expected to weaken IBD relationship over wide regions (Roff 1994). C. pomonella is generally regarded to be a sedentary species, with a sphere of activity extending over 0.5 ha

Genetic structure and diversity of a newly invasive species

(Mani and Wildbolz 1977). The weak flight capacity of the codling moth influences its dispersal and gene flow pattern (Timm et al. 2006; Chen and Dorn 2009). Significant differentiation between codling moth populations occurred across small geographic distances (&1 km) in South Africa (Timm et al. 2006) as well as in Switzerland (\10 km) (Chen and Dorn 2009). Moreover, the distribution areas of codling moth in China occur across very complex topography (huge mountains, large unpopulated areas or large deserts), different climates (temperate, continental arid or frigid climate) and different agricultural landscapes. A high level of population differentiation results from the high potential for adaptation to different environmental conditions (Carde´ and Minks 1995) and high reproductive rates (Myers et al. 2006), which allow these moths to form locally differentiated populations. Their weak flight capacity and humanaided dispersal increase this differentiation. The Bayesian clustering analysis indicated that populations from Ili, Jinghe and Kuytun had a similar proportion of membership coefficients and contributed mainly to cluster 3. Pairwise FST also showed relatively low values in these three populations. This result linking to relatively higher population genetic diversity estimates for the three populations suggests a novel finding where the Ili, Jinghe and Kuytun populations were established by individuals belonging to the same maternal lineage. In northern Xinjiang, two factors may have contributed to this result: the relative proximity of the fruit production regions and the frequent commodities exchange among these locations. This corresponds with indirect evidence obtained in South Africa using AFLPs and finding limited gene flow among C. pomonella populations among local regions due to isolated fruit production areas (Timm et al. 2006). The level of genetic differentiation between populations from northeastern and northwestern China was great. The two populations from Dongning and Mudanjiang contributed mainly to the same cluster 3. There was also strong isolation based on distance. These results add support to the view that the invasion

455

of the codling moth into China came from at least two different pathways. We also found the populations from Heilongjiang had a similar genetic structure with those of western European populations. A phylogenetic study based on a greater number of collection sites is necessary to clarify the relationship among Cydia pomonella populations from European countries before precise conclusions can be drawn, especially the phylogenetic relationships of populations from western Europe and the far eastern region of Russia. The latter is near Heilongjiang and has been documented as a distribution region of codling moth (Willett et al. 2009). In China, management practices against Cydia pomonella have primarily focused on chemical control. Investigation of genetic diversity of C. pomonella populations can provide us a guide for controlling this pest. For example, high inbreeding and low genetic diversity leads to application of pesticides ultimately generating genetic bottlenecks. Furthermore, localized populations with similar genetic structure should be considered the same management unit for most effective control (Ayres et al. 2010). For isolated populations, we should try various management methods, especially a variety of different chemical pesticides. In many cases, immigrating individuals from isolated orchards may result in the sequential introduction of novel resistance genes and promote pesticide resistance in their gene pool (Bouvier et al. 2001). Additional research needs to be done to detail the geographic origins of codling moth in China and its spread through China. Acknowledgments We thank Dr. John Richard Schrock, Emporia State University, Emporia, KS, USA, for language correction of the draft. We are grateful to Anyong Wang, Kun Dong, Chunhan Zheng, Jianqiang Yang, Xiao Zhao and Xinglong Huang for their help in the collection of codling moth samples. This work was supported by grants from the ‘13115’ Sci-Tech Innovation Project of Shaanxi Province (No. 2009ZDKG-06), the Special Fund for Agro-scientific Research in the Public Interest (No. 200903042-03), the National Natural Science Foundation of China (No. 31071687), the International Atomic Energy Agency (No. 16341) and the Talent Recruitment Project of Northwest A & F University.

123

456

Q.-L. Men et al.

Appendix See Table 4. Table 4 Pairwise population differentiation estimates (FST) averaging eight loci between all populations of C. pomonella (below the diagonal) and P values (above diagonal) Zha Zha

Jiu

Dun

Uru

Kum

Kas

Jin

Kor

Ili

Kuy

Mud

Don

Ger

Swi

**

**

**

**

**

**

**

**

**

**

**

**

**

**

**

**

**

**

**

**

**

**

**

**

**

Jiu

0.182

Dun

0.100

0.195

Uru

0.116

0.155

0.146

Kum Kas

0.152 0.144

0.121 0.272

0.142 0.236

0.173 0.205

0.166

Jin

0.244

0.252

0.274

0.209

0.133

0.246

Kor

0.096

0.337

0.306

0.223

0.323

0.332

0.278

Ili

0.154

0.185

0.120

0.124

0.110

0.105

0.083

0.203

Kuy

0.196

0.105

0.199

0.159

0.112

0.210

0.098

0.277

0.080

Mud

0.292

0.303

0.294

0.194

0.316

0.324

0.244

0.231

0.187

0.217

Don

0.277

0.264

0.288

0.204

0.268

0.301

0.196

0.241

0.157

0.161

0.112

Ger

0.264

0.298

0.282

0.154

0.300

0.301

0.214

0.220

0.158

0.198

0.060

0.088

Swi

0.266

0.264

0.266

0.159

0.281

0.289

0.222

0.223

0.157

0.176

0.040

0.110

**

**

**

**

**

**

**

**

**

**

**

**

**

**

**

**

**

**

**

**

**

**

** **

** **

** **

** **

** **

** **

** **

** **

**

**

**

**

**

**

**

**

**

**

**

**

**

**

**

**

**

**

**

**

**

**

**

**

**

**

** *

0.043

* P \ 0.05, ** P \ 0.01; The significances were tested for multi comparisons by the Boferroni method for k = 91, P \ 0.05

References Ahern RG, Hawthorne DJ, Raupp MJ (2009) Founder effects and phenotypic variation in Adelges cooleyi, an insect pest introduced to the eastern United States. Biol Invasions 11:959–971 Ayres RM, Pettigrove VJ, Hoffmann AA (2010) Low diversity and high levels of population genetic structuring in introduced eastern mosquitofish (Gambusia holbrooki) in the greater Melbourne area, Australia. Biol Invasions 12: 3727–3744 Barnes MM (1991) Codling moth occurrence, host race formation and damage. In: Van der Guest LPS, Evenhuis HH (eds) Tortricid pests: their biology, natural enemies and control. Elsevier, Amsterdam, pp 313–328 Bohonak AJ (1999) Dispersal, gene flow and population structure. Q Rev of Biol 74:21–45 Boivin T, Bouvier J-C, Beslay D, Sauphanor B (2004) Variability in diapause propensity within populations of a temperate insect species: interactions between insecticide resistance genes and photoperiodism. Biol J Linn Soc 83:341–351 Bonnet E, Van der Peer Y (2002) ZT: a software tool for simple and partial Mantel tests. J Stat Software 7:1–12 Bouvier J-C, Bue`s R, Boivin T, Boudinhon L, Beslay D, Sauphanor B (2001) Deltamethrin resistance in the codling moth (Lepidoptera: Tortricidae): inheritance and number of genes involved. Heredity 82:456–462

123

Brown JE, Stepien CA (2009) Invasion genetics of the Eurasian round goby in North America: tracing source and spread patterns. Mol Eco 18:64–79 Carde´ RT, Minks AK (1995) Control of moth pests by mating disruption: successes and constraints. Ann Rev Entomol 40:559–585 Cham S (2002) Range expansion of the small red-eyed damselfly Erythromma viridulum (Charp.) in the British Isles. Atropos 15:3–9 Chen MH, Dorn S (2009) Microsatellites reveal genetic differentiation among populations in an insect species with high genetic variability in dispersal, the codling moth, Cydia pomonella (L.) (Lepidoptera: Tortricidae). Bull Entomol Res 100:75–85 Ciosi M, Miller NJ, Kim KS, Giordano R, Estoup A, Guillemaud T (2008) Invasion of Europe by the western corn rootworm, Diabrotica virgifera virgifera: multiple transatlantic introductions with various reductions of genetic diversity. Mol Ecol 17:3614–3627 Cornuet JM, Luikart G (1996) Description and power analysis of two tests for detecting recent population bottlenecks from allele frequency data. Genetics 144:2001–2014 Cristescu R, Sherwin WB, Handasyde K, Cahill V, Cooper DW (2010) Detecting bottlenecks using BOTTLENECK 1.2.02 in wild populations: the importance of the microsatellite structure. Conserv Genet 11:1043–1049 Dakin EE, Avise JC (2004) Microsatellite null alleles in parentage analysis. Heredity 93:504–509

Genetic structure and diversity of a newly invasive species De Barro PJ (2005) Genetic structure of the whitefly Bemisia tabaci in the Asia-Pacific region revealed using microsatellite markers. Mol Ecol 14:3695–3718 Di Rienzo A, Peterson AC, Garza JC, Valdes AM, Slatkin M, Freimer NB (1994) Mutational processes of simplesequence repeat loci in human populations. Proc Nati Acad Sci USA 91:3166–3170 Dieringer D, Schlo¨tterer C (2003) MICROSATELLITE ANALYSER (MSA): a platform independent analysis tool for large microsatellite data sets. Mol Ecol Notes 3:167–169 Dlugosch KM, Parker IM (2008) Founding events in species invasions: genetic variation, adaptive evolution, and the role of multiple introductions. Mol Ecol 17:431–449 Dorn S, Schumacher P, Abivardi C, Meyho¨fer R (1999) Global and regional pest insects and their antagonists in orchards: spatial dynamics. Agric Ecosyst Environ 73:111–118 Espinoza JL, Fuentes-Contreras E, Barros W, Ramı´rez C (2007) Utilizacio´n de microsate´lites para la determinacio´n de la variabilidad gene´tica de la polilla de la manzana Cydia pomonella L. (Lepidoptera: Tortricidae) en Chile Central. Agric Te´c (Chile) 67:244–252 Evanno G, Regnaut S, Goudet J (2005) Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol Ecol 14:2611–2620 Excoffier L, Laval G, Schneider S (2005) Arlequin ver. 3.0: an integrated software package for population genetics data analysis. Evol Bioinform Online 1:47–50 Franck P, Timm AE (2010) Population genetic structure of Cydia pomonella: a review and case study comparing spatiotemporal variation. J Appl Entomol 134:191–200 Franck P, Gue´rin F, Loiseau A, Sauphanor B (2005) Isolation and characterization of microsatellite loci in the codling moth Cydia pomonella L. (Lepidoptera, Tortricidae). Mol Ecol Notes 5:99–102 Franck P, Reyes M, Olivares J, Sauphanor B (2007) Genetic architecture in codling moth populations: comparison between microsatellite and insecticide resistance markers. Mol Ecol 16:3554–3564 Fuentes-Contreras E, Espinoza JL, Lavandero B, Ramı´rez CC (2008) Population genetic structure of codling moth (Lepidoptera: Tortricidae) from apple orchards in central Chile. J Econ Entomol 101:190–198 Ghabooli S, Shiganova TA, Zhan AB, Cristescu ME, EghtesadiAraghi P, MacIsaac HJ (2011) Multiple introductions and invasion pathways for the invasive ctenophore Mnemiopsis leidyi in Eurasia. Biol Invasions 13:679–690 Gillis NK, Walters LJ, Fernandes FC, Hoffman EA (2009) Higher genetic diversity in introduced than in native populations of the mussel Mytella charruana: evidence of population admixture at introduction sites. Divers Distrib 15:784–795 Grapputo A, Boman S, Lindstrom L, Lyytinen A, Mappes J (2005) The voyage of an invasive species across continents: genetic diversity of North American and European Colorado potato beetle populations. Mol Ecol 14:4207– 4219 Herborg L-M, Weetman D, van Oosterhout C, Hanfling B (2007) Genetic population structure and contemporary dispersal patterns of a recent invader, the Chinese mitten crab, Eriocheir sinensis. Mol Ecol 16:231–242

457 Higbee BS, Calkins CO, Temple CA (2001) Overwintering of codling moth (Lepidoptera: Tortricidae) larvae in apple harvest bins and subsequent moth emergence. J Econ Entomol 94:1511–1517 Hufbauer RA, Bogdanowicz SM, Harrison RG (2004) The population genetics of a biological control introduction: mitochondrial DNA and microsatellite variation in native and introduced populations of Aphidus ervis, a parasitoid wasp. Mol Ecol 13:337–348 Lee CE (2002) Evolutionary genetics of invasive species. Trends Ecol Evol 17:386–391 Luikart G, Sherwin WB, Steele BM, Allendorf FW (1998) Usefulness of molecular markers for detecting population bottlenecks via monitoring genetic change. Mol Ecol 7:963–974 Mani E, Wildbolz T (1977) The dispersal of male codling moths (Laspeyresia pomonella L.) in the Upper Rhine Valley. J Appl Entomol 47:39–48 Mantel N (1967) The detection of disease clustering and a generalized regression approach. Canc Res 27:209–220 Margaritopoulos JT, Kasprowicz L, Malloch GL, Fenton B (2009) Tracking the global dispersal of a cosmopolitan insect pest, the peach potato aphid. BMC Ecol 9:13 Megle´cz E, Petenian F, Danchin E, Coeur D’Acier A, Rasplus J-Y, Faure E (2004) High similarity between flanking regions of different microsatellites detected within each of two species of Lepidoptera: Parnassius apollo and Euphydryas aurinia. Mol Ecol 13:1693–1700 Meng XF, Shi M, Chen XX (2008) Population genetic structure of Chilo suppressalis (Walker) (Lepidoptera: Crambidae): strong subdivision in China inferred from microsatellite markers and mtDNA gene sequences. Mol Ecol 17:2880– 2897 Meraner A, Brandsta¨tter A, Thaler R, Aray B, Unterlechner M, Niedersta¨tter H, Parson W, Zelger R, Dalla Via J, Dallinger R (2008) Molecular phylogeny and population structure of the codling moth (Cydia pomonella) in Central Europe: I. Ancient clade splitting revealed by mitochondrial haplotype markers. Mol Phylogenet Evol 48:825–837 Miller NJ, Birley AJ, Overall ADJ, Tatchell GM (2003) Population genetic structure of the lettuce root aphid, Pemphigus bursarius (L.), in relation to geographic distance, gene flow and host plant usage. Heredity 91:217–223 Mozaffarian F, Mardi M, Sarafrazi A, Ganbalani N (2007) Assessment of geographic and host-associated population variations of the carob moth, Ectomyelois ceratoniae, on pomegranate, fig, pistachio and walnut, using AFLP markers. J Insect Sci 8:6 Myers CT, Hull LA, Krawczyk G (2006) Comparative survival rates of oriental fruit moth (Lepidoptera: Tortricidae) larvae on shoots and fruit of apple and peach. J Econ Entomol 99:1299–1309 Nei M, Maruyama T, Chakraborty R (1975) The bottleneck effect and genetic variability in populations. Evolution 29:1–10 Novak SJ, Mack RN (1993) Genetic variation in Bromus tectorum (Poaceae): comparison between native and introduced populations. Heredity 71:167–176 Orsini L, Corander J, Alasentie A, Hanski I (2008) Genetic spatial structure in a butterfly metapopulation correlates

123

458 better with past than present demographic structure. Mol Ecol 17:2629–2642 Parr AJ (2005) Migrant and dispersive dragonflies in Britain during 2004. J Br Dragonfly Soc 21:14–20 Paupy C, Chantha N, Reynes J-M, Failloux A-B (2005) Factors influencing the population structure of Aedes aegypti from the main cities in Cambodia. Heredity 95:144–147 Piry S, Luikart G, Cornuet JM (1999) BOTTLENECK: a computer program for detecting recent reductions in effective population size from allele frequency data. J Hered 90:502–503 Primmer CR, Saino N, Møller AP, Ellegren H (1998) Unraveling the processes of microsatellite evolution through analysis of germ line mutations in barn swallows, Hirundo rustica. Mol Biol Evol 15:1047–1054 Pritchard JK, Stephens M, Donnelly P (2000) Inference of population structure using multilocus genotype data. Genetics 155:945–959 Pruett CI, Winker K (2005) Northwestern song sparrow populations show genetic effects of sequential colonization. Mol Ecol 14:1421–1434 Ramachandran S, Deshpande O, Roseman CC, Rosenberg NA, Feldman MW, Cavalli-Sforza LL (2005) Support from the relationship of genetic and geographic distance in human populations for a serial founder effect originating in Africa. Proc Natl Acad Sci USA 102:15942–15947 Ramstad KM, Woody CA, Sage GK, Allendorf FW (2004) Founding events influence genetic population structure of sockeye salmon (Oncorhynchus nerka) in Lake Clark, Alaska. Mol Ecol 13:277–290 Rice WR (1989) Analyzing tables of statistical tests. Evolution 43:223–225 Roff DA (1994) The evolution of flightlessness: is history important? Evol Ecol 8:639–657 Roman J, Darling JA (2007) Paradox lost: genetic diversity and success of aquatic invasions. Trends Ecol Evol 22:454–464 Rosenberg NA (2004) DISTRUCT: a program for the graphical display of population structure. Ecol Evol Notes 4:137–138 Rousset F (2008) GENEPOP’007: a complete re-implementation of the GENEPOP software for Windows and Linux. Mol Ecol Resour 8:103–106

123

Q.-L. Men et al. Schuelke M (2000) An economic method for the fluorescent labeling of PCR fragments. Nat Biotechnol 18:233–234 Shel‘Deshova GG (1967) Ecological factors determining distribution of the codling moth Laspeyresia pomonella (Lepidoptera: Tortricidae) in the northern and southern hemispheres. Entomol Rev 46:349–361 Slatkin M (1993) Isolation by distance in equilibrium and nonequilibrium populations. Evolution 47:264–279 Szpiech ZA, Jakobsson M, Rosenberg NA (2008) ADZE: a rarefaction approach for counting alleles private to combinations of populations. Bioinformatics 24:2498–2504 Timm AE, Geertsema H, Warnich L (2006) Gene flow among Cydia pomonella (Lepidoptera: Tortricidae) geographic and host populations in South Africa. J Econ Entomol 99:341–348 Torriani MVG, Mazzi D, Hein S, Dorn S (2010) Structured populations of the oriental fruit moth in an agricultural ecosystem. Mol Ecol 19:2651–2660 Van Oosterhout C, Hutchinson WF, Wills DPM, Shipley P (2004) MICRO-CHECKER: software for identifying and correcting genotyping errors in microsatellite data. Mol Ecol Notes 4:535–538 Wan FH, Guo JY, Zhang F (2009) Research on biological invasions in China. Science Press, Beijing Watts PC, Keat S, Thompson DJ (2010) Patterns of spatial genetic structure and diversity at the onset of a rapid range expansion: colonization of the UK by the small red-eyed damselfly Erythromma viridulum. Biol Invasions 12: 3887–3903 Weir BS, Cockerham CC (1984) Estimating F-statistics for the analysis of population-structure. Evolution 38:1358–1370 Willett MJ, Neven L, Miller CE (2009) The occurrence of codling moth in low latitude countries: validation of pest distribution reports. Hort Technology 19:633–637 Wright S (1951) The genetical structure of populations. Ann Eugen 15:323–354 Wright S (1978) Evolution and genetics of populations. University of Chicago, Chicago Zhang XZ (1957) Taxonomic notes on the codling moth, Carpocapsa pomonella L. in Sinkiang. Acta Entomol Sin 7:467–472