A retrospective approach to testing the DNA barcoding method

Integrative species delimitation and five new species of lynx spiders (Araneae, Oxyopidae) in Taiwan

An accurate assessment of species diversity is a cornerstone of biology and conservation. The lynx spiders (Araneae: Oxyopidae) represent one of the most diverse and widespread cursorial spider groups, however their species richness in Asia is highly underestimated. In this study, we revised species diversity with extensive taxon sampling in Taiwan and explored species boundaries based on morphological traits and genetic data using a two-step approach of molecular species delimitation. Firstly, we employed a single COI dataset and applied two genetic distance-based methods: ABGD and ASAP, and two topology-based methods: GMYC and bPTP. Secondly, we further analyzed the lineages that were not consistently delimited, and incorporated H3 to the dataset for a coalescent-based analysis using BPP. A total of eight morphological species were recognized, including five new species, Hamataliwa cordivulva sp. nov., Hamat. leporauris sp. nov., Tapponia auriola sp. nov., T. parva sp. nov. and T. rarobulbus sp. nov., and three newly recorded species, Hamadruas hieroglyphica (Thorell, 1887), Hamat. foveata Tang & Li, 2012 and Peucetia latikae Tikader, 1970. All eight morphological species exhibited reciprocally monophyletic lineages. The results of molecular-based delimitation analyses suggested a variety of species hypotheses that did not fully correspond to the eight morphological species. We found that Hamat. cordivulva sp. nov. and Hamat. foveata showed shallow genetic differentiation in the COI, but they were unequivocally distinguishable according to their genitalia. In contrast, T. parva sp. nov. represented a deep divergent lineage, while differences of genitalia were not detected. This study highlights the need to comprehensively employ multiple evidence and methods to delineate species boundaries and the values of diagnostic morphological characters for taxonomic studies in lynx spiders.

Molecular characterization and genetic diversity of three Stomoxys flies S. bengalensis, S. calcitrans, and S. sitiens (Diptera: Muscidae) in Central Thailand

Stomoxys flies (Diptera: Muscidae) are hematophagous ectoparasites of medical and veterinary importance. In this study, three Stomoxys species, i.e. S. bengalensis, S. calcitrans, and S. sitiens, were collected from three provinces in Central Thailand with the aim of estimating the genetic divergence between species, for species identification, as well as within species, for a genetic diversity study based on the cytochrome c oxidase subunit I (COI) gene. Our results showed that the average intraspecific genetic divergences of Stomoxys flies ranged from 0.11% in S. sitiens to 0.98% in S. calcitrans, whereas the average interspecific genetic divergences ranged from 5.24% between S. sitiens and S. bengalensis to 6.69% between S. calcitrans and S. bengalensis. In addition, there was no overlap between the intraspecific and interspecific genetic divergences. The COI sequence analysis revealed a high haplotype diversity and low nucleotide diversity, reflecting a rapid population expansion after past bottlenecks. Moreover, there was no significant difference (P > 0.05) in the pairwise population differentiation (Fst) among Stomoxys flies in Central Thailand, because of the lack of natural barriers, thus allowing genetic exchange between them. The monitoring of the haplotype network revealed that two lineages of S. calcitrans in Central Thailand were distributed in all study areas, including the Nakhon Pathom, Pathum Thani, and Saraburi Provinces. These findings may improve our understanding of the genetic patterns of these three Stomoxys flies, as well as the underlying biological mechanisms, which is knowledge that can be used for further effective control of these flies.

Evaluation of partial 12S rRNA, 16S rRNA, COI and Cytb gene sequence datasets for potential single DNA barcode for hylids (Anura: Hylidae)

