The effect of geographical scale of sampling on DNA barcoding

"Dark taxonomy": A new protocol for overcoming the taxonomic impediments for dark taxa and broadening the taxon base for biodiversity assessment

We are entering the sixth mass extinction with little data for "dark taxa", although they comprise most species. Much of the neglect is due to the fact that conventional taxonomic methods struggle with handling thousands of specimens belonging to hundreds of species. We thus here propose a new strategy that we call "dark taxonomy". It addresses (i) taxonomic impediments, (ii) the lack of biodiversity baselines and (iii) the low impact of revisionary research. Taxonomic impediments are reduced by carrying out revisions at small geographic scales to keep the number of specimens low. The risk of taxonomic error is reduced by delimiting species based on two types of data. We furthermore show that dark taxonomy can yield important biodiversity baseline data by using samples obtained with biomonitoring traps. Lastly, we argue that the impact of revisionary research can be improved by publishing two papers addressing different readerships. The principles of dark taxonomy are illustrated by our taxonomic treatment of Singapore's fungus gnats (Mycetophilidae) based only on Malaise trap samples. We show that a first batch of specimens (N = 1454) contains 120 species, of which 115 are new to science, thus reducing taxonomic impediments by increasing the number of described Oriental species by 25%. Species delimitation started with using DNA barcodes to estimate the number of Molecular Operational Taxonomic Units (MOTUs) before "LIT" (Large-scale Integrative Taxonomy) was used to obtain the species boundaries for the 120 species by integrating morphological and molecular data. To test the taxonomic completeness of the revision, we next analysed a second batch of 1493 specimens and found that >97% belonged to the 120 species delimited based on the first batch. Indeed, the second batch only contained 18 new and rare MOTUs, i.e. our study suggests that a single revision can simultaneously yield the names for all important species and relevant biodiversity baseline data. Overall, we believe that "dark taxonomy" can quickly ready a large unknown taxon for biomonitoring.

A review of the species groups of the Western Hemisphere Onthophagus Latreille (Coleoptera: Scarabaeidae: Scarabaeinae) using COI barcoding and gene trees

Species groups of Western Hemispheric Onthophagus Latreille (Coleoptera: Scarabaeidae: Scarabaeinae: Onthophagini) are suggested using COI barcoding and gene trees and supported by congruence with external morphology, behavior, ecology, and biogeographic evidence. New species groups, complexes, and taxonomic statuses are offered, and other preexisting proposals are confirmed. No barcoding gap w as found between the intragroup and intergroup genetic distance blocks, but the average intragroup (8.38%) and intergroup (13.88%) Kimura-two-parameter distances are statistically different. The following seven preexisting species groups were supported by the congruence between the mtDNA barcode analysis and other independent evidence: O. chevrolati, O. clypeatus, O. dicranius, O. gazellinus, O. hircus, O. landolti, and O. mexicanus. Eight new species groups are suggested: O. crinitus, O. curvicornis, O. eulophus, O. hecate, O. hoepfneri, O. marginatus, O. nasutus, and O. velutinus. Possible behavioral/ecological adaptations of morphological characters are also discussed. New biogeographic and evolutionary hypotheses are also advanced. An identification key for species groups is presented.

A barcode database for insects associated with the spread of the Cocoa Swollen Shoot Virus Disease in Côte d'Ivoire

Swollen Shoot is a viral disease affecting cocoa trees, transmitted by several species of mealybugs (Insecta, Hemiptera, Sternorrhyncha, Pseudococcidae). These insects maintain trophobiotic relationships with a complex and species-rich assemblage of ants protecting them and natural enemies controlling their populations. Here, we provide a curated DNA barcode database to characterise this insect community. Systematic observation of 7,500 cocoa trees was conducted, coupled with the collection of mealybug colonies and associated insect communities (parasitoids, predators and ants). Natural enemies were reared from mealybug colonies collected from 1,430 cocoa trees. Specimens were identified morphologically and sequenced for fragments of the standard DNA barcode region of the COI. We recovered 17 species of mealybugs from the family Pseudococcidae. Amongst these species, eight are new to the Ivorian cocoa orchard: Dysmicoccusneobrevipes Beardsley, Ferrisiadasylirii (Cockerell), Maconellicoccusugandae (Laing), Paracoccusmarginatus Williams & Granara de Willink, Phenacoccussolenopsis Tinsley, Planococcusminor (Maskell), Pseudococcusconcavocerarii James and Pseudococcusocciduus De Lotto. Three of these species were identified for the first time in cocoa orchards in Africa: D.neobrevipes, Fe.dasylirii and Ph.solenopsis. A total of 54 ant species were identified and represented the first record of these species associated with mealybug colonies in cocoa in Côte d'Ivoire. Amongst the species associated with the mealybugs, 22 primary parasitoids, eight hyperparasitoids, 11 ladybirds beetles (Coccinellidae), seven gall midges (Cecidomyidae), one predatory lepidopteran species and four spider species were identified. Nine species of mealybugs parasitoids are newly recorded in the African cocoa orchards: Acerophagusaff.dysmicocci, Aloencyrtus sp., Anagyruskamali, Anagyrusaff.pseudococci, Aenasiusadvena, Clauseniaaff.corrugata, Gyranusoideaaff.tebygi, Zaplatycerusaff.natalensis (Encyrtidae) and Coccophaguspulvinariae (Aphelinidae) and one hyperparasitoid, Pachyneuronmuscarum (Pteromalidae). For Côte d'Ivoire in particular, besides the previously mentioned nine parasitoids and one hyperparasitoid, five additional species are recorded for the first time, including four primary parasitoids, Blepyrusinsularis (Encyrtidae), Clauseniacorrugata (Encyrtidae), Clausenia sp. (Encyrtidae), and Coccidoctonuspseudococci (Encyrtidae) and one hyperparasitoid, Cheiloneuruscyanonotus (Encyrtidae). These results significantly enhance the knowledge of the diversity of the entomofauna associated with Swollen Shoot disease and pave the way for developing control methods based on the natural regulation of its mealybug (Pseudococcidae) vectors.

Phylogeography of the ambrosia beetle Euwallacea interjectus (Coleoptera: Curculionidae: Scolytinae): an emerging poplar pest and its Fusarium mutualists from poplar plantations in China

Native to Asia, Euwallacea interjectus (Blandford) (Coleoptera: Curculionidae: Scolytinae) is a destructive and invasive pest of live trees, and now it has been found in the United States and Argentina. In recent years, this pest appeared in high densities in poplar monocultures from Eastern China (Jiangsu and Shanghai) and Argentina and caused significant poplar mortality. However, the origin of the pests related to tree damage and the Fusarium mutualists from some poplar zones in China remained unclear. Here, we provided a broader phylogeographic analysis of E. interjectus based on the mitochondrial gene (cytochrome c oxidase I) to determine the global genetic structure of this species. Five mitochondrial lineages were found in the native area. Populations introduced to the United States were originated from 4 localities. The Argentine population was derived from Japan. The species was observed with strikingly high level of cytochrome c oxidase I intraspecific divergence that exceeded interspecific divergence, but the high intraspecific variation was correlated with geographical locations among the native populations. Two nuclear genes (arginine kinase and carbamoyl-phosphate synthetase 2-aspartate transcarbamylase-dihydroorotase) were more conservative, and intraspecific differences were lower than interspecific differences. The mitochondrial genetic variation was probably caused by evolution of lineages among geographically isolated populations. But it is immature to infer the existence of cryptic species based on cytochrome c oxidase I differences. All samples collected from poplar populations were indigenous and formed close relationship with a specimen from eastern and southern China. Surprisingly, pests from poplar populations in Jiangsu and Shanghai showed different haplotypes and mutualists. This suggested that the control strategies should consider the genetic and mutualistic diversity of beetles at different poplar localities.

Performance of DNA metabarcoding, standard barcoding and morphological approaches in the identification of insect biodiversity

