Exploring diversity in cryptorhynchine weevils (Coleoptera) using distance-, character- and tree-based species delineation

Delimiting species, revealing cryptic diversity, and population divergence in Qinghai-Tibet Plateau weevils through DNA barcoding

The Leptomias group represents one of the most diverse taxonomic group of weevils in the Qinghai-Tibet Plateau and its adjacent areas. Despite the potential of hidden diversity, relatively few comprehensive studies have been conducted on species diversity in this taxonomic group. In this study, we performed DNA barcoding analysis for species of the Leptomias group using a comprehensive DNA barcode dataset that included 476 sequences representing 54 morphospecies. Within the dataset, our laboratory contributed 474 sequences, and 390 sequences were newly generated for this study. The average Kimura 2-parameter distances among morphospecies and genera were 0.76% and 19.15%, respectively. In 94.4% of the species, the minimum interspecific distances exceeded the maximum intraspecific distances, indicating the presence of barcode gaps in most species of Leptomias group. The application of Automatic Barcode Gap Discovery, Assemble Species by Automatic Partitioning, Barcode Index Number, Bayesian Poisson tree processes, jMOTU, and Neighbor-joining tree methods revealed 45, 45, 63, 54, and 55 distinct clusters representing single species, respectively. Additionally, a total of four morphospecies, Leptomias kangmarensis, L. midlineatus, L. siahus, and L. sp.9RL, were found to be assigned to multiple subclade each, indicating the geographical divergences and the presence of cryptic diversity. Our findings of this study demonstrate that Qinghai-Tibet Plateau exhibits a higher species diversity of the Leptomias group, and it is imperative to investigate cryptic species within certain morphospecies using integrative taxonomic approaches in future studies. Moreover, the construction of a DNA barcode reference library presented herein establishes a robust foundational dataset to support forthcoming research on weevil taxonomy, phylogenetics, ecology, and evolution.

Molecular Weevil Identification Project: A thoroughly curated barcode release of 1300 Western Palearctic weevil species (Coleoptera, Curculionoidea)

The Molecular Weevil Identification project (MWI) studies the systematics of Western Palearctic weevils (superfamily Curculionoidea) in an integrative taxonomic approach of DNA barcoding, morphology and ecology. This barcode release provides almost 3600 curated CO1 sequences linked to morphological vouchers in about 1300 weevil species. The dataset is presented in statistical distance tables and as a Neighbour-Joining tree. Bayesian Inference trees are computed for the subfamilies Cryptorhynchinae, Apioninae and Ceutorhynchinae. Altogether, 18 unresolved taxonomic issues are discussed. A new barcode primer set is presented. Finally, we establish group-specific genetic distances for many weevil genera to serve as a tool in species delineation. These values are statistically based on distances between "good species" and their congeners. With this morphologically calibrated approach, we could resolve most alpha-taxonomic questions within the MWI project.

Identifying the Genetic Distance Threshold for Entiminae (Coleoptera: Curculionidae) Species Delimitation via COI Barcodes

The subfamily Entiminae is the largest group in the family Curculionidae, and it has long represented a challenge in traditional and molecular classification. Here, we analyzed intra- and interspecific genetic distances of 621 public COI barcode sequences (658bp) from 39 genera and 110 species of Entiminae, to determine parameters most congruent in retaining established species. We found that the mean intraspecific genetic distance (3.07%) was much smaller than the mean interspecific one (21.96%), but there is a wide range of overlap between intra- and interspecific genetic distances (0.77−18.01%), indicating that there is no consistent, universal barcoding gap. Specifically, DNA barcoding gap analysis for morphospecies revealed that 102 of 110 morphospecies had barcoding gaps, and 9.18% was the optimum threshold of genetic distances for 97 species delimitation. We further confirmed this threshold with barcodes from 27 morphologically identified specimens (including 21 newly reported barcodes) sequenced from five genera and seven species. We also identified thresholds to delimit congeneric species within 14 selected genera (species > 2), which varied from 7.42% (Trichalophus) to 13.48% (Barypeithes). We herein present optimal parameters for species identification in the Entiminae. Our study suggests that despite no universal genetic distance threshold value in subfamily Entiminae, 9.18% is optimal for most species. We recommend a wider sampling of geographic populations to better account for intraspecific distance variation, and that genetic distance thresholds for species delimitation should be refined at the genus level.

