Mitochondrial DNA Sequence Variation of the Mushroom Pest Flies, Lycoriella Mali (Diptera: Sciaridae) and Coboldia Fuscipes (Diptera: Scatopsidae), in Korea
Total Page:16
File Type:pdf, Size:1020Kb
Appl. Entomol. Zool. 36 (4): 451–457 (2001) Mitochondrial DNA sequence variation of the mushroom pest flies, Lycoriella mali (Diptera: Sciaridae) and Coboldia fuscipes (Diptera: Scatopsidae), in Korea Jin-Sik Bae, Iksoo Kim, Seong-Ryul Kim, Byung-Rae Jin and Hung-Dae Sohn* College of Natural Resources and Life Science, Dong-A University, Busan 604–714, Korea (Received 12 March 2001; Accepted 14 May 2001) Abstract We analyzed a portion of mitochondrial COI gene sequences (406 bp) to investigate the genetic diversity and geo- graphic variation of the mushroom pest flies Lycoriella mali and Coboldia fuscipes in Korea. L. mali showed minimal sequence divergence (0.2%) in two mtDNA haplotypes, whereas C. fuscipes showed an intermediate level of sequence divergence (1.2% at maximum) compared with other relevant studies. While L. mali was fixed with one haplotype ex- cept for one population, C. fuscipes possessed a total of ten mtDNA haplotypes, and six of these occurred commonly in multiple populations. We ascribed the difference in the level of genetic variation between the two species to a dif- ference in the degree to which they are dependent on cultured mushroom, which is a fluctuating food source. In C. fuscipes, as in other cosmopolitan insect pest species, a high rate of female migration (Nm50.947–infinite) and little 5 2 genetic differentiation (FST 0.345– 0.094) between populations were estimated. Key words: Mushroom fly, Coboldia fuscipes, Lycoriella mali, mitochondrial DNA, genetic diversity Coboldia fuscipes is another major oyster mush- INTRODUCTION room pest in Korea. The population grows in the Understanding the genetic population structure summer season (from July to August), and dimin- of insect pests may enable us to predict increases ishes somewhat in the fall season (from September and decreases in population size, as well as to November), in Korea (Kim et al., 1999; Choi et broaden our understanding of the mode of occur- al., 2000). This species occurs in diverse habitats rence and migration among localities (Roderick, such as decaying plants, animal materials, and 1996). Thus, such knowledge should provide an some flowering grasses (Shiraki, 1981). In Korea, important biological information for the control of mushroom flies recently received attention as pests, pest species. and taxonomy, geographic distribution, life cycle, Mushroom flies cause severe damage to mush- and some physiological aspects have been investi- room: larvae feed on the mycelium and fruitbody gated (Choi et al., 1997, 1999, 2000; Lee et al., of the mushroom and adult flies transport germs 1998, 1999; Kim et al., 1999, 2000). However, such as nematodes, mites and mold spores (Clift, mushroom flies have never been subjected to stud- 1979; Clancy, 1981; Wetzel, 1981; Kim and ies on genetic structure, gene flow, and genetic di- Hwang, 1996). Lycoriella mali is the most abun- versity. dant pest, occurring through at the year in most re- Mitochondrial DNA (mtDNA) has a high evolu- gions of Korea (Kim et al., 1999). This species is tionary rate compared to the functional counterpart distributed in North America, Europe, and Asia of nuclear DNA. In addition, it is inherited mater- (Lee et al., 1999). In North America, crop losses nally, lacking genetic recombination, and is easy to caused by L. mali were averaged to at least 17% handle (Brown et al., 1979, 1982; Cantatore and (Cantelo, 1979). L. mali inhabits mushrooms as Saccone, 1987; Harrison, 1989). These characteris- well as rotting wood, decaying potato, and rotting tics make the mtDNA molecule a particularly ap- vegetable (Stamets and Chilton, 1984; Lee et al., propriate marker for tracing the recent evolutionary 1999). history of animals, including colonization, intro- * To whom correspondence should be addressed. 451 452 J.-S. Bae et al. duction, and population bottlenecks (Wilson et al., 1985). It is known that the cytochrome oxidase subunit I (COI) gene of mtDNA is highly variable in its DNA sequence, especially at the silent sites. We se- lected a portion of the COI gene that includes the membrane-spanning helices M3, M4, and M5, ex- ternal loops E2 and E3, and internal loop I2 (Lunt et al., 1996). This portion has been proven to be useful for the study of intraspecific genetic varia- tion in insects (Simon et al., 1994; Kim et al., 2000a, b, c; Lee et al., 2000). We sequenced a portion of the COI gene from natural populations of L. mali and C. fuscipes in Korea to study the extent and nature of genetic variation, geographic structure, and gene flow among populations. MATERIALS AND METHODS Insects. Adult L. mali and C. fuscipes were col- Fig. 1. Sampling locations of C. fuscipes and L. mali in leted using an aspirator at oyster mushroom houses Korea. General locality names are as follows: 2. Yeongwol- in eight localities in Korea from June 1999 to July gun, Gangwon-do; 3. Gyeongju-si, Gyeonsangbug-do; 4. Hwaseong-gun, Gyeonggi-do; 5. Busan-si; 6. Boseong-gun, 2000 (Fig. 1). One individual of each species (ani- Jeonlabuk-do; 7. Sacheon-si, Gyeongsangnam-do; 8. Chungju- mal numbers JB2 and JB3, respectively) was ob- si, Chungcheongbug-do; 9. Uiseong-gun, Gyeongsangbug-do. tained from indoor-rearing. Samples collected in the field were frozen at 270°C until molecular discovered (LM1 and LM2 in L. mali and CF1 and analysis. CF2 etc., in C. fuscipes). MtDNA. Total DNA was extracted by following Genetic diversity and distance. Haplotype di- the standard Proteinase K method (Kocher et al., versity and nucleotide diversity within local popu- 1989). A part of the COI gene was amplified by lations were estimated for C. fuscipes, using Arle- PCR using primers CI-J-1751 (59-GGAGCTCCT- quin ver. 1.1 (Schneider et al., 1996). Genetic dis- GACATAGCATTCCC-39) and CI-N-2191 (59-CC- tance and female migration rate per-generation CGGTAAAATTAAAATATAAACTTC-39) (Simon were estimated by subroutines in Arlequin ver. 1.1 et al., 1994). PCR conditions were as follows: an (Schneider et al., 1996). Genetic distance was initial denaturation step at 94°C for 5 min, 40 cy- based on pairwise FST indices (F ST) following the cles of 94°C for 30 s, 50°C for 40 s, and 72°C for approach described by Excoffier et al. (1992). Sta- 45 s, and a final extension step at 72°C for 7 min. tistical significance of the difference between a pair To confirm the successful DNA amplification, elec- of populations was tested by permutation (10,000 trophoresis was carried out using 0.5X TAE buffer bootstraps; Excoffier et al., 1992). Pairwise FST in 1% agarose gel. The PCR product was purified values were used to estimate per generation migra- using a PCR purification Kit (QIAGEN, Germany). tion rate, Nm (the product of the effective popula- DNA sequencing was performed using an ABI 377 tion size Ne and migration rate m), based upon the 5 1 Genetic Analyzer (PE Applied Biosystems, USA). equilibrium relationship: FST 1/(2Nm 1). Sequence alignment was performed using an IBI MacVector (ver. 6.5). When homologous se- RESULTS quences from two individuals differed by at least one nucleotide site, the two sequences were consid- Lycoriella mali ered as different haplotypes. Haplotype designa- Sequence analysis of 55 L. mali individuals re- tions were applied to new sequences as they were sulted in only two haplotypes (LM1, LM2) for a Genetic Variation of Mushroom Flies 453 Table 1. A list of trapping localities, animal numbers, sex, mitochondrial COI haplotypes and GenBank accession numbers of C. fuscipes Collecting locality COI GenBank accession Collection date Animal number Sex (no. of individuals) Haplotype number 1. Dong-A Univ., 1998. 12. 17 JB3 F CF1 AF319839 Busan-si (1) 2. Yeongwol-gun, 1999. 7. 28 C2 F CF2 AF319840 Gangwon-do (8) C3 F CF3 AF319841 C5 F CF4 AF319842 C6 M CF5 AF319843 C7 M CF6 AF319844 C8 M CF3 AF319845 C9 M CF6 AF319846 C11 F CF2 AF319847 5. Busan-si (8) 2000. 5. 9 C21 F CF2 AF319848 C22 F CF2 AF319849 C24 F CF4 AF319850 C25 F CF1 AF319851 C27 M CF1 AF319852 C28 M CF1 AF319853 C29 M CF7 AF319854 C30 M CF1 AF319855 6. Boseong-gun, 2000. 5. 12 C33 F CF4 AF319856 Jeonlanam-do (8) C34 F CF8 AF319857 C35 F CF6 AF319858 C36 F CF8 AF319859 C38 M CF4 AF319860 C39 M CF4 AF319861 C41 M CF3 AF319862 C42 M CF9 AF319863 7. Sacheon-si, 2000. 5. 12 C47 F CF1 AF319864 Gyeongsangnam-do (8) C48 F CF6 AF319865 C49 F CF2 AF319866 C50 F CF4 AF319867 C51 F CF1 AF319868 C59 M CF7 AF319869 C60 M CF1 AF319870 C61 M CF4 AF319871 9. Uiseong-gun, 2000. 7. 21 C62 F CF2 AF319872 Gyeongsangbug-do (8) C66 M CF4 AF319873 C67 M CF2 AF319874 C72 M CF6 AF319875 C73 M CF2 AF319876 C76 M CF3 AF320759 C78 M CF6 AF320760 C79 M CF10 AF320761 partial sequencing of the COI gene (406 bp). Se- AF319784–AF319838. quence alignment revealed only one variable nu- cleotide site (position 79), which was the 3rd posi- Coboldia fuscipes tion of the codon and did not substitute an amino A total of ten COI haplotypes (CF1–CF10) was acid. Except for two individuals in Busan (Fig. 1), found from 41 individuals of C. fuscipes (Table 1). all individuals possessed haplotype LM1. GenBank These haplotypes revealed eight variable nu- accession numbers of individual L. mali are cleotide sites, seven of which were the 3rd position 454 J.-S. Bae et al. Table 2. Within-locality diversity estimates of C. fuscipes Locality Na NHb Hc Sd p e 2. Yeongwol 8 5 0.89 5 0.006 5. Busan 8 4 0.75 6 0.005 6.