Universal COI primers for DNA barcoding amphibians

Mountains too high and valleys too deep drive population structuring and demographics in a Qinghai-Tibetan Plateau frog Nanorana pleskei (Dicroglossidae)

Pleistocene glacial–interglacial climatic oscillations greatly shaped the current genetic structure of many species. However, geographic features may influence the impact of climatic cycling. Distinct geographic and environmental characters between northern and southern parts of the eastern Qinghai–Tibetan Plateau (EQTP) facilitate explorations into the impacts of geographic features on species. The northern parts of EQTP contain large areas of marsh, and the environment is rather homogeneous. In contrast, the southern EQTP harbors complex alpine valleys and a much more heterogeneous setting. We evaluate DNA sequence variation from both the mitochondrial and nuclear genomes in Nanorana pleskei, a species endemic to the EQTP. Hypothesis testing on the evolutionary history of N. pleskei indicates that northern populations can disperse freely, but alpine valleys isolate southern populations. Demographic histories between northern and southern populations also differ. Northern populations appear to have experienced population expansions, while southern frogs exhibit a far more stable demographic history. By combining climatic analyses and species' distribution models, our study suggests that geographic and environmental features drive the differences between the northern and southern EQTP.

Evaluating the diversity of Neotropical anurans using DNA barcodes

This study tested the effectiveness of COI barcodes for the discrimination of anuran species from the Amazon basin and other Neotropical regions. Barcodes were determined for a total of 59 species, with a further 58 species being included from GenBank. In most cases, distinguishing species using the barcodes was straightforward. Each species had a distinct COI barcode or codes, with intraspecific distances ranging from 0% to 9.9%. However, relatively high intraspecific divergence (11.4-19.4%) was observed in some species, such as Ranitomeya ventrimaculata, Craugastor fitzingeri, Hypsiboas leptolineatus, Scinax fuscomarginatus and Leptodactylus knudseni, which may reflect errors of identification or the presence of a species complex. Intraspecific distances recorded in species for which samples were obtained from GenBank (Engystomops pustulosus, Atelopus varius, Craugastor podiciferus, and Dendropsophus labialis) were greater than those between many pairs of species. Interspecific distances ranged between 11-39%. Overall, the clear differences observed between most intra- and inter-specific distances indicate that the COI barcode is an effective tool for the identification of Neotropical species in most of the cases analyzed in the present study.

Complete mitochondrial genome of Chinese salamander Batrachuperus pinchonii: genome characterization and phylogenetic analysis

The complete mitogenome sequence of Batrachuperus pinchonii (Caudata: Hynobiidae) has been amplified and sequenced in this study. The overall base composition of B. pinchonii mitogenome is 33.9% for A, 19.7% for C, 13.6% for G, and 32.8% for T and has low GC content of 33.3%. The assembled mitogenome, consisting of 16,381 bp, has 13 protein-coding genes (PCGs), 22 tRNA genes, 2 rRNA genes, and 1 D-loop region. All PCGs, except for ND6 gene, were encoded on H-strand. Phylogenetic analysis with the whole mitogenome sequences revealed a close relationship of B. pinchonii with B. tibetanus.

Improving herpetological surveys in eastern North America using the environmental DNA method

Among vertebrates, herpetofauna has the highest proportion of declining species. Detection of environmental DNA (eDNA) is a promising method towards significantly increasing large-scale herpetological conservation efforts. However, the integration of eDNA results within a management framework requires an evaluation of the efficiency of the method in large natural environments and the calibration of eDNA surveys with the quantitative monitoring tools currently used by conservation biologists. Towards this end, we first developed species-specific primers to detect the wood turtle (Glyptemys insculpta) a species at risk in Canada, by quantitative PCR (qPCR). The rate of eDNA detection obtained by qPCR was also compared to the relative abundance of this species in nine rivers obtained by standardized visual surveys in the Province of Québec (Canada). Second, we developed multi-species primers to detect North American amphibian and reptile species using eDNA metabarcoding analysis. An occurrence index based on the distribution range and habitat type was compared with the eDNA metabarcoding dataset from samples collected in seven lakes and five rivers. Our results empirically support the effectiveness of eDNA metabarcoding to characterize herpetological species distributions. Moreover, detection rates provided similar results to standardized visual surveys currently used to develop conservation strategies for the wood turtle. We conclude that eDNA detection rates may provide an effective semiquantitative survey tool, provided that assay calibration and standardization is performed.