Two new species of Byrrhinus Motschulsky, 1858 (Coleoptera, Limnichidae, Limnichinae) from Negros, Philippines

Two new species of Limnichidae beetles, Byrrhinusnegrosensis sp. nov. and Byrrhinusvillarini sp. nov., are described from the Island of Negros in the Philippines. The adult specimens of the new species can be differentiated by patterns of body punctation, colour and orientation of elytral pubescence, posterolateral angle of pronotum, tarsomere length ratio and aedeagal form. Two clades, representing the two new species, were retrieved in the Maximum Likelihood gene tree using the 3'-end of the COI gene. Maximum genetic divergence within B.negrosensis sp. nov. and B.villarini sp. nov. were recorded to be 2.3% and 1.3%, respectively, while the mean interspecific divergence between the two new species was 19.7%. Morphological descriptions, digital photographs and COI sequences were provided for the two species. The state of knowledge of Byrrhinus is reviewed and an updated Philippine checklist is provided. By coupling morphological and molecular data, this paper provides the first additional new species of Philippine Byrrhinus in the last 28 years.

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.

Cytogenetic and Molecular Characterization of Three Mimetic Species of the Genus Alagoasa Bechyné 1955 (Coleoptera: Alticinae) from the Neotropical Region

Coleoptera is a mega-diverse order, but only about 1% of its species have been analyzed cytogenetically. In this order, the subfamily Alticinae presents many identification problems, mainly due to the occurrence of mimicry. The objective of this work was to cytogenetically characterize 3 very similar species of the genus Alagoasa (A. pantina, A.areata, and A.scissa). We used classical and molecular cytogenetic as well as molecular genetic techniques. All 3 species showed a diploid chromosome number of 2n = 22 (20+X+y), but differences in the morphology of the chromosomes. All had a meiotic formula of 2n = 10II+X+y and an X+y sex determination system with giant, fully asynaptic sex chromosomes, concordant characteristics observed in the subtribe Oedionychina. FISH demonstrated the presence of 18S and 5S rDNA clusters in 1 pair of autosomes, syntenic and colocalizing in the 3 analyzed species. However, in A. areata, heteromorphism between the cistrons was observed. The telomeric (TTAGG)n probe showed signals in all 3 species, with proximal signals in the X and dispersed signals in the y chromosome of A. areata, and 2 proximal signals in the X chromosome of A. scissa. Molecular analysis of the COI gene indicated that they are 3 distinct species, corroborating the observed cytogenetic characteristics.

Taxonomic reassessment of the genus Dichotomius (Coleoptera: Scarabaeinae) through integrative taxonomy

Dung beetles of the subfamily Scarabaeinae are widely recognised as important providers of multiple ecosystem services and are currently experiencing revisions that have improved our understanding of higher-level relationships in the subfamily. However, the study of phylogenetic relationships at the level of genus or species is still lagging behind. In this study we investigated the New World beetle genus Dichotomius, one of the richest within the New World Scarabaeinae, using the most comprehensive molecular and morphological dataset for the genus to date (in terms of number of species and individuals). Besides evaluating phylogenetic relationships, we also assessed species delimitation through a novel Bayesian approach (iBPP) that enables morphological and molecular data to be combined. Our findings support the monophyly of the genus Dichotomius but not that of the subgenera Selenocopris and Dichotomius sensu stricto (s.s). Also, our results do not support the recent synonymy of Selenocopris with Luederwaldtinia. Some species-groups within the genus were recovered, and seem associated with elevational distribution. Our species delimitation analyses were largely congruent irrespective of the set of parameters applied, but the most robust results were obtained when molecular and morphological data were combined. Although our current sampling and analyses were not powerful enough to make definite interpretations on the validity of all species evaluated, we can confidently recognise D. nisus, D. belus and D. mamillatus as valid and well differentiated species. Overall, our study provides new insights into the phylogenetic relationships and classification of dung beetles and has broad implications for their systematics and evolutionary analyses.

