Assembling a DNA barcode reference library for the spiders (Arachnida: Araneae) of Pakistan

Morphological study of 1,795 spiders from sites across Pakistan placed these specimens in 27 families and 202 putative species. COI sequences >400 bp recovered from 1,782 specimens were analyzed using neighbor-joining trees, Bayesian inference, barcode gap, and Barcode Index Numbers (BINs). Specimens of 109 morphological species were assigned to 123 BINs with ten species showing BIN splits, while 93 interim species included representatives of 98 BINs. Maximum conspecific divergences ranged from 0–5.3% while congeneric distances varied from 2.8–23.2%. Excepting one species pair (Oxyopes azhari–Oxyopes oryzae), the maximum intraspecific distance was always less than the nearest-neighbor (NN) distance. Intraspecific divergence values were not significantly correlated with geographic distance. Most (75%) BINs detected in this study were new to science, while those shared with other nations mainly derived from India. The discovery of many new, potentially endemic species and the low level of BIN overlap with other nations highlight the importance of constructing regional DNA barcode reference libraries.

A DNA barcode library for 5,200 German flies and midges (Insecta: Diptera) and its implications for metabarcoding-based biomonitoring

This study summarizes results of a DNA barcoding campaign on German Diptera, involving analysis of 45,040 specimens. The resultant DNA barcode library includes records for 2,453 named species comprising a total of 5,200 barcode index numbers (BINs), including 2,700 COI haplotype clusters without species‐level assignment, so called “dark taxa.” Overall, 88 out of 117 families (75%) recorded from Germany were covered, representing more than 50% of the 9,544 known species of German Diptera. Until now, most of these families, especially the most diverse, have been taxonomically inaccessible. By contrast, within a few years this study provided an intermediate taxonomic system for half of the German Dipteran fauna, which will provide a useful foundation for subsequent detailed, integrative taxonomic studies. Using DNA extracts derived from bulk collections made by Malaise traps, we further demonstrate that species delineation using BINs and operational taxonomic units (OTUs) constitutes an effective method for biodiversity studies using DNA metabarcoding. As the reference libraries continue to grow, and gaps in the species catalogue are filled, BIN lists assembled by metabarcoding will provide greater taxonomic resolution. The present study has three main goals: (a) to provide a DNA barcode library for 5,200 BINs of Diptera; (b) to demonstrate, based on the example of bulk extractions from a Malaise trap experiment, that DNA barcode clusters, labelled with globally unique identifiers (such as OTUs and/or BINs), provide a pragmatic, accurate solution to the “taxonomic impediment”; and (c) to demonstrate that interim names based on BINs and OTUs obtained through metabarcoding provide an effective method for studies on species‐rich groups that are usually neglected in biodiversity research projects because of their unresolved taxonomy.

Using eDNA to biomonitor the fish community in a tropical oligotrophic lake

Environmental DNA (eDNA) is an effective approach for detecting vertebrates and plants, especially in aquatic ecosystems, but prior studies have largely examined eDNA in cool temperate settings. By contrast, this study employs eDNA to survey the fish fauna in tropical Lake Bacalar (Mexico) with the additional goal of assessing the possible presence of invasive fishes, such as Amazon sailfin catfish and tilapia. Sediment and water samples were collected from eight stations in Lake Bacalar on three occasions over a 4-month interval. Each sample was stored in the presence or absence of lysis buffer to compare eDNA recovery. Short fragments (184–187 bp) of the cytochrome c oxidase I (COI) gene were amplified using fusion primers and then sequenced on Ion Torrent PGM or S5 before their source species were determined using a custom reference sequence database constructed on BOLD. In total, eDNA sequences were recovered from 75 species of vertebrates including 47 fishes, 15 birds, 7 mammals, 5 reptiles, and 1 amphibian. Although all species are known from this region, six fish species represent new records for the study area, while two require verification. Sequences for five species (2 birds, 2 mammals, 1 reptile) were only detected from sediments, while sequences from 52 species were only recovered from water. Because DNA from the Amazon sailfin catfish was not detected, we used a mock eDNA experiment to confirm our methods would enable its detection. In summary, we developed protocols that recovered eDNA from tropical oligotrophic aquatic ecosystems and confirmed their effectiveness in detecting fishes and diverse species of vertebrates.

Flower-visitor communities of an arcto-alpine plant-Global patterns in species richness, phylogenetic diversity and ecological functioning

