Species identification in the taxonomically neglected, highly diverse, neotropical parasitoid wasp genus Notiospathius (Braconidae: Doryctinae) based on an integrative molecular and morphological approach

Various DNA sequence-based methods for species delineation have recently been developed to assess the species-richness of highly diverse, neglected invertebrate taxa. These methods, however, need to be tested under a variety of conditions, including the use of different markers and parameters. Here, we explored the species diversity of a species-rich group of braconid parasitoid wasps, the Neotropical genus Notiospathius, including 233 specimens from 10 different countries. We examined sequences of two mitochondrial (mt) (COI, cyt b) and one nuclear (wg) gene fragments. We analysed them separately as well as concatenating the mt data with the general mixed Yule-coalescent (GMYC) model for species delineation using different tree-building methods and parameters for reconstructing ultrametric trees. We evaluated the performance of GMYC analyses by comparing their species delineations with our morphospecies identifications. Reconstructing ultrametric trees with a relaxed lognormal clock rate using the program BEAST gave the most congruent results with morphology for the two mt markers. A tree obtained with wg using the programs MrBayes+Pathd8 had the fewest cases of incongruence with morphology, though the performance of this nuclear marker was considerably lower than that of COI and cyt b. Species delimitation using the coalescent prior to obtain ultrametric trees was morphologically more congruent with COI, whereas the Yule prior was more congruent with cyt b. The analyses concatenating the mt datasets failed to recover some species supported both by morphology and the separate analyses of the mt markers. The highest morphological congruence was obtained with the GMYC analysis on an ultrametric tree reconstructed with cyt b using the relaxed lognormal clock rate and the Yule prior, thus supporting the importance of using alternative markers when the information of the barcoding locus (COI) is not concordant with morphological evidence. Seventy-one species were delimited based on the congruence found among COI, cyt b and morphology. Both mt markers also revealed the existence of seven potential cryptic species. This high species richness from a scattered geographical sampling indicates that there is a remarkable number of Notiospathius species that remains undiscovered.     

Using DNA taxonomy to investigate the ecological determinants of plankton diversity: explaining the occurrence of Synchaeta spp. (Rotifera,Monogononta) in mountain lakes

1. The occurrence of unresolved complexes of cryptic species may hinder the identification of the main ecological drivers of biodiversity when different cryptic taxa have different ecological requirements. 2. We assessed factors influencing the occurrence of Synchaeta species (monogonont rotifers) in 17 waterbodies of the Trentino‐South Tyrol region in the Eastern Alps. To do so, we compared the results of using unresolved complexes of cryptic species, as is common practice in limnological studies based on morphological taxonomy, and having resolved cryptic complexes, made possible by DNA taxonomy. 3. To identify cryptic species, we used the generalised mixed Yule coalescent (GMYC) model. We investigated the relationship between the environment and the occurrence of Synchaeta spp. by multivariate ordination using two definitions of the units of diversity, namely (i) unresolved species complexes (morphospecies) and (ii) putative cryptic species (GMYC entities). Our expectation was that resolving complexes of cryptic species could provide more information than using morphospecies. 4. As expected, DNA taxonomy provided greater taxonomic resolution than morphological taxonomy. Further, environmental‐based multivariate ordination on cryptic species explained a significantly higher proportion of variance than that based on morphospecies. Occurrence of GMYC entities was related to total phosphorus (TP), whereas no relationship could be found between morphospecies and the environment. Moreover, different cryptic species within the same morphospecies showed different, and even opposite, preferences for TP. In addition, the wide geographical distribution of haplotypes and cryptic species indicated the absence of barriers to dispersal in Synchaeta.     

Testing the potential of DNA barcoding in vertebrate radiations: the case of the littoral cichlids (Pisces,Perciformes, Cichlidae) from Lake Tanganyika