Species delimitation and mitogenome phylogenetics in the subterranean genus Pseudoniphargus (Crustacea: Amphipoda)

The amphi-Atlantic distributions exhibited by many thalassoid stygobiont (obligate subterranean) crustaceans have been explained by fragmentation by plate tectonics of ancestral shallow water marine populations. The amphipod stygobiont genus Pseudoniphargus is distributed across the Mediterranean region but also in the North Atlantic archipelagos of Bermuda, Azores, Madeira and the Canaries. We used species delimitation methods and mitogenome phylogenetic analyses to clarify the species diversity and evolutionary relationships within the genus and timing their diversification. Analyses included samples from the Iberian Peninsula, northern Morocco, the Balearic, Canarian, Azores and Madeira archipelagoes plus Bermuda. In most instances, morphological and molecular-based species delimitation analyses yielded consistent results. Notwithstanding, in a few cases either incipient speciation with no involvement of detectable morphological divergence or species crypticism were the most plausible explanations for the disagreement found between morphological and molecular species delimitations. Phylogenetic analyses based on a robust calibrated mitochondrial tree suggested that Pseudoniphargus lineages have a younger age than for other thalassoid amphipods displaying a disjunct distribution embracing both sides of the Atlantic Ocean. A major split within the family was estimated to occur at the Paleocene, when a lineage from Northern Iberian Peninsula diverged from the rest of pseudoniphargids. Species diversification in the peri-Mediterranean area was deduced to occur in early Miocene to Tortonian times, while in the Atlantic islands it started in the Pliocene. Our results show that the current distribution pattern of Pseudoniphargus resulted from a complex admix of relatively ancient vicariance events and several episodes of long- distance dispersal.

Integrative taxonomy by molecular species delimitation: multi-locus data corroborate a new species of Balkan Drusinae micro-endemics

BACKGROUND

Taxonomy offers precise species identification and delimitation and thus provides basic information for biological research, e.g. through assessment of species richness. The importance of molecular taxonomy, i.e., the identification and delimitation of taxa based on molecular markers, has increased in the past decade. Recently developed exploratory tools now allow estimating species-level diversity in multi-locus molecular datasets.

RESULTS

Here we use molecular species delimitation tools that either quantify differences in intra- and interspecific variability of loci, or divergence times within and between species, or perform coalescent species tree inference to estimate species-level entities in molecular genetic datasets. We benchmark results from these methods against 14 morphologically readily differentiable species of a well-defined subgroup of the diverse Drusinae subfamily (Trichoptera, Limnephilidae). Using a 3798 bp (6 loci) molecular data set we aim to corroborate a geographically isolated new species by integrating comparative morphological studies and molecular taxonomy.

CONCLUSIONS

Our results indicate that only multi-locus species delimitation provides taxonomically relevant information. The data further corroborate the new species Drusus zivici sp. nov. We provide differential diagnostic characters and describe the male, female and larva of this new species and discuss diversity patterns of Drusinae in the Balkans. We further discuss potential and significance of molecular species delimitation. Finally we argue that enhancing collaborative integrative taxonomy will accelerate assessment of global diversity and completion of reference libraries for applied fields, e.g., conservation and biomonitoring.

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.

Rarity and Incomplete Sampling in DNA-Based Species Delimitation