For two decades, DNA barcoding and, more recently, DNA metabarcoding have been used for molecular species identification and estimating biodiversity. Despite their growing use, few studies have systematically evaluated these methods. This study aims to evaluate the efficacy of barcoding methods in identifying species and estimating biodiversity, by assessing their consistency with traditional morphological identification and evaluating how assignment consistency is influenced by taxonomic group, sequence similarity thresholds and geographic distance. We first analysed 951 insect specimens across three taxonomic groups: butterflies, bumblebees and parasitic wasps, using both morphological taxonomy and single-specimen COI DNA barcoding. An additional 25,047 butterfly specimens were identified by COI DNA metabarcoding. Finally, we performed a systematic review of 99 studies to assess average consistency between insect species identity assigned via morphology and COI barcoding and to examine the distribution of research effort. Species assignment consistency was influenced by taxonomic group, sequence similarity thresholds and geographic distance. An average assignment consistency of 49% was found across taxonomic groups, with parasitic wasps displaying lower consistency due to taxonomic impediment. The number of missing matches doubled with a 100% sequence similarity threshold and COI intraspecific variation increased with geographic distance. Metabarcoding results aligned well with morphological biodiversity estimates and a strong positive correlation between sequence reads and species abundance was found. The systematic review revealed an 89% average consistency and also indicated taxonomic and geographic biases in research effort. Together, our findings demonstrate that while problems persist, barcoding approaches offer robust alternatives to traditional taxonomy for biodiversity assessment.

Genetic barcodes for species identification and phylogenetic estimation in ghost spiders (Araneae: Anyphaenidae: Amaurobioidinae)

The identification of spider species presents many challenges, since in most cases the characters used are from genital structures that are only fully developed in the adult stage, hence the identification of immatures is most often not possible. Additionally, these structures usually also present some intra-specific variability, which in some cases makes the identification of closely related species difficult. The genetic barcode technique (DNA barcodes), based on sequencing of the mitochondrial marker cytochrome c oxidase subunit I (COI ), has proven a useful, complementary tool to overcome these limitations. In this work, the contribution of DNA barcoding to the taxonomy of the subfamily Amaurobioidinae is explored using the refined single linkage analysis (RESL) algorithm for the delimitation of operational taxonomic units (OTUs), in comparison with the assemble species by automatic partitioning (ASAP) algorithm, and presented in conjunction with an updated molecular phylogenetic analysis of three other markers (28S rRNA, 16S rRNA, Histone H3 ), in addition to COI . Of a total of 97 included species identified by morphology, 82 species were concordant with the operational taxonomic units obtained from RESL, representing an 85% correspondence between the two methods. Similar results were obtained using the ASAP algorithm. Previous observations of morphological variation within the same species are supported, and this technique provides new information on genetic structure and potentially cryptic species. Most of the discrepancies between DNA barcoding and morphological identification are explained by low geographic sampling or by divergent or geographically structured lineages. After the addition of many specimens with only COI data, the multi-marker phylogenetic analysis is consistent with previous results and the support is improved. The markers COI , closely followed by 28S , are the most phylogenetically informative. We conclude that the barcode DNA technique is a valuable source of data for the delimitation of species of Amaurobioidinae, in conjunction with morphological and geographic data, and it is also useful for the detection of cases that require a more detailed and meticulous study.

DNA barcoding of southern African mammal species and construction of a reference library for forensic application

Combating wildlife crimes in South Africa requires accurate identification of traded species and their products. Diagnostic morphological characteristics needed to identify species are often lost when specimens are processed and customs officials lack the expertise to identify species. As a potential solution, DNA barcoding can be used to identify morphologically indistinguishable specimens in forensic cases. However, barcoding is hindered by the reliance on comprehensive, validated DNA barcode reference databases, which are currently limited. To overcome this limitation, we constructed a barcode library of cytochrome c oxidase subunit 1 and cytochrome b sequences for threatened and protected mammals exploited in southern Africa. Additionally, we included closely related or morphologically similar species and assessed the database's ability to identify species accurately. Published southern African sequences were incorporated to estimate intraspecific and interspecific variation. Neighbor-joining trees successfully discriminated 94%-95% of the taxa. However, some widespread species exhibited high intraspecific distances (>2%), suggesting geographic sub-structuring or cryptic speciation. Lack of reliable published data prevented the unambiguous discrimination of certain species. This study highlights the efficacy of DNA barcoding in species identification, particularly for forensic applications. It also highlights the need for a taxonomic re-evaluation of certain widespread species and challenging genera.

Age-dependent variation of aedeagal morphology in Agabusuliginosus and the status of A.lotti (Coleoptera, Dytiscidae)

A doubt has arisen about the taxonomic status of Agabuslotti within the Agabusuliginosus species group due to morphological similarities and lack of molecular data. In this study, a comprehensive morphological and molecular analysis of specimens from Central Europe was conducted, focusing on the Hungarian population. Morphological comparisons of genital structures revealed age-dependent variations, suggesting a gradual transition from A.lotti to A.uliginosus. Molecular analysis of COI sequences further supported this hypothesis, showing minimal genetic differences among most specimens, with only one individual exhibiting distinctiveness. Therefore, A.lotti syn. nov. must be regarded as a junior synonym of A.uliginosus. Our findings also highlight the need for additional multi-marker studies covering a broader geographic range and including both molecular and morphological approaches to elucidate the taxonomic and phylogenetic relationships within this species group. The inclusion of Hungarian samples notably enriched the diversity of haplotypes, emphasizing the importance of expanding sampling efforts in future research.

Genomic, morphological and physiological data support fast ecotypic differentiation and incipient speciation in an alpine diving beetle

An intricate interplay between evolutionary and demographic processes has frequently resulted in complex patterns of genetic and phenotypic diversity in alpine lineages, posing serious challenges to species delimitation and biodiversity conservation planning. Here we integrate genomic data, geometric morphometric analyses and thermal tolerance experiments to explore the role of Pleistocene climatic changes and adaptation to alpine environments on patterns of genomic and phenotypic variation in diving beetles from the taxonomically complex Agabus bipustulatus species group. Genetic structure and phylogenomic analyses revealed the presence of three geographically cohesive lineages, two representing trans-Palearctic and Iberian populations of the elevation-generalist A. bipustulatus and another corresponding to the strictly-alpine A. nevadensis, a narrow-range endemic taxon from the Sierra Nevada mountain range in southeastern Iberia. The best-supported model of lineage divergence, along with the existence of pervasive genetic introgression and admixture in secondary contact zones, is consistent with a scenario of population isolation and connectivity linked to Quaternary climatic oscillations. Our results suggest that A. nevadensis is an alpine ecotype of A. bipustulatus, whose genotypic, morphological and physiological differentiation likely resulted from an interplay between population isolation and local altitudinal adaptation. Remarkably, within the Iberian Peninsula, such ecotypic differentiation is unique to Sierra Nevada populations and has not been replicated in other alpine populations of A. bipustulatus. Collectively, our study supports fast ecotypic differentiation and incipient speciation processes within the study complex and points to Pleistocene glaciations and local adaptation along elevational gradients as key drivers of biodiversity generation in alpine environments.

phyloBARCODER: A Web Tool for Phylogenetic Classification of Eukaryote Metabarcodes Using Custom Reference Databases

