Interspecies introgressive hybridization in spiny frogs Quasipaa (Family Dicroglossidae) revealed by analyses on multiple mitochondrial and nuclear genes

Introgression may lead to discordant patterns of variation among loci and traits. For example, previous phylogeographic studies on the genus Quasipaa detected signs of genetic introgression from genetically and morphologically divergent Quasipaa shini or Quasipaa spinosa. In this study, we used mitochondrial and nuclear DNA sequence data to verify the widespread introgressive hybridization in the closely related species of the genus Quasipaa, evaluate the level of genetic diversity, and reveal the formation mechanism of introgressive hybridization. In Longsheng, Guangxi Province, signs of asymmetrical nuclear introgression were detected between Quasipaa boulengeri and Q. shini. Unidirectional mitochondrial introgression was revealed from Q. spinosa to Q. shini. By contrast, bidirectional mitochondrial gene introgression was detected between Q. spinosa and Q. shini in Lushan, Jiangxi Province. Our study also detected ancient hybridizations between a female Q. spinosa and a male Q. jiulongensis in Zhejiang Province. Analyses on mitochondrial and nuclear genes verified three candidate cryptic species in Q. spinosa, and a cryptic species may also exist in Q. boulengeri. However, no evidence of introgressive hybridization was found between Q. spinosa and Q. boulengeri. Quasipaa exilispinosa from all the sampling localities appeared to be deeply divergent from other communities. Our results suggest widespread introgressive hybridization in closely related species of Quasipaa and provide a fundamental basis for illumination of the forming mechanism of introgressive hybridization, classification of species, and biodiversity assessment in Quasipaa.     

Cryptic genetic diversity in Rattus of the San Francisco Bay region, California

Invasive species can have complex invasion histories, harbor cryptic levels of diversity, and pose taxonomic problems for pest management authorities. Roof rats, Rattus rattus sensu lato, are common invasive pests of the San Francisco Bay Area in California, USA. They are a significant health risk and pest management efforts impose a large financial investment from public institutions and private individuals. Recent molecular genetic and taxonomic studies of black rats in their native range in Asia have shown that the species is a complex of two karyotypic forms and four mitochondrial genetic lineages that may represent four distinct species. We used mtDNA sequences and nuclear microsatellite variation to identify which mitochondrial lineages of the R. rattus group are present in the San Francisco Bay Area and to test for gene flow among them. We recovered specimens with mtDNA sequences representing two of the major mtDNA lineages of the R. rattus group. Microsatellite variation, however, was not structured in concordance with mtDNA lineages, suggesting a more complex history involving hybridization and introgression between these lineages. Although Aplin et al. (2011) and Lack et al. (2012) reported R. rattus Lineage II in North America, this is the first detailed examination of possible gene flow amongst lineages in this region.     

Cryptic variation and the tragedy of unrecognized taxa: the case of international trade in the spiny turtle Heosemys spinosa (Testudines: Geoemydidae)

Loss of habitat and human exploitation have driven many turtle species to the brink of extinction, particularly in many parts of southern Asia. The spiny turtle (Heosemys spinosa) is a terrestrial species distributed throughout the Sundaland region of South‐East Asia. Despite international legislative protection, H. spinosa continues to be illegally collected for the food and traditional medicine markets of China. Given its widespread distribution, taxonomists have reasonably questioned whether H. spinosa truly represents a single evolutionary lineage or multiple undiagnosed species. Recently, a large and illegal shipment of rare, wild‐caught H. spinosa was confiscated in Hong Kong, China, and the turtles were eventually distributed to several zoos and academic collections. Based on analyses of these individuals, along with additional individuals from the pet trade and museum collections, we found concordant genetic and phenotypic variation, indicating that two distinct types of H. spinosa exist in this collection of turtles. Further characterization of this variation will require field surveys and the collection of additional morphological and genetic data from specimens of known geographic provenance. However, our data indicate that this highly exploited, endangered species may contain additional cryptic taxa, and emphasize the critical need for systematic evaluation of species before unrecognized variation is lost forever. © 2012 The Linnean Society of London, Zoological Journal of the Linnean Society, 2012, 164 , 811–824.     

The mitochondrial 16 s rRNA reveals high anthropogenic influence on land snail diversity in a preliminary island survey

