Global Spore Sampling Project: A global, standardized dataset of airborne fungal DNA

Novel methods for sampling and characterizing biodiversity hold great promise for re-evaluating patterns of life across the planet. The sampling of airborne spores with a cyclone sampler, and the sequencing of their DNA, have been suggested as an efficient and well-calibrated tool for surveying fungal diversity across various environments. Here we present data originating from the Global Spore Sampling Project, comprising 2,768 samples collected during two years at 47 outdoor locations across the world. Each sample represents fungal DNA extracted from 24 m3 of air. We applied a conservative bioinformatics pipeline that filtered out sequences that did not show strong evidence of representing a fungal species. The pipeline yielded 27,954 species-level operational taxonomic units (OTUs). Each OTU is accompanied by a probabilistic taxonomic classification, validated through comparison with expert evaluations. To examine the potential of the data for ecological analyses, we partitioned the variation in species distributions into spatial and seasonal components, showing a strong effect of the annual mean temperature on community composition.

Global arthropod beta-diversity is spatially and temporally structured by latitude

Global biodiversity gradients are generally expected to reflect greater species replacement closer to the equator. However, empirical validation of global biodiversity gradients largely relies on vertebrates, plants, and other less diverse taxa. Here we assess the temporal and spatial dynamics of global arthropod biodiversity dynamics using a beta-diversity framework. Sampling includes 129 sampling sites whereby malaise traps are deployed to monitor temporal changes in arthropod communities. Overall, we encountered more than 150,000 unique barcode index numbers (BINs) (i.e. species proxies). We assess between site differences in community diversity using beta-diversity and the partitioned components of species replacement and richness difference. Global total beta-diversity (dissimilarity) increases with decreasing latitude, greater spatial distance and greater temporal distance. Species replacement and richness difference patterns vary across biogeographic regions. Our findings support long-standing, general expectations of global biodiversity patterns. However, we also show that the underlying processes driving patterns may be regionally linked.

A genetic atlas for the butterflies of continental Canada and United States

Multi-locus genetic data for phylogeographic studies is generally limited in geographic and taxonomic scope as most studies only examine a few related species. The strong adoption of DNA barcoding has generated large datasets of mtDNA COI sequences. This work examines the butterfly fauna of Canada and United States based on 13,236 COI barcode records derived from 619 species. It compiles i) geographic maps depicting the spatial distribution of haplotypes, ii) haplotype networks (minimum spanning trees), and iii) standard indices of genetic diversity such as nucleotide diversity (π), haplotype richness (H), and a measure of spatial genetic structure (GST). High intraspecific genetic diversity and marked spatial structure were observed in the northwestern and southern North America, as well as in proximity to mountain chains. While species generally displayed concordance between genetic diversity and spatial structure, some revealed incongruence between these two metrics. Interestingly, most species falling in this category shared their barcode sequences with one at least other species. Aside from revealing large-scale phylogeographic patterns and shedding light on the processes underlying these patterns, this work also exposed cases of potential synonymy and hybridization.

BOLD v4: A Centralized Bioinformatics Platform for DNA-Based Biodiversity Data

BOLD, the Barcode of Life Data System, supports the acquisition, storage, validation, analysis, and publication of DNA barcodes, activities requiring the integration of molecular, morphological, and distributional data. Its pivotal role in curating the reference library of DNA barcodes, coupled with its data management and analysis capabilities, makes it a central resource for biodiversity science. It enables rapid, accurate identification of specimens and also reveals patterns of genetic diversity and evolutionary relationships among taxa.Launched in 2005, BOLD has become an increasingly powerful tool for advancing the understanding of planetary biodiversity. It currently hosts 17 million specimen records and 14 million barcodes that provide coverage for more than a million species from every continent and ocean. The platform has the long-term goal of providing a consistent, accurate system for identifying all species of eukaryotes.BOLD's integrated analytical tools, full data lifecycle support, and secure collaboration framework distinguish it from other biodiversity platforms. BOLD v4 brought enhanced data management and analysis capabilities as well as novel functionality for data dissemination and publication. Its next version will include features to strengthen its utility to the research community, governments, industry, and society-at-large.

Trait-mediated responses of caterpillar communities to spongy moth outbreaks and subsequent tebufenozide treatments

