Appl Entomol Zool (2016) 51:673–676 DOI 10.1007/s13355-016-0442-3

TECHNICAL NOTE

Characterization of eleven microsatellite markers for the tea geometrid obliqua (: Geometridae)

Shengli Jing1,2 · Baoxiao Zheng1 · Li Qiao3 · Kun Zhang1 · Gaoyang Chen1 · Shunqin Guo1 · Xiaoyan Huang1 · Yanping Gao4 · Lili Zhang1,2 · Hongyu Yuan1,2

Received: 23 May 2016 / Accepted: 2 August 2016 / Published online: 29 September 2016 © The Japanese Society of Applied Entomology and Zoology 2016

Abstract The tea geometrid, Ectropis obliqua Prout (Lepi- Keywords Genetic diversity · Genomic library · doptera: Geometridae), is one of the most threatening Geometridae · Lepidoptera · Microsatellite pests of tea ( (L.) O. Kuntze) plants. In this study, we developed eleven new polymorphic micro- satellite markers from tea geometrid by using the micros- Introduction atellite motif enriched library method. Polymorphism of each locus was detected in 69 individuals from three natu- The tea geometrid, Ectropis obliqua Prout, (Lepidop- ral populations. The number of alleles varied from 2 to 8, tera: Geometridae), is one of the most important pests of and the expected and observed heterozygosities ranged tea (Camellia sinensis (L.) O. Kuntze) plants in south- from 0.042 to 0.801 and from 0.042 to 0.875, respectively. ern China, especially in the producing regions of green Six loci had significant deviations from the Hardy–Wein- tea such as Jiangsu, Zhejiang, Anhui and Henan prov- berg equilibrium (HWE). These microsatellite markers inces (Zhang and Tan 2004). These tea-producing areas will be useful for the future study of this agricultural pest are frequently reported to have severe yield losses of tea in genetic diversity, population structure and evolution of production owing to tea geometrid. Currently, these pests resistance. This information will be allowed to establish are mainly controlled with chemical pesticides which new pest management strategies and improve the effective- can cause some problems such as pesticide-resistance, ness of pest control program of this species. environmental pollution and chemical residues on tea (Ye et al. 2014). Moreover, the intensive use of these pesticides can lead to change the genetic variation of tea geometrid population. Therefore, it is important to study * Shengli Jing the genetic diversity, population ecology and evolution [email protected] of E. oblique by molecular tools. Recently, the study of * Lili Zhang genetic analysis of E. oblique using mitochondrial gene [email protected] sequences has suggested deep genetic divergence among * Hongyu Yuan the seventeen geographical populations (Zhang et al. [email protected] 2014). However, there are not any available nuclear molecular markers for E. oblique such as microsatellite 1 College of Life Science, Xinyang Normal University, Xinyang 464000, China which has been used in other (Jing et al. 2012a, b; Kim et al. 2015; Kusumi and Su 2014; Suzuki and Yao 2 Institute for Conservation and Utilization of Agro‑bioresources in Dabie Mountains, Xinyang Normal 2014; Tan et al. 2016). University, Xinyang 464000, China In this study, the eleven new genomic microsatellite 3 Xinyang City Academy of Agricultural Sciences, markers of tea geometrid by a method of enriched genomic Xinyang 464000, China library were developed and used to study the genetic diver- 4 Environmental Monitoring Center Station of Zhumadian, sity in three populations of this species. Zhumadian 463000, China

1 3 674 Appl Entomol Zool (2016) 51:673–676

Materials and methods step of 72 °C for 10 min. PCR amplification products were size-fractionated by electrophoresis on 6 % denaturing The genomic DNA was extracted from tissue of E. obli- polyacrylamide sequencing gels that were run at a constant qua with a CTAB protocol (Tang et al. 2010). The (AC)13 power of 60 W, and then detected by silver staining (You and (AAG)8-enriched partial genomic libraries were con- et al. 2008). Allele sizes were scored by comparing with structed, employing a Fast Isolation by Amplified fragment pBR322 DNA/Msp I DNA size markers (Tiangen Biotech). length polymorphism (AFLP) of Sequences COntaining The level of polymorphism was determined for 69 indi- repeats (FIASCO) protocol (Zane et al. 2002). Fragments viduals from three populations of E. obliqua and collected containing microsatellite repeats were cloned into pUC18- from tea plants growing fields in Guangshan (114°30′E T vector (TaKaRa) and transformed into DH5α cells. 31°59′N), Luoshan (114°51′E 31°55′N) and Maanshan Finally, 128 positive clones with suitable insert length were (114°02′E 32°05′N), Henan Province, China. For each identified by PCR and sequenced on an ABI 3730 DNA locus, the number of alleles (Na), observed heterozygosity sequencer. (Ho), expected heterozygosity (He), tests for linkage dis- A total of 125 sequences were finally obtained, of equilibrium (LD) and deviations from Hardy–Weinberg which 110 sequences contained at least one microsatellite equilibrium (HWE) were calculated by the software Arle- locus. As a result, 93 primer pairs were designed by using quin 3.5 (Excoffier and Lischer 2010). The occurrence of BatchPrimer3 (You et al. 2008). For all PCR amplifica- a null allele was estimated by the software Micro-Checker tions, we used a PTC-100 thermal cycler (MJ Research) (Van Oosterhout et al. 2004). and performed in 10 µL volumes containing 10 ng of tem- plate DNA, 0.3 µM of each of the two primers, 0.2 mM deoxynucleotide triphosphates (dNTPs), 2.5 mM MgCl2, Results and discussion 1 PCR buffer, and 1 unit of Taq DNA polymerase (Fer- × mentas). The PCR cycling program, in each case, was Among the 93 primer pairs, 35 successfully amplified in 94 °C for 5 min, followed by 35 cycles of 94 °C for 15 s, the three tested individuals. And then a total of 11 out of 35 55 °C for 15 s, and 72 °C for 30 s, with a final extension primer pairs assessed were selected because they showed

