PHYLOGENY of Glyphodes Guenee (Lepidoptera: Crambidae
Total Page:16
File Type:pdf, Size:1020Kb
Treubia 2003 33 (1) 35-42 PHYLOGENY OF Glyphodes Guenee (Lepidoptera: Crambidae: Spilomelinae) BASEDON NUCLEOTIDE SEQUENCE VARIATION IN A MITOCHONDRIAL CO I GENE: CONGRUENCE WITH MORPHOLOGICAL DATA Hari Sutrisno Division of Zoology, Research Center for Biology, Indonesian Institute of Sciences, JI. Raya Bogor Km. 46, Cibinong, 16911 Indonesia E-mail: [email protected] Abstract The phylogeny of Glyphodes Cuenee (14 species) and [our outgroup species (Feltia jaculifera, Metallarcha aureodiscalis, Talanga sexpunctalis and Agrioghrypta eurytusalis) was inferred from nucleotide sequence variation across a 686-bp region in the CO I gene. Over the entire 686-bp region, 19.9% sites were informative (3.35% in the 1"-, 0.29% in the 2''''- alld 16.32% in 3"'-codon position). The results also showed that the base composition of this region was high A+T biased (C= 0.258) and the averages of estimated sequence divergence in the comparisons between species within and between groups were 7.1% and 9.0%, respectively. In general, the phylogeny based 011 CO Igene by including all substitutions or any partial data set used in this study was not only able to recover almost the three monophyletic groups within Glyphodes as previously recovered by morphological phylogeny but also showed more clearly the relationships among them: Glyphodes group 2 was branched off first then followed by group 3 and 1. Key words: CO I, Gbjphodes, morphology, phylogeny. Introduction The Glyphodes Cuenee, 1854, is medium-sized moth (15-51 mm) belongs to the subfamily Spilomelinae (Common, 1990; Shaffer et al., 1996). This genus is diverse and widespread in tropical regions with some species penetrating into subtropical and warm temperate areas (Common, 1990;Robinson et al., 1994).The genus consists of 120 species throughout the world, and about 25 and 17 species have been recorded in the Southeast Asia and Australia, respectively (Robinson et al., 1994; Shaffer et al., 1996). As other pyraloids in the Indo-Australian regions, this genus is also poorly defined and has been suggested that it is seriously in need of revision since the current generic classification of this genus was based on superficial similarities only (Robinson et al., 1994;Sutrisno & Horak, 2003). However, the comprehensive study on the taxonomy or phylogeny to assess the monophyly and the relationships among them has never been conducted except for a cladistic analysis based on the Australian Glyphodes and its allied genera using morphological data (Sutrisno, 2002). The cladistic analysis showed that the genus Glyphodes falls into three monophyletic groups (Glyphodes group I, 2 and 3); each of them can be diagnosed by the adult morphology. In general, the result showed that morphological characters used 36 Sutrisno: Phylogeny of Glyphodes Guenee (Lepidoptera: Crambidae: Spilomelinae) in the previous study were able to infer the phylogenetic relationships of this genus. However, the evolutionary change of morphological characters in this genus is extremely complicated that this approach did not produce a clear-cut picture of evolutionary history; i.e., for the relationships among groups within Glyphodes. In addition, the relationships among closely-related species within Glyphodes group 3 seem very difficult to be resolved. Therefore, another possible approach such as molecular phylogeny that has been proved to give more powerful resolution is really needed. Mitochondrial genes have been used intensively to infer the phylogenetic relationships among groups of insects as have been repeatedly reported by numerous authors (reviewed in Simon et al., 1994). CO I gene is one of the most common to be considered infering the relationships among closely-related species in several groups of Lepidoptera, as individual gene or to be combined with other genes (Andrew, 1994; Sperling & Hickey, 1994; Caterino & Sperling, 1999; Landry et al., 1999; and Blum et al., 2003). Therefore, in this issue I used mitochondrial CO I gene to reassess the monophyly of Glyphodes and to estimate the relationships among species within this genus. Material and Methods DNA extraction and Sequencing A total of 17 species representing the four genera (Glyphodes, Talanga, Agrioglypta and Metal/archa) from various localities in Indonesia and Australia were collected using a light trap and preserved in absolute alcohol (Table 1). The DNA extraction method and the PCR amplification conditions used in this study were similar with the procedures described in my previous study (Sutrisno, 2003). I used Mt 0-4 and Mt 0-9 primers to amplify CO I gene for a total of 686-bp and the complete sequences were Mt 0-4: 5'- TACAATTTATCGCCTAAACT TCAG CC-3', and Mt 0-9: 5'-CCCGGTAAAATT AAAATATAAACTTC-3' (http://www.biotech.ubc.ca/services/naps/primers.html). Sequencing was performed using an AB!PRISMDye Terminator Cycle Sequencing Ready Reaction Kit (Perkin-Elmer) on AB! PRISM model 310 Genetic Analyzer (PE Applied Biosystems). Phylogenetic Analysis For the phylogenetic analysis, I used Maximum-Parsimony method using an exhaustive search with different data sets: all substitutions; transversions only; weighting scheme 2:1; ntl + nt3; and nt3 only and Neighbor Joining method using K80 model (Kimura, 1980) in PAUP 4.0b4a (Swofford, 2001); and Bootstrap test (Felsenstein, 1985) to evaluate the statistical confidence of the MP and NJ trees. Four species, Agrioglypta eurytusalis (Walker), Talanga sexpunctalis Moore, Metallarcha aureodiscalis Meyrick and Feltia jaculifera Meyrick (Noctuidae) were included as outgroup taxa in the present analysis. The sequence of CO I-CO II of the last species is available from the gene bank with the accession number U60990. Treubia 2003 33 (1) 37 Table 1. Species examined for the molecular study Species Collection locality Date of Collection Agrioglypta eurytusalis (Walker) G. Pangrango, Java September 2000 Glyphodes G. actorionalis (Walker) Sorong, Papua February, 2002 G. apiospila (Turner) Sorong, Papua February, 2002 G. bicolor (Swainson) Sorong, Papua February, 2002 G. bivitralis Cuenee Patunuang, Sulawesi April, 2001 G. caesalis (Walker) Menado, Sulawesi March,2001 G. conjuncialis Walker Sorong, Papua February, 2002 G. cosrnarcha Meyrick Patunuang, Sulawesi April,2001 G. doleschalii Lederer Sorong,Papua January, 2002 G. jlauizonalis Hampson Sorong, Papua January, 2002 G. margaritaria (CIerck) Patunuang, Sulawesi April, 2001 G. multilinealis Kenrick Bantimurung, Sulawesi April, 2001 G. onychinalis Cuenee Bucasia, Queensland April, 2001 G. puioerulentalis Hampson Sukabumi, Java February, 2002 G. stolalis Cuenee G. Halimun NP, Java March,2002 Metallarcha aureodiscalis Meyrick Bucasia, Queensland April,2001 Talanga sexpuncialis Moore Patunuang, Sulawesi April, 2001 Table 2. Percent variable sites across 14species of Glyphodes and four outgroup species Total lst-codon 2nd-codon 3rtl-codon Constant (%) 67.9 27.25 31.77 8.89 Uninformative (%) 12.1 2.76 1.31 8.01 Informative (%) 19.9 3.35 0.29 16.32 Number of sites 686 229 229 228 Results Base compositions Sequences of 14 species of the Glyphodes, and four species outgroup were aligned with no evidence of insertion and deletion (aligned sequences are available from the author on request). Over the entire 686-bp region, 67.9 % of the nucleotide positions were constant, 12.1 % were uninformative (i.e., any variants were found in single sequence), and 19.9% were informative (Table2). In addition, the most informative 38 Sutrisno: Phylogeny of G/yphodes Guenee (Lepidoptera: Crambidae: Spilomelinae) sites for maximum parsimony analysis lied in the third-codon positions and the least informative site was in the second-codon positions. Table 3 shows the proportion of nucleic acids and bias in CO I. The results showed that the base compositions were strongly biased (C= 0.258) with the A+T proportions ranged from 69.82% to 72.01%. Interspecific variations in base compositions CO I were very low for a total of nucleotides (the Chi-square test, X2= 10.684, df= 51, P=1.00). In addition, the averages of estimated sequence divergence in the comparisons between species within and between groups were 7.1%and 9.0%, respectively. Table 3. Proportion of each nucleotide and bias in CO I Codon position l'l-codon 2nd-codon 3rd-codon Mean A 0.316 0.170 0.485 0.324 C 0.139 0.249 0.042 0.142 G 0.264 0.169 0.097 0.175 T 0.279 0.410 0.462 0.384 Bias 0.258 Bias in base composition is calculated as 4 C = (X)2J, - 0251 /::1 Where C is the composition bias and ci is the frequency of its base. Phylogenetic analyses An exhaustive search in maximum-parsimony analysis using all nucleotides substitutions resulted in two equally MP trees (Fig.l). The MP trees showed that the three groups within the Glyphodes as previously recognized in the morphological study (Sutrisno, 2002) were recovered. However, the relationships among species within Glyphodes group 2 were less-resolved: G. onychinalis was evolved independently. When the weighting scheme 2:1,or transversional substitutions only or the combination between 1"_and 3rd-codon (nt1+nt3) or 3rd-codon (nt3) only were included in the analysis, the relationships among five species of Glyphodes group 2 was also poorly-resolved, except for the relationships between G.flavizonalis and G. apiospila (Figs. 2-5). The descriptive tree statistics for parsimony analyses are presented in Table 4. The NJ tree based on all substitutions revealed that each group was found to be monophyletic (Fig. 6) and the relationships among them were well-resolved: Glyphodes group 2 was branched off first then followed by group 3 and 1. In addition, the relationships among the three species of Glyphodes group 3 were also well-resolved in both MP and NJ analyses: G. doleschalii lied in the basal group, G. actorionalis was the sister group of G. bivitralis. Treubia 2003 33 (1) 39 F.jaculifera Fjaculifera M. aureodiscalis M. aureodiscalis 5 A eurytusalis 58 A eurutusalis T sexpunctalis T sexpunclalis 100 50 Gflavizonalis ] 55 Gflavizonalis ] G apiospila ~ 3 G apiospila ~ I--------G onychinalis .g ,-------+--G conjunctalis -ii 1-----------1iG conjunctalis 2 G bicolor 20 G bicolor 19f--------G onychinalis 51 G caesalis ] G.