Outbreaks of the spongy moth Lymantria dispar can have devastating impacts on forest resources and ecosystems. Lepidoptera-specific insecticides, such as Bacillus thuringiensis var. kurstaki (BTK) and tebufenozide, are often deployed to prevent heavy defoliation of the forest canopy. While it has been suggested that using BTK poses less risk to non-target Lepidoptera than leaving an outbreak untreated, in situ testing of this assumption has been impeded by methodological challenges. The trade-offs between insecticide use and outbreaks have yet to be addressed for tebufenozide, which is believed to have stronger side effects than BTK. We investigated the short-term trade-offs between tebufenozide treatments and no-action strategies for the non-target herbivore community in forest canopies. Over 3 years, Lepidoptera and Symphyta larvae were sampled by canopy fogging in 48 oak stands in southeast Germany during and after a spongy moth outbreak. Half of the sites were treated with tebufenozide and changes in canopy cover were monitored. We contrasted the impacts of tebufenozide and defoliator outbreaks on the abundance, diversity, and functional structure of chewing herbivore communities. Tebufenozide treatments strongly reduced Lepidoptera up to 6 weeks after spraying. Populations gradually converged back to control levels after 2 years. Shelter-building species dominated caterpillar assemblages in treated plots in the post-spray weeks, while flight-dimorphic species were slow to recover and remained underrepresented in treated stands 2 years post-treatment. Spongy moth outbreaks had minor effects on leaf chewer communities. Summer Lepidoptera decreased only when severe defoliation occurred, whereas Symphyta declined 1 year after defoliation. Polyphagous species with only partial host plant overlap with the spongy moth were absent from heavily defoliated sites, suggesting greater sensitivity of generalists to defoliation-induced plant responses. These results demonstrate that both tebufenozide treatments and spongy moth outbreaks alter canopy herbivore communities. Tebufenozide had a stronger and longer lasting impact, but it was restricted to Lepidoptera, whereas the outbreak affected both Lepidoptera and Symphyta. These results are tied to the fact that only half of the outbreak sites experienced severe defoliation. This highlights the limited accuracy of current defoliation forecast methods, which are used as the basis for the decision to spray insecticides.

Two new genera and one new species of the tribe Adeshini (Hymenoptera, Braconidae, Braconinae) from India and South Africa

Two new genera and one new species of the Braconinae tribe Adeshini are described and illustrated: Crenuladesha Ranjith & Quicke, gen. nov., type species Adeshanarendrani Ranjith, 2017, comb. nov. from India, and Protadesha Quicke & Butcher, gen. nov., type species Protadeshaintermedia Quicke & Butcher, sp. nov. from South Africa. The former lacks the mid-longitudinal propodeal carina characteristic of the tribe, and the latter displays less derived fore wing venation with two distinct abscissae of vein 2CU. A molecular phylogenetic analysis is included to confirm their correct placement. Since neither of the two new genera displays all of the characters given in the original diagnosis of the Adeshini a revised diagnosis is provided, as well as an illustrated key to the genera.

Interrogating 1000 insect genomes for NUMTs: A risk assessment for estimates of species richness

The nuclear genomes of most animal species include NUMTs, segments of the mitogenome incorporated into their chromosomes. Although NUMT counts are known to vary greatly among species, there has been no comprehensive study of their frequency/attributes in the most diverse group of terrestrial organisms, insects. This study examines NUMTs derived from a 658 bp 5' segment of the cytochrome c oxidase I (COI) gene, the barcode region for the animal kingdom. This assessment is important because unrecognized NUMTs can elevate estimates of species richness obtained through DNA barcoding and derived approaches (eDNA, metabarcoding). This investigation detected nearly 10,000 COI NUMTs ≥ 100 bp in the genomes of 1,002 insect species (range = 0-443). Variation in nuclear genome size explained 56% of the mitogenome-wide variation in NUMT counts. Although insect orders with the largest genome sizes possessed the highest NUMT counts, there was considerable variation among their component lineages. Two thirds of COI NUMTs possessed an IPSC (indel and/or premature stop codon) allowing their recognition and exclusion from downstream analyses. The remainder can elevate species richness as they showed 10.1% mean divergence from their mitochondrial homologue. The extent of exposure to "ghost species" is strongly impacted by the target amplicon's length. NUMTs can raise apparent species richness by up to 22% when a 658 bp COI amplicon is examined versus a doubling of apparent richness when 150 bp amplicons are targeted. Given these impacts, metabarcoding and eDNA studies should target the longest possible amplicons while also avoiding use of 12S/16S rDNA as they triple NUMT exposure because IPSC screens cannot be employed.