A total of 30 specimens belonging to five species, namely; Cryptozona siamensis, Sarika resplendens and Sarika sp. from the family Ariophantidae as well as Quantula striata and Quantula sp. from the family Dyakiidae were collected from the Langkawi Island in Northern Peninsular Malaysia. All specimens were identified through comparisons of shell morphology and amplification of a 500 bp segment of the 16S rRNA mtDNA gene. To assess phylogenetic insights, the sequences were aligned using ClustalW and phylogenetic trees were constructed. The analyses showed two major lineages in both Maximum Parsimony and Neighbour Joining phylogenetic trees. Each putative taxonomic group formed a monophyletic cluster. Our study revealed low species and intraspecies genetic diversities based on the 16S rRNA gene sequences. Thus, this study has provided an insight of land snail diversity in populations of an island highly influenced by anthropogenic activities through complementary use of shell morphological and molecular data.     

A phylogenetic approach to species delimitation in freshwater Gastrotricha from Sweden

Gastrotricha is a cosmopolitan group of aquatic invertebrates. To date, approximately 765 species have been described. This study is the first to deal with species delimitation and cryptic species of freshwater Gastrotricha. Three commonly encountered species, Heterolepidoderma ocellatum, Lepidochaetus zelinkai, and Lepidodermella squamata, are investigated for cryptic speciation. Most of the material is based on Swedish specimens but closely related species from other parts of the world are also included. Taxonomic revisions are supported by phylogenies based on 18S rRNA, 28S rRNA, and COI mtDNA of freshwater Chaetonotidae from several genera and inferred from Bayesian and maximum likelihood approaches. Heterolepidoderma ocellatum f. sphagnophilum is raised to species level, under the name H. acidophilum n. sp. Moreover, genetic data based on COI indicate large variation between two morphologically very similar groups of Lepidodermella squamata. The extent of cryptic speciation in L. zelinkai appears low. Based on the phylogenetic hypothesis presented in this article, the new species, Lepidodermella intermedia n. sp., from northern Sweden is also described. The phylogenetic hypothesis generated shows that Chaetonotidae is a non-monophyletic group.     

A mirage of cryptic species: Genomics uncover striking mitonuclear discordance in the butterfly Thymelicus sylvestris

Mitochondrial DNA (mtDNA) sequencing has led to an unprecedented rise in the identification of cryptic species. However, it is widely acknowledged that nuclear DNA (nuDNA) sequence data are also necessary to properly define species boundaries. Next generation sequencing techniques provide a wealth of nuclear genomic data, which can be used to ascertain both the evolutionary history and taxonomic status of putative cryptic species. Here, we focus on the intriguing case of the butterfly Thymelicus sylvestris (Lepidoptera: Hesperiidae). We identified six deeply diverged mitochondrial lineages; three distributed all across Europe and found in sympatry, suggesting a potential case of cryptic species. We then sequenced these six lineages using double‐digest restriction‐site associated DNA sequencing (ddRADseq). Nuclear genomic loci contradicted mtDNA patterns and genotypes generally clustered according to geography, i.e., a pattern expected under the assumption of postglacial recolonization from different refugia. Further analyses indicated that this strong mtDNA/nuDNA discrepancy cannot be explained by incomplete lineage sorting, sex‐biased asymmetries, NUMTs, natural selection, introgression or Wolbachia‐mediated genetic sweeps. We suggest that this mitonuclear discordance was caused by long periods of geographic isolation followed by range expansions, homogenizing the nuclear but not the mitochondrial genome. These results highlight T. sylvestris as a potential case of multiple despeciation and/or lineage fusion events. We finally argue, since mtDNA and nuDNA do not necessarily follow the same mechanisms of evolution, their respective evolutionary history reflects complementary aspects of past demographic and biogeographic events.     

Deep sympatric mtDNA divergence in the autumnal moth (Epirrita autumnata)

Deep sympatric intraspecific divergence in mtDNA may reflect cryptic species or formerly distinct lineages in the process of remerging. Preliminary results from DNA barcoding of Scandinavian butterflies and moths showed high intraspecific sequence variation in the autumnal moth, Epirrita autumnata. In this study, specimens from different localities in Norway and some samples from Finland and Scotland, with two congeneric species as outgroups, were sequenced with mitochondrial and nuclear markers to resolve the discrepancy found between mtDNA divergence and present species‐level taxonomy. We found five COI sub‐clades within the E. autumnata complex, most of which were sympatric and with little geographic structure. Nuclear markers (ITS2 and Wingless) showed little variation and gave no indications that E. autumnata comprises more than one species. The samples were screened with primers for Wolbachia outer surface gene (wsp) and 12% of the samples tested positive. Two Wolbachia strains were associated with different mtDNA sub‐clades within E. autumnata, which may indicate indirect selection/selective sweeps on haplotypes. Our results demonstrate that deep mtDNA divergences are not synonymous with cryptic speciation and this has important implications for the use of mtDNA in species delimitation, like in DNA barcoding.     

