The effect of geographical scale of sampling on DNA barcoding

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.

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.

Improving the DNA Barcode Library of Mosquito Species With New Identifications and Discoveries in North-Central Argentina

The family Culicidae is represented by 244 species in Argentina, many of them with epidemiological importance. DNA barcodes are effective tools for identifying mosquito species, for knowing genetic variability, and for establishing phylogenetic relationships. This work aims to explore mosquito diversity employing different species delimitation approaches and to establish formally a DNA barcode library for the Argentinian mosquito fauna. Barcode fragments of 80 specimens of Argentinian mosquitoes of 28 species of the genera Aedeomyia Theobald (Diptera: Culicidae), Anopheles Meigen (Diptera: Culicidae), Coquillettidia Dyar (Diptera: Culicidae), Culex L. (Diptera: Culicidae), Haemagogus Williston (Diptera: Culicidae), Mansonia Blanchard (Diptera: Culicidae), Nyssorhynchus Blanchard (Diptera: Culicidae), Ochlerotatus Lynch-Arribálzaga (Diptera: Culicidae), Psorophora Robinneau-Desvoidy (Diptera: Culicidae) and Uranotaenia Lynch-Arribálzaga (Diptera: Culicidae) were sequenced. Another 82 sequences were obtained from public databases to establish the phylogenetic relationships using Maximum Likelihood and Bayesian Inference, and the species boundaries based on three approaches (ABGD, GMYC, and mPTP). Sixteen of the 28 species sequenced were recovered as monophyletic, of which 12 were also recognized as molecular operational taxonomic units according to the three methodologies. The disparity between morphology and barcode-based identifications could be explained by synonymy, species complexes occurrence, hybridization, incomplete lineage sorting, or the effect of the geographical scale of sampling. Twenty of the 28 sequenced species are new barcodes for Argentina and 11 are the first for science. This increases from 31 to 52 (12.7 to 21.31%) and from six to 10 (28.57 to 47.62%) the number of species and genera, respectively, with barcode sequences in Argentina. New species records are provided.

Species Delimitation of Scavenger Flies in the Valley of Mexico

Identification of species involved in cadaveric decomposition, such as scavenger Diptera, is a fundamental step for the use of entomological evidence in court. Identification based on morphology is widely used in forensic cases; however, taxonomic knowledge of scavenger fauna is poor for many groups and for many countries, particularly Neotropical ones. A number of studies have documented the utility of a DNA barcoding strategy to assist in the identification of poorly known and diverse groups, particularly in cases involving immature states or fragmented organisms. To provide baseline knowledge of the diversity of scavenger Diptera in the Valley of Mexico, we generated a DNA barcode collection comprised of sequences of the cytochrome c oxidase subunit 1 (COI) gene for all families sampled at a nature reserve located in this region. We collected and identified specimens on the basis of morphology and a species delimitation analysis. Our analyses of 339 individuals delineated 42 species distributed across nine families of Diptera. The richest families were Calliphoridae (9 species), Sarcophagidae (7 species), and Phoridae (6 species). We found many of the species previously recorded for the Valley of Mexico, plus 18 new records for the region. Our study highlights the utility of DNA barcoding as a first-step strategy to assess species richness of poorly studied scavenger fly taxa.

Ant crickets and their secrets: Myrmecophilus acervorum is not always parthenogenetic (Insecta: Orthoptera: Myrmecophilidae)

Previously considered as a thelytokous parthenogenetic species, the widespread ant cricket Myrmecophilus acervorum actually turns out to have a mixed reproductive system: our recent surveys in the central part of its distribution area has revealed the presence of both sexes. Detailed morphological and morphometric descriptions of the previously unknown males are here provided. New data on species distribution in south-eastern Europe are presented, including the first records of M. balcanicus in Bulgaria and of M. nonveilleri in Bulgaria and Hungary. Phylogenetic and phylogeographic analyses have revealed several haplotypes of M. acervorum in Europe, with six of them forming a parthenogenetic clade in populations distributed west of the Carpathians. We tested our samples for bacterial infection by Wolbachia and, surprisingly, Wolbachia was identified only in populations with both sexes and no amplification was obtained from parthenogenetic populations. Phylogenetic analyses performed with sequences pertaining to five nominal species related to M. acervorum, yielded topological congruent trees with four well-supported groups: one group with M. acervorum samples, the second group with M. nonveilleri samples, the third group with M. fuscus and M. gallicus samples, and the fourth group with samples of M. balcanicus. We performed species delineation tests on our sequences, which delimited between four to seven putative species.