Pollination is an ecosystem function of global importance. Yet, who visits the flower of specific plants, how the composition of these visitors varies in space and time and how such variation translates into pollination services are hard to establish. The use of DNA barcodes allows us to address ecological patterns involving thousands of taxa that are difficult to identify. To clarify the regional variation in the visitor community of a widespread flower resource, we compared the composition of the arthropod community visiting species in the genus Dryas (mountain avens, family Rosaceae), throughout Arctic and high-alpine areas. At each of 15 sites, we sampled Dryas visitors with 100 sticky flower mimics and identified specimens to Barcode Index Numbers (BINs) using a partial sequence of the mitochondrial COI gene. As a measure of ecosystem functioning, we quantified variation in the seed set of Dryas. To test for an association between phylogenetic and functional diversity, we characterized the structure of local visitor communities with both taxonomic and phylogenetic descriptors. In total, we detected 1,360 different BINs, dominated by Diptera and Hymenoptera. The richness of visitors at each site appeared to be driven by local temperature and precipitation. Phylogeographic structure seemed reflective of geological history and mirrored trans-Arctic patterns detected in plants. Seed set success varied widely among sites, with little variation attributable to pollinator species richness. This pattern suggests idiosyncratic associations, with function dominated by few and potentially different taxa at each site. Taken together, our findings illustrate the role of post-glacial history in the assembly of flower-visitor communities in the Arctic and offer insights for understanding how diversity translates into ecosystem functioning.

Approaches to integrating genetic data into ecological networks

As molecular tools for assessing trophic interactions become common, research is increasingly focused on the construction of interaction networks. Here, we demonstrate three key methods for incorporating DNA data into network ecology and discuss analytical considerations using a model consisting of plants, insects, bats and their parasites from the Costa Rica dry forest. The simplest method involves the use of Sanger sequencing to acquire long sequences to validate or refine field identifications, for example of bats and their parasites, where one specimen yields one sequence and one identification. This method can be fully quantified and resolved and these data resemble traditional ecological networks. For more complex taxonomic identifications, we target multiple DNA loci, for example from a seed or fruit pulp sample in faeces. These networks are also well resolved but gene targets vary in resolution and quantification is difficult. Finally, for mixed templates such as faecal contents of insectivorous bats, we use DNA metabarcoding targeting two sequence lengths (157 and 407 bp) of one gene region and a MOTU, BLAST and BIN association approach to resolve nodes. This network type is complex to generate and analyse, and we discuss the implications of this type of resolution on network analysis. Using these data, we construct the first molecular-based network of networks containing 3,304 interactions between 762 nodes of eight trophic functions and involving parasitic, mutualistic and predatory interactions. We provide a comparison of the relative strengths and weaknesses of these data types in network ecology.

Database establishment for the secondary fungal DNA barcode translational elongation factor 1α (TEF1α) 1

With new or emerging fungal infections, human and animal fungal pathogens are a growing threat worldwide. Current diagnostic tools are slow, non-specific at the species and subspecies levels, and require specific morphological expertise to accurately identify pathogens from pure cultures. DNA barcodes are easily amplified, universal, short species-specific DNA sequences, which enable rapid identification by comparison with a well-curated reference sequence collection. The primary fungal DNA barcode, ITS region, was introduced in 2012 and is now routinely used in diagnostic laboratories. However, the ITS region only accurately identifies around 75% of all medically relevant fungal species, which has prompted the development of a secondary barcode to increase the resolution power and suitability of DNA barcoding for fungal disease diagnostics. The translational elongation factor 1α (TEF1α) was selected in 2015 as a secondary fungal DNA barcode, but it has not been implemented into practice, due to the absence of a reference database. Here, we have established a quality-controlled reference database for the secondary barcode that together with the ISHAM-ITS database, forms the ISHAM barcode database, available online at http://its.mycologylab.org/ . We encourage the mycology community for active contributions.

Large geographic distance versus small DNA barcode divergence: Insights from a comparison of European to South Siberian Lepidoptera

Spanning nearly 13,000 km, the Palearctic region provides an opportunity to examine the level of geographic coverage required for a DNA barcode reference library to be effective in identifying species with broad ranges. This study examines barcode divergences between populations of 102 species of Lepidoptera from Europe and South Siberia, sites roughly 6,000 km apart. While three-quarters of these species showed divergence between their Asian and European populations, these divergence values ranged between 0-1%, distinctly less than the distance to the Nearest-Neighbor species in all but a few cases. Our results suggest that further taxonomic studies may be required for 16 species that showed either extremely low interspecific or high intraspecific variation. For example, seven species pairs showed low or no barcode divergence, but four of these cases are likely to reflect taxonomic over-splitting while the others involve species pairs that are either young or show evidence for introgression. Conversely, some of the nine species with deep intraspecific divergence at varied spatial levels may include overlooked species. Although these 16 cases require further investigation, our overall results indicate that barcode reference libraries based on records from one locality can be very effective in identifying specimens across an extensive geographic area.

