Phylogenetic Status of the Turkish Red Fox (Vulpes Vulpes), Based on Partial Sequences of the Mitochondrial Cytochrome B Gene
Total Page:16
File Type:pdf, Size:1020Kb
64 (2): 273 – 284 © Senckenberg Gesellschaft für Naturforschung, 2014. 25.7.2014 Phylogenetic Status of the Turkish Red Fox (Vulpes vulpes), based on Partial Sequences of the Mitochondrial Cytochrome b Gene Osman İbiş 1, Coşkun Tez 2, * & Servet Özcan 2, 3 1 Graduate School of Natural and Applied Sciences, Erciyes University, Kayseri 38039, Turkey — 2 Department of Biology, Faculty of Sciences, Erciyes University, Kayseri 38039, Turkey — 3 Genome and Stem Cell Center, GENKÖK, Erciyes University, Kayseri 38039, Turkey — * Corresponding author: tezc(at)erciyes.edu.tr Accepted 25.v.2014. Published online at www.senckenberg.de/vertebrate-zoology on 15.vii.2014. Abstract Genetic diversity and multiple mitochondrial phylogroups of the red fox have been revealed from scattered locations in previous studies. There is a still lack of information about the genetic diversity and phylogeographic structure of the red fox in Asia Minor. We investigated the genetic diversity in the Turkish red fox using a part of the cytochrome b mitochondrial gene (375 bp), and attempted to evaluate the phylogeographic structure in various geographic ranges of the species with the use of sequences available from the GenBank from various geographic origins and our data. Bayesian and Network analyses of the cytochrome b sequences from Turkey and GenBank suggested that the red fox is divided into four main phylogroups. They are grouped accordingly: Group 1 (SW Anatolia, Turkey and Hokkaido, Japan), Group 2 (Eurasia and North America), Group 3 (only North America), and Group 4 (Vietnam). The majority of Turkish haplotypes grouped with those of Eurasia. Despite the great distance between the localities, two haplotypes from SW Anatolia, Turkey grouped with previously reported haplotypes from Hokkaido, Japan. The present study shows that the Turkish red fox is nested within two main phylogroups and exhibits a high genetic diversity. Key words Vulpes vulpes, Mitochondrial DNA, Cytochrome b, Turkey. Introduction The red fox Vulpes vulpes is one of the most widely dis- FERNANDES et al. 2008, AUBRY et al. 2009, SACKS et al. tributed canid species throughout Northern Hemisphere 2010, 2011, STATHAM et al. 2011, 2012, TEACHER et al. (also introduced to Australia in the 1800’s – WILSON & 2011, EDWARDS et al. 2012, YU et al. 2012a), control re- REEDER 2005, MACDONALD & REYNOLDS 2008). During gion (INOUE et al. 2007, KIRSCHNING et al. 2007, AUBRY the past two decades, V. vulpes has been one of the most et al. 2009, SACKS et al. 2010, 2011, STATHAM et al. investigated mammals from various regions. There are 2011, 2012, TEACHER et al. 2011, EDWARDS et al. 2012, comprehensive studies dealing with the phylogeograph- KUTSCHERA et al. 2013, GALOV et al. 2014), tRNA-threo- ic structure and pattern of genetic diversity in the red nine and tRNA-proline (INOUE et al. 2007), microsatellite fox populations, using allozyme or isozyme (FRATI et data (LADE et al. 1996, WANDELER et al. 2003, KUKEKOVA al. 1998, SIMONSEN et al. 2003), RAPD (GACHOT-NEVEU et al. 2004, WANDELER & FUNK 2006; SACKS et al. 2010, et al. 2009), mitochondrial DNA cytochrome (Cyt) b 2011, OISHI et al. 2011) and single nucleotide polymor- (FRATI et al. 1998, INOUE et al. 2007, PERRINE et al. 2007, phism (SACKS et al. 2011). ISSN 1864-5755 273 O. İbiş et al.: Phylogenetic status of Turkish red fox Fig. 1. Collection localities and distribution of red foxes (Vulpes vulpes) in Turkey. Map numbers are locality codes (see Table 1). Multiple studies based on mitochondrial DNA of red genetic diversity of Turkish red fox have not been previ- fox emphasize that there has been a considerable impor- ously investigated. Therefore, there is a lack of informa- tance to have a comprehensive picture of genetic varia- tion in this respect. In reference to the red fox dentition tion in the entire range of the species, when new insights study, SZUMA (2008) suggested that the modern red fox are added into phylogeographic and colonization patterns evolved either in Asia Minor (Anatolia) or North Africa of red fox with a wider biogeographic context (FRATI during the early Pleistocene. Consequently, to confirm et al. 1998, INOUE et al. 2007, PERRINE et al. 2007, AUBRY this, the phylogenetic status of Tur kish red fox should be et al. 2009, TEACHER et al. 2011, EDWARDS et al. 2012, YU revealed using different markers. et al. 2012a, KUTSCHERA et al. 2013). Some of the studies In this study, our primary purpose was to determine resulted in the formation of mitochondrial phylogroups, the magnitude of genetic diversity within the Turkish red which were classified as Holoarctic and/or Eurasian, fox by sequencing part of the Cyt b gene. With emphasis Nearctic, and Hokkaido (INOUE et al. 2007, AUBRY et al. on this information we provide new insights into phylo- 2009, YU et al. 2012a, KUTSCHERA et al. 2013). Although geographic processes which may be essential for a better there is a substantial inventory of mitochondrial data understanding of evolution. from Eurasia and North America (FRATI et al. 1998, INOUE et al. 2007, PERRINE et al. 2007, AUBRY et al. 2009, TEACHER et al. 2011, EDWARDS et al. 2012, YU et al. 2012, KUTSCHERA et al. 2013), there is a still lack of information Material and Methods about the genetic diversity and phylogeographic structure of the red fox from Turkey and the Middle East, except Israel. In this sense, more genetic and phylogeographic Localization, tissue sampling and DNA studies are needed to place these new mitochondrial data extraction in a wider context and to elucidate the current pattern of phylogeographic structure throughout the red fox’ range. Various tissues (ear, tail, skin and skeletal muscle) of 51 Therefore, it would not only be beneficial to determine the red foxes, which were road kills, were collected from genetic structure of different red fox populations but also 51 different localities in Turkey (Table 1, Fig. 1) during to focus on regional details, especially those in previously 2002-2010. Tissue samples had been stored in the cor- unsampled areas. Turkey is one of such areas. Because of responding author’s private collection at the Department its continental position, Turkey constitutes a mysterious of Biology, Faculty of Sciences, Erciyes University in and interesting area for mammalian biogeography due Kayseri, Turkey. to its climate, vegetation, heterogeneous topography and Total genomic DNA was extracted by using the geological history; together these unique factors produce DNeasy Blood and Tissue Kit (QIAGEN), following a biodiversity hotspot (MÉDAIL & QUÉZEL 1999, BILGIN the manufacturer’s instructions. A segment of the mito - 2011). Currently, there are at least 150 or more mamma- chondrial cytochrome (Cyt) b (375 base pairs) from each lian species from Turkey (KRYšTUFEK & VOHRALIK 2001, Turkish sample was amplified by polymerase chain reac- 2009, WILSON & REEDER 2005). Although there are a few tion (PCR) using the primer pair studies based on its morphology and distribution (ÖZKURT L14724 (5’ GATATGAAAAACCATCGTTG 3’) and et al. 1998, TEMIZER 2001), the phylogenetic status and H15149 (5’ CAGAATGATATTTGTCCTCA 3’) 274 VERTEBRATE ZOOLOGY — 64 (2) 2014 Table 1. List of Turkish samples and sequences/haplotypes used in the present study (* Map numbers are locality codes in Figure 1). Map Cyt b Tissue Locality Coordinate/Altitude number* haplotype 1 TR1 Tail Özlüce/Boğazlıyan/Yozgat N 39º 29.185’ E 35º 04.298’ H – 1099m 2 TR1 Muscle Boğazkale/Çorum N 39º 56.087’ E 34º 42.963’ H – 1215m 3 TR1 Tail Kemerhisar/Bor/Niğde N 37º 45.509’ E 34º 34.032’ H – 1128m 4 TR1 Tail Karabük N 41º 09.810’ E 32º 38.428’ H – 304m 5 TR1 Tail Bozca/Avanos/Nevşehir N 38º 44.331’ E 35º 02.312’ H – 1048m 6 TR1 Tail Avanos/Nevşehir N 38º 41.409’ E 34º 48.850’ H – 1012m 7 TR1 Tail Konya N 38º 01.340’ E 32º 51.829’ H – 1001m 8 TR1 Tail Konya N 37º 53.613’ E 32º 09.022’ H – 1565m 9 TR1 Tail Yozgat N 39º 48.504’ E 34º 46.330’ H – 1204m 10 TR1 Tail Tokat N 40º 03.833’ E 36º 29.187’ H – 1143m 11 TR1 Ear Çamlıbel/Tokat N 40º 04.261’ E 36º 29.011’ H – 1132m 12 TR1 Tongue Amasya N 40º 41.779’ E 35º 48.765’ H – 408m 13 TR1 Ear Bekdiğin/Havza /Samsun N 40º 59.861’ E 35º 46.125’ H – 725m 14 TR1 Muscle Eğirdir/Isparta N 37º 52.000’ E 30º 15.000’ 15 TR1 Tubercle Acıgöl/Nevşehir N 38º 32.691’ E 34º 29.174’ H – 1264m 16 TR1 Ear Aksaray N 38º 30.969’ E 34º 21.098’ H – 1227m 17 TR1 Ear Kayaönü Mah./Karaman N 37º 14.318’ E 33º 21.128’ H – 1039m 18 TR1 Ear Kirşehir N 39º 05.094’ E 34º 16.287’ H – 1200m 19 TR1 Ear Elmadağ/Ankara N 39º 55.348’ E 33º 15.923’ H – 1072m 20 TR1 Ear Eskişehir N 39º 30.058’ E 31º 13.424’ H – 1028m 21 TR1 Muscle Ilgın/Konya N 38º 18.062’ E 31º 39.356’ H – 1080m 22 TR1 Ear Eşmekaya/Aksaray N 38º 15.152’ E 33º 28.976’ H – 971m 23 TR1 Ear Çalış/Nevşehir N 38º 58.846’ E 34º 54.890’ H – 1316m 24 TR1 Ear Pınarbaşı/Kayseri N 38º 43.000’ E 36º 23.000’ 25 TR2 Tail Kızılırmak/Çankiri N 40º 21.169’ E 33º 58.767’ H – 555m 26 TR2 Tail Yeniyapan/Keskin/Kirikkale N 39º 35.319’ E 33º 42.013’ H – 922m 27 TR3 Tail Kandıra/Kocaeli N 41º 04.907’ E 30º 10.238’ H – 41m 28 TR4 Tail Isparta N 39º 03.203’ E 31º 25.615’ H – 1200m 29 TR4 Tail Dikici Köyü/Dinar/Afyon N 38º 00.140’ E 30º 11.522’ H – 940m 30 TR4 Ear Beyşehir/Konya N 37º 54.691’ E 31º 56.955’ H – 1270m 31 TR5 Tail Bahadin/Sarıkaya/Yozgat N 39º 36.066’ E 35º 15.707’ H – 1123m 32 TR5 Tail Isparta N 37º 52.539’ E 30º 30.730’ H – 935m 33 TR5 Tail Eğirdir/Isparta N 37º 50.753’ E 30º 51.579’ H – 922m 34 TR5 Muscle Arguvan/Malatya N 38º 46.000’ E 38º 16.000’ 35 TR5 Tail Taşköprü/Babaeski/Kirklareli