Molecular Characterization of Codling Moth Cydia pomonella (Linnaeus) (Lepidoptera: Tortricidae) in Swat Valley Pakistan using Randomly Amplified Polymorphic DNA (RAPD) Polymerase Chain Reaction

Hayat Zada1, Ahmad-Ur-Rahman Saljoqi1,*, Ijaz Ali2, Bashir Ahmad1, Abdul Waheed Khan2 and Shabeer Ahmad2

1Department of Plant Protection, The University of Agriculture, Peshawar

2Institute of Biotechnology and Genetic Engineering, The University of Agriculture, Peshawar

ABSTRACT

Codling moth Cydia pomonella (Linnaeus) (Lepidoptera: Tortricidae) is the most important insect pest of apple production in Pakistan. Genomic DNA was extracted from 30 overwintering larvae of each population of Matta, Madyan and Kalam regions. Out of 30 tested primers, 21 amplified 157 polymorphic bands in three populations. The mean gene frequency (f), genetic distance (I) and Shannon’s information index (h) for three populations were 1.336, 0.304 and 0.437, respectively. Nei’s unbiased measures of genetic identity and genetic distance revealed that higher genetic distance was observed among the isolates from Kalam and Madyan (97.87 %) whereas low genetic distance (35.58%) was calculated for C. pomonella isolates from Matta and Madyan. These results revealed that molecular variations among the population of C. pomonella might be attributed to the climatic conditions, geographical location and pest control programs.

Article Information

Received 12 June 2016

Revised 22 March 2018

Accepted 24 July 2018

Available online 28 May 2019

Authors’ Contribution

HZ designed the experiments, collected data. AURS wrote the article. IA provided RAPD primers and helped in lab work. BA helped in sample collection and review of literature. AWK helped in lab work. SA helped in data analysis.

Key words

Cydia pomonella, Molecular characterization, RAPD primers, PCR.

DOI: http://dx.doi.org/10.17582/journal.pjz/2019.51.4.1547.1554

* Corresponding author: [email protected]

0030-9923/2019/0004-1547 $ 9.00/0

Copyright 2019 Zoological Society of Pakistan

Introduction

Codling moth (Cydia pomonella) is a serious insect pest of apple fruit in different parts of the world including Swat Valley in the North West of Khyber Pakhtunkhwa, Pakistan. A high fruit damage (80%) has been observed due to the insect in temperate parts of all major continents (Barnes, 1991). While being an economically important pest of the apple worldwide, the growers have a low tolerance (<1%) for its injury (Dorn et al., 1999). This has necessitated the adoption of different control strategies especially frequent applications of broad-spectrum synthetic insecticides throughout the fruiting period. Besides apple, this pest attacks pear, walnut, quince and some stone fruits causing economic losses in fruit production (Ciglar, 1998). The pest was originated in Eurasia, since last two centuries it dispersed around the world with the spread of the cultivation of apples and pears (Franck et al., 2007). It has achieved a nearly global distribution, being one of the most successful pest insect species known today (Thaler et al., 2008). Currently, Cydia pomonella is present in North America, South America, South Africa, Australia and New Zealand, India, Pakistan and Afghanistan (Franck et al., 2007).

Despite its economic importance, little is known about the genetic structure and patterns of gene flow at the local and regional scale, which are important aspects for establishing an area wide control strategy (Dorn et al., 1999; Calkins and Faust, 2003). Pest management strategies of C. pomonella include regular insecticide treatments, which are known to select for resistance to several insecticide groups (Knight et al., 1994; Sauphanor et al., 1998; Dunley and Welter, 2000). Therefore, a deeper understanding of the population variation of this pest could greatly benefit management decisions for its control (Timm et al., 2006).

It has adapted successfully to different habitats by forming various biotypes and populations, often designated as strains, which differ among each other in several morphological, developmental and physiological features (Meraner et al., 2008). Knowledge of the genetic variation within C. pomonella populations is necessary for their efficient control and management. Molecular technologies provide new ways to study population diversity as well as to differentiate closely related species (Williams et al., 1990; Deverno et al., 1998).