We evaluated the extent of intraspecific and interspecific genetic distances and the effectiveness of predefined threshold values using the main genes for estimates of biodiversity and specimens' identification in anurans. Partial sequences of the mitochondrial genes for small (12S) and large (16S) ribosomal subunits, cytochrome c oxidase subunit I (COI) and Cytochrome b (Cytb) of the family Hylidae were downloaded from GenBank and curated for length, coverage, and potential contaminations. We performed analyses for all sequences of each gene and the same species present in these datasets by distance and tree (monophyly)-based evaluations. We also evaluated the ability to identify specimens using these datasets applying "nearest neighbor" (NN), "best close match" (BCM) and "BOLD ID" tests. Genetic distance thresholds were generated by the function 'threshVal' and "localMinima" from SPIDER package and traditional threshold values (1%, 3%, 6% and 10%) were also evaluated. Coding genes, especially COI, had a better identification capacity than non-coding genes on barcoding gap and monophyly analysis and NN, BCM, BOLD ID tests. Considering the multiple factors involved in global DNA barcoding evaluations, we present a critical assessment of the use of these genes for biodiversity estimation and specimens' identification in anurans (e.g. hylids).

Molecular variation of the cytochrome b DNA and protein sequences in Phytoseiulus macropilis and P. persimilis (Acari: Phytoseiidae) reflect population differentiation

Several phytoseiid mite species are important natural enemies used in biological control strategies. In the present study, Cytb mtDNA sequences of various populations of two species, Phytoseiulus macropolis and P. persimilis, were compared to determine whether the specimens collected in Brazil could belong to P. persimilis as this latter species is reported in South America but not in Brazil. The Cytb marker was used because of its high evolution rate, assumed to capture intraspecific variation. No overlap between intra- and interspecific distances was observed but the distances were quite low for interspecific variation. This can be due to the particular biology of Phytoseiulus species and this shows the difficulty to apply a universal threshold in genetic distances to conclude about the existence of one or several species. Cytb mtDNA sequences were also considered to assess intraspecific variation. The DNA sequences of P. persimilis populations were very similar, probably because they all originated from the West Palearctic region or because of a prevalence of commercialized specimens in natura. For P. macropilis, higher genetic distances were observed and differentiation was noted according to geographic location and, to a smaller extent, pyrethroid resistance. To determine how DNA variation might impact the protein function (CytB fragment considered), the amino acid compositions of the populations studied were compared. No diagnostic mutation was observed between pyrethroid resistant and susceptible populations, whereas four mutations were identified between populations of P. macropilis separated by 1300 km (different climatic conditions). The impact of such mutations is discussed but knowledge is scarce, which makes it difficult to root testable hypotheses. The protein analysis clearly opens new perspectives in Phytoseiidae studies.

Local taxonomic spectra in plants, animals, fungi and terrestrial protists show common mathematical patterns

Taxonomic spectra, i.e. relations between supraspecific taxa by the number of included species, remain poorly understood in aspect of the mathematical properties. We studied taxonomic spectra of plants (Magnoliophyta, Bryophyta), animals (Coleoptera, Aves), fungi (Agaricomycetes) and terrestrial protists (Myxomycetes), found in the Homilsha Forests National Nature Park (North-East of Ukraine), and concluded that they correspond to the hollow-curve distribution at the level of genera, families and orders. The spectra of most taxa, as shown by the Akaike information criterion, are closely approximated by the log-series distribution model at all taxonomic levels. This type of distribution is typical for the species abundance curves, based on collections made from small areas. At the same time, in the genera–families–orders row the similarity to the lognormal distribution increases. The central values and variability vary considerably between different taxonomic groups and ranks, however, without affecting the type of distribution. The number of orders in all taxa except Bryophyta has reached the saturation and coincides with the curve of the estimated number of orders according to the Chao1 coefficient. For families and especially genera the correspondence with estimated number of species is much lower. Our results do not confirm the assumption that hollow-curve distributions of taxonomic spectra result from the artificial fragmentation of taxa. These distributions neither depend on the insufficient knowledge about the species composition at the locality, nor reflect the size of the studied area. The presence of such distributions in both local and global biota of different groups may be explained by the common features of their evolution, especially by the existence of relict orphan groups. The fact that in Homilsha Forests the kurtosis and skewness of distributions decreases in the genera–families–orders row can therefore be explained by the relatively low percentage of the high-rank orphan taxa in the local biota. This may be a common feature of communities studied at small geographical scale, since orphan taxa often demonstrate a high level of endemism. Comparative studies of local communities from different climate zones may help to understand how universal are the patterns, described herein.

