The effect of geographical scale of sampling on DNA barcoding

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.

A re-analysis of the data in Sharkey et al.'s (2021) minimalist revision reveals that BINs do not deserve names, but BOLD Systems needs a stronger commitment to open science

Halting biodiversity decline is one of the most critical challenges for humanity, but monitoring biodiversity is hampered by taxonomic impediments. One impediment is the large number of undescribed species (here called "dark taxon impediment") whereas another is caused by the large number of superficial species descriptions, that can only be resolved by consulting type specimens ("superficial description impediment"). Recently, Sharkey et al. (2021) proposed to address the dark taxon impediment for Costa Rican braconid wasps by describing 403 species based on COI barcode clusters ("BINs") computed by BOLD Systems. More than 99% of the BINs (387 of 390) were converted into species by assigning binominal names (e.g. BIN "BOLD:ACM9419" becomes Bracon federicomatarritai) and adding a minimal diagnosis (consisting only of a consensus barcode for most species). We here show that many of Sharkey et al.'s species are unstable when the underlying data are analyzed using different species delimitation algorithms. Add the insufficiently informative diagnoses, and many of these species will become the next "superficial description impediment" for braconid taxonomy because they will have to be tested and redescribed after obtaining sufficient evidence for confidently delimiting species. We furthermore show that Sharkey et al.'s approach of using consensus barcodes as diagnoses is not functional because it cannot be applied consistently. Lastly, we reiterate that COI alone is not suitable for delimiting and describing species, and voice concerns over Sharkey et al.'s uncritical use of BINs because they are calculated by a proprietary algorithm (RESL) that uses a mixture of public and private data. We urge authors, reviewers and editors to maintain high standards in taxonomy by only publishing new species that are rigorously delimited with open-access tools and supported by publicly available evidence.

Congruence between morphology-based species and Barcode Index Numbers (BINs) in Neotropical Eumaeini (Lycaenidae)

BACKGROUND

With about 1,000 species in the Neotropics, the Eumaeini (Theclinae) are one of the most diverse butterfly tribes. Correct morphology-based identifications are challenging in many genera due to relatively little interspecific differences in wing patterns. Geographic infraspecific variation is sometimes more substantial than variation between species. In this paper we present a large DNA barcode dataset of South American Lycaenidae. We analyze how well DNA barcode BINs match morphologically delimited species.

METHODS

We compare morphology-based species identifications with the clustering of molecular operational taxonomic units (MOTUs) delimitated by the RESL algorithm in BOLD, which assigns Barcode Index Numbers (BINs). We examine intra- and interspecific divergences for genera represented by at least four morphospecies. We discuss the existence of local barcode gaps in a genus by genus analysis. We also note differences in the percentage of species with barcode gaps in groups of lowland and high mountain genera.

RESULTS

We identified 2,213 specimens and obtained 1,839 sequences of 512 species in 90 genera. Overall, the mean intraspecific divergence value of CO1 sequences was 1.20%, while the mean interspecific divergence between nearest congeneric neighbors was 4.89%, demonstrating the presence of a barcode gap. However, the gap seemed to disappear from the entire set when comparing the maximum intraspecific distance (8.40%) with the minimum interspecific distance (0.40%). Clear barcode gaps are present in many genera but absent in others. From the set of specimens that yielded COI fragment lengths of at least 650 bp, 75% of the a priori morphology-based identifications were unambiguously assigned to a single Barcode Index Number (BIN). However, after a taxonomic a posteriori review, the percentage of matched identifications rose to 85%. BIN splitting was observed for 17% of the species and BIN sharing for 9%. We found that genera that contain primarily lowland species show higher percentages of local barcode gaps and congruence between BINs and morphology than genera that contain exclusively high montane species. The divergence values to the nearest neighbors were significantly lower in high Andean species while the intra-specific divergence values were significantly lower in the lowland species. These results raise questions regarding the causes of observed low inter and high intraspecific genetic variation. We discuss incomplete lineage sorting and hybridization as most likely causes of this phenomenon, as the montane species concerned are relatively young and hybridization is probable. The release of our data set represents an essential baseline for a reference library for biological assessment studies of butterflies in mega diverse countries using modern high-throughput technologies an highlights the necessity of taxonomic revisions for various genera combining both molecular and morphological data.

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

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

Improving metabarcoding taxonomic assignment: A case study of fishes in a large marine ecosystem

DNA metabarcoding is an important tool for molecular ecology. However, its effectiveness hinges on the quality of reference sequence databases and classification parameters employed. Here we evaluate the performance of MiFish 12S taxonomic assignments using a case study of California Current Large Marine Ecosystem fishes to determine best practices for metabarcoding. Specifically, we use a taxonomy cross-validation by identity framework to compare classification performance between a global database comprised of all available sequences and a curated database that only includes sequences of fishes from the California Current Large Marine Ecosystem. We demonstrate that the regional database provides higher assignment accuracy than the comprehensive global database. We also document a tradeoff between accuracy and misclassification across a range of taxonomic cutoff scores, highlighting the importance of parameter selection for taxonomic classification. Furthermore, we compared assignment accuracy with and without the inclusion of additionally generated reference sequences. To this end, we sequenced tissue from 597 species using the MiFish 12S primers, adding 252 species to GenBank's existing 550 California Current Large Marine Ecosystem fish sequences. We then compared species and reads identified from seawater environmental DNA samples using global databases with and without our generated references, and the regional database. The addition of new references allowed for the identification of 16 additional native taxa representing 17.0% of total reads from eDNA samples, including species with vast ecological and economic value. Together these results demonstrate the importance of comprehensive and curated reference databases for effective metabarcoding and the need for locus-specific validation efforts.

