Estimating sample sizes for DNA barcoding

Molecular identification and genetic variation of forensically important fly species (Order: Diptera) in Thailand using DNA barcoding

Forensic entomology plays a crucial role in criminal investigations by providing vital insights into minimum postmortem interval (PMImin) and corpse relocation by identifying insect species that colonize in decomposing remains. This study aimed to identify and analyze the genetic variation of forensically significant fly species in Thailand, using DNA barcoding of the mitochondrial cytochrome c oxidase subunit I COI gene. A total of 3,220 fly specimens were collected from 18 provinces across six regions of Thailand from October 2017 to September 2022. These specimens were classified by morphological identification into 21 species among three Dipteran families: Calliphoridae, Muscidae, and Sarcophagidae, with Chrysomya megacephala Diptera: Calliphoridae being the most abundant species. DNA barcoding confirmed the morphological identifications with 100 % accuracy, showing low intraspecific K2P distances0.0 to 1.1 %) and significant interspecific K2P distances 2.5 % to 17.2 %. A Neighbour-Joining (NJ) analysis was conducted to assess the molecular identification capabilities of the barcoding region. This analysis successfully recovered nearly all species as distinct monophyletic groups. The species groupings obtained were generally consistent with both morphological and molecular identifications. These findings underscore the effectiveness of DNA barcoding for precise species identification and contribute to a comprehensive database of forensically important flies in Thailand, thus facilitating improved forensic investigations and biodiversity studies.

Molecular evidence for a new endemic species of Acartia (Copepoda, Calanoida) from the Southeast Pacific coast

The loss of biodiversity in marine populations is one of the consequences of the increased events of extreme environmental conditions in the oceans, which can condition the persistence of populations to future scenarios of climate change. Therefore, it is extremely necessary to explore and monitor the genetic diversity of natural populations. In the Southeast Pacific Ocean (SEPO), specifically on the coast of Chile, the presence of the copepod Acartia tonsa has been indicated solely using morphological evidence, due to the absence of genetic information. In the present work, the genetic diversity, population structure and phylogenetic position within the genus Acartia, of populations identified morphologically as A. tonsa, was evaluated by amplification of the mitochondrial cytochrome c oxidase subunit I and nuclear marker 18 s. Our results showed that the populations identified as A. tonsa correspond to a new monophyletic group endemic to SEPO (GMYC = 1.00; PTP = 0.95). The populations showed moderate to high genetic diversity with an incipient structuring between populations and biogeographic zones. Our results suggest that despite the homogenizing effect of the Humboldt Current, isolation by distance and contrasting environmental conditions at different geographic scales have an important influence on the genetic diversity of zooplankton in the SEPO region.

Unique among high passes: Phylogenetic inferences from DNA barcoding of the butterfly fauna of Ladakh Trans-Himalaya, India

The butterfly assemblage of Ladakh Trans-Himalaya demands a thorough analysis of their population genetic structure owing to their typical biogeographic affinity and their adaptability to extreme cold-desert climates. No such effort has been taken till date, and in this backdrop, we created a barcode reference library of 60 specimens representing 23 species. Barcodes were generated from freshly collected leg samples using the Sanger sequencing method, followed by phylogenetic clade analyses and divergence calculation. Our data represents 22% of Ladakh's Rhopaloceran fauna with the novel barcode submission for six species, including one Schedule II species, Paralasa mani . Contrary to the 3% threshold rule, the interspecific divergence between two species pairs of typical mountain genus Hyponephele and Karanasa was found to be 2.3% and 2.2%, respectively. The addition of conspecific global barcodes revealed that most species showed little increase in divergence value, while a two-fold increase was noted in a few species. Bayesian clade clustering outcomes largely aligned with current morphological classifications, forming monophyletic clades of conspecific barcodes, with only minor exceptions observed for the taxonomically complicated genus Polyommatus and misidentified records of Aulocera in the database. We also observed variations within the same phylogenetic clades forming nested lineages, which may be attributed to the taxonomic intricacies present at the subspecies level globally, mostly among Eurasian species. Overall, our effort not only substantiated the effectiveness of DNA Barcoding for the identification and conservation of this climatically vulnerable assemblage but also highlighted the significance of deciphering the unique genetic composition among this geographically isolated population of Ladakh butterflies.

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.

DNA Barcoding of genus Culicoides biting midges (Diptera: Ceratopogonidae) in the Brazilian Amazon

Culicoides biting midges are capable to transmit Oropouche virus, Bluetongue virus and Mansonella spp. This study aimed to assess the utility of DNA barcode as an alternative method in the Culicoides species identification. The study was conducted in Jamari National Forest. Biting midges were collected using HP light traps during four months, February, April, August and October, 2018. Insects were morphologically identified to the species level, and rest of the body were subjected to DNA extraction and PCR targeting a fragment of the cytochrome c oxidase subunit I (COI) gene, which were analyzed and deposited in GenBank. A phylogenetic gene tree was reconstructed using RAxML software, and the sequences were assigned at Molecular Operational Taxonomic Unit (MOTU) level by species delimitation algorithms. According to morphological approach, 18 species of 2 subgenera and 7 species groups were identified. A total of 191 new COI barcodes from 18 species were generated. Of these, fifteen species have been deposited for the first time in all datasets in the world. These sequences allowed the correct identification of 188 and 187 specimens according to the BM and BCM criteria, respectively. The intraspecific genetic distances ranged from 0 to 16.5%, while the interspecific ones ranged from 2.1 to 27.1%. The nominal species Culicoides glabellus and C. tetrathyris splitted into three and two MOTUs, respectively, except for mPTP, indicating a cryptic diversity in these species. Also, sequences of C. pseudodiabolicus formed two MOTUs using all algorithms, except for PTP and ABGD, suggesting the existence of two potential species. In contrast, some barcodes of C. quasiparaensis and C. paraensis merged into a single MOTU, which can be explained by the complex characteristics of the paraensis group, since these species have similar morphological characters. Here, we provided the first COI barcodes for biting midges in Rondônia and Brazil, and demonstrates that these are sufficient to discriminate between some species.

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.

