Taxonomic Structure and Wing Pattern Evolution in the Parnassius mnemosyne Species Complex (Lepidoptera, Papilionidae)

In our study, using the analysis of DNA barcodes and morphology (wing color, male genitalia, and female sphragis shape), we show that the group of species close to P. mnemosyne comprises the western and eastern phylogenetic lineages. The eastern lineage includes P. stubbendorfii, P. glacialis, and P. hoenei. The western lineage includes three morphologically similar species: P. mnemosyne (Western Eurasia), P. turatii (southwestern Europe), and P. nubilosusstat. nov. (Turkmenistan and NE Iran), as well as the morphologically differentiated P. ariadne (Altai). The latter species differs from the rest of the group in the presence of red spots on the wings. Parnassius mnemosyne s.s. is represented by four differentiated mitochondrial clusters that show clear association with specific geographic regions. We propose to interpret them as subspecies: P. mnemosyne mnemosyne (Central and Eastern Europe, N Caucasus, N Turkey), P. mnemosyne adolphi (the Middle East), P. mnemosyne falsa (Tian Shan), and P. mnemosyne gigantea (Gissar-Alai in Central Asia). We demonstrate that in P. ariadne, the red spots on the wing evolved as a reversion to the ancestral wing pattern. This reversion is observed in Altai, where the distribution areas of the western lineage, represented by P. ariadne, and the eastern lineage, represented by P. stubbendorfii, overlap. These two species hybridize in Altai, and we hypothesize that the color change in P. ariadne is the result of reinforcement of prezygotic isolation in the contact zone. The lectotype of Parnassius mnemosyne var. nubilosus Christoph, 1873, is designated.

Life barcoded by DNA barcodes

The modern concept of DNA-based barcoding for cataloguing biodiversity was proposed in 2003 by first adopting an approximately 600 bp fragment of the mitochondrial COI gene to compare via nucleotide alignments with known sequences from specimens previously identified by taxonomists. Other standardized regions meeting barcoding criteria then are also evolving as DNA barcodes for fast, reliable and inexpensive assessment of species composition across all forms of life, including animals, plants, fungi, bacteria and other microorganisms. Consequently, global DNA barcoding campaigns have resulted in the formation of many online workbenches and databases, such as BOLD system, as barcode references, and facilitated the development of mini-barcodes and metabarcoding strategies as important extensions of barcode techniques. Here we intend to give an overview of the characteristics and features of these barcode markers and major reference libraries existing for barcoding the planet’s life, as well as to address the limitations and opportunities of DNA barcodes to an increasingly broader community of science and society.

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.

Molecular diversity and species delimitation in the family Gasteruptiidae (Hymenoptera: Evanioidea)

Gasteruptiidae Ashmead is an easily recognised family of wasps with ∼589 described species worldwide. Although well characterised by traditional taxonomy, multiple authors have commented on the extreme morphological uniformity of the group, making species-level identification difficult. This problem is enhanced by the lack of molecular data and molecular phylogenetic research for the group. We used 187 cytochrome c oxidase subunit I (COI) barcodes to explore the efficiency of sequence data to delimitate species in Gasteruptiidae. We undertook a graphical and discussion-based comparison of six methods for species delimitation, with the success of methods judged based on known species boundaries and morphology. Both distance-based (ABGD and jMOTU threshold analysis) and tree-based (GMYC and PTP) methods compared across multiple parameters recovered variable molecular operational taxonomic units (MOTUs), ranging from 55 to 123 MOTUs. Tree-based methods tended to split known morphological species less than distance-based methods, with the single-threshold GMYC method the most concordant with known morphospecies. Our results suggest that the incorporation of molecular species delimitation techniques provides a powerful tool to assist in the interpretation of species and help direct informed decisions with taxonomic uncertainty in the family.

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.

Singleton molecular species delimitation based on COI-5P barcode sequences revealed high cryptic/undescribed diversity for Chinese katydids (Orthoptera: Tettigoniidae)

Background

DNA barcoding has been developed as a useful tool for species discrimination. Several sequence-based species delimitation methods, such as Barcode Index Number (BIN), REfined Single Linkage (RESL), Automatic Barcode Gap Discovery (ABGD), a Java program uses an explicit, determinate algorithm to define Molecular Operational Taxonomic Unit (jMOTU), Generalized Mixed Yule Coalescent (GMYC), and Bayesian implementation of the Poisson Tree Processes model (bPTP), were used. Our aim was to estimate Chinese katydid biodiversity using standard DNA barcode cytochrome c oxidase subunit I (COI-5P) sequences.