Polymerase chain reaction (PCR) strategies offer increased sensitivity and speed for the molecular characterization and variation among the populations of insects species (Deverno et al., 1998). Different molecular techniques have been employed, randomly amplified polymorphic DNA (RAPD) is one of them because it is efficient and easy to use for the molecular variation among the populations. This technique has been developed to detect genetic variability by PCR amplification of arbitrary segments of genomic DNA by using short, random primers and thus does not require prior knowledge of a DNA sequence. Its low cost, its efficiency in developing a large number of DNA markers in a short time and the less sophisticated equipment that it requires has made RAPD a valuable technique (Bardakci, 2001; Delaat et al., 2005). RAPD markers are very well suited for use in insect phylogeny areas like the detection of genetic variability among populations as well as the identification of closely related species (Benecke, 1998; Lima et al., 2002).

Since limited studies have been reported regarding the characterization of C. pomonella in this area, the current study was undertaken to ascertain the molecular variation among the population of C. pomonella.

 

Materials and methods

Study area and parameters

Female larvae (In the last instar male larvae well-developed purple-red testes are seen through the skin of the fifth abdominal segment, contrary to the female larvae) trapped and collected by single face cardboard fastened around the apple tree at a distance of 30 cm from the ground following the procedures of Fritsch et al. (2005) with some modifications. In each population, 30 overwintering female larvae were randomly selected for DNA isolation to minimize DNA contamination by endoparasites (Landry et al., 1999) in autumn from Matta (350 55’ 19.11” of latitude North and 720 30’ 37.52” of longitude East) (920.30 meters), Madyan (350 08’ 0.00” of latitude North and 720 32’ 0.00” of longitude East) (1333.84 meters) and Kalam (350 28’ 41.66” of latitude North and 720 34’ 18.61” of longitude East) (2092.30 meters). To eliminate the effect of host association in discrimination of populations, all of the specimens were collected from “Red Delicious” apple orchards. The specimens were washed and stored in 96% alcohol (Ethanol) prior to analysis in the Health laboratory of Institute of Bio-Technology and Genetic Engineering (IBGE), The University of Agriculture, Peshawar.

Genomic DNA extraction

For genomic DNA extraction, Spinklean Genomic DNA Extraction Kit (Thermoscientific® USA) was used and Zimmerman et al. (2000) procedures was followed with some necessary modifications. The specimens were crushed in liquid nitrogen. DNA was extracted using manufacture’s manual. Briefly, TL Buffer (250ul) was added to crushed larvae for tissue lysis. The samples were then vortexed to crash properly mix. Total of 199 units of enzyme was added and mixed thoroughly. Lysis buffer of 220 µl was added to the solutions and mixed thoroughly. A volume of 560 µl Buffer TB was then added to each eppendorf tube and mixed comprehensively through vortex to get a homogeneous solution followed by incubation for 10 minutes at 65oC. Then 200 µl absolute ethanol was added. Sample mixtures were passed through column, assembled in clean collection tube by centrifuging at 8000 xg for 1 min. Column was washed twice with 750 µl wash buffer ‘PS’ and centrifuged at 8000 xg for 1 min. For removing traces of ethanol the column was centrifuged again at 10000 xg for another one minute.

About 200 µl of pre heated TE buffer was added to the column membrane in new tubes, incubated at room temperature for 2 min and centrifuged at 10000xg for 1min to get DNA which was then stored at -20°C.

Polymerase chain reaction

The polymerase chain reaction (PCR) was carried out in a 25 µl volume containing genomic DNA (150 ng), 0.25 mM of RAPD primers (Genlink, USA), 200 µM of each dNTP, 50 mM of KCL, 10 mM Tris, 1.5 mM MgCl2 and 2.5 unit of Taq polymerase (Thermoscientific®, USA). The amplification was carried out in 40 cycles of 50 second at 94°C, 1 min at 28°C, 1 min at 72°C, preceded by an initial denaturation at 94°C for 4 min and final extension at 72°C for 10 min at the end in a gene amp PCR system 2700 programmable thermo cycler. The amplification products were then detected on 2% agarose gel pre-stained with ethidium bromide.