Diet Composition of the Wild Stump-Tailed Macaque (Macaca arctoides) in Perlis State Park, Peninsular Malaysia, Using a Chloroplast tRNL DNA Metabarcoding Approach: A Preliminary Study

Simple Summary

This study investigated plant diet of wild Macaca arctoides in the Malaysia–Thailand border region using a chloroplast tRNL DNA metabarcoding approach. It is a comprehensive molecular technique to assess foods eaten by primates. We have chosen chloroplast tRNL because this region has been widely used for identifying plant species. Chloroplast tRNL DNA was amplified and sequenced using the Illumina MiniSeq platform. Sequences were analyzed using the CLC Genomic Workbench software version 12.0 to check for M. arctoides plant diet. Across these samples, we successfully identified 29 plant orders, 46 families, 124 genera, and 145 species. As the first report in Malaysia, the findings provide an important understanding on diet of wild M. arctoides that only reside in Perlis State Park, Malaysia.

Abstract

Understanding dietary diversity is a fundamental task in the study of stump-tailed macaque, Macaca arctoides in its natural habitat. However, direct feeding observation and morphological identification using fecal samples are not effective and nearly impossible to obtain in natural habitats because this species is sensitive to human presence. As ecological methods are challenging and time-consuming, DNA metabarcoding offers a more powerful assessment of the diet. We used a chloroplast tRNL DNA metabarcoding approach to identify the diversity of plants consumed by free-ranging M. arctoides in the Malaysia–Thailand border region located in Perlis State Park, Peninsular Malaysia. DNA was extracted from three fecal samples, and chloroplast tRNL DNA was amplified and sequenced using the Illumina MiniSeq platform. Sequences were analyzed using the CLC Genomic Workbench software. A total of 145 plant species from 46 families were successfully identified as being consumed by M. arctoides. The most abundant species were yellow saraca, Saraca thaipingensis (11.70%), common fig, Ficus carica (9.33%), aramata, Clathrotropis brachypetala (5.90%), sea fig, Ficus superba (5.44%), and envireira, Malmea dielsiana (1.70%). However, Clathrotropis and Malmea are not considered Malaysian trees because of limited data available from Malaysian plant DNA. Our study is the first to identify plant taxa up to the species level consumed by stump-tailed macaques based on a DNA metabarcoding approach. This result provides an important understanding on diet of wild M. arctoides that only reside in Perlis State Park, Malaysia.

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.

Plant DNA barcoding necessitates marker-specific efforts to establish more comprehensive reference databases

The problem of low species-level identification rates in plants by DNA barcoding is exacerbated by the fact that reference databases are far from being comprehensive. We investigate the impact of increased sampling depth on identification success by analyzing the efficacy of established plant barcode marker sequences (rbcL, matK, trnL-trnF, psbA-trnH, ITS). Adding sequences of the same species to the reference database led to an increase in correct species assignment of +10.9% for rbcL and +19.0% for ITS. Simultaneously, erroneous identification dropped from ∼40% to ∼12.5%. Despite its evolutionary constraints, ITS showed the highest identification rate and identification gain by increased sampling effort, which makes it a very suitable marker in the planning phase of a barcode study. The limited sequence availability of trnL-trnF is problematic for an otherwise very promising plastid plant barcoding marker. Future developments in machine learning algorithms have the potential to give new impetus to plant barcoding, but are dependent on extensive reference databases. We expect that our results will be incorporated into future plans for the development of DNA barcoding reference databases and will lead to these being developed with greater depth and taxonomic coverage.

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.

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.

On the applicability of the Tubulin-Based Polymorphism (TBP) genotyping method: a comprehensive guide illustrated through the application on different genetic resources in the legume family

BACKGROUND

Plant discrimination is of relevance for taxonomic, evolutionary, breeding and nutritional studies. To this purpose, evidence is reported to demonstrate TBP (Tubulin-Based-Polymorphism) as a DNA-based method suitable for assessing plant diversity.

RESULTS

Exploiting one of the most valuable features of TBP, that is the convenient and immediate application of the assay to groups of individuals that may belong to different taxa, we show that the TBP method can successfully discriminate different agricultural species and their crop wild relatives within the Papilionoideae subfamily. Detection of intraspecific variability is demonstrated by the genotyping of 27 different accessions of Phaseolus vulgaris.

CONCLUSIONS

These data illustrate TBP as a useful and versatile tool for plant genotyping. Since its potential has not yet been fully appreciated by the scientific community, we carefully report all the experimental details of a successful TBP protocol, while describing different applications, so that the method can be replicated in other laboratories.

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.

New Data on the Range Expansion of the Thaumetopoea pityocampa (Lepidoptera: Notodontidae) 'ENA clade' in Greece: The Role of Bacterial Endosymbionts

The pine processionary moth, Thaumetopoea pityocampa (Denis and Schiffermüller), is an important insect in the Mediterranean region, as it defoliates pines and its urticating hairs can cause allergic reactions in humans and animals. Moreover, this species exhibits an interesting genetic structure as recently a distinct East-North African mtDNA lineage ('ENA clade') has been described. This clade has been recently detected in Greek populations where it has currently expanded its range by replacing the 'endemic' T. pityocampa lineages. Here, we report new data on the rapid spread of 'ENA clade' in the Greek island Evoia in only a few years. As the underlying mechanisms of the 'ENA clade' range expansion has not been studied so far, we screened T. pityocampa for an infection with the heritable bacterial endosymbionts Wolbachia (Bacteria: Anaplasmataceae), Cardinium (Bacteria: Bacteroidaceae), Rickettsia (Bacteria: Rickettsiaceae) and Spiroplasma (Bacteria: Spiroplasmataceae). These bacteria can manipulate the reproduction of infected hosts, something that could potentially explain the rapid spread of 'ENA clade' lineage. Therefore, we screened 28 individuals that exhibited T. pityocampa 'ENA clade' and 'endemic' T. pityocampa haplotypes from nine populations scattered all over Greece. None of them was infected with any of the four endosymbionts, suggesting that these bacteria do not cause reproductive manipulations in T. pityocampa lineages and, thus, other factors should be explored in future research efforts.

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.

The Prevalence of Single-Specimen/Locality Species in Insect Taxonomy: An Empirical Analysis