Barcoding pest species in a biodiversity hot-spot: the South African polyphagous broad-nosed weevils (Coleoptera, Curculionidae, Entiminae)

Polyphagous broad nosed weevils (Curculionidae: Entiminae) constitute a large and taxonomically challenging subfamily that contains economically significant agricultural pests worldwide. South Africa is a hot-spot for biodiversity and several species of indigenous and endemic genera of Entiminae have shifted on to cultivated plants, with some being phytosanitary pests. The sporadic pest status of many species (where the species has an occasional economic impact on the agricultural industry, but is not encountered often enough that is is readily recognisable by researchers and agricultural extension workers) and the presence of pest complexes and cryptic species represent an identification challenge to non-specialists. Furthermore, no comprehensive identification tools exist to identify immature stages that may be found in crops/soil. In this paper, a curated barcoding database with 70 COI sequences from 41 species (39 Entiminae, 2 Cyclominae) is initiated, to assist with the complexity of identification of species in this group.

Attraction of Trichogramma Wasps to Butterfly Oviposition-Induced Plant Volatiles Depends on Brassica Species, Wasp Strain and Leaf Necrosis

Within the Brassicaceae, wild as well as crop species are challenged by specialist herbivores including cabbage white butterflies (Pieris spp.). The wild crucifer Brassica nigra responds to oviposition by Pieris butterflies by the synergistic expression of two egg-killing traits. Genotypes that express a hypersensitive response (HR)-like necrosis (direct egg-killing) also emit oviposition-induced plant volatiles (OIPVs) attracting Trichogramma egg parasitoids (indirect egg-killing). This so-called double defense line can result in high butterfly egg mortalities. It remains unknown whether this strategy is unique to B. nigra or more common in Brassica species. To test this, we examined the response of different Trichogramma evanescens lines to OIPVs emitted by B. nigra and three close relatives (Brassica napus, Brassica rapa, and Brassica oleracea). Furthermore, we evaluated whether HR-like necrosis played a role in the attraction toward plant volatiles. Our results show a specificity in wasp attraction to different plant species. Three out of four plant species attracted a specific T. evanescens strain, including the crops B. rapa and B. napus. Parasitoid attraction was positively affected by presence of HR-like necrosis in one plant species. Our findings imply that, despite being a true generalist in terms of host range, T. evanescens shows intraspecific variation during host searching, which should be taken into account when selecting parasitoid lines for biocontrol of certain crops. Finally, we conclude that also crop plants within the Brassicaceae family possess egg-killing traits and can exert the double-defense line which may enable effective selection of egg-killing defense traits by cabbage breeders.

Molecular Taxonomy and Diversification of Atlantic Skates (Chondrichthyes, Rajiformes): Adding More Pieces to the Puzzle of Their Evolutionary History

Conservation and long-term management plans of marine species need to be based upon the universally recognized key-feature of species identity. This important assignment is particularly challenging in skates (Rajiformes) in which the phenotypic similarity between some taxa and the individual variability in others, hampers accurate species identification. Here, 432 individual skate samples collected from four major ocean areas of the Atlantic were barcoded and taxonomically analysed. A BOLD project ELASMO ATL was implemented with the aim of establishing a new fully available and well curated barcode library containing both biological and molecular information. The evolutionary histories of the 38 skate taxa were estimated with two concatenated mitochondrial markers (COI and NADH2) through Maximum Likelihood and Bayesian inference. New evolutionary lineages within the genus Raja were discovered off Angola, where paleogeographic history coupled with oceanographic discontinuities could have contributed to the establishment of isolated refugia, playing a fundamental role among skates' speciation events. These data successfully resolved many taxonomic ambiguities, identified cryptic diversity within valid species and demonstrated a highly cohesive monophyletic clustering among the order, laying the background for further inference of evolutionary patterns suitable for addressing management and conservation issues.

Image-based Automated Species Identification: Can Virtual Data Augmentation Overcome Problems of Insufficient Sampling?

Automated species identification and delimitation is challenging, particularly in rare and thus often scarcely sampled species, which do not allow sufficient discrimination of infraspecific versus interspecific variation. Typical problems arising from either low or exaggerated interspecific morphological differentiation are best met by automated methods of machine learning that learn efficient and effective species identification from training samples. However, limited infraspecific sampling remains a key challenge also in machine learning. In this study, we assessed whether a data augmentation approach may help to overcome the problem of scarce training data in automated visual species identification. The stepwise augmentation of data comprised image rotation as well as visual and virtual augmentation. The visual data augmentation applies classic approaches of data augmentation and generation of artificial images using a Generative Adversarial Networks (GAN) approach. Descriptive feature vectors are derived from bottleneck features of a VGG-16 convolutional neural network (CNN) that are then stepwise reduced in dimensionality using Global Average Pooling and PCA to prevent overfitting. Finally, data augmentation employs synthetic additional sampling in feature space by an oversampling algorithm in vector space (SMOTE). Applied on four different image datasets, which include scarab beetle genitalia (Pleophylla, Schizonycha) as well as wing patterns of bees (Osmia) and cattleheart butterflies (Parides), our augmentation approach outperformed a deep learning baseline approach by means of resulting identification accuracy with non-augmented data as well as a traditional 2D morphometric approach (Procrustes analysis of scarab beetle genitalia).

Genetic variation in neotropical butterflies is associated with sampling scale, species distributions, and historical forest dynamics