Statistical analysis

For statistical analysis of randomly amplified polymorphic DNA (RAPD), every scorable band was considered as a single locus/allele. The loci were scored as present (1) or absent (0). Bivariate 1-0 matrix was generated. Genetic distances was calculated using “Unweighted Pair Group of Arithmetic Means” (UPGMA) procedure described by Nei and Lie (1979):

GD = 1 - dxy / dx + dy - dxy

Where, GD is genetic distance between two genotypes, dxy is total No. of common loci (bands) in two genotypes, dx is total No. of loci (bands) in genotype 1 and dy is total No. of loci (bands) in genotypes.

 

Table I.- Name, sequence, size and molecular weight of RAPD primer used for molecular characterization of C. pomonella.

S. No	Primer	Sequence	Size (bp)	Tm	M. wt	% GC
1.	GL Decamer B-12	CCTTGACGCA	10	29.5°C	2987.98	60
2.	GL Decamer D-16	AGGGCGTAAG	10	29.5°C	3117.04	60
3.	GL Decamer C-04	CCGCATCTAC	10	29.5°C	2947.96	60
4.	GL Decamer C-13	AAGCCTCGTC	10	29.5°C	2987.98	60
5.	GL Decamer B-04	GGACTGGAGT	10	29.5°C	3108.04	60
6.	GL Decamer H-02	TGTAGCTGGG	10	29.5°C	3099.04	60
7.	GL Decamer E-09	CTTCACCCGA	10	29.5°C	2947.96	60
8.	GL Decamer F-01	ACGGATCCTG	10	29.5°C	3028.00	60
9.	GL Decamer A-19	CAAACGTCGG	10	29.5°C	3037.00	60
10.	GL Decamer D-08	GTGTGCCCCA	10	33.6°C	3003.99	70
11.	GL Decamer G-11	TCCCCGTCGT	10	33.6°C	2994.99	70
12.	GL Decamer F-07	CCGATATCCC	10	29.5°C	2947.96	60
13.	GL Decamer E-18	GGACTGCAGA	10	29.5°C	3077.02	60
14.	GL Decamer H-13	GACGCCACAC	10	33.6°C	2981.97	70
15.	GL Decamer B-15	GGAGGGTGTT	10	29.5°C	3139.06	60
16.	GL Decamer C-16	CACACTCCAG	10	29.5°C	2956.96	60
17.	GL Decamer C-02	GTGAGGCGTC	10	33.6°C	3084.03	70
18.	GL Decamer H-03	AGACGTCCAC	10	29.5°C	2996.98	60
19.	GL Decamer F-04	GGTGATCAGG	10	29.5°C	3108.04	60
20.	GL Decamer H-13	ACCAGGTTGG	10	29.5°C	3068.02	60
21.	GL Decamer G-02	GGCACTGAGG	10	33.6°C	3093.03	70
22.	GL Decamer* A-06	GGTCCCTGAC	10	33.6°C	3003.99	70
23.	GL Decamer* A-07	GAAACGGGTG	10	29.5°C	3117.04	60
24.	GL Decamer* B-16	TTTGCCCGGA	10	29.5°C	3019.00	60
25.	GL Decamer* D-10	GGTCTACACC	10	29.5°C	2987.98	60
26.	GL Decamer* F-11	TTGGTACCCC	10	29.5°C	2978.98	60
27.	GL Decamer* G-13	CTCTCCGCCA	10	33.6°C	2923.95	70
28.	GL Decamer* G-15	ACTGGGACTC	10	29.5°C	3028.00	60
29.	GL Decamer* H-05	AGTCGTCCCC	10	33.6°C	2963.97	70
30.	GL Decamer* H-10	CCTACGTCAG	10	29.5°C	2987.98	60

* Indicates the RAPD markers giving no results (Bands). With result (Bands), 21; Without result (Bands), 09.

 

Results

RAPD analysis