DNA Barcoding of the Endangered Aquilaria (Thymelaeaceae) and Its Application in Species Authentication of Agarwood Products Traded in the Market

The identification of Aquilaria species from their resinous non-wood product, the agarwood, is challenging as conventional techniques alone are unable to ascertain the species origin. Aquilaria is a highly protected species due to the excessive exploitation of its precious agarwood. Here, we applied the DNA barcoding technique to generate barcode sequences for Aquilaria species and later applied the barcodes to identify the source species of agarwood found in the market. We developed a reference DNA barcode library using eight candidate barcode loci (matK, rbcL, rpoB, rpoC1, psbA-trnH, trnL-trnF, ITS, and ITS2) amplified from 24 leaf accessions of seven Aquilaria species obtained from living trees. Our results indicated that all single barcodes can be easily amplified and sequenced with the selected primers. The combination of trnL-trnF+ITS and trnL-trnF+ITS2 yielded the greatest species resolution using the least number of loci combination, while matK+trnL-trnF+ITS showed potential in detecting the geographical origins of Aquilaria species. We propose trnL-trnF+ITS2 as the best candidate barcode for Aquilaria as ITS2 has a shorter sequence length compared to ITS, which eases PCR amplification especially when using degraded DNA samples such as those extracted from processed agarwood products. A blind test conducted on eight agarwood samples in different forms using the proposed barcode combination proved successful in their identification up to the species level. Such potential of DNA barcoding in identifying the source species of agarwood will contribute to the international timber trade control, by providing an effective method for species identification and product authentication.

Assessing DNA Barcodes for Species Identification in North American Reptiles and Amphibians in Natural History Collections

Background

High rates of species discovery and loss have led to the urgent need for more rapid assessment of species diversity in the herpetofauna. DNA barcoding allows for the preliminary identification of species based on sequence divergence. Prior DNA barcoding work on reptiles and amphibians has revealed higher biodiversity counts than previously estimated due to cases of cryptic and undiscovered species. Past studies have provided DNA barcodes for just 14% of the North American herpetofauna, revealing the need for expanded coverage.

Methodology/Principal Findings

This study extends the DNA barcode reference library for North American herpetofauna, assesses the utility of this approach in aiding species delimitation, and examines the correspondence between current species boundaries and sequence clusters designated by the BIN system. Sequences were obtained from 730 specimens, representing 274 species (43%) from the North American herpetofauna. Mean intraspecific divergences were 1% and 3%, while average congeneric sequence divergences were 16% and 14% in amphibians and reptiles, respectively. BIN assignments corresponded with current species boundaries in 79% of amphibians, 100% of turtles, and 60% of squamates. Deep divergences (>2%) were noted in 35% of squamate and 16% of amphibian species, and low divergences (<2%) occurred in 12% of reptiles and 23% of amphibians, patterns reflected in BIN assignments. Sequence recovery declined with specimen age, and variation in recovery success was noted among collections. Within collections, barcodes effectively flagged seven mislabeled tissues, and barcode fragments were recovered from five formalin-fixed specimens.

Conclusions/Significance

This study demonstrates that DNA barcodes can effectively flag errors in museum collections, while BIN splits and merges reveal taxa belonging to deeply diverged or hybridizing lineages. This study is the first effort to compile a reference library of DNA barcodes for herpetofauna on a continental scale.

A new species of the genus Amolops (Amphibia: Ranidae) from southeastern Tibet, China

A new species of the genus Amolops Cope, 1865 is described from Nyingchi, southeastern Tibet, China, based on morphological and molecular data. The new species, Amolops nyingchiensis sp. nov. is assigned to the Amolops monticola group based on its skin smooth, dorsolateral fold distinct, lateral side of head black, upper lip stripe white extending to the shoulder. Amolops nyingchiensis sp. nov. is distinguished from all other species of Amolops by the following combination of characters: (1) medium body size, SVL 48.5-58.3 mm in males, and 57.6-70.7 mm in females; (2) tympanum distinct, slightly larger than one third of the eye diameter; (3) a small tooth-like projection on anteromedial edge of mandible; (4) the absence of white spine on dorsal surface of body; (5) the presence of circummarginal groove on all fingers; (6) the presence of vomerine teeth; (7) background coloration of dorsal surface brown, lateral body gray with yellow; (8) the presence of transverse bands on the dorsal limbs; (9) the presence of nuptial pad on the first finger in males; (10) the absence of vocal sac in males. Taxonomic status of the populations that were previously identified to A. monticola from Tibet is also discussed.