Table 1 Characteristics of 11 new microsatellite markers developed in Ectropis obliqua

Locus Repeat motif Primer sequence (5′-3′) Size range (bp) GenBank accession number

CH3 (TGA)12 F: TCGCTTAGAAATAAGCTGAGA 190–202 KU175831 R: TTCTGAATGACCTGGAAGTTA CH5 (GCGCA)4 F: ATGCTGAAGCTGTACCAGAG 160–180 KU175833 R: GAGTCAGTGGCACCTAATGT CH6 (TG)11 F: GCGGTCAAGTGTAAACCTAT 134–154 KU175833 R: GTACATGAACGACGAGTCTG CH25 (AAG)6 F: GACCTTCAGCAATGAAGATAA 146–160 KU175835 R: TCGCTTGAAATGGTAGAGATA CH32 (TCT)6 F: TCTTGCCAATTACTCCTATCA 148–190 KU175839 R: ACTTGTTTACAGATGGCTCAA CH48 (AC)6 F: ACGACAACTTGTGCATCTAAT 158–180 KU175846 R: GGTTGCATGAGTGAAAGTAAG CH53 (AG)11 F: AGTTTCCACTCAGAGACTGGT 122–132 KU175850 R: CCTGAGTAACTAAACTGGCAAT CH76 (CA)11 F: TTTATCGTATAGAGCGACCAG 136–180 KU175856 R: CAGTGCTAGGTAATGATTTGG CH77 (AC)9 F: TATCACACACACACACACACA 158–180 KU175856 R: CCAGGGAAATTACAGTGTTTA CH80 (AGA)4 F: CGAGGAGTAATGAGAAGGAAT 138–150 KU175857 R: AGACACACACACAGACAGACA CH81 (TG)10 F: CCAAAAAGTTGACGAAACTC 156–232 KU175858 R: CATAGCTAAGAACGAATCGAA

1 3 Appl Entomol Zool (2016) 51:673–676 675

Table 2 Diversity estimation Locus Population GS (n 24) Population LS (n 24) Population MAS (n 21) in three populations of Ectropis = = = obliqua Na Ho He D F Na Ho He D F Na Ho He D F

CH3 6 0.652 0.757 NS – 6 0.762 0.744 NS – 6 0.789 0.691 NS – CH5 5 0.625 0.752 * – 2 0.042 0.042 NS – 5 0.714 0.688 NS – CH6 8 0.875 0.790 * – 5 0.696 0.647 NS – 5 0.800 0.719 * – CH25 3 0.478 0.617 NS – 3 0.417 0.584 NS 0.10 4 0.350 0.650 * 0.14 CH32 2 0.375 0.439 NS – 4 0.609 0.700 NS – 2 0.619 0.470 NS – CH48 2 0.125 0.120 NS – 8 0.696 0.801 * – 2 0.238 0.215 NS – CH53 5 0.652 0.715 NS – 7 0.833 0.729 NS – 5 0.800 0.700 NS – CH76 7 0.792 0.739 NS – 3 0.458 0.606 NS – 7 0.762 0.678 NS – CH77 5 0.458 0.761 * 0.13 8 0.667 0.791 * 0.19 6 0.619 0.718 NS – CH80 5 0.458 0.595 * – 6 0.375 0.768 * – 4 0.476 0.436 NS – CH81 7 0.708 0.768 NS – 3 0.375 0.577 NS – 6 0.619 0.727 NS – Mean 5 0.564 0.641 – – 5 0.539 0.635 – – 5 0.617 0.608 – –

* significant deviations from Hardy–Weinberg expectations (p < 0.05)