We developed phyloBARCODER (https://github.com/jun-inoue/phyloBARCODER), a new web tool that can identify short DNA sequences to the species level using metabarcoding. phyloBARCODER estimates phylogenetic trees based on the uploaded anonymous DNA sequences and reference sequences from databases. Without such phylogenetic contexts, alternative, similarity-based methods independently identify species names and anonymous sequences of the same group by pairwise comparisons between queries and database sequences, with the caveat that they must match exactly or very closely. By putting metabarcoding sequences into a phylogenetic context, phyloBARCODER accurately identifies (i) species or classification of query sequences and (ii) anonymous sequences associated with the same species or even with populations of query sequences, with clear and accurate explanations. Version 1 of phyloBARCODER stores a database comprising all eukaryotic mitochondrial gene sequences. Moreover, by uploading their own databases, phyloBARCODER users can conduct species identification specialized for sequences obtained from a local geographic region or those of nonmitochondrial genes, e.g. ITS or rbcL.

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.

Towards holistic insect monitoring: species discovery, description, identification and traits for all insects

Holistic insect monitoring needs scalable techniques to overcome taxon biases, determine species abundances, and gather functional traits for all species. This requires that we address taxonomic impediments and the paucity of data on abundance, biomass and functional traits. We here outline how these data deficiencies could be addressed at scale. The workflow starts with large-scale barcoding (megabarcoding) of all specimens from mass samples obtained at biomonitoring sites. The barcodes are then used to group the specimens into molecular operational taxonomic units that are subsequently tested/validated as species with a second data source (e.g. morphology). New species are described using barcodes, images and short diagnoses, and abundance data are collected for both new and described species. The specimen images used for species discovery then become the raw material for training artificial intelligence identification algorithms and collecting trait data such as body size, biomass and feeding modes. Additional trait data can be obtained from vouchers by using genomic tools developed by molecular ecologists. Applying this pipeline to a few samples per site will lead to greatly improved insect monitoring regardless of whether the species composition of a sample is determined with images, metabarcoding or megabarcoding. This article is part of the theme issue 'Towards a toolkit for global insect biodiversity monitoring'.

Utilising public sequence databases to investigate genetic diversity of stoneflies in Medvednica Nature Park

In Medvednica Nature Park, near Croatia's capital Zagreb, urbanisation significantly impacts the fauna. Comprehensive field research has never been conducted in this area, despite the presence of diverse microhabitats and the discovery of several rare species previously unknown in the Croatian fauna. This study provides the Park with first insight into the genetic and morphological diversity of stoneflies, one of the most endangered groups of organisms. Phylogenetic reconstructions and species delineation methods revealed intraspecific haplotype variation in most species (e.g. Brachypteraseticornis, Isoperlagrammatica and Leuctrabraueri), except for Leuctraprima. Additionally, our study has identified isolated populations that merit further in-depth investigation concerning morphology, genetics and ecology.

The genomic diversity of the Eliurus genus in northern Madagascar with a putative new species

Madagascar exhibits extraordinarily high level of species richness and endemism, while being severely threatened by habitat loss and fragmentation (HL&F). In front of these threats to biodiversity, conservation effort can be directed, for instance, in the documentation of species that are still unknown to science, or in investigating how species respond to HL&F. The tufted-tail rats genus (Eliurus spp.) is the most speciose genus of endemic rodents in Madagascar, with 13 described species, which occupy two major habitat types: dry or humid forests. The large species diversity and association to specific habitat types make the Eliurus genus a suitable model for investigating species adaptation to new environments, as well as response to HL&F (dry vs humid). In the present study, we investigated Eliurus spp. genomic diversity across northern Madagascar, a region covered by both dry and humid fragmented forests. From the mitochondrial DNA (mtDNA) and nuclear genomic (RAD-seq) data of 124 Eliurus individuals sampled in poorly studied forests of northern Madagascar, we identified an undescribed Eliurus taxon (Eliurus sp. nova). We tested the hypothesis of a new Eliurus species using several approaches: i) DNA barcoding; ii) phylogenetic inferences; iii) species delimitation tests based on the Multi-Species Coalescent (MSC) model, iv) genealogical divergence index (gdi); v) an ad-hoc test of isolation-by-distance within versus between sister-taxa, vi) comparisons of %GC content patterns and vii) morphological analyses. All analyses support the recognition of the undescribed lineage as a putative distinct species. In addition, we show that Eliurus myoxinus, a species known from the dry forests of western Madagascar, is, surprisingly, found mostly in humid forests in northern Madagascar. In conclusion, we discuss the implications of such findings in the context of Eliurus species evolution and diversification, and use the distribution of northern Eliurus species as a proxy for reconstructing past changes in forest cover and vegetation type in northern Madagascar.

Mitochondrial genome comparison and phylogenetic position of Fannia pusio among the Calyptratae flies

Fannia pusio, the chicken dung fly species, remains unexplored despite its forensic, sanitary, and veterinary importance in the Nearctic and Neotropical regions. In this study, we obtained the complete mitochondrial genome of Fannia pusio for the first time using next-generation sequencing. We compared it with previously published mitogenomes of the genus from the Palearctic region, and its phylogenetic position was studied based on the concatenated protein-coding genes (PCGs) dataset of Calyptratae flies. The circular mitochondrial genome of F. pusio is 16,176 bp in length, with a high A + T content (78.3%), whose gene synteny, codon usage analysis, and amino acid frequency are similar to previously reported Fannia mitogenomes. All PCGs underwent purifying selection except the nad2 gene. Interspecific K2P distances of PCGs of Fannia yielded an average of 12.4% (8.1%-21.1%). The Fannia genus is monophyletic and closely related to Muscidae based on molecular data. Further taxonomic sampling is required to deep into the phylogenetic relationships of the originally proposed species-groups and subgroups within the genus. These results provide a valuable dataset for studying the mitochondrial genome evolution and a resource for the taxonomy and systematics of Fannia.

Cryptic Diversity in Sympatric Migonemyia migonei (Diptera: Psychodidae), Eventual Meaning for Leishmaniasis Transmission

Migonemyia migonei (Franҫa, 1920) (Diptera: Psychodidae) belongs to the subfamily Phlebotominae, of epidemiological importance due to its role as a vector in leishmaniasis transmission cycles and its broad geographic distribution in South America. Few morphometric and genetic studies have demonstrated the existence of variability among geographically distant populations in Brazil. The aim of the study was to estimate the genetic distance within the morphospecies Mg. migonei through the analysis of cytochrome C oxidase subunit I (COI) sequences of specimens captured in Argentina and those available in online databases. The COI sequences from specimens collected in different localities of Argentina and sequences available in online databases were utilized. Genetic distances were analyzed and a median-joining haplotype network was constructed. Finally, phylogenetic reconstruction was performed according to Bayesian inference. The analyses led to the identification of at least two haplogroups: haplogroup I with sequences of specimens from Colombia, Brazil and Argentina, and haplogroup II with sequences of specimens from Argentina. Interestingly, specimens from Argentina whose haplotypes corresponded to both haplogroups, were collected in sympatry. The results suggest that Mg. migonei could be a species complex with at least two distinct members. This hypothesis could explain the known characteristics of adaptability and vector permissiveness of the species, as the putative cryptic species of the complex could differ in traits of epidemiological importance.

Integrative species delimitation uncovers hidden diversity within the Pithecopus hypochondrialis species complex (Hylidae, Phyllomedusinae) and its phylogeography reveals Plio-Pleistocene connectivity among Neotropical savannas

Despite their limited vagility and pronounced habitat heterogeneity in the tropics, many anuran species have unexpectedly extensive geographic ranges. One prominent example of this phenomenon is Pithecopus hypochondrialis, which is found in the Cerrado, Guianan savanna, and Llanos domains, as well as isolated tracts of savanna and open habitat within the Amazon Forest. The present study employs an integrative species delimitation approach to test the hypothesis that P. hypochondrialis is in fact a species complex. We also reconstruct the relationships among the lineages delimited here and other Pithecopus species. In this study, we employ Ecological Niche Modelling (ENM) and spatiotemporal phylogeographic reconstruction approaches to evaluate a multitude of scenarios of connectivity across the Neotropical savannas. We identified three divergent lineages, two of which have been described previously. The lineages were allocated to a lowland Pithecopus clade, although the relationships among these lineages are weakly supported. Both the ENM and the phylogeographic reconstruction highlight the occurrence of periods of connectivity among the Neotropical savannas over the course of the Pliocene and Pleistocene epochs. These processes extended from eastern Amazonia to the northern coast of Brazil. The findings of the present study highlight the presence of hidden diversity within P. hypochondrialis, and reinforce the need for a comprehensive taxonomic review. These findings also indicate intricate and highly dynamic patterns of connectivity across the Neotropical savannas that date back to the Pliocene.

A Measure of the DNA Barcode Gap for Applied and Basic Research

DNA barcoding has largely established itself as a mainstay for rapid molecular taxonomic identification in both academic and applied research. The use of DNA barcoding as a molecular identification method depends on a "DNA barcode gap"-the separation between the maximum within-species difference and the minimum between-species difference. Previous work indicates the presence of a gap hinges on sampling effort for focal taxa and their close relatives. Furthermore, both theory and empirical work indicate a gap may not occur for related pairs of biological species. Here, we present a novel evaluation approach in the form of an easily calculated set of nonparametric metrics to quantify the extent of proportional overlap in inter- and intraspecific distributions of pairwise differences among target species and their conspecifics. The metrics are based on a simple count of the number of overlapping records for a species falling within the bounds of maximum intraspecific distance and minimum interspecific distance. Our approach takes advantage of the asymmetric directionality inherent in pairwise genetic distance distributions, which has not been previously done in the DNA barcoding literature. We apply the metrics to the predatory diving beetle genus Agabus as a case study because this group poses significant identification challenges due to its morphological uniformity despite both relative sampling ease and well-established taxonomy. Results herein show that target species and their nearest neighbor species were found to be tightly clustered and therefore difficult to distinguish. Such findings demonstrate that DNA barcoding can fail to fully resolve species in certain cases. Moving forward, we suggest the implementation of the proposed metrics be integrated into a common framework to be reported in any study that uses DNA barcoding for identification. In so doing, the importance of the DNA barcode gap and its components for the success of DNA-based identification using DNA barcodes can be better appreciated.