The results pertaining to gene frequency of the RAPD primers disclosed that out of 30 RAPD primers (Tables I and II), 21 primers gave the banding pattern as depicted in Table II. Highest gene frequency on the basis of amplification pattern was observed for RAPD primer GLH-13 (3.33) followed by GLH-14 (2.77), GLE-18 (2.66), GLG-11 (2.08), GLF-07 (1.85), GLH-03 and GLB-04 (1.66) each, GLC-16 (1.55), GLC-06 (1.55), GLC-02 (1.46), GLE-09 (1.41), GLF-04 (1.42) and GLD -08 (1.03), whilst the rest of RAPD primers depict the allele frequency below 0.99. The overall mean of the genetic frequency among three population of C. pomonella was 1.33.

On the basis of amplification pattern, the genetic diversity of the three population of the C. pomonella expounded that highest genetic diversity was recorded for the RAPD primers GLC-02, GLE-18 and GLC-16 (0.44 each), followed by GLC-13 (0.38), GLH-14 (0.37), GLA-19 (0.35), GLB-15 (0.34), GLB-12 (0.33) and GLB-04 (0.32), whilst the rest of the RAPD primers explicated below 0.29 genetic diversity among three population of the C. pomonella. The overall mean of the genetic diversity of the three population was 0.29.

The results pertaining to Shannon’s information index for each allele of the primers revealed that maximum Shannon information index was observed for GLC-04, GLE-18 and GLC-16 (0.63), followed by GLC-13 (0.55), GLH-13 (0.53), GLH-0.53) and GLA-19 (0.50), whilst the rest of the primers were less than 0.49 Shannon’s information index values. The overall mean of the Shannon’s information index was 0.44.

 

Table II.- Mean gene frequency, diversity and Shannon information index for RAPD primers used for molecular characterization of C. pomonella at Swat during the year 2012-13.

S. No	RAPD primers used	Total alleles amplified	Range of allele size (bp)	G. F (f)1	G. D (I)2	S. I. I. (h)3
1.	B-12	08	250-10000	0.4166	0.3333	0.4773
2.	D-16	08	250-10000	0.2500	0.2222	0.3182
3.	C-04	06	750-10000	0.5555	0.4444	0.6365
4.	C-13	08	250-10000	0.9166	0.3888	0.5569
5.	B-04	11	250-4500	1.6666	0.3232	0.4629
6.	H-02	10	250-4500	0.9999	0.1333	0.1909
7.	E-09	08	250-10000	1.4166	0.2777	0.3978
8.	F-01	06	250-10000	0.1111	0.1481	0.2121
9.	A-19	10	250-10000	0.8000	0.3555	0.5092
10.	D-08	09	250-3500	1.0370	0.2962	0.4243
11.	G-11	08	250-3500	2.0833	0.2777	0.3978
12.	F-07	09	250-4000	1.8518	0.2468	0.3978
13.	E-18	08	250-3000	2.6666	0.4444	0.6365
14.	H-14	06	250-3000	2.7777	0.3707	0.5304
15.	B-15	09	250-3000	0.5185	0.3456	0.4950
16.	C-16	03	750-1500	1.5555	0.4444	0.6365
17.	C-02	05	250-1500	1.4666	0.2666	0.3819
18.	H-03	04	500-1500	1.6666	0.2222	0.3182
19.	F-04	07	250-2500	1.4284	0.1904	0.2727
20.	H-13	06	250-2000	3.3333	0.4444	0.5304
21.	G-02	08	250-3000	0.5416	0.2222	0.3978
Mean		07	250-10000	1.33618	0.30467	0.4372

1, gene frequency (f); 2, genetic distance (I); 3, Shannon information index (h).

 

Nei’s unbiased measures of genetic identity and genetic distance