Factors affecting the efficiency of molecular species delimitation in a species-rich insect family

In the context of global biodiversity loss, molecular species delimitation approaches can be very useful for accelerating species discovery through DNA taxonomy and inventory through DNA metabarcoding. In this study, the effect of some intrinsic factors on the efficiency of various single-marker species delimitation methods (fixed and variable nucleotide distance thresholds, ABGD, ASAP, GMYC, mPTP) was tested on more than 90 empirical data sets, derived from a set of 7,237 COI sequences attributed to 542 leaf beetles species (Coleoptera: Chrysomelidae). The considered factors were: (i) the number of haplotypes per species (as a proxy for genetic diversity), (ii) the geographic distance among conspecific collection localities (as a proxy of sampling width), (iii) the difficulty related to morphological identification of species, and (iv) the taxonomic rank. Distance-based methods, with on average more than 70% of match with morphological identification, outperformed those relying on phylogenetic trees, with less than 59%. A high number of haplotypes per species was found to have a negative effect on delimitation efficiency, whereas large geographic distances within species had a positive effect. All methods delimitations (except for GMYC) were significantly affected by the presence of species that are difficult to be identified, decreasing their efficiency. Finally, the only method influenced by the taxonomic rank of the data set was GMYC, showing lower efficiency in data sets at the genus than at higher levels. The observed biases we highlighted affecting efficiency could be accounted for when developing input data sets for species delimitation analyses to obtain a more reliable representation of biological diversity.

Barcoding Analysis of Paraguayan Squamata

Paraguay is a key spot in the central region of South America where several ecoregions converge. Its fauna (and specifically its herpetofauna) is getting better studied than years before, but still there is a lack of information regarding molecular genetics, and barcoding analyses have proven to be an excellent tool in this matter. Here, we present results of a barcoding analysis based on 16S rRNA gene sequences, providing valuable data for the scientific community in the region. We based our fieldwork in several areas of Paraguay. We analyzed 249 samples (142 sequenced by us) with a final alignment of 615 bp length. We identified some taxonomic incongruences that can be addressed based on our results. Furthermore, we identify groups, where collecting efforts and research activities should be reinforced. Even though we have some blanks in the geographical coverage of our analysis—and there is still a lot to do towards a better understanding of the taxonomy of the Paraguayan herpetofauna—here, we present the largest genetic dataset for the mitochondrial DNA gene 16S of reptiles (particularly, Squamata) from Paraguay, which can be used to solve taxonomic problems in the region.

Characterization of Fowlpox virus in chickens and bird-biting mosquitoes: a molecular approach to investigating Avipoxvirus transmission

Avian pox is a highly contagious avian disease, yet relatively little is known about the epidemiology and transmission of Avipoxviruses. Using a molecular approach, we report evidence for a potential link between birds and field-caught mosquitoes in the transmission of Fowlpox virus (FWPV) in Singapore. Comparison of fpv167 (P4b), fpv126 (VLTF-1), fpv175-176 (A11R-A12L) and fpv140 (H3L) gene sequences revealed close relatedness between FWPV strains obtained from cutaneous lesions of a chicken and four pools of Culex pseudovishnui, Culex spp. (vishnui group) and Coquellitidea crassipes caught in the vicinity of the study site. Chicken-derived viruses characterized during two separate infections two years later were also identical to those detected in the first event, suggesting repeated transmission of closely related FWPV strains in the locality. Since the study location is home to resident and migratory birds, we postulated that wild birds could be the source of FWPV and that bird-biting mosquitoes could act as bridging mechanical vectors. Therefore, we determined whether the FWPV-positive mosquito pools (n=4) were positive for avian DNA using a polymerase chain reaction-sequencing assay. Our findings confirmed the presence of avian host DNA in all mosquito pools, suggesting a role for Cx. pseudovishnui, Culex spp. (vishnui group) and Cq. crassipes mosquitoes in FWPV transmission. Our study exemplifies the utilization of molecular tools to understand transmission networks of pathogens affecting avian populations, which has important implications for the design of effective control measures to minimize disease burden and economic loss.