Using DNA Barcodes to Identify Road-Killed Animals in Two Atlantic Forest Nature Reserves, Brazil

Road mortality is the leading source of biodiversity loss in the world, especially due to fragmentation of natural habitats and loss of wildlife. The survey of the main species victims of roadkill is of fundamental importance for the better understanding of the problem, being necessary, for this, the correct species identification. The aim of this study was to verify if DNA barcodes can be applied to identify road-killed samples that often cannot be determined morphologically. For this purpose, 222 vertebrate samples were collected in a stretch of the BR-101 highway that crosses two Discovery Coast Atlantic Forest Natural Reserves, the Sooretama Biological Reserve and the Vale Natural Reserve, in Espírito Santo, Brazil. The mitochondrial COI gene was amplified, sequenced and confronted with the BOLD database. It was possible to identify 62.16% of samples, totaling 62 different species, including Pyrrhura cruentata, Chaetomys subspinosus, Puma yagouaroundi and Leopardus wiedii considered Vulnerable in the National Official List of Species of Endangered Wildlife. The most commonly identified animals were a bat (Molossus molossus), an opossum (Didelphis aurita) and a frog (Trachycephalus mesophaeus) species. Only one reptile was identified using the technique, probably due to lack of reference sequences in BOLD. These data may contribute to a better understanding of the impact of roads on species biodiversity loss and to introduce the DNA barcode technique to road ecology scenarios.

DNA barcoding reveals species level divergence between populations of the microhylid frog genus Arcovomer (Anura: Microhylidae) in the Atlantic Rainforest of southeastern Brazil

The microhylid frogs belonging to the genus Arcovomer have been reported from lowland Atlantic Rainforest in the Brazilian states of Espírito Santo, Rio de Janeiro, and São Paulo. Here, we use DNA barcoding to assess levels of genetic divergence between apparently isolated populations in Espírito Santo and Rio de Janeiro. Our mtDNA data consisting of cytochrome oxidase subunit I (COI) nucleotide sequences reveals 13.2% uncorrected and 30.4% TIM2 + I + Γ corrected genetic divergences between these two populations. This level of divergence exceeds the suggested 10% uncorrected divergence threshold for elevating amphibian populations to candidate species using this marker, which implies that the Espírito Santo population is a species distinct from Arcovomer passarellii. Calibration of our model-corrected sequence divergence estimates suggests that the time of population divergence falls between 12 and 29 million years ago.

DNA barcoding and the identification of tree frogs (Amphibia: Anura: Rhacophoridae)

The DNA barcoding gene COI (cytochrome c oxidase subunit I) effectively identifies many species. Herein, we barcoded 172 individuals from 37 species belonging to nine genera in Rhacophoridae to test if the gene serves equally well to identify species of tree frogs. Phenetic neighbor joining and phylogenetic Bayesian inference were used to construct phylogenetic trees, which resolved all nine genera as monophyletic taxa except for Rhacophorus, two new matrilines for Liuixalus, and Polypedates leucomystax species complex. Intraspecific genetic distances ranged from 0.000 to 0.119 and interspecific genetic distances ranged from 0.015 to 0.334. Within Rhacophorus and Kurixalus, the intra- and interspecific genetic distances did not reveal an obvious barcode gap. Notwithstanding, we found that COI sequences unambiguously identified rhacophorid species and helped to discover likely new cryptic species via the synthesis of genealogical relationships and divergence patterns. Our results supported that COI is an effective DNA barcoding marker for Rhacophoridae.

DNA Barcoding Survey of Anurans across the Eastern Cordillera of Colombia and the Impact of the Andes on Cryptic Diversity

Colombia hosts the second highest amphibian species diversity on Earth, yet its fauna remains poorly studied, especially using molecular genetic techniques. We present the results of the first wide-scale DNA barcoding survey of anurans of Colombia, focusing on a transect across the Eastern Cordillera. We surveyed 10 sites between the Magdalena Valley to the west and the eastern foothills of the Eastern Cordillera, sequencing portions of the mitochondrial 16S ribosomal RNA and cytochrome oxidase subunit 1 (CO1) genes for 235 individuals from 52 nominal species. We applied two barcode algorithms, Automatic Barcode Gap Discovery and Refined Single Linkage Analysis, to estimate the number of clusters or “unconfirmed candidate species” supported by DNA barcode data. Our survey included ~7% of the anuran species known from Colombia. While barcoding algorithms differed slightly in the number of clusters identified, between three and ten nominal species may be obscuring candidate species (in some cases, more than one cryptic species per nominal species). Our data suggest that the high elevations of the Eastern Cordillera and the low elevations of the Chicamocha canyon acted as geographic barriers in at least seven nominal species, promoting strong genetic divergences between populations associated with the Eastern Cordillera.