Table III, pertaining to the Nei’s genetic identity (above diagonal) and genetic distance (below diagonal) among the three population of C. pomonella. Higher genetic distance was observed among the isolates from Kalam and Madyan (97.87 %) whereas low genetic distance (35.58%) was calculated from the C. pomonella isolates from Matta and Madyan, which indicates that the population of C. pomonella in both the regions has not so variation/diversity as compared to population in Kalam. Similarly the Nei’s genetic identity revealed that higher genetic similarity (70.06%) was resided by the C. pomonella population at Matta and Madyan while the low level of identity (37.58%) were examined in isolates from Madyan and Kalam (Fig. 3).

 

 

Table III.- Nei’s unbiased measures of genetic identity (above diagonal) and genetic distance (below diagonal) for C. pomonella populations collected from three geographically distant region Swat based on 21 RAPD primers analysis

Population	Matta	Kalam	Madyan
Matta		0.5732	0.7006
Kalam	0.5564		0.3758
Madyan	0.3558	0.9787

 

Discussion

On the basis of amplification pattern, highest gene frequency (f) was evaluated for RAPD primer GLH-13 (f = 3.33) followed by GLH-14 (f = 2.77), GLE-18 (f = 2.66), GLG-11 (f = 2.08), GLF-07 (f = 1.85), GLH-03 and GLB-04 (f = 1.66) each, GLC-16 (f = 1.55), GLC-06 (f = 1.55), GLC-02 (f = 1.46), GLE-09 (f = 1.41), GLF-04 (f = 1.42) and GLD-08 (f = 1.03), whereas the rest of RAPD primers depicted the allele frequency below 0.99. The overall mean of the genetic frequency among three population of C. pomonella observed was 1.33. These results are in close conformity with findings of Lei-Men et al. (2012), who reported that the mean number of alleles per locus ranged from f = 4.3 to f = 12.6 and two populations of C. pomonella from Heilongjiang Province in northeastern China had the largest number of alleles (f = 12.6 and f = 10.6). From populations of northwestern China, one population showed the highest value of mean number of alleles (f = 9.6), followed by the second population (f = 8.6) and third population (f = 8.4). Nevertheless, the gene frequency of null alleles ranged from 0.010 to 0.203, values typical for lepidopterans (Megle’cz et al., 2004; Dakin and Avise, 2004), which further confirmed these results.

On the basis of intensification pattern, the genetic diversity (I) of the three population of the C. pomonella disclosed that highest genetic variation was detected for the RAPD primers GLC-02, GLE-18 and GLC-16 (I = 0.44 each), followed by GLC-13 (I = 0.38), GLH-14 (I = 0.37), GLA-19 (I = 0.35), GLB-15 (I = 0.34), GLB-12 (I = 0.33) and GLB-04 (I = 0.32), while the rest of the RAPD primers elucidated below 0.29 genetic diversity among three population of the C. pomonella. The overall mean of the genetic diversity of the three population was 0.29. These results are in close concordance with findings of Khaghaninia et al. (2009), they find out that by using RAPD primers genetic diversity within population of C. pomonella based on Nie’s gene index ranged from 0.228 to 0.281 at Shabestar and Zunuz populations, respectively. They also observed the maximum (0.14) and minimum (0.04) genetic distances between the population of C. pomonella at different geographical locations in Iran. Contrary to our results regarding genetic diversity among the population of C. pomonella, Bues et al. (1995), observed low genetic differentiation between sampled populations of C. pomonella by using allozyme markers. These results are also in close agreement with findings of Chen and Dorn (2010) reported important genetic differentiation at local geographic scale (even less than 10 km), which they mostly attributed to the sedentary behaviour of C. pomonella through microsatellite study on populations from Switzerland. Bayar et al. (2006) through investigations on population variation of Aeolothrips intermedius Bagnall, found population-specific RAPD primers for their molecular differentiation.

The highest Shannon’s information index (h) for each allele of the primers GLC-04, GLE-18 and GLC-16 were h = 0.63, followed by GLC-13 (h = 0.55), GLH-13 (h = 0.53) and GLA-19 (h = 0.50), whereas the rest of the primers were less than 0.49 Shannon’s information index (h) values. The overall mean of the Shannon’s information index (h) was 0.44. These results are closely corroborated with findings of Timm et al. (2006) who reported that genetic diversity was found to be high within the South African C. pomonella population (h = 0.18). Only three amplified fragment length polymorphism (AFLP) fragments were monomorphic at the 95% level, resulting in 98.60% of the scored loci being polymorphic. They further stated that gene diversity within C. pomonella populations collected from different regions varied, whereas gene diversity within the English and Canadian populations was reckoned as h = 0.046 and 0.052, respectively.