We obtained 398 cytochrome c oxidase subunit I barcodes of 96 morphospecies of Lake Tanganyika (LT) cichlids from the littoral zone. The potential of DNA barcoding in these fishes was tested using both species identification and species delineation methods. The best match (BM) and best close match (BCM) methods were used to evaluate the overall identification success. For this, three libraries were analysed in which the specimens were categorized into Operational Taxonomic Units (OTU) in three alternative ways: (A) morphologically distinct, including undescribed, species, (B) valid species and (C) complexes of morphologically similar or closely related species. For libraries A, B and C, 73, 73 and 96% (BM) and 72, 70 and 94% (BCM) of the specimens were correctly identified. Additionally, the potential of two species delineation methods was tested. The General Mixed Yule Coalescent (GMYC) analysis suggested 70 hypothetical species, while the Automatic Barcode Gap Discovery (ABGD) method revealed 115 putative species. Although the ABGD method had a tendency to oversplit, it outperformed the GMYC analysis in retrieving the species. In most cases where ABGD suggested oversplitting, this was due to intraspecific geographical variation. The failure of the GMYC method to retrieve many species could be attributed to discrepancies between mitochondrial gene trees and the evolutionary histories of LT cichlid species. Littoral LT cichlids have complex evolutionary histories that include instances of hybridization, introgression and rapid speciation. Nevertheless, although the utility of DNA barcoding in identification is restricted to the level of complexes, it has potential for species discovery in cichlid radiations.     

Cryptic diversity and discordance in single‐locus species delimitation methods within horned lizards (Phrynosomatidae: Phrynosoma)

Biodiversity reduction and loss continues to progress at an alarming rate, and thus, there is widespread interest in utilizing rapid and efficient methods for quantifying and delimiting taxonomic diversity. Single‐locus species delimitation methods have become popular, in part due to the adoption of the DNA barcoding paradigm. These techniques can be broadly classified into tree‐based and distance‐based methods depending on whether species are delimited based on a constructed genealogy. Although the relative performance of these methods has been tested repeatedly with simulations, additional studies are needed to assess congruence with empirical data. We compiled a large data set of mitochondrial ND4 sequences from horned lizards (Phrynosoma) to elucidate congruence using four tree‐based (single‐threshold GMYC, multiple‐threshold GMYC, bPTP, mPTP) and one distance‐based (ABGD) species delimitation models. We were particularly interested in cases with highly uneven sampling and/or large differences in intraspecific diversity. Results showed a high degree of discordance among methods, with multiple‐threshold GMYC and bPTP suggesting an unrealistically high number of species (29 and 26 species within the P. douglasii complex alone). The single‐threshold GMYC model was the most conservative, likely a result of difficulty in locating the inflection point in the genealogies. mPTP and ABGD appeared to be the most stable across sampling regimes and suggested the presence of additional cryptic species that warrant further investigation. These results suggest that the mPTP model may be preferable in empirical data sets with highly uneven sampling or large differences in effective population sizes of species.     

Delineating Species with DNA Barcodes: A Case of Taxon Dependent Method Performance in Moths

The accelerating loss of biodiversity has created a need for more effective ways to discover species. Novel algorithmic approaches for analyzing sequence data combined with rapidly expanding DNA barcode libraries provide a potential solution. While several analytical methods are available for the delineation of operational taxonomic units (OTUs), few studies have compared their performance. This study compares the performance of one morphology-based and four DNA-based (BIN, parsimony networks, ABGD, GMYC) methods on two groups of gelechioid moths. It examines 92 species of Finnish Gelechiinae and 103 species of Australian Elachistinae which were delineated by traditional taxonomy. The results reveal a striking difference in performance between the two taxa with all four DNA-based methods. OTU counts in the Elachistinae showed a wider range and a relatively low (ca. 65%) OTU match with reference species while OTU counts were more congruent and performance was higher (ca. 90%) in the Gelechiinae. Performance rose when only monophyletic species were compared, but the taxon-dependence remained. None of the DNA-based methods produced a correct match with non-monophyletic species, but singletons were handled well. A simulated test of morphospecies-grouping performed very poorly in revealing taxon diversity in these small, dull-colored moths. Despite the strong performance of analyses based on DNA barcodes, species delineated using single-locus mtDNA data are best viewed as OTUs that require validation by subsequent integrative taxonomic work.     