Butterfly–parasitoid–hostplant interactions in Western Palaearctic Hesperiidae: a DNA barcoding reference library

The study of ecological interactions between plants, phytophagous insects and their natural enemies is an essential but challenging component for understanding ecosystem dynamics. Molecular methods such as DNA barcoding can help elucidate these interactions. In this study, we employed DNA barcoding to establish hostplant and parasitoid interactions with hesperiid butterflies, using a complete reference library for Hesperiidae of continental Europe and north-western Africa (53 species, 100% of those recorded) based on 2934 sequences from 38 countries. A total of 233 hostplant and parasitoid interactions are presented, some recovered by DNA barcoding larval remains or parasitoid cocoons. Combining DNA barcode results with other lines of evidence allowed 94% species-level identification for Hesperiidae, but success was lower for parasitoids, in part due to unresolved taxonomy. Potential cases of cryptic diversity, both in Hesperiidae and Microgastrinae, are discussed. We briefly analyse the resulting interaction networks. Future DNA barcoding initiatives in this region should focus attention on north-western Africa and on parasitoids, because in these cases barcode reference libraries and taxonomy are less well developed.

Do pseudogenes pose a problem for metabarcoding marine animal communities?

Because DNA metabarcoding typically employs sequence diversity among mitochondrial amplicons to estimate species composition, nuclear mitochondrial pseudogenes (NUMTs) can inflate diversity. This study quantifies the incidence and attributes of NUMTs derived from the 658 bp barcode region of cytochrome c oxidase I (COI) in 156 marine animal genomes. NUMTs were examined to ascertain if they could be recognized by their possession of indels or stop codons. In total, 309 NUMTs ≥ 150 bp were detected, with an average of 1.98 per species (range = 0-33) and a mean length of 391 bp ± 200 bp. Among this total, 75 (24.3 %) lacked indels or stop codons. NUMTs appear to pose the greatest interpretational risk when short (< 313 bp) amplicons are used, such as in eDNA studies, dietary analyses, or processed fish identification. Employing the standard amplicon length (313 bp) for marine metabarcoding, NUMTs could potentially inflate the OTU count by 21% above the true species count while also raising intraspecific variation at COI by 15%. However, when both amplicon length and position are considered, inflation in OTU counts and in barcode variation were just 9% and 10%, respectively, suggesting NUMTs will not seriously distort biodiversity assessments. There was a weak positive correlation between genome size and NUMT count but no variation among phyla or trophic groups. Until bioinformatic advances improve NUMT detection, the best defense involves targeting long amplicons and developing reference databases that include both mitochondrial sequences and their NUMT derivatives.

Culicoides species community composition and feeding preferences in two aquatic ecosystems in northern Spain

Background

Aquatic ecosystems provide breeding sites for blood-sucking insects such as Culicoides biting midges (Diptera: Ceratopogonidae), but factors affecting their distribution and host choice are poorly understood. A study was undertaken at two nature reserves in northern Spain to examine the abundance, species composition, population dynamics and feeding patterns of biting midges between 2018 and 2019.

Methods

Culicoides were captured by light suction traps baited with CO₂ and by sweep netting vegetation. Blood meals and species identification of blood-fed specimens were determined using cytochrome c oxidase I subunit (COI) DNA barcoding. Multivariate generalized linear models were used to evaluate the associations between the abundance of Culicoides, the species richness and other parameters.

Results

