Species diversity can be overestimated by a fixed empirical threshold: insights from DNA barcoding of the genus Cletus (Hemiptera: Coreidae) and the meta-analysis of COI data from previous phylogeographical studies

Genetic and morphometric differentiation between two morphs of Haematobosca sanguinolenta (Diptera: Muscidae) from Thailand

Haematobosca is a genus of biting fly within the subfamily Stomoxyinae of the family Muscidae. It is currently recognized to include 16 species worldwide. These species, acting as ectoparasites, are considered to have significant importance in the veterinary and medical fields. To address the color polymorphism related to the genus Haematobosca in Thailand, herein, we focused on the normal (legs mainly black) and yellow (legs mainly yellow) morphs of Haematobosca sanguinolenta and examined them for genetic differences using three molecular markers: the cytochrome c oxidase subunit 1 (cox1) and cytochrome b (cytb) genes from the mitochondrial genome as well as the internal transcribed spacer 2 (ITS2) region from the nuclear ribosomal DNA. In addition, we analyzed wing differences between the two morphs using geometric morphometrics (GM). The genetic divergences between the two morphs showed that cytb gene showed the greatest divergence, for which the average distance was 5.6%. This was followed by the combination of cox1-cytb-ITS2, exhibiting an average divergence of 4.5%, ITS2 with a divergence of 4.1%, and finally cox1, showing the lowest divergence of 3.5%. Phylogenetic analyses distinctly separated the two morphs of H. sanguinolenta; this separation was supported by high bootstrap values (97-100%). These results were further corroborated by three species delimitation methods, i.e. assemble species by automatic partitioning (ASAP), automated barcode gap discovery (ABGD), and Poisson tree processes (PTP), all of which suggested that the two morphs likely represent separate species. In addition, a GM study identified a statistically significant difference in wing shape between the two morphs of H. sanguinolenta (P < 0.05). This combination of genetic and morphometric results strongly supports the existence of two distinct species within H. sanguinolenta in Thailand.

Comprehensive DNA barcode reference library and optimization of genetic divergence threshold facilitate the exploration of species diversity of green lacewings (Neuroptera: Chrysopidae)

Chrysopidae are a family of Neuroptera of significant importance in biocontrol against agricultural pests because of their predatory larvae. Currently, the taxonomy of Chrysopidae lacks a comprehensive revision, which impedes the exploration of species diversity as well as the selection and the conservation of green lacewings as biocontrol agents. We have established a DNA barcode reference library of the Chinese green lacewings based on an approximately complete sampling (95.63%) in 25 of the 34 provincial regions in China, comprising 1 119 barcodes of 25 genera and 197 species (representing 85% genera and 43.62% species from China). Combining other 1 049 high quality green lacewing DNA barcodes, we first inferred the optimal threshold of interspecific genetic divergence (1.87%) for successful species identification in multiple simulated scenarios based on present data. We further inferred the threshold of genetic divergence (7.77%) among genera with biocontrol significance. The inference and performance of the threshold appears to be mainly associated with the completeness of sampling, the proportion of closely related species, and the analytical approaches. Six new combinations, Apertochrysa platypa (Yang & Yang, 1991) comb. nov., Apertochrysa shennongana (Yang & Wang, 1990) comb. nov., Apertochrysa pictifacialis (Yang, 1988) comb. nov., Apertochrysa helana (Yang, 1993) comb. nov., Plesiochrysa rosulata (Yang & Yang, 2002) comb. nov., and Signochrysa hainana (Yang & Yang, 1991), are proposed according to integrative species delimitation. Our library and optimal threshold will effectively facilitate the exploration of species diversity of green lacewings. Our study also provides a methodological reference in molecular delimitation of other insects.

Mitochondrial introgression and mito-nuclear discordance obscured the closely related species boundaries in Cletus Stål from China (Heteroptera: Coreidae)