Previous studies of butterfly diversification in the Neotropics have focused on Amazonia and the tropical Andes, while southern regions of the continent have received little attention. To address the gap in knowledge about the Lepidoptera of temperate South America, we analysed over 3000 specimens representing nearly 500 species from Argentina for a segment of the mitochondrial cytochrome c oxidase subunit I (COI) gene. Representing 42% of the country's butterfly fauna, collections targeted species from the Atlantic and Andean forests, and biodiversity hotspots that were previously connected but are now isolated. We assessed COI effectiveness for species discrimination and identification and how its performance was affected by geographic distances and taxon coverage. COI data also allowed to study patterns of genetic variation across Argentina, particularly between populations in the Atlantic and Andean forests. Our results show that COI discriminates species well, but that identification success is reduced on average by ~20% as spatial and taxonomic coverage rises. We also found that levels of genetic variation are associated with species' spatial distribution type, a pattern which might reflect differences in their dispersal and colonization abilities. In particular, intraspecific distance between populations in the Atlantic and Andean forests was significantly higher in species with disjunct distributions than in those with a continuous range. All splits between lineages in these forests dated to the Pleistocene, but divergence dates varied considerably, suggesting that historical connections between the Atlantic and Andean forests have differentially affected their shared butterfly fauna. Our study supports the fact that large-scale assessments of mitochondrial DNA variation are a powerful tool for evolutionary studies.

Coverage and quality of DNA barcode references for Central and Northern European Odonata

Background

Dragonflies and damselflies (Odonata) are important components in biomonitoring due to their amphibiotic lifecycle and specific habitat requirements. They are charismatic and popular insects, but can be challenging to identify despite large size and often distinct coloration, especially the immature stages. DNA-based assessment tools rely on validated DNA barcode reference libraries evaluated in a supraregional context to minimize taxonomic incongruence and identification mismatches.

Methods

This study reports on findings from the analysis of the most comprehensive DNA barcode dataset for Central European Odonata to date, with 103 out of 145 recorded European species included and publicly deposited in the Barcode of Life Data System (BOLD). The complete dataset includes 697 specimens (548 adults, 108 larvae) from 274 localities in 16 countries with a geographic emphasis on Central Europe. We used BOLD to generate sequence divergence metrics and to examine the taxonomic composition of the DNA barcode clusters within the dataset and in comparison with all data on BOLD.

Results

Over 88% of the species included can be readily identified using their DNA barcodes and the reference dataset provided. Considering the complete European dataset, unambiguous identification is hampered in 12 species due to weak mitochondrial differentiation and partial haplotype sharing. However, considering the known species distributions only two groups of five species possibly co-occur, leading to an unambiguous identification of more than 95% of the analysed Odonata via DNA barcoding in real applications. The cases of small interspecific genetic distances and the observed deep intraspecific variation in Cordulia aenea (Linnaeus, 1758) are discussed in detail and the corresponding taxa in the public reference database are highlighted. They should be considered in future applications of DNA barcoding and metabarcoding and represent interesting evolutionary biological questions, which call for in depth analyses of the involved taxa throughout their distribution ranges.

A DNA barcode library for the butterflies of North America

Although the butterflies of North America have received considerable taxonomic attention, overlooked species and instances of hybridization continue to be revealed. The present study assembles a DNA barcode reference library for this fauna to identify groups whose patterns of sequence variation suggest the need for further taxonomic study. Based on 14,626 records from 814 species, DNA barcodes were obtained for 96% of the fauna. The maximum intraspecific distance averaged 1/4 the minimum distance to the nearest neighbor, producing a barcode gap in 76% of the species. Most species (80%) were monophyletic, the others were para- or polyphyletic. Although 15% of currently recognized species shared barcodes, the incidence of such taxa was far higher in regions exposed to Pleistocene glaciations than in those that were ice-free. Nearly 10% of species displayed high intraspecific variation (>2.5%), suggesting the need for further investigation to assess potential cryptic diversity. Aside from aiding the identification of all life stages of North American butterflies, the reference library has provided new perspectives on the incidence of both cryptic and potentially over-split species, setting the stage for future studies that can further explore the evolutionary dynamics of this group.

On the Diversity of Phyllodocida (Annelida: Errantia), with a Focus on Glyceridae, Goniadidae, Nephtyidae, Polynoidae, Sphaerodoridae, Syllidae, and the Holoplanktonic Families

Phyllodocida is a clade of errantiate annelids characterized by having ventral sensory palps, anterior enlarged cirri, axial muscular proboscis, compound chaetae (if present) with a single ligament, and of lacking dorsolateral folds. Members of most families date back to the Carboniferous, although the earliest fossil was dated from the Devonian. Phyllodocida holds 27 well-established and morphologically homogenous clades ranked as families, gathering more than 4600 currently accepted nominal species. Among them, Syllidae and Polynoidae are the most specious polychaete groups. Species of Phyllodocida are mainly found in the marine benthos, although a few inhabit freshwater, terrestrial and planktonic environments, and occur from intertidal to deep waters in all oceans. In this review, we (1) explore the current knowledge on species diversity trends (based on traditional species concept and molecular data), phylogeny, ecology, and geographic distribution for the whole group, (2) try to identify the main knowledge gaps, and (3) focus on selected families: Alciopidae, Goniadidae, Glyceridae, Iospilidae, Lopadorrhynchidae, Polynoidae, Pontodoridae, Nephtyidae, Sphaerodoridae, Syllidae, Tomopteridae, Typhloscolecidae, and Yndolaciidae. The highest species richness is concentrated in European, North American, and Australian continental shelves (reflecting a strong sampling bias). While most data come from shallow coastal and surface environments most world oceans are clearly under-studied. The overall trends indicate that new descriptions are constantly added through time and that less than 10% of the known species have molecular barcode information available.

Factors affecting the efficiency of molecular species delimitation in a species-rich insect family