The 4973 identified specimens comprised 28 species of Culicoides. These included two species reported for the first time in northern Spain, thus raising to 54 the number of Culicoides species described in the region. Specimens of all 28 species and 99.6% of the total specimens collected were caught in suction traps, while sweep netting vegetation revealed just 11 species and 0.4% of the total specimens. Midge abundance peaked in June/early July, with five species comprising > 80% of the captures: Culicoides alazanicus (24.9%), Culicoides griseidorsum (20.3%), Culicoides poperinghensis (16.2%), Culicoides kibunensis (10.7%) and Culicoides clastrieri (9.6%). DNA barcode analysis of blood meals from eight Culicoides species revealed that they fed on 17 vertebrate species (3 mammals and 14 birds). Species in the subgenus Avaritia were primarily ornithophilic, except for C. griseidorsum and C. poperinghensis. Host DNA from blood meals was successfully amplified from 75% of blood-fed females. A pictorial blood meal digestion scale is provided to accurately assess the blood-fed status of female Culicoides.

Conclusions

The large number of different blood meal sources identified in the midges captured in this study signals the likely importance of wild birds and mammals (e.g. red deer and wild boar) as reservoir/amplifying hosts for pathogens. Available hosts are more exposed to being bitten by biting midge populations in aquatic ecosystems in late spring and early summer.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-022-05297-5.

Description of Chilearinus Sharkey gen. nov. and status of Nearctic Earinus Wesmael, 1837 (Braconidae, Agathidinae) with the description of new species

The Neotropical members formerly included in Earinus Wesmael, 1837 are transferred to a new genus, Chilearinus Sharkey gen. nov. Presently three Nearctic species of Earinus are recognized, i.e., Earinuserythropoda Cameron, 1887, Earinuslimitaris Say,1835, and Earinuszeirapherae Walley, 1935, and these are retained in Earinus. Earinuschubuquensis Berta, 2000 and Earinusscitus Enderlein, 1920 are transferred to Chilearinus, i.e., C.chubuquensis, and C.scitus, comb. nov. One other species is transferred to Chilearinus, i.e., Microgasterrubricollis Spinola, 1851, Chilearinusrubricollis, comb. nov. Two other Neotropical species, Earinushubrechtae Braet, 2002 and Earinusbourguignoni Braet, 2002 were described under the genus Earinus but are here transferred to Lytopylus, L.hubrechtae, and L.bourguignoni comb. nov. Two new species of Chilearinus are described, C.covidchronos and C.janbert spp. nov. The status of Agathislaevithorax Spinola,1851, Agathisrubricata Spinola,1851, and Agathisareolata Spinola, 1851 is discussed. A neotype is designated for Earinuslimitaris (Say, 1835) and diagnosed with a COI barcode. Earinusaustinbakeri and Earinuswalleyi spp. nov. are described. The status of both Earinus and Chilearinus in the Americas is discussed. A revised key to the genera of Agathidinae of the Americas is presented.

Message in a Bottle-Metabarcoding enables biodiversity comparisons across ecoregions

Background

Traditional biomonitoring approaches have delivered a basic understanding of biodiversity, but they cannot support the large-scale assessments required to manage and protect entire ecosystems. This study used DNA metabarcoding to assess spatial and temporal variation in species richness and diversity in arthropod communities from 52 protected areas spanning 3 Canadian ecoregions.

Results

This study revealed the presence of 26,263 arthropod species in the 3 ecoregions and indicated that at least another 3,000–5,000 await detection. Results further demonstrate that communities are more similar within than between ecoregions, even after controlling for geographical distance. Overall α-diversity declined from east to west, reflecting a gradient in habitat disturbance. Shifts in species composition were high at every site, with turnover greater than nestedness, suggesting the presence of many transient species.

Conclusions

Differences in species composition among their arthropod communities confirm that ecoregions are a useful synoptic for biogeographic patterns and for structuring conservation efforts. The present results also demonstrate that metabarcoding enables large-scale monitoring of shifts in species composition, making it possible to move beyond the biomass measurements that have been the key metric used in prior efforts to track change in arthropod communities.

A DNA barcode survey of insect biodiversity in Pakistan