Although about two million species have been named, our knowledge about the biodiversity of many taxonomic groups remains inadequate and incomplete. There has been increased taxonomic effort for the discovery of more species and their geographical distribution information. During this process, species collected only from a single specimen/locality often appear. However, there are very few empirical data available to understand the current situation of single specimen/locality species in insect taxonomy. In this paper, we collected 1261 articles containing 4811 insect species from ZooKeys between 2009 and 2017, and we extracted data, including the publication date, number of specimens/locality, and DNA usage. Our analyses demonstrated that 21.53% and 21.74% of new species were described from only one specimen and one locality, respectively, and approximately half of all new species were published based on fewer than five specimens. Meanwhile, the rate of single-specimen species in papers with or without DNA data was 15.06% and 23.43%, respectively, which indicates that incorporating DNA data in species descriptions might effectively decrease the occurrence of single-specimen species. We suggest that taxonomists should adopt more beneficial practices, such as increasing specimen diversity, incorporating DNA data, and improving international collaboration, in the description of new species.

Incomplete estimates of genetic diversity within species: Implications for DNA barcoding

DNA barcoding has greatly accelerated the pace of specimen identification to the species level, as well as species delineation. Whereas the application of DNA barcoding to the matching of unknown specimens to known species is straightforward, its use for species delimitation is more controversial, as species discovery hinges critically on present levels of haplotype diversity, as well as patterning of standing genetic variation that exists within and between species. Typical sample sizes for molecular biodiversity assessment using DNA barcodes range from 5 to 10 individuals per species. However, required levels that are necessary to fully gauge haplotype variation at the species level are presumed to be strongly taxon-specific. Importantly, little attention has been paid to determining appropriate specimen sample sizes that are necessary to reveal the majority of intraspecific haplotype variation within any one species. In this paper, we present a brief outline of the current literature and methods on intraspecific sample size estimation for the assessment of COI DNA barcode haplotype sampling completeness. The importance of adequate sample sizes for studies of molecular biodiversity is stressed, with application to a variety of metazoan taxa, through reviewing foundational statistical and population genetic models, with specific application to ray-finned fishes (Chordata: Actinopterygii). Finally, promising avenues for further research in this area are highlighted.

Species identification of horse flies (Diptera: Tabanidae) in Thailand using DNA barcoding

Horse flies (Diptera: Tabanidae) are of medical and veterinary importance because they are known to transmit pathogens. Approximately 80 species of horse flies have been reported in Thailand. Monitoring the distribution of horse fly species is important to control the spread of diseases transmitted by them. Currently, the species identification of horse flies is based on their morphology; this requires considerable skills and taxonomic expertise, and it may be difficult to identify morphologically similar species. DNA-based identification methods are increasingly being developed for rapid and accurate identification of various insect species. In this study, we used mitochondrial cytochrome oxidase subunit I (COI) for species identification of horse flies in Thailand. A 658 bp fragment of COI was amplified from 145 adult horse flies belonging to 48 morphologically distinct species and sequenced. Sequence analysis revealed an intraspecific divergence of 0.0%-4.4% and an interspecific divergence of 0.0%-16.2%. Our results showed that COI barcodes were effective in discriminating a majority of horse flies in Thailand on the basis of the barcoding gap and phylogenetic analyses. However, COI barcodes were unable to distinguish among members of the Tabanus striatus complex and some species within the T. ceylonicus group.

Species identification and connectivity of marine amphipods in Canada's three oceans

Monitoring the distribution of marine biodiversity is a crucial step to better assess the impacts of global changes. Arctic marine fauna is dominated by amphipods in terms of abundance and biomass. These peracarids are an important marine order of crustaceans but the number of species found in the different Canadian oceans is currently unknown. Furthermore, most species are difficult to identify due to poor taxonomic descriptions and morphological convergence. The aim of this study was to assess the species diversity of marine amphipods in the three Canadian oceans using DNA barcoding. To do so, we produced a database of DNA barcodes of amphipods from the three Canadian Oceans publicly available from the BOLD website to which we added 310 new sequences from the Canadian Arctic Archipelago. We first delimited amphipod species based on barcode gap detection techniques and tree based method (bPTP) and then compared the composition of amphipods among the three oceans in order to assess the influence of past transarctic exchanges on Arctic diversity. Our analysis of 2309 sequences which represent more than 250 provisional species revealed a high connectivity between the Atlantic and Arctic Oceans. Our results also suggest that a single threshold to delimitate species is not suitable for amphipods. This study highlights the challenges involved in species delimitation and the need to obtain complete barcoding inventories in marine invertebrates.

Integrating a comprehensive DNA barcode reference library with a global map of yews (Taxus L.) for forensic identification

Rapid and accurate identification of endangered species is a critical component of biosurveillance and conservation management, and potentially policing illegal trades. However, this is often not possible using traditional taxonomy, especially where only small or preprocessed parts of plants are available. Reliable identification can be achieved via a comprehensive DNA barcode reference library, accompanied by precise distribution data. However, these require extensive sampling at spatial and taxonomic scales, which has rarely been achieved for cosmopolitan taxa. Here, we construct a comprehensive DNA barcode reference library and generate distribution maps using species distribution modelling (SDM), for all 15 Taxus species worldwide. We find that trnL-trnF is the ideal barcode for Taxus: It can distinguish all Taxus species and in combination with ITS identify hybrids. Among five analysis methods tested, NJ was the most effective. Among 4,151 individuals screened for trnL-trnF, 73 haplotypes were detected, all species-specific and some population private. Taxonomical, geographical and genetic dimensions of sampling strategy were all found to affect the comprehensiveness of the resulting DNA barcode library. Maps from SDM showed that most species had allopatric distributions, except T. mairei in the Sino-Himalayan region. Using the barcode library and distribution map data, two unknown forensic samples were identified to species (and in one case, population) level and another was determined as a putative interspecific hybrid. This integrated species identification system for Taxus can be used for biosurveillance, conservation management and to monitor and prosecute illegal trade. Similar identification systems are recommended for other IUCN- and CITES-listed taxa.

Molecular phylogeny of Anopheles hyrcanus group (Diptera: Culicidae) based on mtDNA COI

Background