In the context of global biodiversity loss, molecular species delimitation approaches can be very useful for accelerating species discovery through DNA taxonomy and inventory through DNA metabarcoding. In this study, the effect of some intrinsic factors on the efficiency of various single-marker species delimitation methods (fixed and variable nucleotide distance thresholds, ABGD, ASAP, GMYC, mPTP) was tested on more than 90 empirical data sets, derived from a set of 7,237 COI sequences attributed to 542 leaf beetles species (Coleoptera: Chrysomelidae). The considered factors were: (i) the number of haplotypes per species (as a proxy for genetic diversity), (ii) the geographic distance among conspecific collection localities (as a proxy of sampling width), (iii) the difficulty related to morphological identification of species, and (iv) the taxonomic rank. Distance-based methods, with on average more than 70% of match with morphological identification, outperformed those relying on phylogenetic trees, with less than 59%. A high number of haplotypes per species was found to have a negative effect on delimitation efficiency, whereas large geographic distances within species had a positive effect. All methods delimitations (except for GMYC) were significantly affected by the presence of species that are difficult to be identified, decreasing their efficiency. Finally, the only method influenced by the taxonomic rank of the data set was GMYC, showing lower efficiency in data sets at the genus than at higher levels. The observed biases we highlighted affecting efficiency could be accounted for when developing input data sets for species delimitation analyses to obtain a more reliable representation of biological diversity.

Molecular species delimitation of the genera Anodus, Argonectes, Bivibranchia and Micromischodus (Ostariophysi: Characiformes)

ABSTRACT A good taxonomic assessment of specimens is an essential task to many biological studies and DNA data have provided additional sources of information to assist in the disentanglement of taxonomic problems among living organisms, as has been the case of some taxa of the megadiverse Neotropical ichthyofauna. Here we assessed all valid species in the Neotropical freshwater fish genera Anodus, Argonectes, Bivibranchia and Micromischodus of the family Hemiodontidae to establish molecular species boundaries among them. All species delimitation methods defined exactly only one MOTU for Anodus elongatus, Argonectes longiceps, A. robertsi, Bivibranchia bimaculata, B. notata, B. velox, and Micromischodus sugillatus, resulting in total congruence between nominal species and MOTUs for these seven taxa. The three species having discordant results across analyses: Anodus orinocensis, Bivibranchia fowleri, and Bivibranchia simulata, matched more than one MOTU per species in some methods, meaning that cryptic diversity may exist within these taxa. Overall, this great correspondence among morphological and molecular boundaries for thae species analysed seem to be indicative of a reasonably stable taxonomy within these Hemiodontidae genera.

Perlidae (Plecoptera) from the Paranapiacaba Mountains, Atlantic Forest, Brazil: Diversity and implications of the integrative approach and teneral specimens on taxonomy

The study of complementary sources of biological variation (e.g. morphological, molecular) has allowed a better understanding of biodiversity through the construction of an integrative taxonomy. Using this approach, specimens from the Paranapiacaba Mountains, southeastern Brazil, were studied to update the knowledge on the stonefly family Perlidae from the region, characterize the species, and make associations between nymphs and adults using a fragment of the mitochondrial cytochrome c oxidase subunit I (COI) gene. The study also discusses the implications of integrative taxonomy and teneral specimens for the study of South American Perlidae. The molecular data were analyzed using Bayesian inference, Neighbor-joining, and delimiting species methods. Our results revealed that, in general, there was a morphological and molecular congruence between species. In the Paranapiacaba Mountains, three genera and 15 species were recorded: Anacroneuria boraceiensis Froehlich 2004, A. debilis (Pictet 1841) (new record), A. fiorentini De Ribeiro and Froehlich 2007 (new record), A. flintorum Froehlich 2002, A. iporanga Bispo and Froehlich 2004, A. itajaimirim Bispo and Froehlich 2004, A. polita (Burmeister 1913), A. subcostalis Klapálek 1921, A. tupi Bispo and Froehlich 2004 (with a description of the nymph), Kempnyia auberti Froehlich 1996, K. colossica (Navás 1934), K. flava Klapálek 1916, K. neotropica (Jacobson and Bianchi 1905) (including its new junior synonym K. petersorum Froehlich 1996), Kempnyia sp., and Macrogynoplax veneranda Froehlich 1984. COI sequences were obtained for 11 species, five of which had nymphs associated with adults. Among the five associated nymphs, the nymph of A. tupi is described here. The results of this study indicate that the color of adult teneral specimens differs from that of mature specimens. Given this, the synonym of K. neotropica and K. petersorum was proposed since these species have high morphological and molecular similarities and differ only in color patterns. In addition, the previous record of A. petersi Froehlich 2002 from the Paranapiacaba Mountains was invalidated since it was considered a teneral specimen of A. flintorum. These results suggest that the development of an integrative taxonomy is essential to continue advancing the study of Perlidae diversity in South America.

Multiple species delimitation approaches applied to the avian lark genus Alaudala

Species delimitation has advanced from a purely phenotypic exercise to a branch of science that integrates multiple sources of data to identify independently evolving lineages that can be treated as species. We here test species limits in the avian Lesser Short-toed Lark Alaudala rufesens-Sand Lark A. raytal complex, which has an intricate taxonomic history, ranging from a single to three recognised species, with different inclusiveness in different treatments. Our integrative taxonomic approach is based on a combination of DNA sequences, plumage, biometrics, songs, song-flights, geographical distributions, habitat, and bioclimatic data, and using various methods including a species delimitation program (STACEY) based on the multispecies coalescent model. We propose that four species should be recognised: Lesser Short-toed Lark A. rufescens (sensu stricto), Heine's Short-toed Lark A. heinei, Asian Short-toed Lark A. cheleensis and Sand Lark A. raytal. There is also some evidence suggesting lineage separation within A. cheleensis and A. raytal, but additional data are required to evaluate this. The species delimitation based on STACEY agrees well with the non-genetic data. Although computer-based species delimitation programs can be useful in identifying independently evolving lineages, we stress that whenever possible, species hypotheses proposed by these programs should be tested by independent, non-genetic data. Our results highlight the difficulty and subjectivity of delimiting lineages and species, especially at early stages in the speciation process.