Expedited assessment of terrestrial arthropod diversity by coupling Malaise traps with DNA barcoding 1

Monitoring changes in terrestrial arthropod communities over space and time requires a dramatic increase in the speed and accuracy of processing samples that cannot be achieved with morphological approaches. The combination of DNA barcoding and Malaise traps allows expedited, comprehensive inventories of species abundance whose cost will rapidly decline as high-throughput sequencing technologies advance. Aside from detailing protocols from specimen sorting to data release, this paper describes their use in a survey of arthropod diversity in a national park that examined 21 194 specimens representing 2255 species. These protocols can support arthropod monitoring programs at regional, national, and continental scales.

A survey of molecular diversity and population genetic structure in North American clearwing moths (Lepidoptera: Sesiidae) using cytochrome c oxidase I

The phylogeographic structure of insect species in North America is poorly understood. The moth family Sesiidae (Lepidoptera) contains many economically important pests of agriculture and forestry, as well as beneficial species used in biological control. Despite their significance, this study constitutes the first broad-ranging population genetic study on North American sesiids. It probes the population structure of eight species of sesiid moths based on sequence variation in cytochrome c oxidase I (N = 191). Haplotype diversity levels were high in seven of the eight species, while nucleotide diversity varied considerably. Patterns ranged from limited structure and a starburst pattern in the raspberry crown borer Pennisetia marginata to highly geographically structured populations in the peachtree borer Synanthedon exitiosa and the maple callus borer Synanthedon acerni. These varied patterns suggest differing evolutionary histories and dispersal abilities. By elucidating population genetic structure and barriers to dispersal we can begin to devise conservation and management plans.

Use of genetic, climatic, and microbiological data to inform reintroduction of a regionally extinct butterfly

Species reintroductions are increasingly used as means of mitigating biodiversity loss. Besides habitat quality at the site targeted for reintroduction, the choice of source population can be critical for success. The butterfly Melanargia russiae (Esper´s marbled white) was extirpated from Hungary over 100 years ago, and a reintroduction program has recently been approved. We used museum specimens of this butterfly, mitochondrial DNA data (mtDNA), endosymbiont screening, and climatic-similarity analyses to determine which extant populations should be used for its reintroduction. The species displayed 2 main mtDNA lineages across its range: 1 restricted to Iberia and southern France (Iberian lineage) and another found throughout the rest of its range (Eurasian lineage). These 2 lineages possessed highly divergent wsp alleles of the bacterial endosymbiont Wolbachia. The century-old Hungarian specimens represented an endemic haplotype belonging to the Eurasian lineage, differing by one mutation from the Balkan and eastern European populations. The Hungarian populations of M. russiae occurred in areas with a colder and drier climate relative to most sites with extant known populations. Our results suggest the populations used for reintroduction to Hungary should belong to the Eurasian lineage, preferably from eastern Ukraine (genetically close and living in areas with the highest climatic similarity). Materials stored in museum collections can provide unique opportunities to document historical genetic diversity and help direct conservation.

Slippage of degenerate primers can cause variation in amplicon length

It is well understood that homopolymer regions should be avoided for primer binding to prevent off-target amplification. However, in metabarcoding, it is often difficult to avoid primer degeneracy in order to maximize taxa detection. We here investigate primer binding specificity using different primer sets from several invertebrate metabarcoding studies. Our results indicate that primers frequently bound 1-2 bp upstream in taxa where a homopolymer region was present in the amplification direction. Primer binding 1 bp downstream was observed less frequently. This primer slippage leads to taxon-specific length variation in amplicons and subsequent length variation in recovered sequences. Some widely used primer sets were severely affected by this bias, while others were not. While this variation will only have small impacts on the designation of Operational Taxonomic Units (OTUs) by clustering algorithms that ignore terminal gaps, primer sets employed in metabarcoding projects should be evaluated for their sensitivity to slippage. Moreover, steps should be taken to reduce slippage by improving protocols for primer design. For example, the flanking region adjacent to the 3′ end of the primer is not considered by current primer development software although GC clamps in this position could mitigate slippage.

A Sequel to Sanger: amplicon sequencing that scales

BACKGROUND

Although high-throughput sequencers (HTS) have largely displaced their Sanger counterparts, the short read lengths and high error rates of most platforms constrain their utility for amplicon sequencing. The present study tests the capacity of single molecule, real-time (SMRT) sequencing implemented on the SEQUEL platform to overcome these limitations, employing 658 bp amplicons of the mitochondrial cytochrome c oxidase I gene as a model system.