The results pertaining to the Nei’s genetic identity and genetic distance among the three population of C. pomonella was also worked out. Higher genetic distance was resided among the isolates from Kalam and Madyan (97.87%) whereas low genetic distance (35.58%) was observed in the C. pomonella isolates from Matta and Madyan. It indicates that the population of C. pomonella in both of the regions has not showed multiplicity as compared to population at Kalam. These results are in close conformity with findings of Khaghaninia et al. (2009), who observed maximum and minimum genetic distances between the population of C. pomonella at different geographical locations in Iran and significant correlation was noticed between genetic and geographic distance matrices in the population of C. pomonella revealed by Mantel test. It is assumed that due to climate change and frequent insecticide treatments, C. pomonella populations differentiate into many ecotypes with different biological and physiological requirements related to their development (Thaler et al., 2008). These results contrast with those obtained for C. pomonella populations from France and Switzerland, in which allozyme analysis indicated a large degree of genetic similarity between geographic populations (Bue’s and Toubon, 1992). Nonetheless, the molecular variability among and between populations of C. pomonella has been determined in Italy by using RAPD primers (Gomez et al., 2004, 2005).

Likewise, the Nei’s genetic identity was also reckoned. The higher genetic similarity (70.06%) was dwelled by the C. pomonella population at Matta and Madyan whilst the low level of identity (37.58%) was examined in isolates from Madyan and Kalam. These results are supported by the study of Timm et al. (2006), who used AFLP markers and successfully ascertained differences among sampled C. pomonella populations even at small geographic distances. Besides, Timm’s et al. (2006) study was corroborated by Thaler et al. (2008), who also used AFLP markers to study the molecular phylogeny and genetic diversity of C. pomonella population. Deverno et al. (1998) distinguished two closely related species of lepidopterious moths using 17 species-specific RAPD primers. RAPD primers have also been used for studying the population diversity of the Hessian fly, Mayetiola destructor Say, and the wheat stem sawfly, Cephus cinctus Fitch, in Syria and America, respectively (Lou et al., 1998; Naber et al., 2000).

 

Conclusion

RAPD markers are efficient tools for assessing the population variation in insect pests and knowledge of the genetic variation within C. pomonella populations is necessary for their efficient control and management, thus such studies may offer an insight on the possible resistance to insecticides. Higher genetic distances among the populations of C. pomonella could be attributed to climatic conditions of the studied areas, geographical locations and elevations. Nonetheless, further studies are required to assess the molecular variation among the population of this pest.

 

Acknowledgement

This research manuscript is a part of the PhD thesis of the 1st author submitted to The University of Agriculture Peshawar, Khyber Pakhtunkhwa, Pakistan and the authors would like to acknowledge the support and technical help of the Institute of Bio-Technology and Genetic Engineering (IBGE), The University of Agriculture, Peshawar Pakistan and Higher Education Commission (HEC) of Pakistan for their financial support.

 

Statement of conflict of interest

The authors declare no conflict of interest.

 

References

Bardakci, F., 2001. Random amplified polymorphic DNA (RAPD) markers. Turk. J. Biol., 25: 185-196.

Barinova, S., Naiz, A., Barkat, U. and Sarim, F.M., 2013. Ecological adaptation to altitude of algal communities in the Swat Valley (Hindu Cush Mountains, Pakistan). Exp. Opin. environ. Biol., 2: 2. https://doi.org/10.4172/2325-9655.1000104

Barnes, M.M., 1991. Codling moth occurrence, host race formation, and damage. In: Tortricid pests, their biology, natural enemies and control (eds. L.P.S. van der Geest and H.H. Evenhuis). Elsevier, Amsterdam, Netherlands, pp. 313-327.