DNA barcoding for identification of anuran species in the central region of South America

The use of COI barcodes for specimen identification and species discovery has been a useful molecular approach for the study of Anura. Here, we establish a comprehensive amphibian barcode reference database in a central area of South America, in particular for specimens collected in Mato Grosso do Sul state (Brazil), and to evaluate the applicability of the COI gene for species-level identification. Both distance- and tree-based methods were applied for assessing species boundaries and the accuracy of specimen identification was evaluated. A total of 204 mitochondrial COI barcode sequences were evaluated from 22 genera and 59 species (19 newly barcoded species). Our results indicate that morphological and molecular identifications converge for most species, however, some species may present cryptic species due to high intraspecific variation, and there is a high efficiency of specimen identification. Thus, we show that COI sequencing can be used to identify anuran species present in this region.

First DNA sequence reference library for mammals and plants of the Eastern Mediterranean Region

The Mediterranean region is identified as one of the world's 36 biodiversity hotspots, with the Earth's most biologically rich yet threatened areas. Lebanon is a hub for Eastern Mediterranean Region (EMR) biodiversity with 9116 characterized plant and animal species (4486 fauna and 4630 flora). Using DNA barcoding as a tool has become crucial in the accurate identification of species in multiple contexts. It can also complement species morphological descriptions, which will add to our understanding of the biodiversity and richness of ecosystems and benefit conservation projects for endangered and endemic species. In this study, we create the first reference library of standard DNA markers for mammals and plants in the EMR, with a focus on endemic and endangered species. Plant leaves were collected from different nature reserves in Mount Lebanon, and mammal samples were obtained from taxidermized museum specimens or road kills. We generated the 12S rRNA sequences of 18 mammal species from 6 orders and 13 different families. We also obtained the trnL and rbcL barcode sequences of 52 plant species from 24 different families. Twenty-five plant species and two mammal species included in this study were sequenced for the first time using these markers.

Genetic assays for guano-based identification of species and sex in bats of the United States and Canada

Bat guano is a noninvasive, data-rich genetic resource. However, the constituent bat DNA is relatively scant, degraded, and complexed with polymerase chain reaction inhibitors. It also is comingled with a rich pool of nontarget DNA from microbes, parasites, and dietary items. We designed and tested new DNA assays for bat species identification (COX1-Bat) and sex identification (XGXYC) for use with guano and other challenging samples. We reviewed previously published assays that can be used with guano samples to obtain the same species and sex data, and attempted to validate these assays for species in which they had not previously been tested. Our results demonstrate that guano-derived DNA can be used successfully to 1) identify nearly all US and Canadian bats at the species level, or to one of three Myotis species clusters, and 2) identify the sex of at least 23 US and Canadian bat species. Our newly developed assays, and validation of previously published assays, for guano-based identification of species and sex in bats, significantly enhance the power of noninvasive sampling and genetic analysis for bat studies, management, and conservation.

DNA Barcoding and geographical scale effect: The problems of undersampling genetic diversity hotspots

DNA barcoding identification needs a good characterization of intraspecific genetic divergence to establish the limits between species. Yet, the number of barcodes per species is many times low and geographically restricted. A poor coverage of the species distribution range may hamper identification, especially when undersampled areas host genetically distinct lineages. If so, the genetic distance between some query sequences and reference barcodes may exceed the maximum intraspecific threshold for unequivocal species assignation. Taking a group of Quercus herbivores (moths) in Europe as model system, we found that the number of DNA barcodes from southern Europe is proportionally very low in the Barcoding of Life Data Systems. This geographical bias complicates the identification of southern query sequences, due to their high intraspecific genetic distance with respect to barcodes from higher latitudes. Pairwise intraspecific genetic divergence increased along with spatial distance, but was higher when at least one of the sampling sites was in southern Europe. Accordingly, GMYC (General Mixed Yule Coalescent) single-threshold model retrieved clusters constituted exclusively by Iberian haplotypes, some of which could correspond to cryptic species. The number of putative species retrieved was more reliable than that of multiple-threshold GMYC but very similar to results from ABGD and jMOTU. Our results support GMYC as a key resource for species delimitation within poorly inventoried biogeographic regions in Europe, where historical factors (e.g., glaciations) have promoted genetic diversity and singularity. Future European DNA barcoding initiatives should be preferentially performed along latitudinal gradients, with special focus on southern peninsulas.

Italian odonates in the Pandora's box: A comprehensive DNA barcoding inventory shows taxonomic warnings at the Holarctic scale

The Odonata are considered among the most endangered freshwater faunal taxa. Their DNA-based monitoring relies on validated reference data sets that are often lacking or do not cover important biogeographical centres of diversification. This study presents the results of a DNA barcoding campaign on Odonata, based on the standard 658-bp 5' end region of the mitochondrial COI gene, involving the collection of 812 specimens (409 of which barcoded) from peninsular Italy and its main islands (328 localities), belonging to all the 88 species (31 Zygoptera and 57 Anisoptera) known from the country. Additional BOLD and GenBank data from Holarctic samples expanded the data set to 1,294 DNA barcodes. A multi-approach species delimitation analysis involving two distance (OT and ABGD) and four tree-based (PTP, MPTP, GMYC and bGMYC) methods was used to explore these data. Of the 88 investigated morphospecies, 75 (85%) unequivocally corresponded to distinct molecular operational units, whereas the remaining ones were classified as 'warnings' (i.e. showing a mismatch between morphospecies assignment and DNA-based species delimitation). These results are in contrast with other DNA barcoding studies on Odonata showing up to 95% of identification success. The species causing warnings were grouped into three categories depending on if they showed low, high or mixed genetic divergence patterns. The analysis of haplotype networks revealed unexpected intraspecific complexity at the Italian, Palearctic and Holarctic scale, possibly indicating the occurrence of cryptic species. Overall, this study provides new insights into the taxonomy of odonates and a valuable basis for future DNA and eDNA-based monitoring studies.