RESULTS

By examining templates from more than 5000 species and 20,000 specimens, the performance of SMRT sequencing was tested with amplicons showing wide variation in GC composition and varied sequence attributes. SMRT and Sanger sequences were very similar, but SMRT sequencing provided more complete coverage, especially for amplicons with homopolymer tracts. Because it can characterize amplicon pools from 10,000 DNA extracts in a single run, the SEQUEL can reduce greatly reduce sequencing costs in comparison to first (Sanger) and second generation platforms (Illumina, Ion).

CONCLUSIONS

SMRT analysis generates high-fidelity sequences from amplicons with varying GC content and is resilient to homopolymer tracts. Analytical costs are low, substantially less than those for first or second generation sequencers. When implemented on the SEQUEL platform, SMRT analysis enables massive amplicon characterization because each instrument can recover sequences from more than 5 million DNA extracts a year.

Phylogeographic structure in three North American tent caterpillar species (Lepidoptera: Lasiocampidae): Malacosoma americana, M. californica, and M. disstria

While phylogeographic structure has been examined in many North American vertebrate species, insects have received much less attention despite their central ecological roles. The moth genus Malacosoma (Hübner, 1820), is an important group of forestry pests responsible for large-scale defoliation across much of the Nearctic and Palearctic. The present study uses sequence variation in the mitochondrial cytochrome c oxidase 1 (COI) gene to examine the population genetic structure of the three widespread Malacosoma species (M. americana, M. californica, and M. disstria). Populations of all three species showed highest diversity in the south, suggesting that modern populations derived from southern refugia with loss of variation as these lineages dispersed northwards. However, despite similar life histories and dispersal abilities, the extent of regional variation varied among the taxa. M. americana, a species restricted to eastern North America, showed much less genetic structure than the western M. californica or the widespread M. disstria. The regional differentiation in the latter reflects the likely derivation of modern lineages from several refugia, as well as taxonomic uncertainty in M. californica. In these respects, the three species of Malacosoma share phylogeographic patterns similar to those detected in vertebrates which are characterised by greater phylogeographic breaks in the western half of the continent and limited structure in the east.

Using herbarium-derived DNAs to assemble a large-scale DNA barcode library for the vascular plants of Canada

PREMISE OF THE STUDY

Constructing complete, accurate plant DNA barcode reference libraries can be logistically challenging for large-scale floras. Here we demonstrate the promise and challenges of using herbarium collections for building a DNA barcode reference library for the vascular plant flora of Canada.

METHODS

Our study examined 20,816 specimens representing 5076 of 5190 vascular plant species in Canada (98%). For 98% of the specimens, at least one of the DNA barcode regions was recovered from the plastid loci rbcL and matK and from the nuclear ITS2 region. We used beta regression to quantify the effects of age, type of preservation, and taxonomic affiliation (family) on DNA sequence recovery.

RESULTS

Specimen age and method of preservation had significant effects on sequence recovery for all markers, but influenced some families more (e.g., Boraginaceae) than others (e.g., Asteraceae).

DISCUSSION

Our DNA barcode library represents an unparalleled resource for metagenomic and ecological genetic research working on temperate and arctic biomes. An observed decline in sequence recovery with specimen age may be associated with poor primer matches, intragenomic variation (for ITS2), or inhibitory secondary compounds in some taxa.

Counting animal species with DNA barcodes: Canadian insects

Recent estimates suggest that the global insect fauna includes fewer than six million species, but this projection is very uncertain because taxonomic work has been limited on some highly diverse groups. Validation of current estimates minimally requires the investigation of all lineages that are diverse enough to have a substantial impact on the final species count. This study represents a first step in this direction; it employs DNA barcoding to evaluate patterns of species richness in 27 orders of Canadian insects. The analysis of over one million specimens revealed species counts congruent with earlier results for most orders. However, Diptera and Hymenoptera were unexpectedly diverse, representing two-thirds of the 46 937 barcode index numbers (=species) detected. Correspondence checks between known species and barcoded taxa showed that sampling was incomplete, a result confirmed by extrapolations from the barcode results which suggest the occurrence of at least 94 000 species of insects in Canada, a near doubling from the prior estimate of 54 000 species. One dipteran family, the Cecidomyiidae, was extraordinarily diverse with an estimated 16 000 species, a 10-fold increase from its predicted diversity. If Canada possesses about 1% of the global fauna, as it does for known taxa, the results of this study suggest the presence of 10 million insect species with about 1.8 million of these taxa in the Cecidomyiidae. If so, the global species count for this fly family may exceed the combined total for all 142 beetle families. If extended to more geographical regions and to all hyperdiverse groups, DNA barcoding can rapidly resolve the current uncertainty surrounding a species count for the animal kingdom. A newly detailed understanding of species diversity may illuminate processes important in speciation, as suggested by the discovery that the most diverse insect lineages in Canada employ an unusual mode of reproduction, haplodiploidy.