Cryptic speciation in the field vole: a multilocus approach confirms three highly divergent lineages in Eurasia

Species are generally described from morphological features, but there is growing recognition of sister forms that show substantial genetic differentiation without obvious morphological variation and may therefore be considered ‘cryptic species’. Here, we investigate the field vole (Microtus agrestis), a Eurasian mammal with little apparent morphological differentiation but which, on the basis of previous sex‐linked nuclear and mitochondrial DNA (mtDNA) analyses, is subdivided into a Northern and a Southern lineage, sufficiently divergent that they may represent two cryptic species. These earlier studies also provided limited evidence for two major mtDNA lineages within Iberia. In our present study, we extend these findings through a multilocus approach. We sampled 163 individuals from 46 localities, mainly in Iberia, and sequenced seven loci, maternally, paternally and biparentally inherited. Our results show that the mtDNA lineage identified in Portugal is indeed a distinct third lineage on the basis of other markers as well. In fact, multilocus coalescent‐based methods clearly support three separate evolutionary units that may represent cryptic species: Northern, Southern and Portuguese. Divergence among these units was inferred to have occurred during the last glacial period; the Portuguese lineage split occurred first (estimated at c. 70 000 bp ), and the Northern and Southern lineages separated at around the last glacial maximum (estimated at c. 18 500 bp ). Such recent formation of evolutionary units that might be considered species has repercussions in terms of understanding evolutionary processes and the diversity of small mammals in a European context.     

Species boundaries in the Himantura uarnak species complex (Myliobatiformes: Dasyatidae)

Samples of the ‘Himantura uarnak’ species complex (H. leoparda, H. uarnak, H. undulata under their current definitions), mostly from the Coral Triangle, were analyzed using nuclear markers and mitochondrial DNA sequences. Genotypes at five intron loci showed four reproductively isolated clusters of individuals. The COI sequences showed four major mitochondrial lineages, each diagnostic of a cluster as defined by nuclear markers. No mitochondrial introgression was detected. The average Kimura-2 parameter nucleotide distance separating clades was 0.061–0.120 (net: 0.055–0.114), while the distance separating individuals within a clade was 0.002–0.008. Additional, partial cytochrome-b gene sequences were used to link these samples with previously published sequences of reference specimens of the three nominal species. One of the clusters was identified as H. undulata and another one, as H. uarnak, while two cryptic species were uncovered within the recently-described H. leoparda, challenging the current morphology-based taxonomy of species within the H. uarnak species complex.     

Phylogeography, geographic structure, genetic variation, and potential species boundaries in Philippine slender toads

We investigated phylogeography of Philippine slender toads (genus Ansonia) and used a temporal framework for diversification, statistical tests of alternate topologies, and Bayesian approaches to test previous hypotheses concerning dispersal to, and colonization routes within, the southern Philippine island of Mindanao. Two species of Ansonia previously have been documented, with ranges separated by an east–west split corresponding to the approximate boundaries of Mindanao’s paleoisland precursors. We present new mtDNA sequence data (1946 bp from genes encoding ND1, 16S rRNA and tRNA^Leu) for 105 Ansonia specimens sampled from 20 localities on Mindanao Island. Our data suggest that Philippine Ansonia is composed of at least eight, well-supported population lineages, structured into a minimum of four highly divergent mtDNA clades. One clade corresponds to Ansonia mcgregori, a range-restricted species apparently limited to the distal portion of the Zamboanga Peninsula of western Mindanao. Two morphologically indistinguishable, but genetically divergent, lineages possibly are undescribed cryptic species from western Mindanao. We recognize the five remaining lineages as Ansonia muelleri pending data from morphology or bioacoustics that might diagnose separate species among these lineages. Regardless of their species status, the five allopatric lineages of A. muelleri should be viewed as important genetic units for future genetic conservation planning.     