The Anopheles hyrcanus group, which includes at least 25 species, is widely distributed in the Oriental and Palearctic regions. Some group members have been incriminated as vectors of malaria and other mosquito-borne diseases. It is difficult to identify Hyrcanus Group members by morphological features. Thus, molecular phylogeny has been proposed as an important complementary method to traditional morphological taxonomy.

Methods

Based on the GenBank database and our original study data, we used 466 mitochondrial DNA COI sequences belonging to 18 species to reconstruct the molecular phylogeny of the Hyrcanus Group across its worldwide geographic range.

Results

The results are as follows. 1) The average conspecific K2P divergence was 0.008 (range 0.002–0.017), whereas sequence divergence between congroup species averaged 0.064 (range 0.026–0.108). 2) The topology of COI tree of the Hyrcanus Group was generally consistent with classical morphological taxonomy in terms of species classification, but disagreed in subgroup division. In the COI tree, the group was divided into at least three main clusters. The first cluster contained An. nimpe; the second was composed of the Nigerrimus Subgroup and An. argyropus; and the third cluster was comprised of the Lesteri Subgroup and other unassociated species. 3) Phylogenetic analysis of COI indicated that ancient hybridizations probably occurred among the three closely related species, An. sinensis, An. belenrae, and An. kleini. 4) The results supported An. paraliae as a probable synonym of An. lesteri, and it was possible that An. pseudopictus and An. hyrcanus were the same species, as evident from their extremely low interspecific genetic divergence (0.020 and 0.007, respectively) and their phylogenetic positions.

Conclusions

In summary, we reconstructed the molecular phylogeny and analysed genetic divergence of the Hyrcanus Group using mitochondrial COI sequences. Our results suggest that in the future of malaria surveillance, we should not only pay much attention to those known vectors of malaria, but also their closely related species.

Electronic supplementary material

The online version of this article (doi:10.1186/s40249-017-0273-7) contains supplementary material, which is available to authorized users.

Replacing Sanger with Next Generation Sequencing to improve coverage and quality of reference DNA barcodes for plants

We estimate the global BOLD Systems database holds core DNA barcodes (rbcL + matK) for about 15% of land plant species and that comprehensive species coverage is still many decades away. Interim performance of the resource is compromised by variable sequence overlap and modest information content within each barcode. Our model predicts that the proportion of species-unique barcodes reduces as the database grows and that 'false' species-unique barcodes remain >5% until the database is almost complete. We conclude the current rbcL + matK barcode is unfit for purpose. Genome skimming and supplementary barcodes could improve diagnostic power but would slow new barcode acquisition. We therefore present two novel Next Generation Sequencing protocols (with freeware) capable of accurate, massively parallel de novo assembly of high quality DNA barcodes of >1400 bp. We explore how these capabilities could enhance species diagnosis in the coming decades.

A DNA Barcode Library for North American Pyraustinae (Lepidoptera: Pyraloidea: Crambidae)

Although members of the crambid subfamily Pyraustinae are frequently important crop pests, their identification is often difficult because many species lack conspicuous diagnostic morphological characters. DNA barcoding employs sequence diversity in a short standardized gene region to facilitate specimen identifications and species discovery. This study provides a DNA barcode reference library for North American pyraustines based upon the analysis of 1589 sequences recovered from 137 nominal species, 87% of the fauna. Data from 125 species were barcode compliant (>500bp, <1% n), and 99 of these taxa formed a distinct cluster that was assigned to a single BIN. The other 26 species were assigned to 56 BINs, reflecting frequent cases of deep intraspecific sequence divergence and a few instances of barcode sharing, creating a total of 155 BINs. Two systems for OTU designation, ABGD and BIN, were examined to check the correspondence between current taxonomy and sequence clusters. The BIN system performed better than ABGD in delimiting closely related species, while OTU counts with ABGD were influenced by the value employed for relative gap width. Different species with low or no interspecific divergence may represent cases of unrecognized synonymy, whereas those with high intraspecific divergence require further taxonomic scrutiny as they may involve cryptic diversity. The barcode library developed in this study will also help to advance understanding of relationships among species of Pyraustinae.

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

The use of genetic distances to identify species within the framework of DNA barcoding has to some extent improved the development of biodiversity studies. However, using a fixed empirical threshold to delimit species may lead to overestimating species diversity. In this study, we use a new data set of COI sequences for 366 specimens within the genus of Cletus as well as conduct an analysis on the same genetic data for collected morphologically defined species from previous phylogeographical studies, to test whether high intraspecific genetic divergences are common with the premises of comprehensive sampling. The results indicate C. graminis Hsiao & Cheng , is the same species with C. punctiger (Dallas, 1852) and should be synonymized and that the distributional record of C. pugnator (Fabricius, 1787) in China is correct. High intraspecific genetic differentiations (0%-4.35%) were found in C. punctiger. Furthermore, as to the mined data, the maximum intraspecific K2P distances of 186 species (48.44% of 384) exceed 3%, and 101 species (26.30%) can be divided into two or more clusters with a threshold of 3% in cluster analysis. If genetic distance is used to delimit species boundaries, the minimum interspecific K2P distance of the congeneric species should be considered rather than only using the fixed empirical value; otherwise, the species richness may be overestimated in some cases.

Phylogenetics and population genetics of Plotosus canius (Siluriformes: Plotosidae) from Malaysian coastal waters