Bayar, K., Torjak, O. Kiss, E., Gyulai, G. and Heszky, L., 2006. Genetic variation within and among populations of Aeolothrips intermedius. J. Zool., 73: 67-73.

Benecke, M., 1998. Random amplified polymorphic DNA (RAPD) typing of necrophagous insects (Diptera and Coleoptera) in criminal forensic studies: Validation and use in practice. Forensic Sci. Int., 98: 157-168. https://doi.org/10.1016/S0379-0738(98)00150-9

Bues, R. and Toubon, J.F., 1992. Polymorphisme enzymatique dans diffe´ rentes populations de Cydia pomonella L. (Lep., Tortricidae). Acta Oecol., 13: 583-591.

Buès, R., Toubon, J.F. and Poitout, H.S., 1995. Variabilité écophysiologique et enzymatique de Cydia pomonella L. En fonction de l’origine géographique et de la plante hôte. Agronomie, 15: 221-231. https://doi.org/10.1051/agro:19950306

Calkins, C.O. and Faust, R.J., 2003. Overview of area wide programs and the program for suppression of codling moth in the western USA directed by the United States Dep. Agric. Res. Serv. Pest. Manage. Sci., 59: 601-604.

Chen, M.H. and Dorn, S., 2010. 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. https://doi.org/10.1017/S0007485309006786

Ciglar, I., 1998. Integrirana zaštita voćaka i vinove loze. Zrinski, Čakovec.

Dakin, E.E. and Avise, J.C., 2004. Microsatellite null alleles in parentage analysis. Heredity, 93: 504-509. https://doi.org/10.1038/sj.hdy.6800545

Delaat, D.M., Carvalho, M.R.S., Acedo, M.D.P. and de Fonseca, C.G., 2005. Applicability of RAPD markers on silver-stained polyacrylamide gels to ascertain genetic diversity in Peripatus acacioi (Onychophora, Peripatidae). Genet. Mol. Res., 4: 716-725.

Deverno, L.L., Smith, G.A. and Harrison, K.J., 1998. Randomly amplified polymorphic DNA evidence of introgression in two closely related sympatric species of Coniferophagous choristoneura (Lepidoptera: Tortricidae) in Atlantic Canada. Annls entomol. Soc. Am., 91: 248-259. https://doi.org/10.1093/aesa/91.3.248

Dorn, S., Schumacher, P., Abivardi, C. and Meyhofer, R., 1999. Global and regional pest insects and their antagonists in orchards: Spatial dynamics. Agric. Ecosyst. Environ., 73: 111-118. https://doi.org/10.1016/S0167-8809(99)00020-1

Dunley, J.E. and Welter, S.C., 2000. Correlated insecticide cross-resistance in azinphosmethyl resistant codling moth (Lepidoptera: Tortricidae). J. econ. Ent., 93: 955-962. https://doi.org/10.1603/0022-0493-93.3.955

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

FAOSTAT, 1998. FAO database. Available at: http://faostat.fao.org/site/339/default.aspx

Fritsch, E., Undorf-Spahn, K., Kienzle, J., Zebitz, C.P.W. and Huber, J., 2005. Codling moth granulovirus: First indications of variations in the susceptibility of local codling moth populations. Nachrichtenbl. Deut. Pflanzenschutz., 57: 29-34.

Gomez, D.R., de Silva, J.J. Costa, F. and Binneck, E., 2005. Population structure of the Brazilian Southern green stink bug, Nezara viridula. J. Insect Sci., 36: 216-227.

Gomez, D.R., Delpin, K.E., Almeida, M.R. and Hirose, E., 2004. Genetic differentiation among Brazilian populations of Euschistus heros (Fab.) (Heteroptera: Pentatomidae) based on RAPD analysis. Neotrop. Ent., 33: 179-187. https://doi.org/10.1590/S1519-566X2004000200009