An exhaustive phylogeny of the combtooth blenny genus Salaria (Pisces, Blenniidae) shows introgressive hybridization and lack of reciprocal mtDNA monophyly between the marine species Salaria basilisca and Salaria pavo

A comprehensive phylogeny of the genus Salaria based on mitochondrial and nuclear markers grouped the extant species of the genus in well-characterised marine and freshwater clades, thus rejecting the hypothesis of a polytypic origin of the freshwater Salaria populations and supporting the occurrence of a single invasion event of the inland waters by the genus. Based on both mitochondrial and nuclear DNA datasets, the Salaria species of the freshwater clade proved to be vicariant taxa originating from a common ancestor which could possibly spread throughout the circum-Mediterranean inland waters during the late Miocene Messinian salinity crisis, then experiencing a process of allopatric differentiation after the re-flooding of the Mediterranean basin. Within the marine clade, although the nuDNA datasets showed the existence of well-supported subclades in accordance to the morphological identification of the studied specimens, one of the two subclades obtained in the phylogenetic tree based on the mtDNA dataset included both S. basilisca and S. pavo specimens, thus failing to find the two species as reciprocally monophyletic. Such a mito-nuclear discordance is here ascribed to multiple mtDNA unidirectional introgression events from S. basilisca to S. pavo, and the molecular diversity pattern of the marine Salaria species is here ascribed to a Pleistocene speciation event nowadays partly concealed by the occurrence of introgressive hybridization phenomena between the two taxa. Our results urge for prudence when implementing DNA barcoding approaches since, in the presence of mito-nuclear discordance phenomena, single-marker mtDNA-only analyses might lead to significant misidentifications.

A citizen science model for implementing statewide educational DNA barcoding

Our aim was to develop a widely available educational program in which students conducted authentic research that met the expectations of both the scientific and educational communities. This paper describes the development and implementation of a citizen science project based on DNA barcoding of reptile specimens obtained from the Museums Victoria frozen tissue collection. The student program was run by the Gene Technology Access Centre (GTAC) and was delivered as a "one day plus one lesson" format incorporating a one-day wet laboratory workshop followed by a single lesson at school utilising online bioinformatics tools. The project leveraged the complementary resources and expertise of the research and educational partners to generate robust scientific data that could be analysed with confidence, meet the requirements of the Victorian state education curriculum, and provide participating students with an enhanced learning experience. During two 1-week stints in 2013 and 2014, 406 students mentored by 44 postgraduate university students participated in the project. Students worked mainly in pairs to process ~200 tissue samples cut from 53 curated reptile specimens representing 17 species. A total of 27 novel Cytochrome Oxidase subunit 1 (CO1) sequences were ultimately generated for 8 south-east Australian reptile species of the families Scincidae and Agamidae.

Evaluating the Feasibility of Five Candidate DNA Barcoding Loci for Philippine Lasianthus Jack (Lasiantheae: Rubiaceae)

Introduction

The pantropical genus Lasianthus Jack is identified for high phenotypic plasticity making traditional taxonomic identification difficult. Having some members with important medicinal properties, a precise complimentary identification through DNA barcoding is needed for species delineation.

Materials and Methods

In this study, 12 samples representing six Philippine Lasianthus species were used to determine the most efficient barcoding loci among the cpDNA markers (matK, rbcL, rps16, and trnT-F) and nrDNA (ITS) based on the criteria of universality, discriminatory power, and resolution of species.

Results

The results revealed that ITS has the recommended primer universality, greatest interspecific divergences, and average resolution of species. Among the cpDNA markers, matK and rbcL are recommended but with minimal resolution of species. While trnT-F showed moderate interspecific variations and resolution of Lasianthus species, rps16 has the lowest interspecific divergence and resolution of species.

Conclusion

Consequently, ITS is the potential ideal DNA barcode for Lasianthus species.

SUMMARY

ITS, matK, and rps16 markers have the excellent amplification and sequence qualityITS marker has the highest interspecific divergence with the maximum values, followed by matK, rbcL, trnT-F, and rps16, respectivelyAll markers except rps16 yielded average resolution to Lasianthus speciesITS marker is the most ideal locus in terms of excellent universality, high interspecific discriminatory ability, and average species resolution. Abbreviations used: ITS: Internal Transcribe Spacer, matK: maturase K, rbcL: ribulose-1,5-biphospahte-carboxylase, rps16: ribosomal protein 16 small subunit gene.