Sharpening the DNA barcoding tool through a posteriori taxonomic validation: The case of Longitarsus flea beetles (Coleoptera: Chrysomelidae)

The accuracy of the DNA barcoding tool depends on the existence of a comprehensive archived library of sequences reliably determined at species level by expert taxonomists. However, misidentifications are not infrequent, especially following large-scale DNA barcoding campaigns on diverse and taxonomically complex groups. In this study we used the species-rich flea beetle genus Longitarsus, that requires a high level of expertise for morphological species identification, as a case study to assess the accuracy of the DNA barcoding tool following several optimization procedures. We built a cox1 reference database of 1502 sequences representing 78 Longitarsus species, among which 117 sequences (32 species) were newly generated using a non-invasive DNA extraction method that allows keeping reference voucher specimens. Within this dataset we identified 69 taxonomic inconsistencies using barcoding gap analysis and tree topology methods. Threshold optimisation and a posteriori taxonomic revision based on newly generated reference sequences and metadata allowed resolving 44 sequences with ambiguous and incorrect identification and provided a significant improvement of the DNA barcoding accuracy and identification efficacy. Unresolved taxonomic uncertainties, due to overlapping intra- and inter-specific levels of divergences, mainly regards the Longitarsus pratensis species complex and polyphyletic groups L. melanocephalus, L. nigrofasciatus and L. erro. Such type of errors indicates either poorly established taxonomy or any biological processes that make mtDNA groups poorly predictive of species boundaries (e.g. recent speciation or interspecific hybridisation), thus providing directions for further integrative taxonomic and evolutionary studies. Overall, this study underlines the importance of reference vouchers and high-quality metadata associated to sequences in reference databases and corroborates, once again, the key role of taxonomists in any step of the DNA barcoding pipeline in order to generate and maintain a correct and functional reference library.

Highly variable COI haplotype diversity between three species of invasive pest fruit fly reflects remarkably incongruent demographic histories

Distance decay principles predict that species with larger geographic ranges would have greater intraspecific genetic diversity than more restricted species. However, invasive pest species may not follow this prediction, with confounding implications for tracking phenomena including original ranges, invasion pathways and source populations. We sequenced an 815 base-pair section of the COI gene for 441 specimens of Bactrocera correcta, 214 B. zonata and 372 Zeugodacus cucurbitae; three invasive pest fruit fly species with overlapping hostplants. For each species, we explored how many individuals would need to be included in a study to sample the majority of their haplotype diversity. We also tested for phylogeographic signal and used demographic estimators as a proxy for invasion potency. We find contrasting patterns of haplotype diversity amongst the species, where B. zonata has the highest diversity but most haplotypes were represented by singletons; B. correcta has ~7 dominant haplotypes more evenly distributed; Z. cucurbitae has a single dominant haplotype with closely related singletons in a 'star-shape' surrounding it. We discuss how these differing patterns relate to their invasion histories. None of the species showed meaningful phylogeographic patterns, possibly due to gene-flow between areas across their distributions, obscuring or eliminating substructure.

New insights into the haplotype diversity of the cosmopolitan cat flea Ctenocephalides felis (Siphonaptera: Pulicidae)

The present study investigated the genetic profile of the cosmopolitan cat flea, Ctenocephalides felis (Siphonaptera: Pulicidae) from Malaysia and the reference data available in the National Center for Biotechnology Information (NCBI) GenBank. A set of sequences of 100 Malaysian samples aligned as 550 characters of the cytochrome c oxidase subunit I (cox1) and 706 characters of the II (cox2) genes revealed ten haplotypes (A1-A10) and eight haplotypes (B1-B8), respectively. The concatenated sequences of cox1 and cox2 genes with a total of 1256 characters revealed 15 haplotypes (AB1-AB15). Analyses indicated that haplotype AB1 was the most frequent and the most widespread haplotype in Malaysia. Overall haplotype and nucleotide diversities of the concatenated sequences were 0.52909 and 0.00424, respectively, with moderate genetic differentiation (F_ST = 0.17522) and high gene flow (N_m = 1.18). The western population presented the highest genetic diversity (H_d = 0.78333, P_i = 0.01269, N_h = 9), whereas the southern population demonstrated the lowest diversity (H_d = 0.15667, P_i = 0.00019, N_h = 3). The concatenated sequences showed genetic distances ranged from 0.08 % to 4.39 %. There were three aberrant haplotypes in cox2 sequences that highly divergent, suggesting the presence of cryptic species or occurrence of introgression. In the global point of view, the aligned sequences of C. felis revealed 65 haplotypes (AA1-AA65) by the cox1 gene (n = 586), and 27 haplotypes (BB1-BB27) by the cox2 gene (n = 204). Mapping of the haplotype network showed that Malaysian C. felis possesses seven unique haplotypes in both genes with the common haplotypes demonstrated genetic affinity with C. felis from Southeast Asia for cox1 and South America for cox2. The topologies of cox1 and cox2 phylogenetic trees were concordant with relevant grouping pattern of haplotypes in the network but revealed two major lineages by which Malaysian haplotypes were closely related with haplotypes from the tropical region.

Twenty-six new species of Hoploscopa (Lepidoptera, Crambidae) from South-East Asia revealed by morphology and DNA barcoding