Accurate taxonomy and delimitation are of great importance for pest control strategies and management programs. Here, we focus on Cletus (Insecta: Hemiptera: Coreidae), which includes many crop pests. The species boundaries still conflict and only cytochrome c oxidase subunit I (COI) barcoding has been previously used for molecular studies. We generated new mitochondrial genome and nuclear genome-wide SNPs to explore the species boundaries of 46 Cletus samples from China using multiple species delimitation approaches. All results recovered a monophyly with high support, except for two closely related species in clade I - C. punctiger and C. graminis. Mitochondrial data demonstrated admixture in clade I, while genome-wide SNPs unambiguously identified two separate species, which were confirmed by morphological classification. Inconsistent nuclear and mitochondrial data indicated mito-nuclear discordance. Mitochondrial introgression is the most likely explanation, and more extensive sampling and more comprehensive data are needed to ascertain a pattern. Accurate species delimitation will shed light on species status; thus, an accurate taxonomy is of particular concern, as there is a pressing need to implement precise control of agricultural pests and to perform further research on diversification.

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.

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

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

Integrative species delimitation based on COI, ITS, and morphological evidence illustrates a unique evolutionary history of the genus Paracercion (Odonata: Coenagrionidae)

Paracercion are common ‘blue and black’ colored damselflies. We explore the species boundaries of Paracercion (Odonata: Coenagrionidae) using ABGD, bPTP, GMYC and Distance-based clustering. We finally got the molecular data of all nine species of Paracercion. P. hieroglyphicum and P. melanotum were combined into one putative species based on cytochrome c oxidase I (COI). However, they were separated into two putative species based on the nuclear segment including ITS1-5.8S-ITS2 (ITS). This suggests the introgression of mtDNA in Paracercion. Paracercion barbatum and Paracercion melanotum can be separated into two species based on COI, whereas they were combined into one putative species based on ITS, which suggests a hybridization event between them. The lower interspecific divergence (COI: 0.49%) between P. barbatum and Paracercion v-nigrum indicates a recent speciation event in Paracercion. Paracercion sieboldii and P. v-nigrum can be separated into two putative species based on COI, while they were frequently merged into the same putative species based on ITS. This can be explained by incomplete lineage sorting in nDNA. Besides, P. pendulum and P. malayanum were synonymized as junior synonyms of P. melanotum. P. luzonicum was confirmed not to belong to Paracercion. The possibility of introgression, hybridization, recent speciation and incomplete lineage sorting makes species delimitation, based on molecular data, difficult and complicates understanding of the evolutionary history of Paracercion. The discordance in COI and ITS also indicates the value of using markers from different sources in species delimitation studies.

Identifying error and accurately interpreting environmental DNA metabarcoding results: A case study to detect vertebrates at arid zone waterholes

Environmental DNA (eDNA) metabarcoding surveys enable rapid, noninvasive identification of taxa from trace samples with wide-ranging applications from characterizing local biodiversity to identifying food-web interactions. However, the technique is prone to error from two major sources: (a) contamination through foreign DNA entering the workflow, and (b) misidentification of DNA within the workflow. Both types of error have the potential to obscure true taxon presence or to increase taxonomic richness by incorrectly identifying taxa as present at sample sites, but multiple error sources can remain unaccounted for in metabarcoding studies. Here, we use data from an eDNA metabarcoding study designed to detect vertebrate species at waterholes in Australia's arid zone to illustrate where and how in the workflow errors can arise, and how to mitigate those errors. We detected the DNA of 36 taxa spanning 34 families, 19 orders and five vertebrate classes in water samples from waterholes, demonstrating the potential for eDNA metabarcoding surveys to provide rapid, noninvasive detection in remote locations, and to widely sample taxonomic diversity from aquatic through to terrestrial taxa. However, we initially identified 152 taxa in the samples, meaning there were many false positive detections. We identified the sources of these errors, allowing us to design a stepwise process to detect and remove error, and provide a template to minimize similar errors that are likely to arise in other metabarcoding studies. Our findings suggest eDNA metabarcoding surveys need to be carefully conducted and screened for errors to ensure their accuracy.

An integrative DNA barcoding framework of ladybird beetles (Coleoptera: Coccinellidae)