Plotosus canius (Hamilton, 1822) is a significant marine species in Malaysia from nutritional and commercial perspectives. Despite numerous fundamental research on biological characteristics of P. canius, there are various concerns on the level of population differentiation, genomic structure, and the level of genetic variability among their populations due to deficiency of genetic-based studies. Deficiency on basic contexts such as stock identification, phylogenetic relationship and population genetic structure would negatively impact their sustainable conservation. Hence, this study was conducted to characterize the genetic structure of P. canius for the first time through the application of mitochondrial Cytochrome Oxidase I (COI) gene, cross amplification of Tandanus tandanus microsatellites, and a total of 117 collected specimens across five selected populations of Malaysia. The experimental results of the mitochondrial analysis revealed that the haplotype diversity and nucleotide diversity varied from 0.395–0.771 and 0.033–0.65 respectively. Moreover, the statistical analysis of microsatellites addressed a considerable heterozygote insufficiency in all populations, with average observed heterozygosity (H_o) value of 0.2168, which was lower than the standard heterozygosity in marine populations (H_o = 0.79). This alongside the high Fis values estimation, high pairwise differentiation among populations and low within population variations are supposed to be associated with small sample size, and inbreeding system. Besides, the significant finding of this study was the sharing of common haplotype KR086940, which reflects a historical genetic connectivity between Peninsular Malaysia and Borneo populations due to the geological history of Southeast Asia during Pleistocene era. Demographic analyses showed that all populations were in an equilibrium state with no significant evidence of population expansion. To put it briefly, the current study has managed to provide an initial genomic database toward understanding of the genetic characterization, phylogenetic, molecular diversification and population structure in P. canius, and should be necessary highlighted for appropriate management and conservation of species. Further studies must be carried out involving more geographical and sampling sites, larger population size per site, and utilization of species specific microsatellites loci.

Evolutionary histories determine DNA barcoding success in vascular plants: seven case studies using intraspecific broad sampling of closely related species

Background

Four plastid regions, rpoB, rpoC1, matK, and trnH-psbA, have been recommended as DNA barcodes for plants. Their success in delimiting species boundaries depends on the existence of a clear-cut difference between inter- and intraspecific variability. We tested the ability of these regions to discriminate among closely related species in seven genera of flowering plants with different generation times (trees, perennials, and annuals). To ensure a maximum coverage of intraspecific diversity, and therefore to better evaluate the resolution power of each barcode, we applied a population genetics approach by sampling three to 45 individuals per species over a wide geographical range.

Results

All possible combinations between loci were analysed, which showed that using more than one locus does not always improve the resolution power. The trnH-psbA locus was most effective at discriminating among closely related species (Acer, Lonicera, Geranium, and Veronica), singly or in combination. For Salix, Adenostyles, and Gentiana, the best results were obtained with the combination of matK, rpoB, and trnH-psbA. No barcoding gap was found within six genera analysed, excepting Lonicera. This is due to shared polymorphisms among species, combined with very divergent sequences within species. These genetic patterns reflect incomplete lineage sorting and hybridization events followed by chloroplast capture.

Conclusions

Our results strongly suggest that adding trnH-psbA to the two obligate DNA barcodes proposed by the CBOL plant-working group (matK and rbcL) should be mandatory for closely related species. In our sampling, generation time had no influence on DNA barcoding success, as the best and worst identification successes were found for the two tree genera (Acer, 64 % success and Salix, 86 % failure). Evolutionary histories are the main factor influencing DNA barcoding success in the studied genera.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0678-0) contains supplementary material, which is available to authorized users.

Molecular systematics of the freshwater stingrays (myliobatiformes: potamotrygonidae) of the Amazon, Orinoco, Magdalena, Esequibo, Caribbean, and Maracaibo basins (Colombia - Venezuela): evidence from three mitochondrial genes

Lack of adequate information about the taxonomic and evolutionary relationships, ecology, biology, and distribution of several species belonging to the family Potamotrygonidae makes these species vulnerable to anthropic activities, including commercial overexploitation for the ornamental fish market. The aim of this study was to investigate the systematic relationships among genera and species belonging to this family by analyses of three mitochondrial gene regions. Samples were collected from the main river basins in Colombia and Venezuela for four genera and seven species of the family, as well as for what appear to be unidentified species. Three mitochondrial molecular markers COI, Cytb, and ATP6 were amplified and sequenced. Maximum likelihood and Bayesian inference analysis were performed to obtain topologies for each marker and for a concatenated dataset including the three genes. Small dataset may compromise some methods estimations of sequence divergence in the ATP6 marker. Monophyly of the four genera in Potamotrygonidae was confirmed and phylogenetic relationships among members of the Potamotrygon genus were not clearly resolved. However, results obtained with the molecular marker Cytb appear to offer a good starting point to differentiate among genera and species as a tool that could be used for barcoding. The application of this gene as a barcode could be applied for management and regulation of extraction practices for these genera. Sequencing complete mitochondrial genomes would be the next step for testing evolutionary hypothesis among these genera. Population structure analyses should be undertaken for Paratrygon, Potamotrygon magdalenae and motoro.

A simulation study of sample size for DNA barcoding

For some groups of organisms, DNA barcoding can provide a useful tool in taxonomy, evolutionary biology, and biodiversity assessment. However, the efficacy of DNA barcoding depends on the degree of sampling per species, because a large enough sample size is needed to provide a reliable estimate of genetic polymorphism and for delimiting species. We used a simulation approach to examine the effects of sample size on four estimators of genetic polymorphism related to DNA barcoding: mismatch distribution, nucleotide diversity, the number of haplotypes, and maximum pairwise distance. Our results showed that mismatch distributions derived from subsamples of ≥20 individuals usually bore a close resemblance to that of the full dataset. Estimates of nucleotide diversity from subsamples of ≥20 individuals tended to be bell‐shaped around that of the full dataset, whereas estimates from smaller subsamples were not. As expected, greater sampling generally led to an increase in the number of haplotypes. We also found that subsamples of ≥20 individuals allowed a good estimate of the maximum pairwise distance of the full dataset, while smaller ones were associated with a high probability of underestimation. Overall, our study confirms the expectation that larger samples are beneficial for the efficacy of DNA barcoding and suggests that a minimum sample size of 20 individuals is needed in practice for each population.

A molecular framework for the identification of planthopper vectors (Hemiptera: Delphacidae) of central Argentina