Characterization of a Plasmopara species on Ambrosia artemisiifolia, and notes on P. halstedii, based on morphology and multiple gene phylogenies

Common ragweed (Ambrosia artemisiifolia) is an invasive and highly allergenic plant species, on which two species, Plasmopara halstedii and Plasmopara angustiterminalis, have been recognized to cause downy mildew disease. In this study, morphological and molecular patterns of seven Plasmopara specimens collected from A. artemisiifolia in Canada, Hungary, and USA were compared with those of P. halstedii and P. angustiterminalis from Helianthus and Xanthium, respectively. Analyses of partial sequences of three genes, namely those for the large subunit (28S) of rDNA, cytochrome c oxidase subunit II (COX2), and NADH dehydrogenase subunit I (ND1) of mtDNA, were carried out to examine the phylogenetic relationships among these specimens using both Bayesian and maximum parsimony methods. All the phylogenetic analyses revealed that the downy mildew pathogens infecting A. artemisiifolia in Hungary and North America clearly represent a lineage distinct from other Plasmopara taxa investigated. The shape of sporangia and the width of trunks and branches also allowed the separation of the specimens parasitic to A. artemisiifolia from P. halstedii on Helianthus annuus and P. angustiterminalis on Xanthium strumarium. Surprisingly, the Hungarian and the Canadian specimens were more closely related to each other than to those from the USA based on COX2 and ND1 mtDNA data, although the D1/D2/D3 sequences of 28S rDNA were identical in all these Plasmopara specimens. The regional distribution of the mtDNA haplotypes seen in this study suggests a transatlantic migration has occurred and would be interesting to follow up with a more detailed sampling. To investigate the diversity within P. halstedii sensu lato, infecting different host plant species, specimens from six asteraceous genera, Ambrosia, Flaveria, Helianthus, Siegesbeckia, Solidago, and Xanthium, were also included in molecular analyses. These represented six distinct lineages according to the host plant genera. These findings might serve as a basis for a taxonomical reassessment of the P. halstedii complex and also for the delimitation of several well-defined species within this complex.     

Molecular phylogeny of Diplostomum,Tylodelphys, Austrodiplostomum and Paralaria (Digenea: Diplostomidae) necessitates systematic changes and reveals a history of evolutionary host switching events

The Diplostomidae Poirier, 1886 is a large, globally distributed family of digeneans parasitic in intestines of their definitive hosts. Diplostomum and Tylodelphys spp. are broadly distributed, commonly reported, and the most often sequenced diplostomid genera. The majority of published DNA sequences from these genera originated from larval stages only, which typically cannot be identified to the species level based on morphology alone. We generated partial large ribosomal subunit (28S) rRNA and cytochrome c oxidase subunit 1 (cox1) mtDNA gene sequences from 14 species/species-level lineages of Diplostomum, six species/species-level lineages of Tylodelphys, two species/species-level lineages of Austrodiplostomum, one species previously assigned to Paralaria, two species/species-level lineages of Dolichorchis and one unknown diplostomid. Our DNA sequences of 11 species/species-level lineages of Diplostomum (all identified to species), four species/species-level lineages of Tylodelphys (all identified to species), Austrodiplostomum compactum, Paralaria alarioides and Dolichorchis lacombeensis originated from adult specimens. 28S sequences were used for phylogenetic inference to demonstrate the position of Paralaria alarioides and Dolichorchis spp. within the Diplostomoidea and study the interrelationships of Diplostomum, Tylodelphys and Austrodiplostomum. Our results demonstrate that two diplostomids from the North American river otter (P. alarioides and a likely undescribed taxon) belong within Diplostomum. Further, our results demonstrate the non-monophyly of Tylodelphys due to the position of Austrodiplostomum spp., based on our phylogenetic analyses and morphology. Furthermore, the results of phylogenetic analysis of 28S confirmed the status of Dolichorchis as a separate genus. The phylogenies suggest multiple definitive host-switching events (birds to otters and among major avian groups) and a New World origin of Diplostomum and Tylodelphys spp. Our DNA sequences from adult digeneans revealed identities of 10 previously published lineages of Diplostomum and Tylodelphys, which were previously identified to genus only. The novel DNA data from this work provide opportunities for future comparisons of larval diplostomines collected in ecological studies.     

Vicariance and Its Impact on the Molecular Ecology of a Chinese Ranid Frog Species-Complex (Odorrana schmackeri,Ranidae)