Exploring the Leaf Beetle Fauna (Coleoptera: Chrysomelidae) of an Ecuadorian Mountain Forest Using DNA Barcoding

Background

Tropical mountain forests are hotspots of biodiversity hosting a huge but little known diversity of insects that is endangered by habitat destruction and climate change. Therefore, rapid assessment approaches of insect diversity are urgently needed to complement slower traditional taxonomic approaches. We empirically compare different DNA-based species delimitation approaches for a rapid biodiversity assessment of hyperdiverse leaf beetle assemblages along an elevational gradient in southern Ecuador and explore their effect on species richness estimates.

Methodology/Principal Findings

Based on a COI barcode data set of 674 leaf beetle specimens (Coleoptera: Chrysomelidae) of 266 morphospecies from three sample sites in the Podocarpus National Park, we employed statistical parsimony analysis, distance-based clustering, GMYC- and PTP-modelling to delimit species-like units and compared them to morphology-based (parataxonomic) species identifications. The four different approaches for DNA-based species delimitation revealed highly similar numbers of molecular operational taxonomic units (MOTUs) (n = 284–289). Estimated total species richness was considerably higher than the sampled amount, 414 for morphospecies (Chao2) and 469–481 for the different MOTU types. Assemblages at different elevational levels (1000 vs. 2000 m) had similar species numbers but a very distinct species composition for all delimitation methods. Most species were found only at one elevation while this turnover pattern was even more pronounced for DNA-based delimitation.

Conclusions/Significance

Given the high congruence of DNA-based delimitation results, probably due to the sampling structure, our study suggests that when applied to species communities on a regionally limited level with high amount of rare species (i.e. ~50% singletons), the choice of species delimitation method can be of minor relevance for assessing species numbers and turnover in tropical insect communities. Therefore, DNA-based species delimitation is confirmed as a valuable tool for evaluating biodiversity of hyperdiverse insect communities, especially when exact taxonomic identifications are missing.     

Multiple species delimitation approaches with COI barcodes poorly fit each other and morphospecies – An integrative taxonomy case of Sri Lankan Sericini chafers (Coleoptera: Scarabaeidae)

DNA taxonomy including barcoding and metabarcoding is widely used to explore the diversity in biodiversity hotspots. In most of these hotspot areas, chafers are represented by a multitude of species, which are well defined by the complex shape of male genitalia. Here, we explore how well COI barcode data reflect morphological species entities and thus their usability for accelerated species inventorization. We conducted dedicated field surveys in Sri Lanka to collect the species‐rich and highly endemic Sericini chafers (Coleoptera: Scarabaeidae). Congruence among results of a series of protocols for de novo species delimitation and with morphology‐based species identifications was investigated. Different delimitation methods, such as the Poisson tree processes (PTP) model, Statistical Parsimony Analysis (TCS), Automatic Barcode Gap Discovery (ABGD), Assemble Species by Automatic Partitioning (ASAP), and Barcode Index Number (BIN) assignments, resulted in different numbers of molecular operational taxonomic units (MOTUs). All methods showed both over‐splitting and lumping of morphologically identified species. Only 18 of the observed 45 morphospecies perfectly matched MOTUs from all methods. The congruence of delimitation between MOTUs and morphospecies expressed by the match ratio was low, ranging from 0.57 to 0.67. TCS and multirate PTP (mPTP) showed the highest match ratio, while (BIN) assignment resulted in the lowest match ratio and most splitting events. mPTP lumped more species than any other method. Principal coordinate analysis (PCoA) on a match ratio‐based distance matrix revealed incongruent outcomes of multiple DNA delimitation methods, although applied to the same data. Our results confirm that COI barcode data alone are unlikely to correctly delimit all species, in particular, when using only a single delimitation approach. We encourage the integration of various approaches and data, particularly morphology, to validate species boundaries.     

Mitochondrial cox1 sequence data reliably uncover patterns of insect diversity but suffer from high lineage-idiosyncratic error rates