Khaghaninia, S., Mohamadi, S.A., Sarafrazi, A.M., Iraninejad, H. and Ebrahimi, K.E., 2009. Analysis of population structure in codling moth, Cydia pomonella (Lepidoptera: Tortricidae) from north west of Iran using RAPD markers. Türkiye III. Bitki Koruma Kongresi, 15-18 Temmuz, Van.

Knight, A.L., Brunner, J.F. and Alston, D., 1994. Survey of azinphosmethyl resistance in codling moth (Lepidoptera: Tortricidae) in Washington and Utah. J. econ. Ent., 87: 285-292. https://doi.org/10.1093/jee/87.2.285

Landry, B., Powell, J.A. and Sperling, F.A.H., 1999. Systematics of the Argyrotaenia franciscana (Lepidoptera: Tortricidae) species group: Evidence from mitochondrial DNA. Annls entomol. Soc. Am., 92: 40-46. https://doi.org/10.1093/aesa/92.1.40

Lei-Men, Q., Chen, M.H., Zhang, Y.L. and Feng, J.N., 2012. Genetic structure and diversity of a newly invasive species, the codling moth, Cydia pomonella (L.) (Lepidoptera: Tortricidae) in China. Springer Science Business Media B.V. Biol. Invasions, 15: 447-458.

Lima, L.H.C., Campos, L., Moretzsohn, M.C. and de Oliveira, M.R.V., 2002. Genetic diversity of Bemesia tabaci (Genn.) populations in Brazil revealed by RAPD markers. Genet. mol. Biol., 25: 217-223. https://doi.org/10.1590/S1415-47572002000200016

Lou, K.F., Weiss, M.J., Bruckner, P.L., Talbert, L. and Martin, J.M., 1998. RAPD variation within and among geographic populations of wheat stem sawfly (Cephus cinctus Norton). J. Hered., 89: 329-335. https://doi.org/10.1093/jhered/89.4.329

Megle´cz, E., Petenian, F. Danchin, E., D’Acier, A.C., Rasplus, J.Y. and 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. https://doi.org/10.1111/j.1365-294X.2004.02163.x

Meraner, A., Brandstatter, A., Thaler, R., Aray, B., Unterlechner, M., Niederstatter, H., Parson, W., Zelger, R., Dalla, J. and 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. https://doi.org/10.1016/j.ympev.2008.05.026

Naber, N., El-Bouhssini, M., Labhilili, M., Udupa, S.M. and Nachit, M.M., 2000. Genetic variation among populations of the Hessian fly Mayetiola destructor (Diptera: Cecidomyiidae) in Morocco and Syria. Bull. entomol. Res., 90: 245-252. https://doi.org/10.1017/S0007485300000365

Nei, M. and Li, W.H., 1979. Mathematical model for studying genetic variation in terms of restriction endonucleases. Proc. natl. Acad. Sci., 76: 5269-5273. https://doi.org/10.1073/pnas.76.10.5269

Sauphanor, B., Brosse, V., Monier, C. and Bouvier, J.C., 1998. Differential ovicidal and larvicidal resistance to benzoylureas in the codling moth, Cydia pomonella. Ent. Exp. Appl., 88: 247-253. https://doi.org/10.1046/j.1570-7458.1998.00369.x

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

Timm, A.E., Geertsema, H. and Warnich, L., 2006. Gene flow among Cydia pomonella (Lepidoptera: Tortricidae) geographic and host populations in South Africa. J. econ. Ent., 99: 341-348. https://doi.org/10.1603/0022-0493-99.2.341

Williams, J.G.K., Kubelik, A.R., Livak, K.J., Rafalski, J.A. and Tingey, S.V., 1990. DNA polymorphisms amplified by arbitrary primers are useful as genetic markers. Nucl. Acids Res., 18: 6531-6535. https://doi.org/10.1093/nar/18.22.6531

Zimmermann, M., Wahlberg, N. and Descimon, H., 2000. Phylogeny of Euphydryas Checkerspot butterflies (Lepidoptera: Nymphalidae) based on mitochondrial DNA sequence data. Annls entomol. Soc. Am., 93: 347-355.