Planthoppers are important worldwide crop pests as well as vectors of numerous diseases. Different species transmit Mal de Río Cuarto virus, which causes the most economically important corn disease in central Argentina. Epidemiological studies rely on the accurate identification of the species present in the field. Presently, morphological identification of planthoppers requires taxonomic expertise and there are no taxonomic keys for females and nymphs. Nevertheless, no molecular protocols are available for accurate species identification of most frequent delphacid species from central Argentina. In this context, the aim of this study was to evaluate the utility of the cytochrome oxidase I gene (COI) as a DNA barcode and its digestion with restriction enzymes (Restriction Fragment Length Polymorphism, RFLP) for the identification of the most common species of planthoppers in central Argentina. We amplified and sequenced a 843 bp fragment of the COI gene of taxonomically identified specimens and evaluated its use as a DNA barcode. Restriction enzymes were also selected for digesting the COI fragment via RFLP. The high interspecific variability (20.79%; ± 2.32%) and low intraspecific divergence (0.12%; ± 0.17%) observed in the studied species, demonstrate the effectiveness of the COI gene for species identification of major vector delphacids affecting corn crops in Argentina. Moreover, the digestion of this COI gene fragment with Bfa I and Apo I enzymes allows a fast and cost-effective species identification method when numerous specimens need to be processed. Both molecular techniques developed here, allow the accurate identification of planthopper species at regional scale. These new tools would assist traditional identification of these insects, especially for aiding non-experts in morphological taxonomy.

Culicoides Latreille (Diptera: Ceratopogonidae) taxonomy: current challenges and future directions

Culicoides Latreille biting midges (Diptera: Ceratopogonidae) cause a significant biting nuisance to humans, livestock and equines, and are the biological vectors of a range of internationally important pathogens of both veterinary and medical importance. Despite their economic significance, the delimitation and identification of species and evolutionary relationships between species within this genus remains at best problematic. To date no phylogenetic study has attempted to validate the subgeneric classification of the genus and the monophyly of many of the subgenera remains doubtful. Many informal species groupings are also known to exist but few are adequately described, further complicating accurate identification. Recent contributions to Culicoides taxonomy at the species level have revealed a high correlation between morphological and molecular analyses although molecular analyses are revealing the existence of cryptic species. This review considers the methods for studying the systematics of Culicoides using both morphological and genetic techniques, with a view to understanding the factors limiting our current understanding of Culicoides biology and hence arbovirus epidemiology. In addition, we examine the global status of Culicoides identification, highlighting areas that are poorly addressed, including the potential implementation of emerging technologies.

An exploration of sufficient sampling effort to describe intraspecific DNA barcode haplotype diversity: examples from the ray-finned fishes (Chordata: Actinopterygii)

Estimating appropriate sample sizes to measure species abundance and richness is a fundamental problem for most biodiversity research. In this study, we explore a method to measure sampling sufficiency based on haplotype diversity in the ray-finned fishes (Animalia: Chordata: Actinopterygii). To do this, we use linear regression and hypothesis testing methods on haplotype accumulation curves from DNA barcodes for 18 species of fishes, in the statistics platform R. We use a simple mathematical model to estimate sampling sufficiency from a sample-number based prediction of intraspecific haplotype diversity, given an assumption of equal haplotype frequencies. Our model finds that haplotype diversity for most of the 18 fish species remains largely unsampled, and this appears to be a result of small sample sizes. Lastly, we discuss how our overly simple model may be a useful starting point to develop future estimators for intraspecific sampling sufficiency in studies using DNA barcodes.

DNA barcoding and the associated PhylAphidB@se website for the identification of European aphids (Insecta: Hemiptera: Aphididae)

Aphids constitute a diverse group of plant-feeding insects and are among the most important crop pests in temperate regions. Their morphological identification is time-consuming and requires specific knowledge, training and skills that may take years to acquire. We assessed the advantages and limits of DNA barcoding with the standard COI barcode fragment for the identification of European aphids. We constructed a large reference dataset of barcodes from 1020 specimens belonging to 274 species and 87 genera sampled throughout Europe and set up a database-driven website allowing species identification from query sequences.

RESULTS

In this unbiased sampling of the taxonomic diversity of European aphids, intraspecific divergence ranged from 0.0% to 3.9%, with a mean value of 0.29%, whereas mean congeneric divergence was 6.4%, ranging from 0.0% to 15%. Neighbor-joining analysis generated a tree in which most species clustered in distinct genetic units. Most of the species with undifferentiated or overlapping barcodes belonged to the genus Aphis or, to a lesser extent, the genera Brachycaudus, Dysaphis and Macrosiphum. The taxa involved were always morphologically similar or closely related and belonged to species groups known to present taxonomic difficulties.

CONCLUSIONS

These data confirm that COI barcoding is a useful identification tool for aphids. Barcode identification is straightforward and reliable for 80% of species, including some difficult to distinguish on the basis of morphological characters alone. Unsurprisingly, barcodes often failed to distinguish between species from groups for which classical taxonomy has also reached its limits, leading to endless revisions and discussions about species and subspecies definitions. In such cases, the development of an effective procedure for the accurate identification of aphid specimens continues to pose a difficult challenge.

Identifying species of moths (Lepidoptera) from Baihua Mountain, Beijing, China, using DNA barcodes

DNA barcoding has become a promising means for the identification of organisms of all life-history stages. Currently, distance-based and tree-based methods are most widely used to define species boundaries and uncover cryptic species. However, there is no universal threshold of genetic distance values that can be used to distinguish taxonomic groups. Alternatively, DNA barcoding can deploy a "character-based" method, whereby species are identified through the discrete nucleotide substitutions. Our research focuses on the delimitation of moth species using DNA-barcoding methods. We analyzed 393 Lepidopteran specimens belonging to 80 morphologically recognized species with a standard cytochrome c oxidase subunit I (COI) sequencing approach, and deployed tree-based, distance-based, and diagnostic character-based methods to identify the taxa. The tree-based method divided the 393 specimens into 79 taxa (species), and the distance-based method divided them into 84 taxa (species). Although the diagnostic character-based method found only 39 so-identifiable species in the 80 species, with a reduction in sample size the accuracy rate substantially improved. For example, in the Arctiidae subset, all 12 species had diagnostics characteristics. Compared with traditional morphological method, molecular taxonomy performed well. All three methods enable the rapid delimitation of species, although they have different characteristics and different strengths. The tree-based and distance-based methods can be used for accurate species identification and biodiversity studies in large data sets, while the character-based method performs well in small data sets and can also be used as the foundation of species-specific biochips.

DNA barcoding for the identification of sand fly species (Diptera, Psychodidae, Phlebotominae) in Colombia