Even though ladybirds are well known as economically important biological control agents, an integrative framework of DNA barcoding research was not available for the family so far. We designed and present a set of efficient mini-barcoding primers to recover full DNA barcoding sequences for Coccinellidae, even for specimens collected 40 years ago. Based on these mini-barcoding primers, we obtained 104 full DNA barcode sequences for 104 species of Coccinellidae, in which 101 barcodes were newly reported for the first time. We also downloaded 870 COI barcode sequences (658 bp) from GenBank and BOLD database, belonging to 108 species within 46 genera, to assess the optimum genetic distance threshold and compare four methods of species delimitation (GMYC, bPTP, BIN and ABGD) to determine the most accurate approach for the family. The results suggested the existence of a 'barcode gap' and that 3% is likely an appropriate genetic distance threshold to delimit species of Coccinellidae using DNA barcodes. Species delimitation analyses confirm ABGD as an accurate and efficient approach, more suitable than the other three methods. Our research provides an integrative framework for DNA barcoding and descriptions of new taxa in Coccinellidae. Our results enrich DNA barcoding public reference libraries, including data for Chinese coccinellids. This will facilitate taxonomic identification and biodiversity monitoring of ladybirds using metabarcoding.

Molecular and Morphological Identifications Reveal Species Composition of Lygus (Hemiptera: Miridae) Bugs in Potatoes Fields in the Lower Columbia Basin of the United States

Lygus bugs are highly polyphagous insect pests. In recent years, Lygus bugs have become more conspicuous on potato, Solanum tuberosum L., fields in the Pacific Northwest, particularly in the Lower Columbia Basin. There are concerns that direct feeding damage or potential pathogen transmission can reduce yield. Lygus species on potatoes in the region are collectively identified as 'Lygus bugs'. Overlapping physical traits and the fact that the same species exhibit morphological variations across a geographic range makes it difficult to identify Lygus to species level. Thus, in this study we used DNA barcodes in combination with morphological characters to identify Lygus species on potatoes. Three species were identified in the Lower Columbia Basin: Lygus hesperus (Knight) and Lygus elisus L. were the most common, whereas Lygus keltoni L. was the least common. Interspecific genetic distances among Lygus species were relatively low, ranging from 0.013 to 0.004. Neighbor-joining (NJ) tree clustered L. hesperus and L. elisus into two major clades, with L. keltoni forming a subclade within L. hesperus clade. Statistical parsimony analysis corroborated findings from phylogenetic analysis with L. keltoni and L. hesperus sharing one haplotype. Our study demonstrates the utility of integrating morphology and molecular markers in identifying morphologically similar species such as Lygus bugs. The study also serves as a fundamental step in contributing to developing suitable management strategies against Lygus bugs on potato.

A two-step DNA barcoding approach for delimiting moth species: moths of Dongling Mountain (Beijing, China) as a case study

DNA barcoding, based on a fragment of cytochrome c oxidase I (COI) mtDNA, is as an effective molecular tool for identification, discovery, and biodiversity assessment for most animals. However, multiple gene markers coupled with more sophisticated analytical approaches may be necessary to clarify species boundaries in cases of cryptic diversity or morphological plasticity. Using 339 moths collected from mountains surrounding Beijing, China, we tested a pipeline consisting of two steps: (1) rapid morphospecies sorting and screening of the investigated fauna with standard COI barcoding approaches; (2) additional analyses with multiple molecular markers for those specimens whose morphospecies and COI barcode grouping were incongruent. In step 1, 124 morphospecies were delimited into 116 barcode units, with 90% of the conflicts being associated with specimens identified to the genus Hypena. In step 2, 55 individuals representing all 12 Hypena morphospecies were analysed using COI, COII, 28S, EF-1a, Wgl sequences or their combinations with the BPP (Bayesian Phylogenetics and Phylogeography) multigene species delimitation method. The multigene analyses supported the delimitation of 5 species, consistent with the COI analysis. We conclude that a two-step barcoding analysis pipeline is able to rapidly characterize insect biodiversity and help to elucidate species boundaries for taxonomic complexes without jeopardizing overall project efficiency by substantially increasing analytical costs.