Results

Detection of a barcoding gap by similarity-based analyses and clustering-base analyses indicated that 131 identified morphological species (morphospecies) were assigned to 196 BINs and were divided into four categories: (i) MATCH (83/131 = 64.89%), morphospecies were a perfect match between morphospecies and BINs (including 61 concordant BINs and 22 singleton BINs); (ii) MERGE (14/131 = 10.69%), morphospecies shared its unique BIN with other species; (iii) SPLIT (33/131 = 25.19%, when 22 singleton species were excluded, it rose to 33/109 = 30.28%), morphospecies were placed in more than one BIN; (iv) MIXTURE (4/131 = 5.34%), morphospecies showed a more complex partition involving both a merge and a split. Neighbor-joining (NJ) analyses showed that nearly all BINs and most morphospecies formed monophyletic cluster with little variation. The molecular operational taxonomic units (MOTUs) were defined considering only the more inclusive clades found by at least four of seven species delimitation methods. Our results robustly supported 61 of 109 (55.96%) morphospecies represented by more than one specimen, 159 of 213 (74.65%) concordant BINs, and 3 of 8 (37.5%) discordant BINs.

Conclusions

Molecular species delimitation analyses generated a larger number of MOTUs compared with morphospecies. If these MOTU splits are proven to be true, Chinese katydids probably contain a seemingly large proportion of cryptic/undescribed taxa. Future amplification of additional molecular markers, particularly from the nuclear DNA, may be especially useful for specimens that were identified here as problematic taxa.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1404-5) contains supplementary material, which is available to authorized users.

Metabarcoding using multiplexed markers increases species detection in complex zooplankton communities

Metabarcoding combines DNA barcoding with high-throughput sequencing, often using one genetic marker to understand complex and taxonomically diverse samples. However, species-level identification depends heavily on the choice of marker and the selected primer pair, often with a trade-off between successful species amplification and taxonomic resolution. We present a versatile metabarcoding protocol for biomonitoring that involves the use of two barcode markers (COI and 18S) and four primer pairs in a single high-throughput sequencing run, via sample multiplexing. We validate the protocol using a series of 24 mock zooplanktonic communities incorporating various levels of genetic variation. With the use of a single marker and single primer pair, the highest species recovery was 77%. With all three COI fragments, we detected 62%-83% of species across the mock communities, while the use of the 18S fragment alone resulted in the detection of 73%-75% of species. The species detection level was significantly improved to 89%-93% when both markers were used. Furthermore, multiplexing did not have a negative impact on the proportion of reads assigned to each species and the total number of species detected was similar to when markers were sequenced alone. Overall, our metabarcoding approach utilizing two barcode markers and multiple primer pairs per barcode improved species detection rates over a single marker/primer pair by 14% to 35%, making it an attractive and relatively cost-effective method for biomonitoring natural zooplankton communities. We strongly recommend combining evolutionary independent markers and, when necessary, multiple primer pairs per marker to increase species detection (i.e., reduce false negatives) in metabarcoding studies.

eDNA metabarcoding as a new surveillance approach for coastal Arctic biodiversity

Because significant global changes are currently underway in the Arctic, creating a large-scale standardized database for Arctic marine biodiversity is particularly pressing. This study evaluates the potential of aquatic environmental DNA (eDNA) metabarcoding to detect Arctic coastal biodiversity changes and characterizes the local spatio-temporal distribution of eDNA in two locations. We extracted and amplified eDNA using two COI primer pairs from ~80 water samples that were collected across two Canadian Arctic ports, Churchill and Iqaluit, based on optimized sampling and preservation methods for remote regions surveys. Results demonstrate that aquatic eDNA surveys have the potential to document large-scale Arctic biodiversity change by providing a rapid overview of coastal metazoan biodiversity, detecting nonindigenous species, and allowing sampling in both open water and under the ice cover by local northern-based communities. We show that DNA sequences of ~50% of known Canadian Arctic species and potential invaders are currently present in public databases. A similar proportion of operational taxonomic units was identified at the species level with eDNA metabarcoding, for a total of 181 species identified at both sites. Despite the cold and well-mixed coastal environment, species composition was vertically heterogeneous, in part due to river inflow in the estuarine ecosystem, and differed between the water column and tide pools. Thus, COI-based eDNA metabarcoding may quickly improve large-scale Arctic biomonitoring using eDNA, but we caution that aquatic eDNA sampling needs to be standardized over space and time to accurately evaluate community structure changes.

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.