Background

The demand for scientific biodiversity data is increasing, but taxonomic expertise is often limited or not available. DNA sequencing is a potential remedy to overcome this taxonomic impediment. Mitochondrial DNA is most commonly used, e.g., for species identification (“DNA barcoding”). Here, we present the first study in arthropods based on a near-complete species sampling of a family-level taxon from the entire Australian region. We aimed to assess how reliably mtDNA data can capture species diversity when many sister species pairs are included. Then, we contrasted phylogenetic subsampling with the hitherto more commonly applied geographical subsampling, where sister species are not necessarily captured.

Methodology/Principal Findings

We sequenced 800 bp cox1 for 1,439 individuals including 260 Australian species (78% species coverage). We used clustering with thresholds of 1 to 10% and general mixed Yule Coalescent (GMYC) analysis for the estimation of species richness. The performance metrics used were taxonomic accuracy and agreement between the morphological and molecular species richness estimation. Clustering (at the 3% level) and GMYC reliably estimated species diversity for single or multiple geographic regions, with an error for larger clades of lower than 10%, thus outperforming parataxonomy. However, the rates of error were higher for some individual genera, with values of up to 45% when very recent species formed nonmonophyletic clusters. Taxonomic accuracy was always lower, with error rates above 20% and a larger variation at the genus level (0 to 70%). Sørensen similarity indices calculated for morphospecies, 3% clusters and GMYC entities for different pairs of localities was consistent among methods and showed expected decrease over distance.

Conclusion/Significance

Cox1 sequence data are a powerful tool for large-scale species richness estimation, with a great potential for use in ecology and β-diversity studies and for setting conservation priorities. However, error rates can be high in individual lineages.     

Toward a DNA taxonomy of Alpine Rhithrogena (Ephemeroptera: Heptageniidae) using a mixed Yule-coalescent analysis of mitochondrial and nuclear DNA