Hoploscopa Meyrick (Lepidoptera: Crambidae) is a fern-feeding genus found in montane areas of South-East Asia and Melanesia, eastwards up to the Samoan Islands. It includes sixteen described species, with at least 70 further undescribed species known from scientific collections. An iterative approach including morphological and molecular characters was used in order to explore the diversity of Hoploscopa. The hitherto described species are revised, and descriptions authored by T. Léger and M. Nuss are provided for an additional 26 new species: H.agtuuganonensissp. nov., H.albipunctasp. nov., H.albomaculatasp. nov., H.anacanthasp. nov., H.boletasp. nov., H.cynodontasp. nov., H.danaoensissp. nov., H.gombongisp. nov., H.gracilissp. nov., H.ignitamaculaesp. nov., H.isarogensissp. nov., H.jubatasp. nov., H.kelamasp. nov., H.kinabaluensissp. nov., H.mallyisp. nov., H.marijoweissaesp. nov., H.matheaesp. nov., H.niveofasciasp. nov., H.pangrangoensissp. nov., H.parvimaculasp. nov., H.pseudometacrossasp. nov., H.sepanggisp. nov., H.sumatrensissp. nov., H.titikasp. nov., H.tonsepisp. nov., H.ypsilonsp. nov. Using a protocol specific for the amplification of DNA from old museum specimens, we recovered 101 COI barcodes for all but one of the newly described species, with 76 being barcode compliant (>487 bp). Species delimitation analyses suggest cryptic diversity, with six cases reflecting allopatric divergence, and two further cases found in sympatry.

DNA barcode approaches to reveal interspecies genetic variation of Indian ungulates

In the past two decades, identification of species from noninvasive sampling has turned out to be an important tool for wildlife conservation. In this study a total 93 specimens representing 22 species of ungulates were analyzed from partial sequences of mtDNA COI and Cytb genes. All the species showed unique clades, and sequences divergence within species was between 0.01–3.9% in COI and 0.01–13.7 in Cytb, whereas divergence between species ranged from 2.2 to 29.5% in COI and 2.3 to 28.8% in Cytb. Highest intraspecific divergence was observed within the Ovis aries in COI and Porcula salvania in Cytb. Bayesian (BA) phylogeny analysis of both genes combined distinguishes all the studied species as monophyletic criteria. The Indian rhinoceros (Rhinoceros unicornis) exhibited closer relation to horse (Equus caballus). No barcode gap was observed between species in COI. This study demonstrates that even short fragments of COI and Cytb generated from fecal pellets can efficiently identify the Indian ungulates, thus demonstrating its high potential for use in wildlife conservation activities.

HACSim: an R package to estimate intraspecific sample sizes for genetic diversity assessment using haplotype accumulation curves

Assessing levels of standing genetic variation within species requires a robust sampling for the purpose of accurate specimen identification using molecular techniques such as DNA barcoding; however, statistical estimators for what constitutes a robust sample are currently lacking. Moreover, such estimates are needed because most species are currently represented by only one or a few sequences in existing databases, which can safely be assumed to be undersampled. Unfortunately, sample sizes of 5-10 specimens per species typically seen in DNA barcoding studies are often insufficient to adequately capture within-species genetic diversity. Here, we introduce a novel iterative extrapolation simulation algorithm of haplotype accumulation curves, called HACSim (Haplotype Accumulation Curve Simulator) that can be employed to calculate likely sample sizes needed to observe the full range of DNA barcode haplotype variation that exists for a species. Using uniform haplotype and non-uniform haplotype frequency distributions, the notion of sampling sufficiency (the sample size at which sampling accuracy is maximized and above which no new sampling information is likely to be gained) can be gleaned. HACSim can be employed in two primary ways to estimate specimen sample sizes: (1) to simulate haplotype sampling in hypothetical species, and (2) to simulate haplotype sampling in real species mined from public reference sequence databases like the Barcode of Life Data Systems (BOLD) or GenBank for any genomic marker of interest. While our algorithm is globally convergent, runtime is heavily dependent on initial sample sizes and skewness of the corresponding haplotype frequency distribution.

DNA Barcoding Subtropical Aphids and Implications for Population Differentiation

DNA barcoding has proven its worth in species identification, discovering cryptic diversity, and inferring genetic divergence. However, reliable DNA barcode reference libraries that these applications depend on are not available for many taxonomic groups and geographical regions. Aphids are a group of plant sap sucking insects, including many notorious pests in agriculture and forestry. The aphid fauna of the subtropical region has been understudied. In this study, based on extensive sampling effort across main subtropical areas, we sequenced 1581 aphid specimens of 143 morphospecies, representing 75 genera, and 13 subfamilies, to build the first comprehensive DNA barcode library for subtropical aphids. We examined the utility of DNA barcodes in identifying aphid species and population differentiation and evaluated the ability of different species delimitation methods (automatic barcode gap discovery (ABGD), generalized mixed Yule-coalescent (GMYC), and Bayesian Poisson tree processes (bPTP)). We found that most aphid species demonstrated barcode gaps and that a threshold value of 2% genetic distance is suitable for distinguishing most species. Our results indicated that ten morphospecies may have species divergence related to factors such as host plant or geography. By using two pest species Aphis spiraecola and A. gossypii as examples, we also discussed the effect of the sampling scale of host plants on the results and reliability of DNA barcoding of phytophagous insects. This DNA barcode library will be valuable for future studies and applications.

BIN overlap confirms transcontinental distribution of pest aphids (Hemiptera: Aphididae)