Although Pakistan has rich biodiversity, many groups are poorly known, particularly insects. To address this gap, we employed DNA barcoding to survey its insect diversity. Specimens obtained through diverse collecting methods at 1,858 sites across Pakistan from 2010-2019 were examined for sequence variation in the 658 bp barcode region of the cytochrome c oxidase 1 (COI) gene. Sequences from nearly 49,000 specimens were assigned to 6,590 Barcode Index Numbers (BINs), a proxy for species, and most (88%) also possessed a representative image on the Barcode of Life Data System (BOLD). By coupling morphological inspections with barcode matches on BOLD, every BIN was assigned to an order (19) and most (99.8%) were placed to a family (362). However, just 40% of the BINs were assigned to a genus (1,375) and 21% to a species (1,364). Five orders (Coleoptera, Diptera, Hemiptera, Hymenoptera, Lepidoptera) accounted for 92% of the specimens and BINs. More than half of the BINs (59%) are so far only known from Pakistan, but others have also been reported from Bangladesh (13%), India (12%), and China (8%). Representing the first DNA barcode survey of the insect fauna in any South Asian country, this study provides the foundation for a complete inventory of the insect fauna in Pakistan while also contributing to the global DNA barcode reference library.

DNA barcoding uncovers cryptic diversity in minute herbivorous mites (Acari, Eriophyoidea)

Eriophyoid mites (Acari: Eriophyoidea) are among the smallest of terrestrial arthropods and the most species-rich group of herbivorous mites with a high host specificity. However, knowledge of their species diversity has been impeded by the difficulty of their morphological differentiation. This study assembles a DNA barcode reference library that includes 1,850 mitochondrial COI sequences which provides coverage for 45% of the 930 species of eriophyoid mites known from China, and for 37 North American species. Sequence analysis showed a clear barcode gap in nearly all species, reflecting the fact that intraspecific divergences averaged 0.97% versus a mean of 18.51% for interspecific divergences (minimum nearest-neighbor distances) in taxa belonging to three families. Based on these results, we used DNA barcoding to explore the species diversity of eriophyoid mites as well as their host interactions. The 1,850 sequences were assigned to 531 Barcode Index Numbers (BINs). Analyses examining the correspondence between these BINs and species identifications based on morphology revealed that members of 45 species were assigned to two or more BINs, resulting in 1.16 times more BINs than morphospecies. Richness projections suggest that over 2,345 BINs occurred at the sampled locations. Host plant analysis showed that 89% of these mites (BINs) attack only one or two congeneric host species, but the others have several hosts. Furthermore, host-mite network analyses demonstrate that eriophyoid mites are high host-specific, and modularity is high in plant-mite networks. By creating a highly effective identification system for eriophyoid mites in BOLD, DNA barcoding will advance our understanding of the diversity of eriophyoid mites and their host interactions.

Unearthing soil arthropod diversity through DNA metabarcoding

DNA metabarcoding has the potential to greatly advance understanding of soil biodiversity, but this approach has seen limited application for the most abundant and species-rich group of soil fauna-the arthropods. This study begins to address this gap by comparing information on species composition recovered from metabarcoding two types of bulk samples (specimens, soil) from a temperate zone site and from bulk soil samples collected at eight sites in the Arctic. Analysis of 22 samples (3 specimen, 19 soil) revealed 410 arthropod OTUs belonging to 112 families, 25 orders, and nine classes. Studies at the temperate zone site revealed little overlap in species composition between soil and specimen samples, but more overlap at higher taxonomic levels (families, orders) and congruent patterns of α- and β-diversity. Expansion of soil analyses to the Arctic revealed locally rich, highly dissimilar, and spatially structured assemblages compatible with dispersal limited and environmentally driven assembly. The current study demonstrates that DNA metabarcoding of bulk soil enables rapid, large-scale assessments of soil arthropod diversity. However, deep sequence coverage is required to adequately capture the species present in these samples, and expansion of the DNA barcode reference library is necessary to improve taxonomic resolution of the sequences recovered through this approach.

Peering into the Darkness: DNA Barcoding Reveals Surprisingly High Diversity of Unknown Species of Diptera (Insecta) in Germany

Determining the size of the German insect fauna requires better knowledge of several megadiverse families of Diptera and Hymenoptera that are taxonomically challenging. This study takes the first step in assessing these “dark taxa” families and provides species estimates for four challenging groups of Diptera (Cecidomyiidae, Chironomidae, Phoridae, and Sciaridae). These estimates are based on more than 48,000 DNA barcodes (COI) from Diptera collected by Malaise traps that were deployed in southern Germany. We assessed the fraction of German species belonging to 11 fly families with well-studied taxonomy in these samples. The resultant ratios were then used to estimate the species richness of the four “dark taxa” families (DT families hereafter). Our results suggest a surprisingly high proportion of undetected biodiversity in a supposedly well-investigated country: at least 1800–2200 species await discovery in Germany in these four families. As this estimate is based on collections from one region of Germany, the species count will likely increase with expanded geographic sampling.