How many species and under what names? Using DNA barcoding and GenBank data for west Central African amphibian conservation

Development projects in west Central Africa are proceeding at an unprecedented rate, often with little concern for their effects on biodiversity. In an attempt to better understand potential impacts of a road development project on the anuran amphibian community, we conducted a biodiversity assessment employing multiple methodologies (visual encounter transects, auditory surveys, leaf litter plots and pitfall traps) to inventory species prior to construction of a new road within the buffer zone of Moukalaba-Doudou National Park, Gabon. Because of difficulties in morphological identification and taxonomic uncertainty of amphibian species observed in the area, we integrated a DNA barcoding analysis into the project to improve the overall quality and accuracy of the species inventory. Based on morphology alone, 48 species were recognized in the field and voucher specimens of each were collected. We used tissue samples from specimens collected at our field site, material available from amphibians collected in other parts of Gabon and the Republic of Congo to initiate a DNA barcode library for west Central African amphibians. We then compared our sequences with material in GenBank for the genera recorded at the study site to assist in identifications. The resulting COI and 16S barcode library allowed us to update the number of species documented at the study site to 28, thereby providing a more accurate assessment of diversity and distributions. We caution that because sequence data maintained in GenBank are often poorly curated by the original submitters and cannot be amended by third-parties, these data have limited utility for identification purposes. Nevertheless, the use of DNA barcoding is likely to benefit biodiversity inventories and long-term monitoring, particularly for taxa that can be difficult to identify based on morphology alone; likewise, inventory and monitoring programs can contribute invaluable data to the DNA barcode library and the taxonomy of complex groups. Our methods provide an example of how non-taxonomists and parataxonomists working in understudied parts of the world with limited geographic sampling and comparative morphological material can use DNA barcoding and publicly available sequence data (GenBank) to rapidly identify the number of species and assign tentative names to aid in urgent conservation management actions and contribute to taxonomic resolution.

Diet Assessment Based on Rumen Contents: A Comparison between DNA Metabarcoding and Macroscopy

Dietary choices are central to our understanding of ecology and evolution. Still, many aspects of food choice have been hampered by time consuming procedures and methodological problems. Faster and cheaper methods, such as DNA metabarcoding, have therefore been widely adopted. However, there is still very little empirical support that this new method is better and more accurate compared to the classic methods. Here, we compare DNA metabarcoding to macroscopic identifications of rumen contents in two species of wild free-ranging ungulates: roe deer and fallow deer. We found that the methods were comparable, but they did not completely overlap. Sometimes the DNA method failed to identify food items that were found macroscopically, and the opposite was also true. However, the total number of taxa identified increased using DNA compared to the macroscopic analysis. Moreover, the taxonomic precision of metabarcoding was substantially higher, with on average 90% of DNA-sequences being identified to genus or species level compared to 75% of plant fragments using macroscopy. In niche overlap analyses, presence/absence data showed that both methods came to very similar conclusions. When using the sequence count data and macroscopic weight, niche overlap was lower than when using presence-absence data yet tended to increase when using DNA compared to macroscopy. Nevertheless, the significant positive correlation between macroscopic quantity and number of DNA sequences counted from the same plant group give support for the use of metabarcoding to quantify plants in the rumen. This study thus shows that there is much to be gained by using metabarcoding to quantitatively assess diet composition compared to macroscopic analysis, including higher taxonomic precision, sensitivity and cost efficiency.

Hybridization relics complicate barcode-based identification of species in earthworms