Species Diagnosis and DNA Taxonomy

The use of DNA has helped to improve and speed up species identification and delimitation. However, it also provides new challenges to taxonomists. Incongruence of outcome from various markers and delimitation methods, bias from sampling and skewed species distribution, implemented models, and the choice of methods/priors may mislead results and also may, in conclusion, increase elements of subjectivity in species taxonomy. The lack of direct diagnostic outcome from most contemporary molecular delimitation approaches and the need for a reference to existing and best sampled trait reference systems reveal the need for refining the criteria of species diagnosis and diagnosability in the current framework of nomenclature codes and good practices to avoid nomenclatorial instability, parallel taxonomies, and consequently more and new taxonomic impediment.

Delimiting Species with Single-Locus DNA Sequences

DNA sequences are increasingly used for large-scale biodiversity inventories. Because these genetic data avoid the time-consuming initial sorting of specimens based on their phenotypic attributes, they have been recently incorporated into taxonomic workflows for overlooked and diverse taxa. Major statistical developments have accompanied this new practice, and several models have been proposed to delimit species with single-locus DNA sequences. However, proposed approaches to date make different assumptions regarding taxon lineage history, leading to strong discordance whenever comparisons are made among methods. Distance-based methods, such as Automatic Barcode Gap Discovery (ABGD) and Assemble Species by Automatic Partitioning (ASAP), rely on the detection of a barcode gap (i.e., the lack of overlap in the distributions of intraspecific and interspecific genetic distances) and the associated threshold in genetic distances. Network-based methods, as exemplified by the REfined Single Linkage (RESL) algorithm for the generation of Barcode Index Numbers (BINs), use connectivity statistics to hierarchically cluster-related haplotypes into molecular operational taxonomic units (MOTUs) which serve as species proxies. Tree-based methods, including Poisson Tree Processes (PTP) and the General Mixed Yule Coalescent (GMYC), fit statistical models to phylogenetic trees by maximum likelihood or Bayesian frameworks.Multiple webservers and stand-alone versions of these methods are now available, complicating decision-making regarding the most appropriate approach to use for a given taxon of interest. For instance, tree-based methods require an initial phylogenetic reconstruction, and multiple options are now available for this purpose such as RAxML and BEAST. Across all examined species delimitation methods, judicious parameter setting is paramount, as different model parameterizations can lead to differing conclusions. The objective of this chapter is to guide users step-by-step through all the procedures involved for each of these methods, while aggregating all necessary information required to conduct these analyses. The "Materials" section details how to prepare and format input files, including options to align sequences and conduct tree reconstruction with Maximum Likelihood and Bayesian inference. The Methods section presents the procedure and options available to conduct species delimitation analyses, including distance-, network-, and tree-based models. Finally, limits and future developments are discussed in the Notes section. Most importantly, species delimitation methods discussed herein are categorized based on five indicators: reliability, availability, scalability, understandability, and usability, all of which are fundamental properties needed for any approach to gain unanimous adoption within the DNA barcoding community moving forward.

Species Delimitation in a Polyploid Group of Iberian Jasione (Campanulaceae) Unveils Coherence between Cryptic Speciation and Biogeographical Regionalization

Groups with morphological stasis are an interesting framework to address putative cryptic species that may be hidden behind traditional taxonomic treatments, particularly when distribution ranges suggest disjunct and environmentally heterogeneous biogeographic patterns. New hypotheses of delimitation of evolutionary independent units can lead to the identification of different biogeographic processes, laying the foundation to investigate their historical and ecological significance. Jasione is a plant genus with a distribution centered in the Mediterranean basin, characterized by significant morphological stasis. Within the western Mediterranean J. gr. crispa species complex, J. sessiliflora s.l. and allied taxa form a distinct group, occupying environmentally diverse regions. At least two ploidy levels, diploid and tetraploid, are known to occur in the group. The internal variability is assessed with phylogenetic tools, viz. GMYC and ASAP, for species delimitation. The results are compared with other lines of evidence, including morphology and cytology. The fitting of distribution patterns of the inferred entities to chorological subprovinces is also used as a biogeographical and environmental framework to test the species hypothesis. Despite the scarcity of diagnostic morphological characters in the group, phylogenetic delimitation supports the description of at least one cryptic species, a narrow endemic in the NE Iberian Peninsula. Moreover, the results support the segregation of a thermophilic group of populations in eastern Iberia from J. sessiliflora. Ploidy variation from a wide geographical survey supports the systematic rearrangement suggested by species delimitation. Taxonomic reorganization in J. sessiliflora s.l. would allow ecological interpretations of distribution patterns in great accordance with biogeographical regionalization at the subprovince level, supporting geobotanical boundaries as a framework to interpret species ecological coherence of cryptic lineages. These results suggest that species differentiation, together with geographic isolation and polyploidization, is associated with adaptation to different environments, shifting from more to less thermophilic conditions. Thus, the recognition of concealed evolutionary entities is essential to correctly interpret biogeographical patterns in regions with a complex geologic and evolutionary history, such as the Mediterranean basin, and biogeographical units emerge as biologically sound frameworks to test the species hypothesis.

taxalogue: a toolkit to create comprehensive CO1 reference databases

Background

Taxonomic identification through DNA barcodes gained considerable traction through the invention of next-generation sequencing and DNA metabarcoding. Metabarcoding allows for the simultaneous identification of thousands of organisms from bulk samples with high taxonomic resolution. However, reliable identifications can only be achieved with comprehensive and curated reference databases. Therefore, custom reference databases are often created to meet the needs of specific research questions. Due to taxonomic inconsistencies, formatting issues, and technical difficulties, building a custom reference database requires tremendous effort. Here, we present taxalogue, an easy-to-use software for creating comprehensive and customized reference databases that provide clean and taxonomically harmonized records. In combination with extensive geographical filtering options, taxalogue opens up new possibilities for generating and testing evolutionary hypotheses.

Methods

taxalogue collects DNA sequences from several online sources and combines them into a reference database. Taxonomic incongruencies between the different data sources can be harmonized according to available taxonomies. Dereplication and various filtering options are available regarding sequence quality or metadata information. taxalogue is implemented in the open-source Ruby programming language, and the source code is available at https://github.com/nwnoll/taxalogue. We benchmark four reference databases by sequence identity against eight queries from different localities and trapping devices. Subsamples from each reference database were used to compare how well another one is covered.

Results

taxalogue produces reference databases with the best coverage at high identities for most tested queries, enabling more accurate, reliable predictions with higher certainty than the other benchmarked reference databases. Additionally, the performance of taxalogue is more consistent while providing good coverage for a variety of habitats, regions, and sampling methods. taxalogue simplifies the creation of reference databases and makes the process reproducible and transparent. Multiple available output formats for commonly used downstream applications facilitate the easy adoption of taxalogue in many different software pipelines. The resulting reference databases improve the taxonomic classification accuracy through high coverage of the query sequences at high identities.

DNA barcoding the ichthyofauna of the Beibu Gulf: Implications for fisheries management in a seafood market hub

The Beibu Gulf in China is situated in the tropics, in the western Pacific Ocean. It is an emblematic region combining proximity to a marine biodiversity hotspot and a major seafood hub. Intensification of marine fishing and ocean warming led to a drastic decline in fish populations in the Beibu Gulf during the last decades. This situation urges the development of molecular resources of the Beibu Gulf fish fauna in order to enable automated molecular identifications at the species level for next-generation monitoring. With this objective, we present the results of a large-scale campaign to DNA barcode fishes of the Beibu Gulf. We successfully generated 789 new DNA barcodes corresponding to 263 species which, together with 291 sequences mined from Genbank and BOLD, resulted in a reference library of 1080 sequences from 285 species. Based on the use of four DNA-based species delimitation methods (BIN, ASAP, mPTP, mGMYC), a total of 285 Molecular Operational Taxonomical Units (MOTUs). A single case of cryptic diversity was detected in Scomberomorus guttatus and a single species pair was not captured by delimitation methods. Intraspecific K2P genetic distances averaged 0.36% among sequences within species, whereas K2P genetic distances among species within genera averaged 6.96%. The most speciose families in open water trawling differ from those at fish market, and discrepancies with historical data are discussed in the light of recently documented stock collapses.