A molecular-based identification resource for the arthropods of Finland

To associate specimens identified by molecular characters to other biological knowledge, we need reference sequences annotated by Linnaean taxonomy. In this study, we (1) report the creation of a comprehensive reference library of DNA barcodes for the arthropods of an entire country (Finland), (2) publish this library, and (3) deliver a new identification tool for insects and spiders, as based on this resource. The reference library contains mtDNA COI barcodes for 11,275 (43%) of 26,437 arthropod species known from Finland, including 10,811 (45%) of 23,956 insect species. To quantify the improvement in identification accuracy enabled by the current reference library, we ran 1000 Finnish insect and spider species through the Barcode of Life Data system (BOLD) identification engine. Of these, 91% were correctly assigned to a unique species when compared to the new reference library alone, 85% were correctly identified when compared to BOLD with the new material included, and 75% with the new material excluded. To capitalize on this resource, we used the new reference material to train a probabilistic taxonomic assignment tool, FinPROTAX, scoring high success. For the full-length barcode region, the accuracy of taxonomic assignments at the level of classes, orders, families, subfamilies, tribes, genera, and species reached 99.9%, 99.9%, 99.8%, 99.7%, 99.4%, 96.8%, and 88.5%, respectively. The FinBOL arthropod reference library and FinPROTAX are available through the Finnish Biodiversity Information Facility (www.laji.fi) at https://laji.fi/en/theme/protax. Overall, the FinBOL investment represents a massive capacity-transfer from the taxonomic community of Finland to all sectors of society.

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

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

Assessment of current taxonomic assignment strategies for metabarcoding eukaryotes

The effective use of metabarcoding in biodiversity science has brought important analytical challenges due to the need to generate accurate taxonomic assignments. The assignment of sequences to genus or species level is critical for biodiversity surveys and biomonitoring, but it is particularly challenging as researchers must select the approach that best recovers information on species composition. This study evaluates the performance and accuracy of seven methods in recovering the species composition of mock communities by using COI barcode fragments. The mock communities varied in species number and specimen abundance, while upstream molecular and bioinformatic variables were held constant, and using a set of COI fragments. We evaluated the impact of parameter optimization on the quality of the predictions. Our results indicate that BLAST top hit competes well with more complex approaches if optimized for the mock community under study. For example, the two machine learning methods that were benchmarked proved more sensitive to reference database heterogeneity and completeness than methods based on sequence similarity. The accuracy of assignments was impacted by both species and specimen counts (query compositional heterogeneity) which ultimately influence the selection of appropriate software. We urge researchers to: (i) use realistic mock communities to allow optimization of parameters, regardless of the taxonomic assignment method employed; (ii) carefully choose and curate the reference databases including completeness; and (iii) use QIIME, BLAST or LCA methods, in conjunction with parameter tuning to better assign taxonomy to diverse communities, especially when information on species diversity is lacking for the area under study.

A DNA barcode library for the butterflies of North America

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

Debar: A sequence-by-sequence denoiser for COI-5P DNA barcode data

DNA barcoding and metabarcoding are now widely used to advance species discovery and biodiversity assessments. High-throughput sequencing (HTS) has expanded the volume and scope of these analyses, but elevated error rates introduce noise into sequence records that can inflate estimates of biodiversity. Denoising -the separation of biological signal from instrument (technical) noise-of barcode and metabarcode data currently employs abundance-based methods which do not capitalize on the highly conserved structure of the cytochrome c oxidase subunit I (COI) region employed as the animal barcode. This manuscript introduces debar, an R package that utilizes a profile hidden Markov model to denoise indel errors in COI sequences introduced by instrument error. In silico studies demonstrated that debar recognized 95% of artificially introduced indels in COI sequences. When applied to real-world data, debar reduced indel errors in circular consensus sequences obtained with the Sequel platform by 75%, and those generated on the Ion Torrent S5 by 94%. The false correction rate was less than 0.1%, indicating that debar is receptive to the majority of true COI variation in the animal kingdom. In conclusion, the debar package improves DNA barcode and metabarcode workflows by aiding the generation of more accurate sequences aiding the characterization of species diversity.