Paleogeological events and Pleistocene climatic fluctuations have had profound influences on the genetic patterns and phylogeographic structure of species in southern China. In this study, we investigated the population genetic structure and Phylogeography of the Odorrana schmackeri species complex, mountain stream-dwelling odorous frogs, endemic to southern China. We obtained mitochondrial sequences (1,151bp) of the complete ND2 gene and two flanking tRNAs of 511 individuals from 25 sites for phylogeographic analyses. Phylogenetic reconstruction revealed seven divergent evolutionary lineages, with mean pairwise (K2P) sequence distances from 7.8% to 21.1%, except for a closer ND2 distance (3.4%). The complex geological history of southern China drove matrilineal divergence in the O. schmackeri species complex into highly structured geographical units. The first divergence between lineage A+B and other lineages (C-G) had likely been influenced by the uplift of coastal mountains of Southeast China during the Mio-Pliocene period. The subsequent divergences between the lineages C-G may have followed the formation of the Three Gorges and the intensification of the East Asian summer monsoon during the late Pliocene and early Pleistocene. Demographic analyses indicated that major lineages A and C have been experienced recent population expansion (c. 0.045–0.245 Ma) from multiple refugia prior to the Last Glacial Maximum (LGM). Molecular analysis suggest that these seven lineages may represent seven different species, three described species and four cryptic species and should at least be separated into seven management units corresponding to these seven geographic lineages for conservation.     

Cryptic speciation in the Anadenobolus excisus millipede species complex on the Island of Jamaica

Anadenobolus excisus is a large species of millipede endemic to the Caribbean Island of Jamaica. Initial detailed morphological studies showed little or no discrete variation across this species' distribution in somatic or, in particular. genitalic morphology. However, a molecular survey based on approximately 1000 base pairs of the mitochondrial (mtDNA) 16S rRNA gene that examines 242 individuals sampled from 54 localities reveals three highly divergent mtDNA lineages. A lack of discrete morphological differentiation suggests that genetic and morphological divergence may be decoupled, a pattern inconsistent with a number of evolutionary models. In contrast to minimal morphological divergence, size variation among mtDNA lineages suggests that character displacement has occurred and that these lineages are cohesive in sympatry. We conclude that A. excisus is actually a complex of three cryptic species and that morphological approaches to delineating millipede species may sometimes underestimate evolutionary diversity.     

The Role of Integrative Taxonomy in the Conservation Management of Cryptic Species: The Taxonomic Status of Endangered Earless Dragons (Agamidae: Tympanocryptis) in the Grasslands of Queensland,Australia

Molecular phylogenetics is increasingly highlighting the prevalence of cryptic species, where morphologically similar organisms have long independent evolutionary histories. When such cryptic species are known to be declining in numbers and are at risk of extinction due to a range of threatening processes, the disjunction between molecular systematics research and conservation policy becomes a significant problem. We investigate the taxonomic status of Tympanocryptis populations in Queensland, which have previously been assigned to T. tetraporophora, using three species delimitation approaches. The taxonomic uncertainties in this species-group are of particular importance in the Darling Downs Earless Dragon (T. cf. tetraporophora), which is ranked as an endangered ‘species’ of high priority for conservation by the Queensland Department of Environment and Heritage Protection. We undertook a morphological study, integrated with a comprehensive genetic study and species delimitation analyses, to investigate the species status of populations in the region. Phylogenetic analyses of two gene regions (mtDNA: ND2; nuclear: RAG1) revealed high levels of genetic divergence between populations, indicating isolation over long evolutionary time frames, and strongly supporting two independent evolutionary lineages in southeastern Queensland, from the Darling Downs, and a third in the Gulf Region of northern Queensland. Of the three species delimitation protocols used, we found integrative taxonomy the most applicable to this cryptic species complex. Our study demonstrates the utility of integrative taxonomy as a species delimitation approach in cryptic complexes of species with conservation significance, where limited numbers of specimens are available.     

Evidence for Cryptic Speciation in Directly Transmitted Gyrodactylid Parasites of Trinidadian Guppies

