Cytochrome C Oxidase Subunit 1 (COI) Profile of the Philippine Helicostylinae (Gastropoda: Stylommatophora: Camaenidae)
Total Page:16
File Type:pdf, Size:1020Kb
Philippine Journal of Science 148 (S1): 1-13, Special Issue on Genomics ISSN 0031 - 7683 Date Received: 31 Jan 2019 Cytochrome C Oxidase Subunit 1 (COI) Profile of the Philippine Helicostylinae (Gastropoda: Stylommatophora: Camaenidae) Gizelle A. Batomalaque1,4,*, Gerard Clinton L. Que1, Tyrill Adolf B. Itong5, Anna Regina L. Masanga1, Emmanuel Ryan C. de Chavez3, and Ian Kendrich C. Fontanilla1,2 1Insitute of Biology, College of Science, University of the Philippines Diliman, Quezon City 1101 Philippines 2Natural Sciences Research Institute, University of the Philippines Diliman, Quezon City 1101 Philippines 3Institute of Biological Sciences, College of Arts and Sciences, University of the Philippines Los Baños 4031 Laguna, Philippines 4Department of Biodiversity, Earth and Environmental Sciences, College of Arts and Sciences, Drexel University, Philadelphia, PA 19104 USA 5College of Science, University of the Philippines Cebu, Cebu City 6000 Philippines The Philippines is the center of radiation of the land snail subfamily Helicostylinae, with around 253 recognized species. Despite their morphological diversity, research on their biology and taxonomy is lacking. We present here the first mitochondrial COI profiles of 32 species of Philippine helicostyline land snails. With the addition of sequences downloaded from GenBank, we tested the utility of the COI for species identification. Relative distributions of intraspecific and interspecific distances overlapped; hence, no barcoding gap was observed. However, 90% of uncorrected interspecific comparisons can distinguish species at 14% genetic distance or lower. Furthermore, the COI barcodes could not discriminate several co-distributed species that have similar conchological features, which should be flagged for taxonomic re-evaluation. Keywords: DNA barcoding, Helicostylinae, mitochondrial COI, Philippine land snails INTRODUCTION exhibit a range in shell forms from discoidal, depressed and keeled, globose, to elongated conical forms (Parkinson The Helicostylinae, a subfamily under family Camaenidae et al. 1987). Within the Philippines, different helicostyline (sensu Bouchet et al. 2017) and order Stylommatophora, species vary in distributions, with most occurring in are hermaphroditic ground and tree snails whose center single islands (e.g., Anixa siquijorensis in Siquijor Is. of diversity is the Philippine Islands (Parkinson et al. and Helicostyla (Calocochlea) chrysocheila in Luzon 1987, Abbott 1989, de Chavez et al. 2015) and whose Is.) and some occurring in multiple adjacent islands (e.g., distribution extends to Taiwan, the Moluccas, and the Leytia fragilis in Samar and Leyte islands and Trachystyla smaller islands off the coast of Borneo (Schileyko 2004, cryptica in the islands of Samar, Leyte, and Mindanao). Schilthuizen et al. 2013). Members of this subfamily There are about 245 species (Batomalaque n/p, Faustino 1930, Richardson 1983, Abbott 1989) belonging to *Corresponding author: [email protected] 1 Special Issue on Genomics Batomalaque et al.: Philippine Helicostylinae DNA Barcodes 23 genera (Schileyko 2004, Bouchet et al. 2017). The amplification using Taq DNA Polymerase and dNTPack current taxonomy of the helicostylines is based on shell (Roche, USA), and universal (forward LCO 1490 morphology, although the reproductive anatomy for GGTCAACAAATCATAAACATATTGG, reverse HCO some species has been described (Schileyko 2004). No 2198 TAAACTTCAGGGTGACCAAAAAATCA; Folmer molecular work has been done to evaluate their current et al. 1994) and stylommatophoran-specific (forward STY_ classification, and phylogenetic relationships among the LCOii ACGAATCATAAGGATATTGGTAC, reverse species are unknown. STY_HCO GAATTAAAATATATACTTCTGGGTG; Fontanilla et al. 2017) primers for the mitochondrial The mitochondrial cytochrome c oxidase subunit 1 (COI) cytochrome c oxidase subunit I (COI) gene. A 50 μL PCR gene has been the gene of choice for DNA barcoding in mix consisted of the following components: 5 μL of 10X animals (Hebert et al. 2003, Meyer and Paulay, 2005, PCR buffer; 1 μL of 10 mM dNTP; 2.5 μL each of 10 mM Park et al. 2011, Siddall et al. 2012, Perez et al. 2014). forward and reverse primers; 22.75 μL distilled water; However, its utility in species discrimination in low- 0.25 μL Taq-polymerase (5 units/μL); 10 μL Q-buffer vagility species (Davison et al. 2009, Virgilio et al. 2010) (Qiagen, USA); 2 μL of 15 mM MgCl2; and 4 μL of 10 appears to be fraught with high error rates due to lack of mM DNA. Amplification protocol consisted of 2 min at baseline differences established through morphology and 94°C followed by 38 cycles of 30 sec at 94°C, 30 sec at a DNA sequence database. In land snails, several studies 45°C, 60 sec at 65°C, and a final extension of 5 min at have shown high levels of mtDNA sequence divergence 72°C. PCR products were visualized in a 1% agarose gel in intraspecific populations (Watanabe and Chiba 2001, using EtBr UV illumination. Pfenniger and Posada 2002, Davison et al. 2009). Amplified PCR products were extracted using QIAquick® In this paper, we present the first COI profiles of the Gel Extraction Kit (Qiagen, USA). Purified samples were Philippine helicostyline land snails, and we test the sent to 1stBASE, Malaysia for sequencing. utility of 463-bp COI barcodes in distinguishing among morphological species. Sequence Alignment and DNA Barcoding Analysis Sequences were assembled using the STADEN package v.1.5.3 (Staden et al. 2000), and aligned using BioEdit MATERIALS AND METHODS v.5.0.5 (Hall 1999). Only unambiguously aligned nucleotide sites were included in the analyses. All sequences were deposited in GenBank (Table 1). Taxon Sampling and Identification A total of 134 specimens attributed to 35 camaenid An additional 292 stylommatophoran COI sequences species were collected from 27 localities across the from GenBank were analyzed together with those from Philippines (Table 1), representing approximately 15% this study for a total of 423 sequences to test their utility of the total nominal species of Helicostylinae. The 41 in species discrimination. Sequences of species under camaenid species comprised 35 species under subfamily family Polygyridae were used as outgroups, following the Helicostylinae and three under subfamily Bradybaeninae. recent molecular helicoid phylogeny of Sei et al. (2017). Cuttings of foot tissue were preserved in 95% (v/v) Sequences were viewed and trimmed to 463 nucleotides ethanol, while the vouchers were preserved in 70% (v/v) common to all taxa using BioEdit v.7.0.9 (Hall 1999) ethanol. Identification of species was based on shell and aligned using the ClustalW (Thompson et al. 1994) morphology, using literature (Springsteen and Leobrera accessory program. Haplotypes were counted using 1986, Abbott 1989), and by examination of reference DnaSP v. 6.12.01 (Rozas et al. 2017). collections in the University of the Philippines Diliman The substitution model was determined using ModelTest- Invertebrate Museum (UPDIM), Quezon City, Philippines. NG v.0.1.5 (Darriba et al. 2015), and the model with the best log-likelihood score was chosen using the Akaike DNA Extraction and Sequencing Information Criterion (Akaike 1973, 1974; Hurvich and DNA was extracted using the relatively rapid and Tsai 1993). The Xia Test (Xia et al. 2003, Xia and Lemey inexpensive modified NaOH-lysis method (Fontanilla et 2009) for substitution saturation was also performed in al. 2017). In this method, tissue slices were ground using DAMBE v. 6.4.81 (Xia 2013, 2017). Pairwise comparisons glass beads with 200 μL of 0.1 N NaOH and centrifuged of likelihood scores were performed on the dataset using with 300 μL chloroform-isoamyl alcohol (24:1). The uncorrected p-distances and GTR+Γ corrected distances upper phase was then collected and centrifuged with generated by PAUP* v4.10b (Swofford 2003), with ~300 μL isopropanol. The pellets were washed with GTR+Γ (Tavare 1986) determined as an optimal model ethanol, air-dried, and finally re-suspended in 150 by ModelTest-NG. Comparisons of pairwise distances μL TE buffer. DNA extracts were subjected to PCR 2 Special Issue on Genomics Batomalaque et al.: Philippine Helicostylinae DNA Barcodes Table 1. Species of Philippine camaenids collected with their corresponding GenBank accession numbers. Subfamily Species name Locality Collector/s GenBank accession no. Helicostylinae Chloraea amoena (Pfeiffer, 1845) Alaminos, I.K.C. Fontanilla, G.A. KM279469 Pangasinan, Luzon Batomalaque, E. de Vera KM279470 KM279471 KM279472 KM279473 Chloraea hennigiana Moellendorff, 1893 Alaminos, I.K.C. Fontanilla, G.A. KM279464 Pangasinan, Luzon Batomalaque, E. de Vera KM279465 KM279466 KM279467 KM279468 Chloraea fibula Not Uploaded to GenBank Chrysallis chrysalidiformis (Sowerby, 1833) Puerto Galera, A.U. Luczon KM056693 Mindoro KM056694 KM056695 Cochlostyla bicolorata (Lea, 1840) Laguna, Luzon C.P. Española KM056744 Cochlostyla daphnis (Broderip, 1841) Borbon, Cebu R.J.C. Canoy KM056706 Santander, Cebu P. Olvis KM056707 Cochlostyla fauna (Broderip, 1841) Bantayan Is., Cebu P. Olvis KM056713 KM056714 KM056715 Cochlostyla imperator (Pfeiffer, 1848) Sibulan Is., Polillo E.R.C. de Chavez KM056725 KM056726 Cochlostyla intermedia (Quadras and Moellendorff, 1896) Benguet, Luzon D. Constantino-Santos, KM056728 I.K.C. Fontanilla, A.U. Luczon KM056729 KM056730 KM056731 Cochlostyla marinduquensis (Hidalgo, 1887) Gasan, Marinduque B.O. Sosa III, R.D.C. KM279486 Pedales KM279487 KM279488 KM279489 KM279490 KM279491 KM279492 KM279493 KM279494