COI-based barcoding of Chinese vipers (Reptilia: Squamata: Viperidae)

DNA barcoding seeks to assemble a standardized reference library for rapid and unambiguous identification of species, and can be used to screen for potentially cryptic species. The 5′ region of cytochrome oxidase subunit I (COI), which is a mitochondrial DNA (mtDNA) gene fragment, has been proposed as a universal marker for this purpose among animals. However, DNA barcoding of reptiles is still supported only by few datasets compared with other groups. We investigated the utilization of COI to discriminate 34 putative species of vipers, representing almost 92% of the recorded species in China. Based on a total of 241 sequences, our results indicated that the average degree of intraspecific variability (0.0198) tends to be one-sixth the average of interspecific divergence (0.0931), but no barcoding gap was detected between them. The threshold method, BLOG analyses and tree-based methods all can identify species with a high success rate. These results consistently suggested the usefulness and reliability of the DNA barcoding approach in Chinese vipers.

DNA barcodes and species distribution models evaluate threats of global climate changes to genetic diversity: a case study from Nanorana parkeri (Anura: Dicroglossidae)

Anthropogenic global climate changes are one of the greatest threats to biodiversity. Distribution modeling can predict the effects of climate changes and potentially their effects on genetic diversity. DNA barcoding quickly identifies patterns of genetic diversity. As a case study, we use DNA barcodes and distribution models to predict threats under climate changes in the frog Nanorana parkeri, which is endemic to the Qinghai-Tibetan Plateau. Barcoding identifies major lineages W and E. Lineage W has a single origin in a refugium and Lineage E derives from three refugia. All refugia locate in river valleys and each greatly contributes to the current level of intraspecific genetic diversity. Species distribution models suggest that global climate changes will greatly influence N. parkeri, especially in the level of genetic diversity, because two former refugia will fail to provide suitable habitat. Our pipeline provides a novel application of DNA barcoding and has important implications for the conservation of biodiversity in southern areas of the Qinghai-Tibetan Plateau.

Current and future directions of DNA in wildlife forensic science

Wildlife forensic science may not have attained the profile of human identification, yet the scale of criminal activity related to wildlife is extensive by any measure. Service delivery in the arena of wildlife forensic science is often ad hoc, unco-ordinated and unregulated, yet many of those currently dedicated to wildlife conservation and the protection of endangered species are striving to ensure that the highest standards are met. The genetic markers and software used to evaluate data in wildlife forensic science are more varied than those in human forensic identification and are rarely standardised between species. The time and resources required to characterise and validate each genetic maker is considerable and in some cases prohibitive. Further, issues are regularly encountered in the construction of allelic databases and allelic ladders; essential in human identification studies, but also applicable to wildlife criminal investigations. Accreditation and certification are essential in human identification and are currently being strived for in the forensic wildlife community. Examples are provided as to how best practice can be demonstrated in all areas of wildlife crime analysis and ensure that this field of forensic science gains and maintains the respect it deserves. This review is aimed at those conducting human identification to illustrate how research concepts in wildlife forensic science can be used in the criminal justice system, as well as describing the real importance of this type of forensic analysis.

DNA barcode reference data for the Korean herpetofauna and their applications

Recently, amphibians and reptiles have drawn attention because of declines in species and populations caused mainly by habitat loss, overexploitation and climate change. This study constructed a DNA barcode database for the Korean herpetofauna, including all the recorded amphibians and 68% of the recorded reptiles, to provide a useful, standardized tool for species identification in monitoring and management. A total of 103 individuals from 18 amphibian and 17 reptile species were used to generate barcode sequences using partial sequences of the mitochondrial cytochrome c oxidase subunit I (COI) gene and to compare it with other suggested barcode loci. Comparing 16S rRNA, cytochrome b (Cytb) and COI for amphibians and 12S rRNA, Cytb and COI for reptiles, our results revealed that COI is better than the other markers in terms of a high level of sequence variation without length variation and moderate amplification success. Although the COI marker had no clear barcoding gap because of the high level of intraspecific variation, all of the analysed individuals from the same species clustered together in a neighbour-joining tree. High intraspecific variation suggests the possibility of cryptic species. Finally, using this database, confiscated snakes were identified as Elaphe schrenckii, designated as endangered in Korea and a food contaminant was identified as the lizard Takydromus amurensis.