Genetic Diversity and Population Structure of the Asian Tiger Mosquito (Aedes albopictus) in Vietnam: Evidence for Genetic Differentiation by Climate Region

Aedes albopictus is a native mosquito to Southeast Asia with a high potential for disease transmission. Understanding how Ae. albopictus populations that develop in the species' native range is useful for planning future control strategies and for identifying the sources of invasive ranges. The present study aims to investigate the genetic diversity and population structure of Ae. albopictus across various climatic regions of Vietnam. We analyzed mitochondrial cytochrome oxidase I (COI) gene sequences from specimens collected from 16 localities, and we used distance-based redundancy analysis to evaluate the amount of variation in the genetic distance that could be explained by both geographic distance and climatic factors. High levels of genetic polymorphism were detected, and the haplotypes were similar to those sequences from both temperate and tropical regions worldwide. Of note, these haplotype groups were geographically distributed, resulting in a distinct population structure in which northeastern populations and the remaining populations were genetically differentiated. Notably, genetic variation among the Ae. albopictus populations was driven primarily by climatic factors (64.55%) and to a lesser extent was also influenced by geographic distance (33.73%). These findings fill important gaps in the current understanding of the population genetics of Ae. albopictus in Vietnam, especially with respect to providing data to track the origin of the invaded regions worldwide.

DROP: Molecular voucher database for identification of Drosophila parasitoids

Molecular identification is increasingly used to speed up biodiversity surveys and laboratory experiments. However, many groups of organisms cannot be reliably identified using standard databases such as GenBank or BOLD due to lack of sequenced voucher specimens identified by experts. Sometimes a large number of sequences are available, but with too many errors to allow identification. Here, we address this problem for parasitoids of Drosophila by introducing a curated open-access molecular reference database, DROP (Drosophila parasitoids). Identifying Drosophila parasitoids is challenging and poses a major impediment to realize the full potential of this model system in studies ranging from molecular mechanisms to food webs, and in biological control of Drosophila suzukii. In DROP, genetic data are linked to voucher specimens and, where possible, the voucher specimens are identified by taxonomists and vetted through direct comparison with primary type material. To initiate DROP, we curated 154 laboratory strains, 856 vouchers, 554 DNA sequences, 16 genomes, 14 transcriptomes, and six proteomes drawn from a total of 183 operational taxonomic units (OTUs): 114 described Drosophila parasitoid species and 69 provisional species. We found species richness of Drosophila parasitoids to be heavily underestimated and provide an updated taxonomic catalogue for the community. DROP offers accurate molecular identification and improves cross-referencing between individual studies that we hope will catalyse research on this diverse and fascinating model system. Our effort should also serve as an example for researchers facing similar molecular identification problems in other groups of organisms.

Meta-Fish-Lib: A generalised, dynamic DNA reference library pipeline for metabarcoding of fishes

The accuracy and reliability of DNA metabarcoding analyses depend on the breadth and quality of the reference libraries that underpin them. However, there are limited options available to obtain and curate the huge volumes of sequence data that are available on public repositories such as NCBI and BOLD. Here, we provide a pipeline to download, clean and annotate mitochondrial DNA sequence data for a given list of fish species. Features of this pipeline include (a) support for multiple metabarcode markers; (b) searches on species synonyms and taxonomic name validation; (c) phylogeny assisted quality control for identification and removal of misannotated sequences; (d) automatically generated coverage reports for each new GenBank release update; and (e) citable, versioned DOIs. As an example we provide a ready-to-use curated reference library for the marine and freshwater fishes of the U.K. To augment this reference library for environmental DNA metabarcoding specifically, we generated 241 new MiFish-12S sequences for 88 U.K. marine species, and make available new primer sets useful for sequencing these. This brings the coverage of common U.K. species for the MiFish-12S fragment to 93%, opening new avenues for scaling up fish metabarcoding across wide spatial gradients. The Meta-Fish-Lib reference library and pipeline is hosted at https://github.com/genner-lab/meta-fish-lib.