DNA-based species delimitation may be compromised by limited sampling effort and species rarity, including "singleton" representatives of species, which hampers estimates of intra- versus interspecies evolutionary processes. In a case study of southern African chafers (beetles in the family Scarabaeidae), many species and subclades were poorly represented and 48.5% of species were singletons. Using cox1 sequences from >500 specimens and ∼100 species, the Generalized Mixed Yule Coalescent (GMYC) analysis as well as various other approaches for DNA-based species delimitation (Automatic Barcode Gap Discovery (ABGD), Poisson tree processes (PTP), Species Identifier, Statistical Parsimony), frequently produced poor results if analyzing a narrow target group only, but the performance improved when several subclades were combined. Hence, low sampling may be compensated for by "clade addition" of lineages outside of the focal group. Similar findings were obtained in reanalysis of published data sets of taxonomically poorly known species assemblages of insects from Madagascar. The low performance of undersampled trees is not due to high proportions of singletons per se, as shown in simulations (with 13%, 40% and 52% singletons). However, the GMYC method was highly sensitive to variable effective population size ([Formula: see text]), which was exacerbated by variable species abundances in the simulations. Hence, low sampling success and rarity of species affect the power of the GMYC method only if they reflect great differences in [Formula: see text] among species. Potential negative effects of skewed species abundances and prevalence of singletons are ultimately an issue about the variation in [Formula: see text] and the degree to which this is correlated with the census population size and sampling success. Clade addition beyond a limited study group can overcome poor sampling for the GMYC method in particular under variable [Formula: see text] This effect was less pronounced for methods of species delimitation not based on coalescent models.

Comparison of five methods for delimitating species in Ophion Fabricius, a diverse genus of parasitoid wasps (Hymenoptera, Ichneumonidae)