An Integrated Molecular Approach to Untangling Host-Vector-Pathogen Interactions in Mosquitoes (Diptera: Culicidae) From Sylvan Communities in Mexico

There are ~240 species of Culicidae in Mexico, of which some are vectors of arthropod-borne viruses such as Zika virus, dengue virus, chikungunya virus, and West Nile virus. Thus, the identification of mosquito feeding preferences is paramount to understanding of vector–host–pathogen interactions that, in turn, can aid the control of disease outbreaks. Typically, DNA and RNA are extracted separately for animal (insects and blood meal hosts) and viral identification, but this study demonstrates that multiple organisms can be analyzed from a single RNA extract. For the first time, residual DNA present in standard RNA extracts was analyzed by DNA barcoding in concert with Sanger and next-generation sequencing (NGS) to identify both the mosquito species and the source of their meals in blood-fed females caught in seven sylvan communities in Chiapas State, Mexico. While mosquito molecular identification involved standard barcoding methods, the sensitivity of blood meal identification was maximized by employing short primers with NGS. In total, we collected 1,634 specimens belonging to 14 genera, 25 subgenera, and 61 morphospecies of mosquitoes. Of these, four species were new records for Mexico (Aedes guatemala, Ae. insolitus, Limatus asulleptus, Trichoprosopon pallidiventer), and nine were new records for Chiapas State. DNA barcode sequences for >300 bp of the COI gene were obtained from 291 specimens, whereas 130 bp sequences were recovered from another 179 specimens. High intraspecific divergence values (>2%) suggesting cryptic species complexes were observed in nine taxa: Anopheles eiseni (5.39%), An. pseudopunctipennis (2.79%), Ae. podographicus (4.05%), Culex eastor (4.88%), Cx. erraticus (2.28%), Toxorhynchites haemorrhoidalis (4.30%), Tr. pallidiventer (4.95%), Wyeomyia adelpha/Wy. guatemala (7.30%), and Wy. pseudopecten (4.04%). The study increased the number of mosquito species known from 128 species to 138 species for Chiapas State, and 239 for Mexico as a whole. Blood meal analysis showed that Aedes angustivittatus fed on ducks and chicken, whereas Psorophora albipes fed on humans. Culex quinquefasciatus fed on diverse hosts including chicken, human, turkey, and Mexican grackle. No arbovirus RNA was detected by reverse transcriptase–polymerase chain reaction in the surveyed specimens. This study demonstrated, for the first time, that residual DNA present in RNA blood meal extracts can be used to identify host vectors, highlighting the important role of molecular approaches in both vector identification and revealing host–vector–pathogen interactions.

High resolution DNA barcode library for European butterflies reveals continental patterns of mitochondrial genetic diversity

The study of global biodiversity will greatly benefit from access to comprehensive DNA barcode libraries at continental scale, but such datasets are still very rare. Here, we assemble the first high-resolution reference library for European butterflies that provides 97% taxon coverage (459 species) and 22,306 COI sequences. We estimate that we captured 62% of the total haplotype diversity and show that most species possess a few very common haplotypes and many rare ones. Specimens in the dataset have an average 95.3% probability of being correctly identified. Mitochondrial diversity displayed elevated haplotype richness in southern European refugia, establishing the generality of this key biogeographic pattern for an entire taxonomic group. Fifteen percent of the species are involved in barcode sharing, but two thirds of these cases may reflect the need for further taxonomic research. This dataset provides a unique resource for conservation and for studying evolutionary processes, cryptic species, phylogeography, and ecology.

A SMRT approach for targeted amplicon sequencing of museum specimens (Lepidoptera)-patterns of nucleotide misincorporation