Insights into the divergent evolution of the oceanic squid Sthenoteuthis oualaniensis (Cephalopoda: Ommastrephidae) from the Indian Ocean

Sthenoteuthis oualaniensis is known for its complex population structure with three major transoceanic forms (viz. middle-sized, dwarf, and giant forms) whose taxonomic status has been disputed for decades. This integrated taxonomic study examines these prevenient morphotypes gathered on cruises in the Indian Ocean to ascertain their status in the evolutionary history of the species. Molecular analyses employing mitochondrial (COI, ND2) and nuclear (H3) markers revealed four genetically distinct and novel lineages of the species in the Indian Ocean, representing three morphotypes from the Arabian Sea and one from the Southern Indian Ocean. The mitochondrial-based phylograms revealed two distinct clades in the species: "dwarf forms + giant form" and "middle-sized forms," which further branch into geographically structured evolutionary units. Species delimitation analyses recovered five distinct clades, namely, the Arabian Sea giant and dwarf forms, Equatorial, Eastern Typical, and Other Middle-sized forms, representing the consensus molecular operational taxonomic units. H3 being heterozygous could not resolve the phylogeny. Haplotype network and AMOVA analysis of mtDNA genes indicated explicit phylogeographic structuring of haplotypes, whereas these outputs and PCA results were incongruent with the morphological grouping. Phenetic features distinguishing the morphotypes were sometimes plastic and mismatched with the genotypes. The giant form was genetically close to the dwarf forms, contradicting the earlier notion that it descended from the middle-sized form. It may be assumed that the dwarf form evolved following sympatric speciation and adaptation to warm equatorial waters, while the focal features of the Western Arabian Sea guide toward allopatric speciation of the giant form.

Penicillium rhizophilum, a novel species in the section Exilicaulis isolated from the rhizosphere of sugarcane in Southwest Iran

During a survey of species diversity of Penicillium and Talaromyces in sugarcane (Saccharum officinarum) rhizosphere in the Khuzestan province of Iran [1], 195 strains were examined, from which 187 belonged to Penicillium (11 species) and eight to Talaromyces (one species). In the present study, three strains of Penicillium belonging to section Exilicaulis series Restricta, identified as P. restrictum by Ansari et al. [1], were subjected to a phylogenetic study. The multilocus phylogeny of partial β-tubulin, calmodulin and RNA polymerase II second largest subunit genes enabled the recognition of one new phylogenetic species that is here formally described as Penicillium rhizophilum sp. nov. This species is phylogenetically distinct in series Restricta, but it does not show significant morphological differences from other species previously classified in the series. Therefore, we here placed bias on the phylogenetic species concept. The holotype of Penicillium rhizophilum sp. nov. is IRAN 18169F and the ex-type culture is LA30T (=IRAN 4042CT=CBS 149737T).

How we study cryptic species and their biological implications: A case study from marine shelled gastropods

Methodological and biological considerations are intertwined when studying cryptic species. A potentially large component of modern biodiversity, the frequency of cryptic species among taxonomic groups is not well documented. The term "cryptic species" is imprecisely used in scientific literature, causing ambiguity when interpreting their evolutionary and ecological significance. This study reviews how cryptic species have been defined, discussing implications for taxonomy and biology, and explores these implications with a case study based on recently published literature on extant shelled marine gastropods. Reviewed gastropods were recorded by species. Records of cryptic gastropods were presented by authors with variable levels of confidence but were difficult to disentangle from inherent biases in the study effort. These complexities notwithstanding, most gastropod species discussed were not cryptic. To the degree that this review's sample represents extinct taxa, the results suggest that a high proportion of shelled marine gastropod species are identifiable for study in the fossil record. Much additional work is needed to provide a more adequate understanding of the relative frequency of cryptic species in shelled marine gastropods, which should start with more explicit definitions and targeted case studies.

Molecular and morphological approaches redefine the limits among polymorphic species in the Neotropical longhorn beetle genus, Myzomorphus Sallé (Coleoptera: Cerambycidae: Prioninae)

Myzomorphus Sallé is a charismatic genus of prionine longhorn beetles (Cerambycidae) composed of nine species. Myzomorphus species are found from Costa Rica to southern Brazil, but only two species have wide distributions across this range: M. scutellatus Sallé from Costa Rica to northern Brazil, and M. quadripunctatus (Gray) from Colombia to southern Brazil. These species are highly polymorphic and their limits are difficult to determine due to their strong morphological similarities--males are only distinguishable by subtle size variations and females by color patterns. Here, we used mitochondrial DNA (cox1 and 12S) to reconstruct the first phylogeny of Myzomorphus and, in combination with morphological data, assess the taxonomic limits between M. scutellatus and M. quadripunctatus. Our phylogenetic results confirm the monophyly of Myzomorphus and reveal a close relationship among M. birai, M. quadripunctatus and M. scutellatus. Using pairwise distance estimations, we found that the intraspecific variation of M. quadripunctatus is remarkably high (K2P: 0-11.7%; p-distances: 0-9.7%) and the interspecific distances of M. quadripunctatus in relation to M. birai and M. scutellatus (K2P: 14.8-20.1%; p-distances: 12-15%) are close to the intraspecific distances of M. quadripunctatus. We further analyzed the diagnostic characters of these species and found that their morphological intraspecific variations largely overlap. Altogether, our results demonstrate that the variability of M. birai, M. scutellatus and M. quadripunctatus represent polymorphisms of a single species. We thus argue for the synonymy of M. birai and M. scutellatus under M. quadripunctatus (syn. nov.) and highlight the need for multiple lines of evidence to solve the taxonomic problems in polymorphic species of Cerambycidae.

Molecular Species Delimitation Using COI Barcodes of Mealybugs (Hemiptera: Pseudococcidae) from Coffee Plants in Espírito Santo, Brazil

Mealybugs are insects belonging to the family Pseudococcidae. This family includes many plant-pest species with similar morphologies, which may lead to errors in mealybug identification and delimitation. In the present study, we employed molecular-species-delimitation approaches based on distance (ASAP) and coalescence (GMYC and mPTP) methods to identify mealybugs collected from coffee and other plant hosts in the states of Espírito Santo, Bahia, Minas Gerais, and Pernambuco, Brazil. We obtained 171 new COI sequences, and 565 from the BOLD Systems database, representing 26 candidate species of Pseudococcidae. The MOTUs estimated were not congruent across different methods (ASAP-25; GMYC-30; mPTP-22). Misidentifications were revealed in the sequences from the BOLD Systems database involving Phenacoccus solani × Ph. solenopsis, Ph. tucumanus × Ph. baccharidis, and Planacoccus citri × Pl. minor species. Ten mealybug species were collected from coffee plants in Espírito Santo. Due to the incorrect labeling of the species sequences, the COI barcode library of the dataset from the database needs to be carefully analyzed to avoid the misidentification of species. The systematics and taxonomy of mealybugs may be improved by integrative taxonomy which may facilitate the integrated pest management of these pests.

DNA Barcode Library of Megadiverse Lepidoptera in an Alpine Nature Park (Italy) Reveals Unexpected Species Diversity

Species inventories are a prerequisite for biodiversity monitoring and conservation, particularly in protected areas. However, the possibilities of a standardized survey of species diversity using DNA barcoding have so far hardly been implemented, especially in species-rich groups. A first-time molecular-based and nearly complete inventory of the megadiverse insect order Lepidoptera in a protected area in the Alps (Cottian Alps, Italy) was intended to test the possibilities and reliability of DNA-based identifications. From voucher material collected between 2019 and 2022, we successfully sequenced 1213 morphospecies that grouped into 1204 BINs (barcode index numbers), whereas DNA barcoding failed for another 18 species. A total of 35 species shared a BIN with one or more taxa, but a majority of 19 species could still be discriminated by divergent sequences. A total of 12 morphospecies split into two BINs. These species and a further 22 taxa with unique BINs and barcode divergences >2% to the nearest neighbor require taxonomic re-assessment. Two additional cryptic species from the study area were described recently. Finally, 16 species are newly recorded for Italy. Our study, therefore, demonstrates the importance of DNA barcoding for both faunistics and the discovery of cryptic diversity, even in apparently well-studied protected areas.

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.

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).

From DNA barcodes to ecology: Meta-analysis of central European beetles reveal link with species ecology but also to data pattern and gaps