This article is part of the themed issue ‘From DNA barcodes to biomes’.

Barcoding the butterflies of southern South America: Species delimitation efficacy, cryptic diversity and geographic patterns of divergence

Because the tropical regions of America harbor the highest concentration of butterfly species, its fauna has attracted considerable attention. Much less is known about the butterflies of southern South America, particularly Argentina, where over 1,200 species occur. To advance understanding of this fauna, we assembled a DNA barcode reference library for 417 butterfly species of Argentina, focusing on the Atlantic Forest, a biodiversity hotspot. We tested the efficacy of this library for specimen identification, used it to assess the frequency of cryptic species, and examined geographic patterns of genetic variation, making this study the first large-scale genetic assessment of the butterflies of southern South America. The average sequence divergence to the nearest neighbor (i.e. minimum interspecific distance) was 6.91%, ten times larger than the mean distance to the furthest conspecific (0.69%), with a clear barcode gap present in all but four of the species represented by two or more specimens. As a consequence, the DNA barcode library was extremely effective in the discrimination of these species, allowing a correct identification in more than 95% of the cases. Singletons (i.e. species represented by a single sequence) were also distinguishable in the gene trees since they all had unique DNA barcodes, divergent from those of the closest non-conspecific. The clustering algorithms implemented recognized from 416 to 444 barcode clusters, suggesting that the actual diversity of butterflies in Argentina is 3%-9% higher than currently recognized. Furthermore, our survey added three new records of butterflies for the country (Eurema agave, Mithras hannelore, Melanis hillapana). In summary, this study not only supported the utility of DNA barcoding for the identification of the butterfly species of Argentina, but also highlighted several cases of both deep intraspecific and shallow interspecific divergence that should be studied in more detail.

Escaping introns in COI through cDNA barcoding of mushrooms: Pleurotus as a test case

DNA barcoding involves the use of one or more short, standardized DNA fragments for the rapid identification of species. A 648-bp segment near the 5' terminus of the mitochondrial cytochrome c oxidase subunit I (COI) gene has been adopted as the universal DNA barcode for members of the animal kingdom, but its utility in mushrooms is complicated by the frequent occurrence of large introns. As a consequence, ITS has been adopted as the standard DNA barcode marker for mushrooms despite several shortcomings. This study employed newly designed primers coupled with cDNA analysis to examine COI sequence diversity in six species of Pleurotus and compared these results with those for ITS. The ability of the COI gene to discriminate six species of Pleurotus, the commonly cultivated oyster mushroom, was examined by analysis of cDNA. The amplification success, sequence variation within and among species, and the ability to design effective primers was tested. We compared ITS sequences to their COI cDNA counterparts for all isolates. ITS discriminated between all six species, but some sequence results were uninterpretable, because of length variation among ITS copies. By comparison, a complete COI sequences were recovered from all but three individuals of Pleurotus giganteus where only the 5' region was obtained. The COI sequences permitted the resolution of all species when partial data was excluded for P. giganteus. Our results suggest that COI can be a useful barcode marker for mushrooms when cDNA analysis is adopted, permitting identifications in cases where ITS cannot be recovered or where it offers higher resolution when fresh tissue is. The suitability of this approach remains to be confirmed for other mushrooms.

Probing planetary biodiversity with DNA barcodes: The Noctuoidea of North America

This study reports the assembly of a DNA barcode reference library for species in the lepidopteran superfamily Noctuoidea from Canada and the USA. Based on the analysis of 69,378 specimens, the library provides coverage for 97.3% of the noctuoid fauna (3565 of 3664 species). In addition to verifying the strong performance of DNA barcodes in the discrimination of these species, the results indicate close congruence between the number of species analyzed (3565) and the number of sequence clusters (3816) recognized by the Barcode Index Number (BIN) system. Distributional patterns across 12 North American ecoregions are examined for the 3251 species that have GPS data while BIN analysis is used to quantify overlap between the noctuoid faunas of North America and other zoogeographic regions. This analysis reveals that 90% of North American noctuoids are endemic and that just 7.5% and 1.8% of BINs are shared with the Neotropics and with the Palearctic, respectively. One third (29) of the latter species are recent introductions and, as expected, they possess low intraspecific divergences.