Introgressive hybridization results in mito-nuclear discordance which could obscure the delimitation of closely related taxa. Although such events are increasingly reported, they have been poorly studied in earthworms. Here, we propose a method for investigating the degree of introgressive hybridization between three taxa of the Allolobophora chlorotica aggregate within two field populations (N = 67 and N = 105) using a reference data set including published DNA barcoding and microsatellite data of all known A. chlorotica lineages (N = 85). For this, we used both molecular phylogenetic and population genetic approaches. The test of correspondence between mitochondrial cytochrome c oxidase I (COI) lineages and clusters of nuclear microsatellite genotypes allowed individuals to be sorted in three categories (matching, admixed and nonmatching) and additional markers (mitochondrial NADH dehydrogenase subunit 1, nuclear Histone 3 and Internal transcribed Spacer Region 2) were used for phylogenetic reconstructions in order to check assignments. Although 15 admixed individuals were observed, no early-generation hybrids were detected within the two populations. Interestingly, 14 nonmatching individuals (i.e. with a mtDNA haplotype that did not correspond to their nuclear cluster) were detected, a pattern that would result after multiple generations of unidirectional hybridization of female from one taxon to male of the other taxon. Because earthworms are simultaneous hermaphrodites, these events of unidirectional hybridization suggest sterility of the male function in several crosses and highlight that some individuals can be misidentified if reliance is placed on COI barcodes alone. These findings could improve the use of these barcodes in earthworms for species delineation.

Unexpectedly High Levels of Cryptic Diversity Uncovered by a Complete DNA Barcoding of Reptiles of the Socotra Archipelago

Few DNA barcoding studies of squamate reptiles have been conducted. Due to the significance of the Socotra Archipelago (a UNESCO Natural World Heritage site and a biodiversity hotspot) and the conservation interest of its reptile fauna (94% endemics), we performed the most comprehensive DNA barcoding study on an island group to date to test its applicability to specimen identification and species discovery. Reptiles constitute Socotra's most important vertebrate fauna, yet their taxonomy remains under-studied. We successfully DNA-barcoded 380 individuals of all 31 presently recognized species. The specimen identification success rate is moderate to high, and almost all species presented local barcoding gaps. The unexpected high levels of intra-specific variability found within some species suggest cryptic diversity. Species richness may be under-estimated by 13.8-54.4%. This has implications in the species' ranges and conservation status that should be considered for conservation planning. Other phylogenetic studies using mitochondrial and nuclear markers are congruent with our results. We conclude that, despite its reduced length (663 base pairs), cytochrome c oxidase 1, COI, is very useful for specimen identification and for detecting intra-specific diversity, and has a good phylogenetic signal. We recommend DNA barcoding to be applied to other biodiversity hotspots for quickly and cost-efficiently flagging species discovery, preferentially incorporated into an integrative taxonomic framework.

Does a global DNA barcoding gap exist in Annelida?

Accurate identification of unknown specimens by means of DNA barcoding is contingent on the presence of a DNA barcoding gap, among other factors, as its absence may result in dubious specimen identifications - false negatives or positives. Whereas the utility of DNA barcoding would be greatly reduced in the absence of a distinct and sufficiently sized barcoding gap, the limits of intraspecific and interspecific distances are seldom thoroughly inspected across comprehensive sampling. The present study aims to illuminate this aspect of barcoding in a comprehensive manner for the animal phylum Annelida. All cytochrome c oxidase subunit I sequences (cox1 gene; the chosen region for zoological DNA barcoding) present in GenBank for Annelida, as well as for "Polychaeta", "Oligochaeta", and Hirudinea separately, were downloaded and curated for length, coverage and potential contaminations. The final datasets consisted of 9782 (Annelida), 5545 ("Polychaeta"), 3639 ("Oligochaeta"), and 598 (Hirudinea) cox1 sequences and these were either (i) used as is in an automated global barcoding gap detection analysis or (ii) further analyzed for genetic distances, separated into bins containing intraspecific and interspecific comparisons and plotted in a graph to visualize any potential global barcoding gap. Over 70 million pairwise genetic comparisons were made and results suggest that although there is a tendency towards separation, no distinct or sufficiently sized global barcoding gap exists in either of the datasets rendering future barcoding efforts at risk of erroneous specimen identifications (but local barcoding gaps may still exist allowing for the identification of specimens at lower taxonomic ranks). This seems to be especially true for earthworm taxa, which account for fully 35% of the total number of interspecific comparisons that show 0% divergence.