DNA barcoding has been used worldwide to identify biological specimens and to delimit species. It represents a cost-effective, fast, and efficient way to assess biodiversity with help of the public Barcode of Life Database (BOLD) accounting for more than 236,000 animal species and more than 10 million barcode sequences. Here, we performed a meta-analysis of available barcode data of central European Coleoptera to detect intraspecific genetic patterns among ecological groups in relation to geographic distance with the aim to investigate a possible link between infraspecific variation and species ecology. We collected information regarding feeding style, body size, as well as habitat and biotope preferences. Mantel tests and two variants of Procrustes analysis, both involving the Principal Coordinates Neighborhood Matrices (PCNM) approach, were applied on genetic and geographic distance matrices. However, significance levels were too low to further use the outcome for further trait investigation: these were in mean for all ecological guilds only 7.5, 9.4, or 15.6% for PCNM + PCA, NMDS + PCA, and Mantel test, respectively, or at best 28% for a single guild. Our study confirmed that certain ecological traits were associated with higher species diversity and foster stronger genetic differentiation. Results suggest that increased numbers of species, sampling localities, and specimens for a chosen area of interest may give new insights to explore barcode data and species ecology for the scope of conservation on a larger scale. We performed a meta-analysis of available barcode data of central European beetles to detect intraspecific genetic patterns among ecological groups in relation to geographic distance, regarding feeding style, body size, as well as habitat and biotope preferences. Our study confirmed that certain ecological traits were associated with higher species diversity and foster stronger genetic differentiation. However, significance levels were too low to further use the outcome for further trait investigation.

DNA Barcoding of Culicoides Latreille (Diptera: Ceratopogonidae) From Thailand Reveals Taxonomic Inconsistencies and Novel Diversity Among Reported Sequences

Recent focus on Culicoides species diversity in Thailand was prompted by a need to identify vectors responsible for the transmission of African Horse Sickness in that country. To assist rapid genetic identification of species, we sampled mitochondrial cytochrome c oxidase subunit I (COI) DNA barcodes (N = 78) from 40 species of Culicoides biting midge from Thailand, including 17 species for which DNA barcodes were previously unavailable. The DNA barcodes were assigned to 39 Barcode Identification Numbers (BINs) representing terminal genetic clusters at the Barcode of Life Data systems (BOLD). BINs assisted with comparisons to published conspecific DNA barcodes and allowed partial barcodes obtained from seven specimens to be associated with BINs by their similarity. Some taxonomic issues were revealed and attributed to the possible misidentification of earlier reported specimens as well as a potential synonymy of C. elbeli Wirth & Hubert and C. menglaensis Chu & Liu. Comparison with published BINs also revealed genetic evidence of divergent population processes and or potentially cryptic species in 16 described taxa, flagged by their high levels of COI sequence difference among conspecifics. We recommend the BOLD BIN system to researchers preparing DNA barcodes of vouchered species for public release. This will alert them to taxonomic incongruencies between their records and publicly released DNA barcodes, and also flag genetically deep and potentially novel diversity in described species.

A manager's guide to using eDNA metabarcoding in marine ecosystems

Environmental DNA (eDNA) metabarcoding is a powerful tool that can enhance marine ecosystem/biodiversity monitoring programs. Here we outline five important steps managers and researchers should consider when developing eDNA monitoring program: (1) select genes and primers to target taxa; (2) assemble or develop comprehensive barcode reference databases; (3) apply rigorous site occupancy based decontamination pipelines; (4) conduct pilot studies to define spatial and temporal variance of eDNA; and (5) archive samples, extracts, and raw sequence data. We demonstrate the importance of each of these considerations using a case study of eDNA metabarcoding in the Ports of Los Angeles and Long Beach. eDNA metabarcoding approaches detected 94.1% (16/17) of species observed in paired trawl surveys while identifying an additional 55 native fishes, providing more comprehensive biodiversity inventories. Rigorous benchmarking of eDNA metabarcoding results improved ecological interpretation and confidence in species detections while providing archived genetic resources for future analyses. Well designed and validated eDNA metabarcoding approaches are ideally suited for biomonitoring applications that rely on the detection of species, including mapping invasive species fronts and endangered species habitats as well as tracking range shifts in response to climate change. Incorporating these considerations will enhance the utility and efficacy of eDNA metabarcoding for routine biomonitoring applications.

Between allopatry and secondary contact: differentiation and hybridization among three sympatric Gentiana species in the Qinghai-Tibet Plateau

BACKGROUND

Mountains of the world host a significant portion of all terrestrial biodiversity, and the region of the Qinghai-Tibet Plateau (QTP) stands as one of the most remarkable mountain regions on Earth.  Because many explosive radiations occurred there, the QTP is a natural laboratory which is ideal to investigate patterns and processes linked to speciation and diversification. Indeed, understanding how closely related and sympatric species diverged is vital to explore drivers fostering speciation, a topic only rarely investigated in the QTP. By combining genomic and environmental data, we explored the speciation process among three closely related and sympatric species, Gentiana hexaphylla, G. lawrencei and G. veitchiorum in the QTP region.

RESULTS

Combining genome sizes and cytological data, our results showed that G. hexaphylla and G. veitchiorum are diploid, whereas G. lawrencei is tetraploid. Genetic clustering and phylogenetic reconstruction based on genomic SNPs indicated a clear divergence among the three species. Bayesian clustering, migrant, and D-statistic analyses all showed an obvious signature of hybridization among the three species, in particular between G. lawrencei and both G. hexaphylla and G. veitchiorum in almost all populations. Environmental variables related to precipitation and particularly temperature showed significant differences among the three gentians, and in fact a redundancy analysis confirmed that temperature and precipitation were the major climatic factors explaining the genetic differentiation among the three species.

CONCLUSION

Our study suggested that ancient hybridization, polyploidization, geological isolation and the evolution of different climatic preferences were all likely to be involved in the divergence of the three Gentiana species, as may be the case for many other taxa in the QTP region.

Luminescent characteristics and mitochondrial COI barcodes of nine cohabitated Taiwanese fireflies

Background

Over 50 Taiwanese firefly species have been discovered, but scientists lack information regarding most of their genetics, bioluminescent features, and cohabitating phenomena. In this study, we focus on morphological species identification and phylogeny reconstructed by COI barcoding, as well as luminescent characteristics of cohabited Taiwanese firefly species to determine the key factors that influenced how distinct bioluminescent species evolved to coexist and proliferate within the same habitat.

Methods

In this study, 366 specimens from nine species were collected in northern Taiwan from April to August, 2016-2019. First, the species and sex of the specimens were morphologically and genetically identified. Then, their luminescent spectra and intensities were recorded using a spectrometer and a power meter, respectively. The habitat temperature, relative humidity, and environmental light intensity were also measured. The cytochrome oxidase I (COI) gene sequence was used as a DNA barcode to reveal the phylogenetic relationships of cohabitated species.

Results

Nine species-eight adult species (Abscondita chinensis, Abscondita cerata, Aquatica ficta, Luciola curtithorax, Luciola kagiana, Luciola filiformis, Curtos sauteri, and Curtos costipennis) and one larval Pyrocoelia praetexta-were morphologically identified. The nine species could be found in April-August. Six of the eight adult species shared an overlap occurrence period in May. Luminescent spectra analysis revealed that the λ _max of studied species ranged from 552-572 nm (yellow-green to orange-yellow). The average luminescent intensity range of these species was about 1.2-14 lux (182.1-2,048 nW/cm²) for males and 0.8-5.8 lux (122.8-850 nW/cm²) for females, and the maximum luminescent intensity of males was 1.01-7.26-fold higher than that of females. Compared with previous studies, this study demonstrates that different λ _max, species-specific flash patterns, microhabitat choices, nocturnal activity time, and/or an isolated mating season are key factors that may lead to the species-specific courtship of cohabitated fireflies. Moreover, we estimated that the fireflies start flashing or flying when the environmental light intensity decreased to 6.49-28.1 lux. Thus, based on a rough theoretical calculation, the sensing distance between male and female fireflies might be 1.8-2.7 m apart in the dark. In addition, the mitochondrial COI barcode identified species with high resolution and suggested that most of the studied species have been placed correctly with congeners in previous phylogenies. Several cryptic species were revealed by the COI barcode with 3.27%-12.3% variation. This study renews the idea that fireflies' luminescence color originated from the green color of a Lampyridae ancestor, then red-shifted to yellow-green in Luciolinae, and further changed to orange-yellow color in some derived species.

Analysis of the Holarctic Dictyoptera aurora Complex (Coleoptera, Lycidae) Reveals Hidden Diversity and Geographic Structure in Müllerian Mimicry Ring