Natural history collections are a valuable resource for molecular taxonomic studies and for examining patterns of evolutionary diversification, particularly in the case of rare or extinct species. However, the recovery of sequence information is often complicated by DNA degradation. This article describes use of the Sequel platform (Pacific Biosciences) to recover the 658 bp barcode region of the mitochondrial cytochrome c oxidase I (COI) gene from 380 butterflies with an average age of 50 years. Nested multiplex PCR was employed for library preparation to facilitate sequence recovery from extracts with low concentrations of highly degraded DNA. By employing circular consensus sequencing (CCS) of short amplicons (circa 150 bp), full-length barcodes could be assembled without a reference sequence, an important advance from earlier protocols which required reference sequences to guide contig assembly. The Sequel protocol recovered COI sequences (499 bp on average) from 318 of 380 specimens (84%), much higher than for Sanger sequencing (26%). Because each read derives from a single molecule, it was also possible to quantify the incidence of substitutions arising from DNA damage. In agreement with past work on sequence changes induced by DNA degradation, the transition C/G → T/A was the most prevalent category of change, but its rate of occurrence (4.58E−4) was so low that it did not impede the recovery of reliable sequences. Because the current protocol recovers COI sequence from most museum specimens, and because sequence fidelity is unaffected by nucleotide misincorporations, large-scale sequence characterization of museum specimens is feasible.

Vertebrate-Aedes aegypti and Culex quinquefasciatus (Diptera)-arbovirus transmission networks: Non-human feeding revealed by meta-barcoding and next-generation sequencing

Background

Aedes aegypti mosquito-borne viruses including Zika (ZIKV), dengue (DENV), yellow fever (YFV), and chikungunya (CHIKV) have emerged and re-emerged globally, resulting in an elevated burden of human disease. Aedes aegypti is found worldwide in tropical, sub-tropical, and temperate areas. The characterization of mosquito blood meals is essential to understand the transmission dynamics of mosquito-vectored pathogens.

Methodology/principal findings

Here, we report Ae. aegypti and Culex quinquefasciatus host feeding patterns and arbovirus transmission in Northern Mexico using a metabarcoding-like approach with next-generation deep sequencing technology. A total of 145 Ae. aegypti yielded a blood meal analysis result with 107 (73.8%) for a single vertebrate species and 38 (26.2%) for two or more. Among the single host blood meals for Ae. aegypti, 28.0% were from humans, 54.2% from dogs, 16.8% from cats, and 1.0% from tortoises. Among those with more than one species present, 65.9% were from humans and dogs. For Cx. quinquefasciatus, 388 individuals yielded information with 326 (84%) being from a single host and 63 (16.2%) being from two or more hosts. Of the single species blood meals, 77.9% were from dogs, 6.1% from chickens, 3.1% from house sparrows, 2.4% from humans, while the remaining 10.5% derived from other 12 host species. Among those which had fed on more than one species, 11% were from dogs and humans, and 89% of other host species combinations. Forage ratio analysis revealed dog as the most over-utilized host by Ae. aegypti (= 4.3) and Cx. quinquefasciatus (= 5.6) and the human blood index at 39% and 4%, respectively. A total of 2,941 host-seeking female Ae. aegypti and 3,536 Cx. quinquefasciatus mosquitoes were collected in the surveyed area. Of these, 118 Ae. aegypti pools and 37 Cx. quinquefasciatus pools were screened for seven arboviruses (ZIKV, DENV 1–4, CHIKV, and West Nile virus (WNV)) using qRT-PCR and none were positive (point prevalence = 0%). The 95%-exact upper limit confidence interval was 0.07% and 0.17% for Ae. aegypti and Cx. quinquefasciatus, respectively

Conclusions/significance

The low human blood feeding rate in Ae. aegypti, high rate of feeding on mammals by Cx. quinquefasciatus, and the potential risk to transmission dynamics of arboviruses in highly urbanized areas of Northern Mexico is discussed.

Elucidating arbovirus-vector-host contact networks is critical to understand and control mosquito-borne virus transmission, including pathogens such as ZIKV, DENV 1–4, CHIKV, and WNV. Here, we report the results of metabarcoding of blood meals of two primary pathogen mosquito vectors, Aedes aegypti and Culex quinquefasciatus. We found limited human blood feeding by Ae. aegypti and high preference for feeding on mammals by Cx. quinquefasciatus. Interestingly, blood meal analysis revealed dogs as the most utilized host for both vector species suggesting the potential for zooprophylaxis for human-amplified urban arboviruses. Pools of these vector species were tested for seven arboviruses and all were negative. We calculated vectorial capacity to discuss the potential risk and transmission dynamics of pathogens transmitted by these two important vectors in an urban location in Northern Mexico.