A molecular-based identification resource for the arthropods of Finland

To associate specimens identified by molecular characters to other biological knowledge, we need reference sequences annotated by Linnaean taxonomy. In this study, we (1) report the creation of a comprehensive reference library of DNA barcodes for the arthropods of an entire country (Finland), (2) publish this library, and (3) deliver a new identification tool for insects and spiders, as based on this resource. The reference library contains mtDNA COI barcodes for 11,275 (43%) of 26,437 arthropod species known from Finland, including 10,811 (45%) of 23,956 insect species. To quantify the improvement in identification accuracy enabled by the current reference library, we ran 1000 Finnish insect and spider species through the Barcode of Life Data system (BOLD) identification engine. Of these, 91% were correctly assigned to a unique species when compared to the new reference library alone, 85% were correctly identified when compared to BOLD with the new material included, and 75% with the new material excluded. To capitalize on this resource, we used the new reference material to train a probabilistic taxonomic assignment tool, FinPROTAX, scoring high success. For the full-length barcode region, the accuracy of taxonomic assignments at the level of classes, orders, families, subfamilies, tribes, genera, and species reached 99.9%, 99.9%, 99.8%, 99.7%, 99.4%, 96.8%, and 88.5%, respectively. The FinBOL arthropod reference library and FinPROTAX are available through the Finnish Biodiversity Information Facility (www.laji.fi) at https://laji.fi/en/theme/protax. Overall, the FinBOL investment represents a massive capacity-transfer from the taxonomic community of Finland to all sectors of society.

Molecular phylogeny of the Anopheles hyrcanus group (Diptera: Culicidae) based on rDNA-ITS2 and mtDNA-COII

BACKGROUND

The Anopheles hyrcanus group, which includes 25 species, is widely distributed in the Oriental and Palaearctic regions. Given the difficulty in identifying cryptic or sibling species based on their morphological characteristics, molecular identification is regarded as an important complementary approach to traditional morphological taxonomy. The aim of this study was to reconstruct the phylogeny of the Hyrcanus group using DNA barcoding markers in order to determine the phylogenetic correlations of closely related taxa and to compare these markers in terms of identification efficiency and genetic divergence among species.

METHODS

Based on data extracted from the GenBank database and data from the present study, we used 399 rDNA-ITS2 sequences of 19 species and 392 mtDNA-COII sequences of 14 species to reconstruct the molecular phylogeny of the Hyrcanus group across its worldwide range. We also compared the performance of rDNA-ITS2 against that of mtDNA-COII to assess the genetic divergence of closely related species within the Hyrcanus group.

RESULTS

Average interspecific divergence for the rDNA-ITS2 sequence (0.376) was 125-fold higher than the average intraspecies divergence (0.003), and average interspecific divergence for the mtDNA-COII sequence (0.055) was eightfold higher than the average intraspecies divergence (0.007). The barcoding gap ranged from 0.015 to 0.073 for rDNA-ITS2, and from 0.017 to 0.025 for mtDNA-COII. Two sets of closely related species, namely, Anophels lesteri and An. paraliae, and An. sinensis, An. belenrae and An. kleini, were resolved by rDNA-ITS2. In contrast, the relationship of An. sinensis/An. belenrae/An. kleini was poorly defined in the COII tree. The neutrality test and mismatch distribution revealed that An. peditaeniatus, An. hyrcanus, An. sinensis and An. lesteri were likely to undergo hitchhiking or population expansion in accordance with both markers. In addition, the population of an important vivax malaria vector, An. sinensis, has experienced an expansion after a bottleneck in northern and southern Laos.

CONCLUSIONS

The topology of the Hyrcanus group rDNA-ITS2 and mtDNA-COII trees conformed to the morphology-based taxonomy for species classification rather than for that for subgroup division. rDNA-ITS2 is considered to be a more reliable diagnostic tool than mtDNA-COII in terms of investigating the phylogenetic correlation between closely related mosquito species in the Hyrcanus group. Moreover, the population expansion of an important vivax malaria vector, An. sinensis, has underlined a potential risk of malaria transmission in northern and southern Laos. This study contributes to the molecular identification of the Anopheles hyrcanus group in vector surveillance.