N population sample size, Ho observed heterozygosity, He expected heterozygosity, Na number of alleles, D deviation from Hardy–Weinberg equilibrium, F frequencies of null allele, NS not significant, Population GS Guangshan population, Population LS Luoshan population, Population MAS Maanshan population polymorphisms among 69 individuals of the three popula- To the best of our knowledge, this study is the first report tions which were collected from the tea producing areas of on the development of polymorphic microsatellite markers China. The sequences of these 11 microsatellite loci were in this pest. These new loci will provide powerful molecu- deposited in NCBI database (Genbank accession numbers lar tools for investigating the genetic divergence among showed in Table 1). different geographic populations and studying the genetic It was found that the degree of polymorphism among diversity and population structure of this threatening agri- three populations was similar. The average number of cultural pest. alleles per locus for each population was 5, and the mean observed heterozygosity was small different. The num- Acknowledgments This work was supported by the National Natural bers of detected alleles per locus in 69 individuals ranged Science Foundation of China (31401732) and the Project of Henan Provincial Science and Technology (Grant No. 162102110089). from 2 to 8, with an average of 5 alleles. The observed heterozygosity ranged from 0.042 to 0.875 (mean 0.573) and expected heterozygosity from 0.042 to 0.801 (mean 0.628) (Table 2). Six loci (CH5, CH6, CH25, CH76, CH77 References and CH80) deviated significantly from Hardy–Weinberg equilibrium (HWE) (p < 0.05) at least in one population Excoffier L, Lischer HE (2010) Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux due to heterozygote deficiency or the existence of the null and Windows. Mol Ecol Resour 10:564–567 alleles (Table 2). Among of these 6 loci, the microsatellite Jing S, Liu B, Peng L, Peng X, Zhu L, Fu Q, He G (2012a) Develop- marker CH6 displayed significant deviations from HWE in ment and use of EST-SSR markers for assessing genetic diversity all three popupulation. And the results of null alleles test in the brown planthopper (Nilaparvata lugens Stål). Bull Ento- mol Res 102:113–122 showed that the markers CH77 (in Guangshan and Luoshan Jing S, Zhou X, Yu H, Liu B, Zhang C, Wang S, Peng X, Zhu L, Ding populations) and CH25 (in Maanshan population) indicated Y, He G (2012b) Isolation and characterization of microsatellite the existence of the null allele by using Micro-Checker markers in brown planthopper (Nilaparvata lugens Stål). Int J software. Within the populations, significant linkage dis- Mol Sci 13:9527–9533 Kim H, Rodriguez-Saona C, Kwon DH, Park S, Kang T-J, Kim S-J, equilibrium (p < 0.05) was found in 11 pairs of loci: four Hong K-J, Lee H-S (2015) Development and characterization of pairs (CH3 and CH5; CH5 and CH6; CH25 and CH77; 12 microsatellite loci from the blueberry gall midge Dasineura and CH76 and CH77) in Guangshan population, four pairs oxycoccana (Diptera: Cecidomyiidae). Appl Entomol Zool (CH6 and CH53; CH25 and CH53; CH5 and CH77; CH76 50:415–418 Kusumi J, Su Z-H (2014) Isolation and characterization of 15 poly- and CH77) in Luoshan population, and three pairs (CH3 morphic microsatellite markers for the fig-pollinating wasp, and CH5; CH6 and CH80; CH 76 and CH77) in Luoshan Blastophaga nipponica (Hymenoptera: Agaonidae). Appl Ento- population (Table 2). mol Zool 49:487–491

1 3 676 Appl Entomol Zool (2016) 51:673–676

Suzuki S, Yao I (2014) Isolation of nine polymorphic microsatel- You FM, Huo N, Gu YQ, Luo M-c, Ma Y, Hane D, Lazo GR, Dvorak lite loci from the burying beetle, Nicrophorus quadripunctatus J, Anderson OD (2008) BatchPrimer3: a high throughput web (Coleoptera: Silphidae). Appl Entomol Zool 49:493–497 application for PCR and sequencing primer design. BMC Bio- Tan L, Deng J-Y, Fu R-H, Chen Y (2016) Isolation and characteri- inform 9:253 zation of 13 polymorphic microsatellite loci for the fig wasp, Zane L, Bargelloni L, Patarnello T (2002) Strategies for microsatellite Ceratosolen sp. (Hymenoptera: Agaonidae). Appl Entomol Zool isolation: a review. Mol Ecol 11:1–16 51:317–320 Zhang H, Tan J (2004) Tea pests and non-pollution management in Tang M, Lv L, Jing SL, Zhu LL, He GC (2010) Bacterial symbionts China. Anhui Science and Technology Publishing House, Hefei, of the brown planthopper, Nilaparvata lugens (Homoptera: Del- pp 3–9 phacidae). Appl Environ Microbiol 76:1740–1745 Zhang GH, Yuan ZJ, Zhang CX, Yin KS, Tang MJ, Guo HW, Fu JY, Van Oosterhout C, Hutchinson WF, Wills DP, Shipley P (2004) Xiao Q (2014) Detecting deep divergence in seventeen popula- MICRO-CHECKER: software for identifying and correct- tions of tea geometrid (Ectropis obliqua Prout) in China by COI ing genotyping errors in microsatellite data. Mol Ecol Notes mtDNA and cross-breeding. PLoS One 9:e99373 4:535–538 Ye G-Y, Xiao Q, Chen M, Chen X-x, Yuan Z-j, Stanley DW, Hu C (2014) Tea: biological control of and mite pests in China. Biol Control 68:73–91

1 3