Aquatic larvae of many Rhithrogena mayflies (Ephemeroptera) inhabit sensitive Alpine environments. A number of species are on the IUCN Red List and many recognized species have restricted distributions and are of conservation interest. Despite their ecological and conservation importance, ambiguous morphological differences among closely related species suggest that the current taxonomy may not accurately reflect the evolutionary diversity of the group. Here we examined the species status of nearly 50% of European Rhithrogena diversity using a widespread sampling scheme of Alpine species that included 22 type localities, general mixed Yule-coalescent (GMYC) model analysis of one standard mtDNA marker and one newly developed nDNA marker, and morphological identification where possible. Using sequences from 533 individuals from 144 sampling localities, we observed significant clustering of the mitochondrial (cox1) marker into 31 GMYC species. Twenty-one of these could be identified based on the presence of topotypes (expertly identified specimens from the species' type locality) or unambiguous morphology. These results strongly suggest the presence of both cryptic diversity and taxonomic oversplitting in Rhithrogena. Significant clustering was not detected with protein-coding nuclear PEPCK, although nine GMYC species were congruent with well supported terminal clusters of nDNA. Lack of greater congruence in the two data sets may be the result of incomplete sorting of ancestral polymorphism. Bayesian phylogenetic analyses of both gene regions recovered four of the six recognized Rhithrogena species groups in our samples as monophyletic. Future development of more nuclear markers would facilitate multi-locus analysis of unresolved, closely related species pairs. The DNA taxonomy developed here lays the groundwork for a future revision of the important but cryptic Rhithrogena genus in Europe.     

Using Different Methods to Access the Difficult Task of Delimiting Species in a Complex Neotropical Hyperdiverse Group

The genus Rineloricaria is a Neotropical freshwater fish group with a long and problematic taxonomic history, attributed to the large number of species and the pronounced similarity among them. In the present work, taxonomic information and different molecular approaches were used to identify species boundaries and characterize independent evolutionary units. We analyzed 228 samples assembled in 53 distinct morphospecies. A general mixed yule-coalescent (GMYC) analysis indicated the existence of 70 entities, while BOLD system analyses showed the existence of 56 distinct BINs. When we used a new proposed integrative taxonomy approach, mixing the results obtained by each analysis, we identified 73 OTUs. We suggest that Rineloricaria probably has some complexity in the known species and several species not formally described yet. Our data suggested that other hyperdiverse fish groups with wide distributions can be further split into many new evolutionary taxonomic units.     

A population genetic assessment of taxonomic species: The case of Lake Malawi cichlid fishes

Organisms sampled for population‐level research are typically assigned to species by morphological criteria. However, if those criteria are limited to one sex or life stage, or the organisms come from a complex of closely related forms, the species assignments may misdirect analyses. The impact of such sampling can be assessed from the correspondence of genetic clusters, identified only from patterns of genetic variation, to the species identified using only phenotypic criteria. We undertook this protocol with the rock‐dwelling mbuna cichlids of Lake Malawi, for which species within genera are usually identified using adult male coloration patterns. Given high local endemism of male colour patterns, and considerable allele sharing among species, there persists considerable taxonomic uncertainty in these fishes. Over 700 individuals from a single transect were photographed, genotyped and separately assigned: (a) to morphospecies using photographs; and (b) to genetic clusters using five widely used methods. Overall, the correspondence between clustering methods was strong for larger clusters, but methods varied widely in estimated number of clusters. The correspondence between morphospecies and genetic clusters was also strong for larger clusters, as well as some smaller clusters for some methods. These analyses generally affirm (a) adult male‐limited sampling and (b) the taxonomic status of Lake Malawi mbuna, as the species in our study largely appear to be well‐demarcated genetic entities. More generally, our analyses highlight the challenges for clustering methods when the number of populations is unknown, especially in cases of highly uneven sample sizes.     

Evolutionary criteria outperform operational approaches in producing ecologically relevant fungal species inventories

Analyses of the structure and function of microbial communities are highly constrained by the diversity of organisms present within most environmental samples. A common approach is to rely almost entirely on DNA sequence data for estimates of microbial diversity, but to date there is no objective method of clustering sequences into groups that is grounded in evolutionary theory of what constitutes a biological lineage. The general mixed Yule-coalescent (GMYC) model uses a likelihood-based approach to distinguish population-level processes within lineages from processes associated with speciation and extinction, thus identifying a distinct point where extant lineages became independent. Using two independent surveys of DNA sequences associated with a group of ubiquitous plant-symbiotic fungi, we compared estimates of species richness derived using the GMYC model to those based on operational taxonomic units (OTUs) defined by fixed levels of sequence similarity. The model predicted lower species richness in these surveys than did traditional methods of sequence similarity. Here, we show for the first time that groups delineated by the GMYC model better explained variation in the distribution of fungi in relation to putative niche-based variables associated with host species identity, edaphic factors, and aspects of how the sampled ecosystems were managed. Our results suggest the coalescent-based GMYC model successfully groups environmental sequences of fungi into clusters that are ecologically more meaningful than more arbitrary approaches for estimating species richness.     

Excluding spatial sampling bias does not eliminate oversplitting in DNA‐based species delimitation analyses

DNA barcoding and DNA‐based species delimitation are major tools in DNA taxonomy. Sampling has been a central debate in this context, because the geographical composition of samples affects the accuracy and performance of DNA barcoding. Performance of complex DNA‐based species delimitation is to be tested under simpler conditions in absence of geographic sampling bias. Here, we present an empirical dataset sampled from a single locality in a Southeast‐Asian biodiversity hotspot (Laos: Phou Pan mountain). We investigate the performance of various species delimitation approaches on a megadiverse assemblage of herbivorous chafer beetles (Coleoptera: Scarabaeidae) to infer whether species delimitation suffers in the same way from exaggerate infraspecific variation despite the lack of geographic genetic variation that led to inconsistencies between entities from DNA‐based and morphology‐based species inference in previous studies. For this purpose, a 658 bp fragment of the mitochondrial cytochrome c oxidase subunit 1 (cox1) was analyzed for a total of 186 individuals of 56 morphospecies. Tree‐based and distance‐based species delimitation methods were used. All approaches showed a rather limited match ratio (max. 77%) with morphospecies. Poisson tree process (PTP) and statistical parsimony network analysis (TCS) prevailingly over‐splitted morphospecies, while 3% clustering and Automatic Barcode Gap Discovery (ABGD) also lumped several species into one entity. ABGD revealed the highest congruence between molecular operational taxonomic units (MOTUs) and morphospecies. Disagreements between morphospecies and MOTUs have to be explained by historically acquired geographic genetic differentiation, incomplete lineage sorting, and hybridization. The study once again highlights how important morphology still is in order to correctly interpret the results of molecular species delimitation.     

Testing the reliability of genetic methods of species identification via simulation

Although genetic methods of species identification, especially DNA barcoding, are strongly debated, tests of these methods have been restricted to a few empirical cases for pragmatic reasons. Here we use simulation to test the performance of methods based on sequence comparison (BLAST and genetic distance) and tree topology over a wide range of evolutionary scenarios. Sequences were simulated on a range of gene trees spanning almost three orders of magnitude in tree depth and in coalescent depth; that is, deep or shallow trees with deep or shallow coalescences. When the query's conspecific sequences were included in the reference alignment, the rate of positive identification was related to the degree to which different species were genetically differentiated. The BLAST, distance, and liberal tree-based methods returned higher rates of correct identification than did the strict tree-based requirement that the query was within, but not sister to, a single-species clade. Under this more conservative approach, ambiguous outcomes occurred in inverse proportion to the number of reference sequences per species. When the query's conspecific sequences were not in the reference alignment, only the strict tree-based approach was relatively immune to making false-positive identifications. Thresholds affected the rates at which false-positive identifications were made when the query's species was unrepresented in the reference alignment but did not otherwise influence outcomes. A conservative approach using the strict tree-based method should be used initially in large-scale identification systems, with effort made to maximize sequence sampling within species. Once the genetic variation within a taxonomic group is well characterized and the taxonomy resolved, then the choice of method used should be dictated by considerations of computational efficiency. The requirement for extensive genetic sampling may render these techniques inappropriate in some circumstances.     

DNA barcoding for community ecology--how to tackle a hyperdiverse, mostly undescribed Melanesian fauna

Background

Trigonopterus weevils are widely distributed throughout Melanesia and hyperdiverse in New Guinea. They are a dominant feature in natural forests, with narrow altitudinal zonation. Their use in community ecology has been precluded by the “taxonomic impediment”.

Methodology/Principal Findings

We sampled >6,500 specimens from seven areas across New Guinea; 1,002 specimens assigned to 270 morphospecies were DNA sequenced. Objective clustering of a refined dataset (excluding nine cryptic species) at 3% threshold revealed 324 genetic clusters (DNA group count relative to number of morphospecies = 20.0% overestimation of species diversity, or 120.0% agreement) and 85.6% taxonomic accuracy (the proportion of DNA groups that “perfectly” agree with morphology-based species hypotheses). Agreement and accuracy were best at an 8% threshold. GMYC analysis revealed 328 entities (21.5% overestimation) with 227 perfect GMYC entities (84.1% taxonomic accuracy). Both methods outperform the parataxonomist (19% underestimation; 31.6% taxonomic accuracy). The number of species found in more than one sampling area was highest in the Eastern Highlands and Huon (Sørensen similarity index 0.07, 4 shared species); ⅓ of all areas had no species overlap. Success rates of DNA barcoding methods were lowest when species showed a pronounced geographical structure. In general, Trigonopterus show high α and β-diversity across New Guinea.

Conclusions/Significance

DNA barcoding is an excellent tool for biodiversity surveys but success rates might drop when closer localities are included. Hyperdiverse Trigonopterus are a useful taxon for evaluating forest remnants in Melanesia, allowing finer-grained analyses than would be possible with vertebrate taxa commonly used to date. Our protocol should help establish other groups of hyperdiverse fauna as target taxa for community ecology. Sequencing delivers objective data on taxa of incredible diversity but mostly without a solid taxonomic foundation and should help pave the road for the eventual formal naming of new species.     

DNA barcode‐based delineation of putative species: efficient start for taxonomic workflows

The analysis of DNA barcode sequences with varying techniques for cluster recognition provides an efficient approach for recognizing putative species (operational taxonomic units, OTUs). This approach accelerates and improves taxonomic workflows by exposing cryptic species and decreasing the risk of synonymy. This study tested the congruence of OTUs resulting from the application of three analytical methods (ABGD, BIN, GMYC) to sequence data for Australian hypertrophine moths. OTUs supported by all three approaches were viewed as robust, but 20% of the OTUs were only recognized by one or two of the methods. These OTUs were examined for three criteria to clarify their status. Monophyly and diagnostic nucleotides were both uninformative, but information on ranges was useful as sympatric sister OTUs were viewed as distinct, while allopatric OTUs were merged. This approach revealed 124 OTUs of Hypertrophinae, a more than twofold increase from the currently recognized 51 species. Because this analytical protocol is both fast and repeatable, it provides a valuable tool for establishing a basic understanding of species boundaries that can be validated with subsequent studies.     

The influence of spatial scale on the congruence of classifications circumscribing morphological units of biodiversity

Aim The ‘taxonomic impediment’ has led to a growing trend in ecology and conservation biology to use operational surrogates for species within the context of a particular research project. Because such ‘parataxonomic’ classifications are typically spatially limited in scope, we examined the influence of increasing spatial scale on the congruence of two such approaches with a more traditional taxonomic classification. Location Sturt National Park, north‐western New South Wales, Australia. Methods Specimens of two ant genera, Camponotus and Rhytidoponera, were classified by three independent methods. The ‘parataxonomic’ classification assigned specimens to morphospecies without specialist taxonomic expertise; the ‘taxonomic’ classification assigned specimens to either described species or, where this was not possible, to operational taxonomic units (OTUs) using specialist taxonomic expertise; the ‘phenetic’ classification assigned specimens to putative species using a K‐means partitioning algorithm on basic morphometric data. Specimens of each genus were pooled into ‘assemblages’, which were defined at multiple spatial scales using a nested sampling design. Congruence in the interspecimen relationships of the different classifications was tested for each assemblage using pair‐wise Mantel correlations. Results Classification congruence tended to decrease with increasing spatial scale. There were, however, clear differences between the genera. Parataxonomic–taxonomic congruence was consistently greater for Camponotus, while phenetic–taxonomic congruence showed the opposite pattern. Conclusions Observed patterns in classification congruence are attributed to two principal causes: (i) within‐species morphological variation, including ecotypic variation in Rhytidoponera and caste polymorphism in Camponotus; and (ii) a limit to the morphological similarity of potentially competing species at small spatial scales. Regardless of cause, the decline in agreement as the spatial scale of observation is increased has important implications for the measurement of biodiversity, particularly when comparing samples over regional, continental, and global scales.     

Spatial heterogeneity in the Mediterranean Biodiversity Hotspot affects barcoding accuracy of its freshwater fishes

Incomplete knowledge of biodiversity remains a stumbling block for conservation planning and even occurs within globally important Biodiversity Hotspots (BH). Although technical advances have boosted the power of molecular biodiversity assessments, the link between DNA sequences and species and the analytics to discriminate entities remain crucial. Here, we present an analysis of the first DNA barcode library for the freshwater fish fauna of the Mediterranean BH (526 spp.), with virtually complete species coverage (498 spp., 98% extant species). In order to build an identification system supporting conservation, we compared species determination by taxonomists to multiple clustering analyses of DNA barcodes for 3165 specimens. The congruence of barcode clusters with morphological determination was strongly dependent on the method of cluster delineation, but was highest with the general mixed Yule‐coalescent (GMYC) model‐based approach (83% of all species recovered as GMYC entity). Overall, genetic morphological discontinuities suggest the existence of up to 64 previously unrecognized candidate species. We found reduced identification accuracy when using the entire DNA‐barcode database, compared with analyses on databases for individual river catchments. This scale effect has important implications for barcoding assessments and suggests that fairly simple identification pipelines provide sufficient resolution in local applications. We calculated Evolutionarily Distinct and Globally Endangered scores in order to identify candidate species for conservation priority and argue that the evolutionary content of barcode data can be used to detect priority species for future IUCN assessments. We show that large‐scale barcoding inventories of complex biotas are feasible and contribute directly to the evaluation of conservation priorities.     

A Good Compromise: Rapid and Robust Species Proxies for Inventorying Biodiversity Hotspots Using the Terebridae (Gastropoda: Conoidea)

Devising a reproducible approach for species delimitation of hyperdiverse groups is an ongoing challenge in evolutionary biology. Speciation processes combine modes of passive and adaptive trait divergence requiring an integrative taxonomy approach to accurately generate robust species hypotheses. However, in light of the rapid decline of diversity on Earth, complete integrative approaches may not be practical in certain species-rich environments. As an alternative, we applied a two-step strategy combining ABGD (Automated Barcode Gap Discovery) and Klee diagrams, to balance speed and accuracy in producing primary species hypotheses (PSHs). Specifically, an ABGD/Klee approach was used for species delimitation in the Terebridae, a neurotoxin-producing marine snail family included in the Conoidea. Delimitation of species boundaries is problematic in the Conoidea, as traditional taxonomic approaches are hampered by the high levels of variation, convergence and morphological plasticity of shell characters. We used ABGD to analyze gaps in the distribution of pairwise distances of 454 COI sequences attributed to 87 morphospecies and obtained 98 to 125 Primary Species Hypotheses (PSHs). The PSH partitions were subsequently visualized as a Klee diagram color map, allowing easy detection of the incongruences that were further evaluated individually with two other species delimitation models, General Mixed Yule Coalescent (GMYC) and Poisson Tree Processes (PTP). GMYC and PTP results confirmed the presence of 17 putative cryptic terebrid species in our dataset. The consensus of GMYC, PTP, and ABGD/Klee findings suggest the combination of ABGD and Klee diagrams is an effective approach for rapidly proposing primary species proxies in hyperdiverse groups and a reliable first step for macroscopic biodiversity assessment.     

Fast Census of Moth Diversity in the Neotropics: A Comparison of Field-Assigned Morphospecies and DNA Barcoding in Tiger Moths

The morphological species delimitations (i.e. morphospecies) have long been the best way to avoid the taxonomic impediment and compare insect taxa biodiversity in highly diverse tropical and subtropical regions. The development of DNA barcoding, however, has shown great potential to replace (or at least complement) the morphospecies approach, with the advantage of relying on automated methods implemented in computer programs or even online rather than in often subjective morphological features. We sampled moths extensively for two years using light traps in a patch of the highly endangered Atlantic Forest of Brazil to produce a nearly complete census of arctiines (Noctuoidea: Erebidae), whose species richness was compared using different morphological and molecular approaches (DNA barcoding). A total of 1,075 barcode sequences of 286 morphospecies were analyzed. Based on the clustering method Barcode Index Number (BIN) we found a taxonomic bias of approximately 30% in our initial morphological assessment. However, a morphological reassessment revealed that the correspondence between morphospecies and molecular operational taxonomic units (MOTUs) can be up to 94% if differences in genitalia morphology are evaluated in individuals of different MOTUs originated from the same morphospecies (putative cases of cryptic species), and by recording if individuals of different genders in different morphospecies merge together in the same MOTU (putative cases of sexual dimorphism). The results of two other clustering methods (i.e. Automatic Barcode Gap Discovery and 2% threshold) were very similar to those of the BIN approach. Using empirical data we have shown that DNA barcoding performed substantially better than the morphospecies approach, based on superficial morphology, to delimit species of a highly diverse moth taxon, and thus should be used in species inventories.