Sand flies include a group of insects that are of medical importance and that vary in geographic distribution, ecology, and pathogen transmission. Approximately 163 species of sand flies have been reported in Colombia. Surveillance of the presence of sand fly species and the actualization of species distribution are important for predicting risks for and monitoring the expansion of diseases which sand flies can transmit. Currently, the identification of phlebotomine sand flies is based on morphological characters. However, morphological identification requires considerable skills and taxonomic expertise. In addition, significant morphological similarity between some species, especially among females, may cause difficulties during the identification process. DNA-based approaches have become increasingly useful and promising tools for estimating sand fly diversity and for ensuring the rapid and accurate identification of species. A partial sequence of the mitochondrial cytochrome oxidase gene subunit I (COI) is currently being used to differentiate species in different animal taxa, including insects, and it is referred as a barcoding sequence. The present study explored the utility of the DNA barcode approach for the identification of phlebotomine sand flies in Colombia. We sequenced 700 bp of the COI gene from 36 species collected from different geographic localities. The COI barcode sequence divergence within a single species was <2% in most cases, whereas this divergence ranged from 9% to 26.6% among different species. These results indicated that the barcoding gene correctly discriminated among the previously morphologically identified species with an efficacy of nearly 100%. Analyses of the generated sequences indicated that the observed species groupings were consistent with the morphological identifications. In conclusion, the barcoding gene was useful for species discrimination in sand flies from Colombia.

Testing mitochondrial marker efficacy for DNA barcoding in spiders: a test case using the dwarf spider genus Oedothorax (Araneae : Linyphiidae : Erigoninae)

The present study focusses on comparatively assessing the efficacy for DNA barcoding of the two most commonly used mitochondrial markers (cox1 and 16S) in a genus of erigonine spiders. In total, 53 specimens representing five species, including four multi-sampled species, were sampled from several European localities. Initial evaluation of species monophyly was performed through parsimony and Bayesian phylogenetic analyses. Efficacy of mitochondrial markers was tested using operational (including distance-, tree-based measures and Barcode Gap) and evolutionary criteria (using the General Mixed Yule-coalescent Model) for species delimitation. We propose that the cox1 marker can potentially overestimate analyses of biodiversity and thus might not be the preferred marker for DNA species identification and delimitation methods in Oedothorax. Instead, our results suggest that the 16S marker appears to be a promising candidate for such endeavour. Evaluating the contribution and suitability of markers to the re-identification of species, measured by their recovery of well established morphological species, is critical for future studies and for reliable results in species identification in spiders.

DNA barcoding and species boundary delimitation of selected species of Chinese Acridoidea (Orthoptera: Caelifera)

We tested the performance of DNA barcoding in Acridoidea and attempted to solve species boundary delimitation problems in selected groups using COI barcodes. Three analysis methods were applied to reconstruct the phylogeny. K2P distances were used to assess the overlap range between intraspecific variation and interspecific divergence. "Best match (BM)", "best close match (BCM)", "all species barcodes (ASB)" and "back-propagation neural networks (BP-based method)" were utilized to test the success rate of species identification. Phylogenetic species concept and network analysis were employed to delimitate the species boundary in eight selected species groups. The results demonstrated that the COI barcode region performed better in phylogenetic reconstruction at genus and species levels than at higher-levels, but showed a little improvement in resolving the higher-level relationships when the third base data or both first and third base data were excluded. Most overlaps and incorrect identifications may be due to imperfect taxonomy, indicating the critical role of taxonomic revision in DNA barcoding study. Species boundary delimitation confirmed the presence of oversplitting in six species groups and suggested that each group should be treated as a single species.

A multi-locus approach to barcoding in the Anopheles strodei subgroup (Diptera: Culicidae)

Background

The ability to successfully identify and incriminate pathogen vectors is fundamental to effective pathogen control and management. This task is confounded by the existence of cryptic species complexes. Molecular markers can offer a highly effective means of species identification in such complexes and are routinely employed in the study of medical entomology. Here we evaluate a multi-locus system for the identification of potential malaria vectors in the Anopheles strodei subgroup.

Methods

Larvae, pupae and adult mosquitoes (n = 61) from the An. strodei subgroup were collected from 21 localities in nine Brazilian states and sequenced for the COI, ITS2 and white gene. A Bayesian phylogenetic approach was used to describe the relationships in the Strodei Subgroup and the utility of COI and ITS2 barcodes was assessed using the neighbor joining tree and “best close match” approaches.

Results

Bayesian phylogenetic analysis of the COI, ITS2 and white gene found support for seven clades in the An. strodei subgroup. The COI and ITS2 barcodes were individually unsuccessful at resolving and identifying some species in the Subgroup. The COI barcode failed to resolve An. albertoi and An. strodei but successfully identified approximately 92% of all species queries, while the ITS2 barcode failed to resolve An. arthuri and successfully identified approximately 60% of all species queries. A multi-locus COI-ITS2 barcode, however, resolved all species in a neighbor joining tree and successfully identified all species queries using the “best close match” approach.

Conclusions

Our study corroborates the existence of An. albertoi, An. CP Form and An. strodei in the An. strodei subgroup and identifies four species under An. arthuri informally named A-D herein. The use of a multi-locus barcode is proposed for species identification, which has potentially important utility for vector incrimination. Individuals previously found naturally infected with Plasmodium vivax in the southern Amazon basin and reported as An. strodei are likely to have been from An. arthuri C identified in this study.

Reexamination of the species assignment of Diacavolinia pteropods using DNA barcoding