A de novo case of floating chromosomal polymorphisms by translocation in Quasipaa boulengeri (Anura, Dicroglossidae)

Very few natural polymorphisms involving interchromosomal reciprocal translocations are known in amphibians even in vertebrates. In this study, thirty three populations, including 471 individuals of the spiny frog Quasipaa boulengeri, were karyotypically examined using Giemsa stain or FISH. Five different karyomorphs were observed. The observed heteromorphism was autosomal but not sex-related, as the same heteromorphic chromosomes were found both in males and females. Our results indicated that the variant karyotypes resulted from a mutual interchange occurring between chromosomes 1 and 6. The occurrence of a nearly whole-arm translocation between chromosome no. 1 and no. 6 gave rise to a high frequency of alternate segregation and probably resulted in the maintenance of the translocation polymorphisms in a few populations. The translocation polymorphism is explained by different frequencies of segregation modes of the translocation heterozygote during meiosis. Theoretically, nine karyomorphs should be investigated, however, four expected karyotypes were not found. The absent karyomorphs may result from recessive lethal mutations, position effects, duplications and deficiencies. The phylogenetic inference proved that all populations of Q. boulengeri grouped into a monophyletic clade. The mutual translocation likely evolved just once in this species and the dispersal of the one karyomorph (type IV) can explain the chromosomal variations among populations.

First large-scale DNA barcoding assessment of reptiles in the biodiversity hotspot of Madagascar, based on newly designed COI primers

BACKGROUND

DNA barcoding of non-avian reptiles based on the cytochrome oxidase subunit I (COI) gene is still in a very early stage, mainly due to technical problems. Using a newly developed set of reptile-specific primers for COI we present the first comprehensive study targeting the entire reptile fauna of the fourth-largest island in the world, the biodiversity hotspot of Madagascar.

METHODOLOGY/PRINCIPAL FINDINGS

Representatives of the majority of Madagascan non-avian reptile species (including Squamata and Testudines) were sampled and successfully DNA barcoded. The new primer pair achieved a constantly high success rate (72.7-100%) for most squamates. More than 250 species of reptiles (out of the 393 described ones; representing around 64% of the known diversity of species) were barcoded. The average interspecific genetic distance within families ranged from a low of 13.4% in the Boidae to a high of 29.8% in the Gekkonidae. Using the average genetic divergence between sister species as a threshold, 41-48 new candidate (undescribed) species were identified. Simulations were used to evaluate the performance of DNA barcoding as a function of completeness of taxon sampling and fragment length. Compared with available multi-gene phylogenies, DNA barcoding correctly assigned most samples to species, genus and family with high confidence and the analysis of fewer taxa resulted in an increased number of well supported lineages. Shorter marker-lengths generally decreased the number of well supported nodes, but even mini-barcodes of 100 bp correctly assigned many samples to genus and family.

CONCLUSIONS/SIGNIFICANCE

The new protocols might help to promote DNA barcoding of reptiles and the established library of reference DNA barcodes will facilitate the molecular identification of Madagascan reptiles. Our results might be useful to easily recognize undescribed diversity (i.e. novel taxa), to resolve taxonomic problems, and to monitor the international pet trade without specialized expert knowledge.

DNA barcoding amphibians and reptiles

Only a few major research programs are currently targeting COI barcoding of amphibians and reptiles (including chelonians and crocodiles), two major groups of tetrapods. Amphibian and reptile species are typically old, strongly divergent, and contain deep conspecific lineages which might lead to problems in species assignment with incomplete reference databases. As far as known, there is no single pair of COI primers that will guarantee a sufficient rate of success across all amphibian and reptile taxa, or within major subclades of amphibians and reptiles, which means that the PCR amplification strategy needs to be adjusted depending on the specific research question. In general, many more amphibian and reptile taxa have been sequenced for 16S rDNA, which for some purposes may be a suitable complementary marker, at least until a more comprehensive COI reference database becomes available. DNA barcoding has successfully been used to identify amphibian larval stages (tadpoles) in species-rich tropical assemblages. Tissue sampling, DNA extraction, and amplification of COI is straightforward in amphibians and reptiles. Single primer pairs are likely to have a failure rate between 5 and 50% if taxa of a wide taxonomic range are targeted; in such cases the use of primer cocktails or subsequent hierarchical usage of different primer pairs is necessary. If the target group is taxonomically limited, many studies have followed a strategy of designing specific primers which then allow an easy and reliable amplification of all samples.