DNA taxonomy has been proposed as a method to quickly assess diversity and species limits in highly diverse, understudied taxa. Here we use five methods for species delimitation and two genetic markers (COI and ITS2) to assess species diversity within the parasitoid genus, Ophion. We searched for compensatory base changes (CBC's) in ITS2, and determined that they are too rare to be of practical use in delimiting species in this genus. The other four methods used both COI and ITS2, and included distance-based (threshold analysis and ABGD) and tree-based (GMYC and PTP) models. We compared the results of these analyses to each other under various parameters and tested their performance with respect to 11 Nearctic species/morphospecies and 15 described Palearctic species. We also computed barcode accumulation curves of COI sequences to assess the completeness of sampling. The species count was highly variable depending on the method and parameters used, ranging from 47 to 168 species, with more conservative estimates of 89-121 species. Despite this range, many of the Nearctic test species were fairly robust with respect to method. We concluded that while there was often good congruence between methods, GMYC and PTP were less reliant on arbitrary parameters than the other two methods and more easily applied to genetic markers other than COI. However, PTP was less successful at delimiting test species than was GMYC. All methods, as well as the barcode accumulation curves, indicate that several Palearctic species remain undescribed and that we have scarcely begun to appreciate the Nearctic diversity within this genus.

Phylogenetic molecular species delimitations unravel potential new species in the pest genus Spodoptera Guenée, 1852 (Lepidoptera, Noctuidae)

Nowadays molecular species delimitation methods promote the identification of species boundaries within complex taxonomic groups by adopting innovative species concepts and theories (e.g. branching patterns, coalescence). As some of them can efficiently deal with large single-locus datasets, they could speed up the process of species discovery compared to more time consuming molecular methods, and benefit from the existence of large public datasets; these methods can also particularly favour scientific research and actions dealing with threatened or economically important taxa. In this study we aim to investigate and clarify the status of economically important moths species belonging to the genus Spodoptera (Lepidoptera, Noctuidae), a complex group in which previous phylogenetic analyses and integrative approaches already suggested the possible occurrence of cryptic species and taxonomic ambiguities. In this work, the effectiveness of innovative (and faster) species delimitation approaches to infer putative species boundaries has been successfully tested in Spodoptera, by processing the most comprehensive dataset (in terms of number of species and specimens) ever achieved; results are congruent and reliable, irrespective of the set of parameters and phylogenetic models applied. Our analyses confirm the existence of three potential new species clusters (for S. exigua (Hübner, 1808), S. frugiperda (J.E. Smith, 1797) and S. mauritia (Boisduval, 1833)) and support the synonymy of S. marima (Schaus, 1904) with S. ornithogalli (Guenée, 1852). They also highlight the ambiguity of the status of S. cosmiodes (Walker, 1858) and S. descoinsi Lalanne-Cassou & Silvain, 1994. This case study highlights the interest of molecular species delimitation methods as valuable tools for species discovery and to emphasize taxonomic ambiguities.

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.

Species delimitation in Neoplecostomus (Siluriformes: Loricariidae) using morphologic and genetic approaches

In the present study, we compare the results of alpha taxonomy (based on morphology), DNA Barcoding method with a 2% genetic divergence threshold, and the GMYC (General Model Yule Coalescent) technique to identify species clusters in Neoplecostomus. We used partial sequences of mitochondrial COI (Cytochrome oxidase subunit I) for 59 specimens representing 13 valid species (Neoplecostomus bandeirante, N. jaguari, N. langeanii, N. paranensis, N. yapo, N. botucatu, N. selenae, N. doceensis, N. corumba, N. ribeirensis, N. microps, N. espiritosantensis and N. franciscoensis) of Neoplecostomus collected in all its distribution area. For the analysis we used Bayesian inference of phylogeny with relaxed clock methods on an arbitrary timescale, using BEAST. The ultrametrics genes trees obtained for each tree prior evaluated (Yule, Birth-death and coalescent Population Size) were used in the GMYC analysis to identify a time in the tree when the branching rate shifts (in forward time) from Yule and Birth-Death (species) to a coalescent (population) process. We found that the GMYC model using the Yule prior identified 11 groups, closer to the current taxonomy (13 species). GMYC analyses using other ultrametric gene trees estimated under alternative prior such as Birth-Death and Coalescent Population Size, identified 9 groups, results observed in the traditional 2% genetic distance threshold, resulting in a low number of species recognized compared to the number of species identified with current taxonomy. Based on these results we conclude that the molecular approaches are helpful to distinguish species of Neoplecostomus, nevertheless it is important to combine molecular methodologies with current taxonomy in order to identify correctly species that recently originated.

Building-up of a DNA barcode library for true bugs (insecta: hemiptera: heteroptera) of Germany reveals taxonomic uncertainties and surprises

During the last few years, DNA barcoding has become an efficient method for the identification of species. In the case of insects, most published DNA barcoding studies focus on species of the Ephemeroptera, Trichoptera, Hymenoptera and especially Lepidoptera. In this study we test the efficiency of DNA barcoding for true bugs (Hemiptera: Heteroptera), an ecological and economical highly important as well as morphologically diverse insect taxon. As part of our study we analyzed DNA barcodes for 1742 specimens of 457 species, comprising 39 families of the Heteroptera. We found low nucleotide distances with a minimum pairwise K2P distance <2.2% within 21 species pairs (39 species). For ten of these species pairs (18 species), minimum pairwise distances were zero. In contrast to this, deep intraspecific sequence divergences with maximum pairwise distances >2.2% were detected for 16 traditionally recognized and valid species. With a successful identification rate of 91.5% (418 species) our study emphasizes the use of DNA barcodes for the identification of true bugs and represents an important step in building-up a comprehensive barcode library for true bugs in Germany and Central Europe as well. Our study also highlights the urgent necessity of taxonomic revisions for various taxa of the Heteroptera, with a special focus on various species of the Miridae. In this context we found evidence for on-going hybridization events within various taxonomically challenging genera (e.g. Nabis Latreille, 1802 (Nabidae), Lygus Hahn, 1833 (Miridae), Phytocoris Fallén, 1814 (Miridae)) as well as the putative existence of cryptic species (e.g. Aneurus avenius (Duffour, 1833) (Aradidae) or Orius niger (Wolff, 1811) (Anthocoridae)).

Phylogeny, historical biogeography and divergence time estimates of the genus Colophon Gray (Coleoptera : Lucanidae)

The flightless Cape High-mountain stag beetle genus Colophon (Coleoptera: Lucanidae) is studied. All species are endemic to the Cape Floristic Region of the Western Cape, South Africa. The study aimed to determine the specific and phylogenetic status of the described species of Colophon and to determine the main factors driving their evolution by testing the hypothesis of a lowland origin. This was achieved by analysing the mitochondrial COI, 16S rRNA and nuclear CAD genes, using a maximum likelihood and Bayesian approach. Timing of key biogeographical events in the diversification of Colophon was estimated in BEAST. The combined molecular dataset supports the described species. High genetic divergence was found between the described taxa. The mean estimated divergence of the genus was mid-Cretaceous, with a split into two lineages during the early Paleocene to mid Eocene. Species divergence was shown to have occurred during mid to late Miocene. Most species of Colophon showed an allopatric distribution, although contact zones between geographically adjacent species are likely. Climate seems to be the main driving factor behind Colophon evolution and the hypothesis of a lowland origin appears to be supported.

How to describe a cryptic species? Practical challenges of molecular taxonomy

BACKGROUND

Molecular methods of species delineation are rapidly developing and widely considered as fast and efficient means to discover species and face the 'taxonomic impediment' in times of biodiversity crisis. So far, however, this form of DNA taxonomy frequently remains incomplete, lacking the final step of formal species description, thus enhancing rather than reducing impediments in taxonomy. DNA sequence information contributes valuable diagnostic characters and -at least for cryptic species - could even serve as the backbone of a taxonomic description. To this end solutions for a number of practical problems must be found, including a way in which molecular data can be presented to fulfill the formal requirements every description must meet. Multi-gene barcoding and a combined molecular species delineation approach recently revealed a radiation of at least 12 more or less cryptic species in the marine meiofaunal slug genus Pontohedyle (Acochlidia, Heterobranchia). All identified candidate species are well delimited by a consensus across different methods based on mitochondrial and nuclear markers.

RESULTS

The detailed microanatomical redescription of Pontohedyle verrucosa provided in the present paper does not reveal reliable characters for diagnosing even the two major clades identified within the genus on molecular data. We thus characterize three previously valid Pontohedyle species based on four genetic markers (mitochondrial cytochrome c oxidase subunit I, 16S rRNA, nuclear 28S and 18S rRNA) and formally describe nine cryptic new species (P. kepii sp. nov., P. joni sp. nov., P. neridae sp. nov., P. liliae sp. nov., P. wiggi sp. nov., P. wenzli sp. nov., P. peteryalli sp. nov., P. martynovi sp. nov., P. yurihookeri sp. nov.) applying molecular taxonomy, based on diagnostic nucleotides in DNA sequences of the four markers. Due to the minute size of the animals, entire specimens were used for extraction, consequently the holotype is a voucher of extracted DNA ('DNA-type'). We used the Character Attribute Organization System (CAOS) to determine diagnostic nucleotides, explore the dependence on input data and data processing, and aim for maximum traceability in our diagnoses for future research. Challenges, pitfalls and necessary considerations for applied DNA taxonomy are critically evaluated.

CONCLUSIONS

To describe cryptic species traditional lines of evidence in taxonomy need to be modified. DNA sequence information, for example, could even serve as the backbone of a taxonomic description. The present contribution demonstrates that few adaptations are needed to integrate into traditional taxonomy novel diagnoses based on molecular data. The taxonomic community is encouraged to join the discussion and develop a quality standard for molecular taxonomy, ideally in the form of an automated final step in molecular species delineation procedures.

Integrative analysis of DNA phylogeography and morphology of the European rose chafer (Cetonia aurata) to infer species taxonomy and patterns of postglacial colonisation in Europe

Integrative taxonomy has been proposed as a framework to unify new conceptual and methodological developments in quantitative assessment of trait variation used in species delimitation, but empirical studies in this young branch of systematics are rare. Here we use standard phylogenetic and parsimony network analyses on nuclear and mitochondrial DNA (Cox1, ITS1) of 230 individuals from 65 European sampling sites in order to deduce population structure of Cetonia beetles from geno- and haplotypes. Statistical measures of population differentiation are inferred on genealogical and geographical scales to test hypotheses about species limits and population history. By combining results of phylogenetic structure with features of morphology, including genital shape morphometrics and discrete external body characters, as well as with measures of population genetics, we attempt to integrate the results as a test of the validity of species limits, in particular of currently recognised subspecies. Despite high Cox1 divergence between some haplotype lineages, even some sympatric lineages (9%, e.g. N2 vs. N4), nDNA and morphology, as well as pattern of geographical and genealogical divergence measured by AMOVA analysis did not support the hypothesis of separate species. Highest divergence in nuclear markers was found among Italian populations of C. aurata pisana and C. a. sicula, and moderately high fixation indices along measurable morphological divergence suggest the correctness of their status as 'subspecies'. Divergence time estimates of the lineages suggest a glacial divergence in different refugia between the major haplogroups, while population differentiation in mtDNA among these was primarily attributable to restricted gene flow caused by geographic isolation.

Barcoding against a paradox? Combined molecular species delineations reveal multiple cryptic lineages in elusive meiofaunal sea slugs

BACKGROUND

Many marine meiofaunal species are reported to have wide distributions, which creates a paradox considering their hypothesized low dispersal abilities. Correlated with this paradox is an especially high taxonomic deficit for meiofauna, partly related to a lower taxonomic effort and partly to a high number of putative cryptic species. Molecular-based species delineation and barcoding approaches have been advocated for meiofaunal biodiversity assessments to speed up description processes and uncover cryptic lineages. However, these approaches show sensitivity to sampling coverage (taxonomic and geographic) and the success rate has never been explored on mesopsammic Mollusca.

RESULTS

We collected the meiofaunal sea-slug Pontohedyle (Acochlidia, Heterobranchia) from 28 localities worldwide. With a traditional morphological approach, all specimens fall into two morphospecies. However, with a multi-marker genetic approach, we reveal multiple lineages that are reciprocally monophyletic on single and concatenated gene trees in phylogenetic analyses. These lineages are largely concordant with geographical and oceanographic parameters, leading to our primary species hypothesis (PSH). In parallel, we apply four independent methods of molecular based species delineation: General Mixed Yule Coalescent model (GMYC), statistical parsimony, Bayesian Species Delineation (BPP) and Automatic Barcode Gap Discovery (ABGD). The secondary species hypothesis (SSH) is gained by relying only on uncontradicted results of the different approaches ('minimum consensus approach'), resulting in the discovery of a radiation of (at least) 12 mainly cryptic species, 9 of them new to science, some sympatric and some allopatric with respect to ocean boundaries. However, the meiofaunal paradox still persists in some Pontohedyle species identified here with wide coastal and trans-archipelago distributions.

CONCLUSIONS

Our study confirms extensive, morphologically cryptic diversity among meiofauna and accentuates the taxonomic deficit that characterizes meiofauna research. We observe for Pontohedyle slugs a high degree of morphological simplicity and uniformity, which we expect might be a general rule for meiofauna. To tackle cryptic diversity in little explored and hard-to-sample invertebrate taxa, at present, a combined approach seems most promising, such as multi-marker-barcoding (i.e., molecular systematics using mitochondrial and nuclear markers and the criterion of reciprocal monophyly) combined with a minimum consensus approach across independent methods of molecular species delineation to define candidate species.