DNA barcoding is highly effective for identifying specimens once a reference sequence library is available for the species assemblage targeted for analysis. Despite the great need for an improved capacity to identify the insect pests of crops, the use of DNA barcoding is constrained by the lack of a well-parameterized reference library. The current study begins to address this limitation by developing a DNA barcode reference library for the pest aphids of Pakistan. It also examines the affinities of these species with conspecific populations from other geographic regions based on both conventional taxonomy and Barcode Index Numbers (BINs). A total of 809 aphids were collected from a range of plant species at sites across Pakistan. Morphological study and DNA barcoding allowed 774 specimens to be identified to one of 42 species while the others were placed to a genus or subfamily. Sequences obtained from these specimens were assigned to 52 BINs whose monophyly were supported by neighbor-joining (NJ) clustering and Bayesian inference. The 42 species were assigned to 41 BINs with 38 showing BIN concordance. These species were represented on BOLD by 7,870 records from 69 countries. Combining these records with those from Pakistan produced 60 BINs with 12 species showing a BIN split and three a BIN merger. Geo-distance correlations showed that intraspecific divergence values for 49% of the species were not affected by the distance between populations. Forty four of the 52 BINs from Pakistan had counterparts in 73 countries across six continents, documenting the broad distributions of pest aphids.

Molecular phylogeny and species delimitation of the genus Dicerapanorpa (Mecoptera: Panorpidae)

Given that species is the fundamental unit in systematic biology, rigorous species delimitation is crucial for taxonomic studies, yet routine species delimitation remains an ongoing challenge in the taxonomic practice of insects. The two-horned scorpionfly Dicerapanorpa is a small genus in Panorpidae (Mecoptera) endemic to the Qinling-Bashan and Hengduan mountains, a biodiversity hotspot. However, species of Dicerapanorpa are difficult to delineate owing to marked intraspecific variation and interspecific similarity. Here, we investigate the diversity and species boundaries of Dicerapanorpa using an integrative approach based on DNA barcoding, morphological, geometric morphometric and molecular phylogenetic analyses. This integrative analyses confirmed the 13 described species of Dicerapanorpa and revealed three new species: Dicerapanorpa lativalva sp. nov., Dicerapanorpa hualongshana sp. nov. and Dicerapanorpa minshana sp. nov. Most molecular operational taxonomic units are in congruence with morphological clusters. Possible reasons for several discordances in Dicerapanorpa are tentatively discussed.

DNA barcoding of large oak-living cerambycids: diagnostic tool, phylogenetic insights and natural hybridization between Cerambyx cerdo and Cerambyx welensii (Coleoptera: Cerambycidae)

Three large saproxylic cerambycids with different pest/legal status co-occur in the Iberian oak woodlands, Cerambyx welensii (Cw), Cerambyx cerdo (Cc) and Prinobius myardi (Pm): Cw is an emerging pest, Cc is a protected but sometimes harmful species and Pm is a secondary/minor pest. A precise taxonomic diagnosis is necessary for research, management or protection purposes, but may be problematic mainly because Cw and Cc larvae are morphologically indistinguishable. To resolve this constraint, we genotyped adults, larvae and eggs collected over a wide geographical range using the mitochondrial barcoding of the cytochrome c oxidase subunit I (COI). A Neighbour-Joining tree phylogram revealed three distinct clusters corresponding to Cw, Cc and Pm. We further first sequenced for Cw and Cc two mitochondrial (12S rRNA and 16S rRNA) and one nuclear (28S rRNA) gene fragments. For the first two genes, interspecific divergence was lower than in COI, and for the 28S (lower mutation rate), the two species shared identical haplotypes. Two approaches for species delimitation (General Mixed Yule Coalescent (GMYC), Barcode Index Number (BIN)) confirmed the species distinctiveness of Cc and Cw. The Bayesian COI gene tree showed a remarkable genetic divergence between Cc populations from Iberia and the rest of Europe. Such divergence has relevant taxonomic connotations and stresses the importance of a wide geographical scale sampling for accurate DNA barcoding species identification. Incongruities between morphology/lineage and COI barcodes in some individuals revealed natural hybridization between Cw and Cc. Natural hybridization is important from a phylogenetic/evolutionary perspective in these cerambycids, but the prevalence of (and the behavioural/ecological factors involved in) interspecific cross-breeding remain to be investigated.

The impact of genetic diversity on the accuracy of DNA barcoding to identify species: A study on the genus Phellodendron

DNA barcoding is widely used in species identification, but there is considerable controversy regarding the extent of sampling in research methods. Some scholars have proposed that this small sample size underestimates the intraspecific genetic diversity, which would impact on the accuracy of DNA barcoding to identify species. In study, we selected all Phellodendron species (including P. amurense Rupr., P. chinense Schneid., and P. chinense var. glabriusculum Schneid.) as the materials, collected 59 P. amurense samples from 35 populations greatly to represent the genetic diversity, and analyzed the haplotype, genetic distance, barcoding gap, and Neighbor-Joining (NJ) trees based on psbA-trnH and internal transcribed spacer gene sequences. Additionally, a sampling simulation was conducted to assess the correlation between genetic diversity and the number of populations. Finally, analysis of critical geographical populations was performed. Based on analysis of haplotype, genetic distance, barcoding gap, and NJ trees, we found that eight P. amurense samples impacted on the effectiveness of DNA barcoding, which genetic information were very important to identify Phellodendron species. Moreover, the result of the NJ tree analysis performed the small-scale P. amurense sample size did not completely match the objective phylogenetic relationship in Phellodendron. In simulation sampling analysis, the data showed the genetic diversity indexes at the same population level gradually decreased and stabilized as the number of simulation sampling populations increased. We found that 1-2 samples from over 24 populations based on uniform geographical distribution could represent 80% of the genetic diversity of P. amurense and ensure authenticity and reliability of DNA barcoding. Thus, we proposed it is particularly important adequately samples to cover infraspecific genetic diversity in order to ensure identification accuracy of DNA barcoding.