The elateroid family Lycidae is known for limited dispersal propensity and high species-level endemism. The red net-winged beetle, Dictyoptera aurora (Herbst, 1874), differs from all relatives by the range comprising almost the entire Holarctic region. Based on a five-marker phylogeny and 67 barcode entries (cox1-5' mtDNA) from the whole range, we recovered two genetically distinct species within traditionally defined D. aurora and resurrected the name D. coccinata (Say, 1835) as the oldest available synonym for Nearctic populations. Yet, no reliable morphological trait distinguishes these species except for minute differences in the male genitalia. D. coccinata is a monophylum resulting from a single Miocene dispersal event, ~15.8 million years ago, and genetic divergence implies long-term isolation by the Bering Strait. Far East Asian and west European populations are also genetically distinct, although to a lower extent. Two independent colonization events established the Fennoscandian populations after the last glacial maximum. Besides intrinsic factors, the high morphological similarity might result from stabilizing selection for shared aposematic signals. The rapidly accumulating barcode data provide valuable information on the evolutionary history and the origins of regional faunas.

DNA barcoding and species delimitation of the genus Oxynoemacheilus (Teleostei: Nemacheilidae) in Anatolia

The DNA barcoding approach was used for the determination of evolutionary relationships and species delimitation of the genus Oxynoemacheilus (Teleostei: Nemacheilidae). The COI barcode region (615 bp amplicon) was used to barcode 444 individuals from 64 morphologically identified species in the genus Oxynoemacheilus and 189 haplotypes were identified. The average of the interspecific p distance (9.59%) was about 21-fold higher than the average intraspecific distance (0.44%). A general genetic threshold of 1.46% sequence divergence was defined for species delimitation. The multiple species delimitation methods (BCM, GMYC, bPTP and TCS) revealed a total of 62 molecular operational taxonomic units for 64 morphospecies with a new loach species from the BuyukMelen River. Neighbour-joining, maximum likelihood and Bayesian inference analyses indicated that all haplotypes were clustered into 62 clades, which corresponded to Oxynoemacheilus species, with strong bootstrap support (≥95%). Furthermore, all samples grouped in concurrence with the taxonomic status of the species except for species groups (O. germencicus-O. cinicus-O. mesudae and O. leontinae-O. namiri) that were showed intraspecific overlap in genetic diversity for COI-based barcodes. In conclusion, our analyses indicate that COI-based barcodes provide reliable species discrimination. Therefore, we currently recommend COI barcodes as the suitable barcode for genus Oxynoemacheilus.

An integrative insight into the diversity, distribution, and biogeography of the freshwater endemic clade of the Ponticola syrman group (Teleostei: Gobiidae) in the Caucasus biodiversity hotspot

Freshwater habitats of the Caucasus biodiversity hotspot represent a center of endemism for the gobiid genus Ponticola Iljin, 1927. Hitherto, large-scale molecular studies, owing to restricted taxon and geographical sampling, have failed to give an elaborate picture of diversity and evolutionary history of these species. Here, to contribute to filling this gap, we assessed taxonomic diversity, phylogeography and evolutionary history for the south Caspian populations of Ponticola presently classified as P. iranicus and P. patimari, using an integrative taxonomic approach comprising an entire geographic range sampling, and analyses of mitochondrial DNA haplotypes, the head lateral line system, otolith shape, and meristic and morphometric variation. All freshwater samples of the P. syrman group belong to a monophyletic clade with two main subclades: a small subclade confined to the upper Sefidroud sub-basin including the type locality of P. iranicus and a large subclade with three geographically constrained haplogroups (Hg1, Hg2, and Hg3), comprising the rest of the distribution. Hg1 showed an eastern distribution including the type locality of P. patimari, while Hg2 and Hg3 are sister groups with central and western-central distributions, respectively. The freshwater clade diverged from P. syrman during the Tyurkyanian low stand (~150 m b.s.l. lasting ~0.1 Myr), while the divergence of P. iranicus and P. patimari and radiations within P. patimari took place during the Bakunian high stand (up to 50 m a.s.l. lasting ~378-480 kya). Species delimitation analyses indicated two distinct species, corresponding to each main subclade. Although the otolith shape and lateral line analyses did not reflect with phylogeographic pattern, PCA and DFA plots of meristic and morphometric data showed a clear separation of the two major subclades corresponding to P. iranicus and P. patimari, suggesting the presence of significant morphological variation meriting formal taxonomic recognition. Overall, our findings (i) reveal the presence of two freshwater endemic species in the P. syrman group, and pending further investigation, hypothesize the presence of a third cryptic species; (ii) revise and document a narrow distributional range and low diversity for P. iranicus, in contrast to a wider distributional range and high diversity for P. patimari; (iii) suggest that the climatic oscillations of the Pleistocene were associated with the cladogenesis within the P. syrman group; and (iv) allowed for the recognition of conservation units and proposition of management measures.

DNA barcoding of a lesser-known catfish, Clupisoma bastari (Actinopterygii: Ailiidae) from Deccan Peninsula, India

DNA barcoding substantiates species identification, and simultaneously indicates the misnomer taxa. Based on the morphological descriptions, we identified a lesser-known catfish, Clupisoma bastari, from Godavari River basin, and contributed novel DNA barcode data to the GenBank. The Kimura 2 parameter genetic divergence between species, and the neighbour-joining phylogeny clearly depicted a distinct clade of C. bastari in the studied dataset. Clupisoma bastari maintained sufficient K2P genetic divergence (8.3% to 11.2%) with other congeners, and branched as a sister-species of C. garua. The present study highlights possible existence of a few misnomer taxa in the GenBank. We encourage further extensive sampling of different congeners of Clupisoma from a wide range of habitats to explore the species diversity and phylogenetic relationship. 

Butterfly–parasitoid–hostplant interactions in Western Palaearctic Hesperiidae: a DNA barcoding reference library

The study of ecological interactions between plants, phytophagous insects and their natural enemies is an essential but challenging component for understanding ecosystem dynamics. Molecular methods such as DNA barcoding can help elucidate these interactions. In this study, we employed DNA barcoding to establish hostplant and parasitoid interactions with hesperiid butterflies, using a complete reference library for Hesperiidae of continental Europe and north-western Africa (53 species, 100% of those recorded) based on 2934 sequences from 38 countries. A total of 233 hostplant and parasitoid interactions are presented, some recovered by DNA barcoding larval remains or parasitoid cocoons. Combining DNA barcode results with other lines of evidence allowed 94% species-level identification for Hesperiidae, but success was lower for parasitoids, in part due to unresolved taxonomy. Potential cases of cryptic diversity, both in Hesperiidae and Microgastrinae, are discussed. We briefly analyse the resulting interaction networks. Future DNA barcoding initiatives in this region should focus attention on north-western Africa and on parasitoids, because in these cases barcode reference libraries and taxonomy are less well developed.

Large scale DNA barcoding of the subfamily Culterinae (Cypriniformes: Xenocyprididae) in East Asia unveils geographic scale effect, taxonomic warnings and cryptic diversity

Geographical scale might be expected to impact significantly the efficiency of DNA barcoding as spatially comprehensive sampling provides opportunities to uncover intricate relationships among closely related species and to detect cryptic diversity for widespread taxa. Here, we present a DNA barcoding study on a Xencyprididae subfamily (Culterinae) involving the production of 998 newly generated DNA barcodes from East Asian drainages (80 localities). Together with 513 barcodes mined from BOLD and GenBank, a reference library consisting of 1511 DNA barcodes (116 localities) for 42 species was assembled, accounting for 66% of known Culterinae species. Intraspecific genetic distances are positively correlated to geographical scale, while a negative correlation is detected between interspecific genetic distances and geographical scale. The present study demonstrates that geographical scale influences the efficiency of DNA barcoding by narrowing the width of the barcoding gap. DNA‐based species delimitation analyses delimited 44 molecular operational taxonomic units (MOTUs). Rampant cryptic diversity is detected within eight species with multiple MOTUs, whereas 25 species present mismatch between morphological and molecular delimitations. A total of 18 species are lumped into nine MOTUs due to low interspecific divergence and/or mixed lineages. Several MOTU divergences are hypothesized to relate to known biogeographical barriers and geological events during the Pliocene and Pleistocene. This study provides new insights into the taxonomy and phylogeography of the subfamily Culterinae.

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.

Towards Large-scale Integrative Taxonomy (LIT): resolving the data conundrum for dark taxa

New, rapid, accurate, scalable, and cost-effective species discovery and delimitation methods are needed for tackling “dark taxa,” here defined as groups for which \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$<$\end{document}10\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$\%$\end{document} of all species are described and the estimated diversity exceeds 1,000 species. Species delimitation for these taxa should be based on multiple data sources (“integrative taxonomy”) but collecting multiple types of data risks impeding a discovery process that is already too slow. We here develop large-scale integrative taxonomy (LIT), an explicit method where preliminary species hypotheses are generated based on inexpensive data that can be obtained quickly and cost-effectively. These hypotheses are then evaluated based on a more expensive type of “validation data” that is only obtained for specimens selected based on objective criteria applied to the preliminary species hypotheses. We here use this approach to sort 18,000 scuttle flies (Diptera: Phoridae) into 315 preliminary species hypotheses based on next-generation sequencing barcode (313 bp) clusters (using objective clustering [OC] with a 3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$\%$\end{document} threshold). These clusters are then evaluated with morphology as the validation data. We develop quantitative indicators for predicting which barcode clusters are likely to be incongruent with morphospecies by randomly selecting 100 clusters for in-depth validation with morphology. A linear model demonstrates that the best predictors for incongruence between barcode clusters and morphology are maximum p-distance within the cluster and a newly proposed index that measures cluster stability across different clustering thresholds. A test of these indicators using the 215 remaining clusters reveals that these predictors correctly identify all clusters that are incongruent with morphology. In our study, all morphospecies are true or disjoint subsets of the initial barcode clusters so that all incongruence can be eliminated by varying clustering thresholds. This leads to a discussion of when a third data source is needed to resolve incongruent grouping statements. The morphological validation step in our study involved 1,039 specimens (5.8\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$\%$\end{document} of the total). The formal LIT protocol we propose would only have required the study of 915 (5.1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$\%$\end{document}: 2.5 specimens per species), as we show that clusters without signatures of incongruence can be validated by only studying two specimens representing the most divergent haplotypes. To test the generality of our results across different barcode clustering techniques, we establish that the levels of incongruence are similar across OC, Automatic Barcode Gap Discovery (ABGD), Poisson Tree Processes (PTP), and Refined Single Linkage (RESL) (used by Barcode of Life Data System to assign Barcode Index Numbers [BINs]). OC and ABGD achieved a maximum congruence score with the morphology of 89\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$\%$\end{document} while PTP was slightly less effective (84\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$\%$\end{document}). RESL could only be tested for a subset of the specimens because the algorithm is not public. BINs based on 277 of the original 1,714 haplotypes were 86\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$\%$\end{document} congruent with morphology while the values were 89\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$\%$\end{document} for OC, 74\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$\%$\end{document} for PTP, and 72\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$\%$\end{document} for ABGD. [Biodiversity discovery; dark taxa; DNA barcodes; integrative taxonomy.]

The Probability of Joint Monophyly of Samples of Gene Lineages for All Species in an Arbitrary Species Tree

Monophyly is a feature of a set of genetic lineages in which every lineage in the set is more closely related to all other members of the set than it is to any lineage outside the set. Multiple sets of lineages that are separately monophyletic are said to be reciprocally monophyletic, or jointly monophyletic. The prevalence of reciprocal monophyly, or joint monophyly (JM), has been used to evaluate phylogenetic and phylogeographic hypotheses, as well as to delimit species. These applications often make use of a probability of JM under models of gene lineage evolution. Studies in coalescent theory have computed this JM probability for small numbers of separate groups in arbitrary species trees and for arbitrary numbers of separate groups in trivial species trees. In this study, generalizing existing results on monophyly probabilities under the multispecies coalescent, we derive the probability of JM for arbitrary numbers of separate groups in arbitrary species trees. We illustrate how our result collapses to previously examined cases. We also study the effect of tree height, sample size, and number of species on the probability of JM. We obtain relatively simple lower and upper bounds on the JM probability. Our results expand the scope of JM calculations beyond small numbers of species, subsuming past formulas that have been used in simpler cases.

Exploring Large-Scale Patterns of Genetic Variation in the COI Gene among Insecta: Implications for DNA Barcoding and Threshold-Based Species Delimitation Studies

The genetic variation in the COI gene has had a great effect on the final results of species delimitation studies. However, little research has comprehensively investigated the genetic divergence in COI among Insecta. The fast-growing COI data in BOLD provide an opportunity for the comprehensive appraisal of the genetic variation in COI among Insecta. We calculated the K2P distance of 64,414 insect species downloaded from BOLD. The match ratios of the clustering analysis, based on different thresholds, were also compared among 4288 genera (35,068 species). The results indicate that approximately one-quarter of the species of Insecta showed high intraspecific genetic variation (>3%), and a conservative estimate of this proportion ranges from 12.05% to 22.58%. The application of empirical thresholds (e.g., 2% and 3%) in the clustering analysis may result in the overestimation of the species diversity. If the minimum interspecific genetic distance of the congeneric species is greater than or equal to 2%, it is possible to avoid overestimating the species diversity on the basis of the empirical thresholds. In comparison to the fixed thresholds, the “threshOpt” and “localMinima” algorithms are recommended for the provision of a reference threshold for threshold-based species delimitation studies.

Lack of Statistical Rigor in DNA Barcoding Likely Invalidates the Presence of a True Species' Barcode Gap

DNA barcoding has been largely successful in satisfactorily exposing levels of standing genetic diversity for a wide range of taxonomic groups through the employment of only one or a few universal gene markers. However, sufficient coverage of geographically-broad intra-specific haplotype variation within genomic databases like the Barcode of Life Data Systems (BOLD) and GenBank remains relatively sparse. As reference sequence libraries continue to grow exponentially in size, there is now the need to identify novel ways of meaningfully analyzing vast amounts of available DNA barcode data. This is an important issue to address promptly for the routine tasks of specimen identification and species discovery, which have seen broad adoption in areas as diverse as regulatory forensics and resource conservation. Here, it is demonstrated that the interpretation of DNA barcoding data is lacking in statistical rigor. To highlight this, focus is set specifically on one key concept that has become a household name in the field: the DNA barcode gap. Arguments outlined herein specifically center on DNA barcoding in animal taxa and stem from three angles: (1) the improper allocation of specimen sampling effort necessary to capture adequate levels of within-species genetic variation, (2) failing to properly visualize intra-specific and interspecific genetic distances, and (3) the inconsistent, inappropriate use, or absence of statistical inferential procedures in DNA barcoding gap analyses. Furthermore, simple statistical solutions are outlined which can greatly propel the use of DNA barcoding as a tool to irrefutably match unknowns to knowns on the basis of the barcoding gap with a high degree of confidence. Proposed methods examined herein are illustrated through application to DNA barcode sequence data from Canadian Pacific fish species as a case study.

SPEDE‐sampler: an R Shiny application to assess how methodological choices and taxon‐sampling can affect Generalised Mixed Yule Coalescent (GMYC) output and interpretation

Species delimitation tools are vital to taxonomy and the discovery of new species. These tools can make use of genetic data to estimate species boundaries, where one of the most widely used methods is the Generalized Mixed Yule Coalescent (GMYC) model. Despite its popularity, a number of factors are known to influence the performance and resulting inferences of the GMYC. Moreover, the few studies that have assessed model performance to date have been predominantly based on simulated data sets, where model assumptions are not violated. Here, we present a user‐friendly R Shiny application, ‘SPEDE‐sampler’ (SPEcies DElimitation sampler), that assesses the effect of computational and methodological choices, in combination with sampling effects, on the GMYC model. Output phylogenies are used to test the effect that (1) sample size, (2) BEAST and GMYC parameters (e.g. prior settings, single vs multiple threshold, clock model), and (3) singletons have on GMYC output. Optional predefined grouping information (e.g. morphospecies/ecotypes) can be uploaded in order to compare it with GMYC species and estimate percentage match scores. Additionally, predefined groups that contribute to inflated species richness estimates are identified by SPEDE‐sampler, allowing for the further investigation of potential cryptic species or geographical substructuring in those groups. Merging by the GMYC is also recorded to identify where traditional taxonomy has overestimated species numbers. Four worked examples are provided to illustrate the functionality of the program's workflow, and the variation that can arise when applying the GMYC model to empirical data sets. The R Shiny program is available for download at https://github.com/clarkevansteenderen/spede_sampler_R.