Cryptic species complexes are common among parasites, which tend to have large populations and are subject to rapid evolution. Such complexes may arise through host-parasite co-evolution and/or host switching. For parasites that reproduce directly on their host, there might be increased opportunities for sympatric speciation, either by exploiting different hosts or different micro-habitats within the same host. The genus Gyrodactylus is a specious group of viviparous monogeneans. These ectoparasites transfer between teleosts during social contact and cause significant host mortality. Their impact on the guppy (Poecilia reticulata), an iconic evolutionary and ecological model species, is well established and yet the population genetics and phylogenetics of these parasites remains understudied. Using mtDNA sequencing of the host and its parasites, we provide evidence of cryptic speciation in Gyrodactylus bullatarudis, G. poeciliae and G. turnbulli. For the COII gene, genetic divergence of lineages within each parasite species ranged between 5.7 and 17.2%, which is typical of the divergence observed between described species in this genus. Different lineages of G. turnbulli and G. poeciliae appear geographically isolated, which could imply allopatric speciation. In addition, for G. poeciliae, co-evolution with a different host species cannot be discarded due to its host range. This parasite was originally described on P. caucana, but for the first time here it is also recorded on the guppy. The two cryptic lineages of G. bullatarudis showed considerable geographic overlap. G. bullatarudis has a known wide host range and it can also utilize a killifish (Anablepsoides hartii) as a temporary host. This killifish is capable of migrating overland and it could act as a transmission vector between otherwise isolated populations. Additional genetic markers are needed to confirm the presence of these cryptic Gyrodactylus species complexes, potentially leading to more in-depth genetic, ecological and evolutionary analyses on this multi-host-parasite system.     

Mitochondrial DNA diversity and taxa delineation in the land snails of the Iberus gualtieranus (Pulmonata,Helicidae) complex

Iberus gualtieranus is a species complex of land snails that is endemic to the Iberian Peninsula. The species taxonomy of the group is based merely on the basis of shell morphology, but validity of the existing taxonomy is uncertain. Using mitochondrial DNA (mtDNA) data (cytochrome oxidase I and 16S rRNA sequences) we were able to validate the observed phylogenetic taxa within the I. gualtieranus s.l. complex by means of the analysis of specimens of the different morphospecies, together with the study of topotypes. Strong incongruences were obtained between morphology and molecular data. The Iberus alonensis morphospecies comprised several genetically divergent but morphologically cryptic lineages. Considering (1) the allopatric distribution of the operational taxonomic units (OTUs), (2) the morphological differentiation, (3) the possible occurrence of hybridization among the different lineages, and (4) the strong differentiation of the mtDNA phylogroups, we suggest the main lineages obtained, for the time being, may be treated as evolutionary species. The robust phylogenetic reconstruction obtained allows us to consider I. alonensis s.s., Iberus campesinus, Iberus carthaginiensis, and Iberus gualtieranus s.s. as valid species. Two additional unnominated taxa of the alonensis shell type have also been identified. Further subdivisions are also considered, including Iberus gualtieranus mariae and Iberus gualtieranus ornatissimus. © 2008 The Linnean Society of London, Zoological Journal of the Linnean Society, 2008, 154 , 722–737.     

Implications of Hybridization,NUMTs, and Overlooked Diversity for DNA Barcoding of Eurasian Ground Squirrels

The utility of DNA Barcoding for species identification and discovery has catalyzed a concerted effort to build the global reference library; however, many animal groups of economical or conservational importance remain poorly represented. This study aims to contribute DNA barcode records for all ground squirrel species (Xerinae, Sciuridae, Rodentia) inhabiting Eurasia and to test efficiency of this approach for species discrimination. Cytochrome c oxidase subunit 1 (COI) gene sequences were obtained for 97 individuals representing 16 ground squirrel species of which 12 were correctly identified. Taxonomic allocation of some specimens within four species was complicated by geographically restricted mtDNA introgression. Exclusion of individuals with introgressed mtDNA allowed reaching a 91.6% identification success rate. Significant COI divergence (3.5–4.4%) was observed within the most widespread ground squirrel species (Spermophilus erythrogenys, S. pygmaeus, S. suslicus, Urocitellus undulatus), suggesting the presence of cryptic species. A single putative NUMT (nuclear mitochondrial pseudogene) sequence was recovered during molecular analysis; mitochondrial COI from this sample was amplified following re-extraction of DNA. Our data show high discrimination ability of 100 bp COI fragments for Eurasian ground squirrels (84.3%) with no incorrect assessments, underscoring the potential utility of the existing reference librariy for the development of diagnostic ‘mini-barcodes’.