Thecosome pteropods (Mollusca, Gastropoda) are an ecologically important, diverse, and ubiquitous group of holoplanktonic animals that are the focus of intense research interest due to their external aragonite shell and vulnerability to ocean acidification. Characterizing the response of these animals to low pH and other environmental stressors has been hampered by continued uncertainty in their taxonomic identification. An example of this confusion in species assignment is found in the genus Diacavolinia. All members of this genus were originally indentified as a single species, Cavolinia longirostris, but over the past fifty years the taxonomy has been revisited multiple times; currently the genus comprises 22 different species. This study examines five species of Diacavolinia, including four sampled in the Northeast Atlantic (78 individuals) and one from the Eastern tropical North Pacific (15 individuals). Diacavolina were identified to species based on morphological characteristics according to the current taxonomy, photographed, and then used to determine the sequence of the "DNA barcoding" region of the cytochrome c oxidase subunit I (COI). Specimens from the Atlantic, despite distinct differences in shell morphology, showed polyphyly and a genetic divergence of <3% (K2P distance) whereas the Pacific and Atlantic samples were more distant (≈ 19%). Comparisons of Diacavolinia spp. with other Cavolinia spp. reveal larger distances (≈ 24%). These results indicate that specimens from the Atlantic comprise a single monophyletic species and suggest possible species-level divergence between Atlantic and Pacific populations. The findings support the maintenance of Diacavolinia as a separate genus, yet emphasize the inadequacy of our current taxonomic understanding of pteropods. They highlight the need for accurate species identifications to support estimates of biodiversity, range extent and natural exposure of these planktonic calcifiers to environmental variability; furthermore, the apparent variation of the pteropods shell may have implications for our understanding of the species' sensitivity to ocean acidification.

The effect of geographical scale of sampling on DNA barcoding

Eight years after DNA barcoding was formally proposed on a large scale, CO1 sequences are rapidly accumulating from around the world. While studies to date have mostly targeted local or regional species assemblages, the recent launch of the global iBOL project (International Barcode of Life), highlights the need to understand the effects of geographical scale on Barcoding's goals. Sampling has been central in the debate on DNA Barcoding, but the effect of the geographical scale of sampling has not yet been thoroughly and explicitly tested with empirical data. Here, we present a CO1 data set of aquatic predaceous diving beetles of the tribe Agabini, sampled throughout Europe, and use it to investigate how the geographic scale of sampling affects 1) the estimated intraspecific variation of species, 2) the genetic distance to the most closely related heterospecific, 3) the ratio of intraspecific and interspecific variation, 4) the frequency of taxonomically recognized species found to be monophyletic, and 5) query identification performance based on 6 different species assignment methods. Intraspecific variation was significantly correlated with the geographical scale of sampling (R-square = 0.7), and more than half of the species with 10 or more sampled individuals (N = 29) showed higher intraspecific variation than 1% sequence divergence. In contrast, the distance to the closest heterospecific showed a significant decrease with increasing geographical scale of sampling. The average genetic distance dropped from > 7% for samples within 1 km, to < 3.5% for samples up to > 6000 km apart. Over a third of the species were not monophyletic, and the proportion increased through locally, nationally, regionally, and continentally restricted subsets of the data. The success of identifying queries decreased with increasing spatial scale of sampling; liberal methods declined from 100% to around 90%, whereas strict methods dropped to below 50% at continental scales. The proportion of query identifications considered uncertain (more than one species < 1% distance from query) escalated from zero at local, to 50% at continental scale. Finally, by resampling the most widely sampled species we show that even if samples are collected to maximize the geographical coverage, up to 70 individuals are required to sample 95% of intraspecific variation. The results show that the geographical scale of sampling has a critical impact on the global application of DNA barcoding. Scale-effects result from the relative importance of different processes determining the composition of regional species assemblages (dispersal and ecological assembly) and global clades (demography, speciation, and extinction). The incorporation of geographical information, where available, will be required to obtain identification rates at global scales equivalent to those in regional barcoding studies. Our result hence provides an impetus for both smarter barcoding tools and sprouting national barcoding initiatives-smaller geographical scales deliver higher accuracy.

Sampling strategy and potential utility of indels for DNA barcoding of closely related plant species: a case study in taxus

Although DNA barcoding has become a useful tool for species identification and biodiversity surveys in plant sciences, there remains little consensus concerning appropriate sampling strategies and the treatment of indels. To address these two issues, we sampled 39 populations for nine Taxus species across their entire ranges, with two to three individuals per population randomly sampled. We sequenced one core DNA barcode (matK) and three supplementary regions (trnH-psbA, trnL-trnF and ITS) for all samples to test the effects of sampling design and the utility of indels. Our results suggested that increasing sampling within-population did not change the clustering of individuals, and that meant within-population P-distances were zero for most populations in all regions. Based on the markers tested here, comparison of methods either including or excluding indels indicated that discrimination and nodal support of monophyletic groups were significantly increased when indels were included. Thus we concluded that one individual per population was adequate to represent the within-population variation in these species for DNA barcoding, and that intra-specific sampling was best focused on representing the entire ranges of certain taxa. We also found that indels occurring in the chloroplast trnL-trnF and trnH-psbA regions were informative to differentiate among for closely related taxa barcoding, and we proposed that indel-coding methods should be considered for use in future for closed related plant species DNA barcoding projects on or below generic level.

Mitochondrial barcodes are diagnostic of shared refugia but not species in hybridizing oak gallwasps

Mitochondrial DNA barcodes provide a simple taxonomic tool for systematic and ecological research, with particular benefit for poorly studied or species-rich taxa. Barcoding assumes genetic diversity follows species boundaries; however, many processes disrupt species-level monophyly of barcodes leading to incorrect classifications. Spatial population structure, particularly when shared across closely related and potentially hybridizing taxa, can invalidate barcoding approaches yet few data exist to examine its impacts. We test how shared population structure across hybridizing species impacts upon mitochondrial barcodes by sequencing the cytochrome b gene for 518 individuals of four well-delimited Western Palaearctic gallwasp species within the Andricus quercuscalicis species group. Mitochondrial barcodes clustered individuals into mixed-species clades corresponding to refugia, with no difference in within- and between-species divergence. Four nuclear genes were also sequenced from 4 to 11 individuals per refugial population of each species. Multi-locus analyses of these data supported established species, with no support for the refugial clustering across species seen in mitochondrial barcodes. This pattern is consistent with mitochondrial introgression among populations of species sharing the same glacial refugium, such that mitochondrial barcodes identify a shared history of population structure rather than species. Many taxa show phylogeographic structure across glacial refugia, suggesting that mitochondrial barcoding may fail when applied to other sets of co-distributed, closely related species. Robust barcoding approaches must sample extensively across population structure to disentangle spatial from species-level variation. Methods incorporating multiple unlinked loci are also essential to accommodate coalescent variation among genes and provide power to resolve recently diverged species.