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.

Dynamics of a host-parasitoid interaction clarified by modelling and DNA sequencing.

It has been hypothesised that the 2-year oscillations in abundance of Xestia moths are mediated by interactions with 1-year Ophion parasitoid wasps. We tested this hypothesis by modelling a 35-year time series of Xestia and Ophion from Northern Finland. Additionally, we used DNA barcoding to ascertain the species diversity of Ophion and targeted amplicon sequencing of their gut contents to confirm their larval hosts. Modelling of the time-series data strongly supported the hypothesised host-parasitoid dynamics and that periodic occurrence of Xestia moths is mediated by Ophion. DNA barcodes revealed that Ophion included five species rather than just one while targeted amplicon sequencing verified that Ophion does parasitise Xestia. At least one Ophion species employs 1-year Syngrapha interrogationis as an alternate host, but it did not detectably affect Xestia-Ophion dynamics. We also demonstrate the previously unrecognised complexity of this system due to cryptic parasitoid diversity.

The School Malaise Trap Program: Coupling educational outreach with scientific discovery.

The School Malaise Trap Program (SMTP) provides a technologically sophisticated and scientifically relevant educational experience that exposes students to the diversity of life, enhancing their understanding of biodiversity while promoting environmental stewardship. Since 2013, the SMTP has allowed 15,000 students at 350 primary and secondary schools to explore insect diversity in Canadian schoolyards. Students at each school collected hundreds of insects for an analysis of DNA sequence variation that enabled their rapid identification to a species. Through this hands-on approach, they participated in a learning exercise that conveys a real sense of scientific discovery. As well, the students contributed valuable data to the largest biodiversity genomics initiative ever undertaken: the International Barcode of Life project. To date, the SMTP has sequenced over 80,000 insect specimens, which includes representatives of 7,990 different species, nearly a tenth of the Canadian fauna. Both surprisingly and importantly, the collections generated the first DNA barcode records for 1,288 Canadian species.

Using DNA-barcoded Malaise trap samples to measure impact of a geothermal energy project on the biodiversity of a Costa Rican old-growth rain forest.

We report one year (2013-2014) of biomonitoring an insect community in a tropical old-growth rain forest, during construction of an industrial-level geothermal electricity project. This is the first-year reaction by the species-rich insect biodiversity; six subsequent years are being analyzed now. The site is on the margin of a UNESCO Natural World Heritage Site, Área de Conservación Guanacaste (ACG), in northwestern Costa Rica. This biomonitoring is part of Costa Rica's ongoing efforts to sustainably retain its wild biodiversity through biodevelopmental integration with its societies. Essential tools are geothermal engineering needs, entomological knowledge, insect species-rich forest, government-NGO integration, common sense, DNA barcoding for species-level identification, and Malaise traps. This research is tailored for integration with its society at the product level. We combine an academic view with on-site engineering decisions. This biomonitoring requires alpha-level DNA barcoding combined with centuries of morphology-based entomological taxonomy and ecology. Not all desired insect community analyses are performed; they are for data from subsequent years combined with this year. We provide enough analysis to be used by both guilds now. This biomonitoring has shown, for the first year, that the geothermal project impacts only the biodiversity within a zone less than 50 m from the project margin.

Nuclear genomes distinguish cryptic species suggested by their DNA barcodes and ecology.

DNA sequencing brings another dimension to exploration of biodiversity, and large-scale mitochondrial DNA cytochrome oxidase I barcoding has exposed many potential new cryptic species. Here, we add complete nuclear genome sequencing to DNA barcoding, ecological distribution, natural history, and subtleties of adult color pattern and size to show that a widespread neotropical skipper butterfly known as Udranomia kikkawai (Weeks) comprises three different species in Costa Rica. Full-length barcodes obtained from all three century-old Venezuelan syntypes of U. kikkawai show that it is a rainforest species occurring from Costa Rica to Brazil. The two new species are Udranomia sallydaleyae Burns, a dry forest denizen occurring from Costa Rica to Mexico, and Udranomia tomdaleyi Burns, which occupies the junction between the rainforest and dry forest and currently is known only from Costa Rica. Whereas the three species are cryptic, differing but slightly in appearance, their complete nuclear genomes totaling 15 million aligned positions reveal significant differences consistent with their 0.00065-Mbp (million base pair) mitochondrial barcodes and their ecological diversification. DNA barcoding of tropical insects reared by a massive inventory suggests that the presence of cryptic species is a widespread phenomenon and that further studies will substantially increase current estimates of insect species richness.

DNA barcodes expose unexpected diversity in Canadian mites.

Mites (Arachnida: Acariformes, Parasitiformes) are the most abundant and species-rich group of arthropods in soil, but are also diverse in freshwater habitats, on plants, and as symbionts of larger animals. However, assessment of their diversity has been impeded by their small size and often cryptic morphology. As a consequence, published estimates of their species richness span more than two orders of magnitude (0.4-114 million). In this study we employ DNA barcoding and the Barcode Index Number (BIN) system to investigate mite diversity at over 1,800 sites across Canada, primarily from soil and litter habitats with smaller contributions from freshwater, plants, and animal hosts. Barcodes from 73,394 specimens revealed 7,077 BINs with representatives from all four orders (Ixodida, Mesostigmata, Sarcoptiformes, Trombidiformes) and 60% (186) of the known families. The BIN total is 2.4 times the number of species previously recorded from Canada (2,999), reflecting the unexpectedly high richness of several families. Richness projections suggest that more than 28,000 BINs occur at the sampled locations, indicating that the Canadian mite fauna almost certainly includes more than 30,000 species-a total similar to that for the most diverse insect order in Canada, Diptera. This unexpected diversity was partitioned into highly dissimilar, spatially-structured assemblages that likely reflect dispersal limitation and environmental heterogeneity. Further sampling of a greater diversity of habitats will refine understanding of mite diversity in Canada, but similar analyses in other geographic regions will be essential to ascertain their diversity at a global scale.

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.

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.

Recalibrating the molecular clock for Arctic marine invertebrates based on DNA barcodes 1.

Divergence times for species assemblages of Arctic marine invertebrates have often been estimated using a standard rate (1.4%/MY) of molecular evolution calibrated using a single sister pair of tropical crustaceans. Because rates of molecular evolution vary among taxa and environments, it is essential to obtain clock calibrations from northern lineages. The recurrent opening and closure of the Bering Strait provide an exceptional opportunity for clock calibration. Here, we apply the iterative calibration approach to investigate patterns of molecular divergence among lineages of northern marine molluscs and arthropods using publicly available sequences of the cytochrome c oxidase subunit I (COI) gene and compare these results with previous estimates of trans-Bering divergences for echinoderms and polychaetes. The wide range of Kimura two-parameter (K2P) divergences among 73 trans-Bering sister pairs (0.12%-16.89%) supports multiple pulses of migration through the Strait. Overall, the results indicate a rate of K2P divergence of 3.2%/MY in molluscs, 5%-5.2%/MY in arthropods, and 3.5%-4.7%/MY in polychaetes. While these rates are considerably higher than the often-adopted 1.4%/MY rate, they are similar to calibrations (3%-5%/MY) in several other studies of marine invertebrates. This upward revision in rates means there is a need both to reevaluate the evolutionary history of marine lineages and to reexamine the impact of prior climatic changes upon the diversification of marine life.

Characterization and comparison of poorly known moth communities through DNA barcoding in two Afrotropical environments in Gabon 1.

Biodiversity research in tropical ecosystems-popularized as the most biodiverse habitats on Earth-often neglects invertebrates, yet invertebrates represent the bulk of local species richness. Insect communities in particular remain strongly impeded by both Linnaean and Wallacean shortfalls, and identifying species often remains a formidable challenge inhibiting the use of these organisms as indicators for ecological and conservation studies. Here we use DNA barcoding as an alternative to the traditional taxonomic approach for characterizing and comparing the diversity of moth communities in two different ecosystems in Gabon. Though sampling remains very incomplete, as evidenced by the high proportion (59%) of species represented by singletons, our results reveal an outstanding diversity. With about 3500 specimens sequenced and representing 1385 BINs (Barcode Index Numbers, used as a proxy to species) in 23 families, the diversity of moths in the two sites sampled is higher than the current number of species listed for the entire country, highlighting the huge gap in biodiversity knowledge for this country. Both seasonal and spatial turnovers are strikingly high (18.3% of BINs shared between seasons, and 13.3% between sites) and draw attention to the need to account for these when running regional surveys. Our results also highlight the richness and singularity of savannah environments and emphasize the status of Central African ecosystems as hotspots of biodiversity.

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.

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 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.

Stable baselines of temporal turnover underlie high beta diversity in tropical arthropod communities.

While the high species diversity of tropical arthropod communities has often been linked to marked spatial heterogeneity, their temporal dynamics have received little attention. This study addresses this gap by examining spatio-temporal variation in the arthropod communities of a tropical montane forest in Honduras. By employing DNA barcode analysis and Malaise trap sampling across 4 years and five sites, 51,596 specimens were assigned to 8,193 presumptive species. High beta diversity was linked more strongly to elevation than geographic distance, decreasing by 12% when only the dominant species were considered. When sampling effort was increased by deploying more traps at a site, beta diversity only decreased by 2%, but extending sampling across years decreased beta diversity by 27%. Species inconsistently detected among years, likely transients from other settings, drove the low similarity in species composition among traps only a few metres apart. The dominant, temporally persistent species substantially influenced the cyclic pattern of change in community composition among years. This pattern likely results from divergence-convergence dynamics, suggesting a stable baseline of temporal turnover in each community. The overall results establish that large sample sizes are necessary to reveal species richness, but are not essential for quantifying beta diversity. This study further highlights the need for standardized methods of sampling and species identification to generate the comparative data required to evaluate biodiversity change in space and time.

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.

Species-Level Para- and Polyphyly in DNA Barcode Gene Trees: Strong Operational Bias in European Lepidoptera.

The proliferation of DNA data is revolutionizing all fields of systematic research. DNA barcode sequences, now available for millions of specimens and several hundred thousand species, are increasingly used in algorithmic species delimitations. This is complicated by occasional incongruences between species and gene genealogies, as indicated by situations where conspecific individuals do not form a monophyletic cluster in a gene tree. In two previous reviews, non-monophyly has been reported as being common in mitochondrial DNA gene trees. We developed a novel web service "Monophylizer" to detect non-monophyly in phylogenetic trees and used it to ascertain the incidence of species non-monophyly in COI (a.k.a. cox1) barcode sequence data from 4977 species and 41,583 specimens of European Lepidoptera, the largest data set of DNA barcodes analyzed from this regard. Particular attention was paid to accurate species identification to ensure data integrity. We investigated the effects of tree-building method, sampling effort, and other methodological issues, all of which can influence estimates of non-monophyly. We found a 12% incidence of non-monophyly, a value significantly lower than that observed in previous studies. Neighbor joining (NJ) and maximum likelihood (ML) methods yielded almost equal numbers of non-monophyletic species, but 24.1% of these cases of non-monophyly were only found by one of these methods. Non-monophyletic species tend to show either low genetic distances to their nearest neighbors or exceptionally high levels of intraspecific variability. Cases of polyphyly in COI trees arising as a result of deep intraspecific divergence are negligible, as the detected cases reflected misidentifications or methodological errors. Taking into consideration variation in sampling effort, we estimate that the true incidence of non-monophyly is ∼23%, but with operational factors still being included. Within the operational factors, we separately assessed the frequency of taxonomic limitations (presence of overlooked cryptic and oversplit species) and identification uncertainties. We observed that operational factors are potentially present in more than half (58.6%) of the detected cases of non-monophyly. Furthermore, we observed that in about 20% of non-monophyletic species and entangled species, the lineages involved are either allopatric or parapatric-conditions where species delimitation is inherently subjective and particularly dependent on the species concept that has been adopted. These observations suggest that species-level non-monophyly in COI gene trees is less common than previously supposed, with many cases reflecting misidentifications, the subjectivity of species delimitation or other operational factors.

DNA barcoding as an aid for species identification in austral black flies (Insecta: Diptera: Simuliidae).

In this paper, the utility of a partial sequence of the COI gene, the DNA barcoding region, for the identification of species of black flies in the austral region was assessed. Twenty-eight morphospecies were analyzed: eight of the genus Austrosimulium (four species in the subgenus Austrosimulium s. str., three species in the subgenus Novaustrosimulium, and one species unassigned to subgenus), two of the genus Cnesia, eight of Gigantodax, three of Paracnephia, one of Paraustrosimulium, and six of Simulium (subgenera Morops, Nevermannia, and Pternaspatha). The neighbour-joining tree derived from the DNA barcode sequences grouped most specimens according to species or species groups recognized by morphotaxonomic studies. Intraspecific sequence divergences within morphologically distinct species ranged from 0% to 1.8%, while higher divergences (2%-4.2%) in certain species suggested the presence of cryptic diversity. The existence of well-defined groups within S. simile revealed the likely inclusion of cryptic diversity. DNA barcodes also showed that specimens identified as C. dissimilis, C. nr. pussilla, and C. ornata might be conspecific, suggesting possible synonymy. DNA barcoding combined with a sound morphotaxonomic framework would provide an effective approach for the identification of black flies in the region.

Complementary molecular information changes our perception of food web structure.

How networks of ecological interactions are structured has a major impact on their functioning. However, accurately resolving both the nodes of the webs and the links between them is fraught with difficulties. We ask whether the new resolution conferred by molecular information changes perceptions of network structure. To probe a network of antagonistic interactions in the High Arctic, we use two complementary sources of molecular data: parasitoid DNA sequenced from the tissues of their hosts and host DNA sequenced from the gut of adult parasitoids. The information added by molecular analysis radically changes the properties of interaction structure. Overall, three times as many interaction types were revealed by combining molecular information from parasitoids and hosts with rearing data, versus rearing data alone. At the species level, our results alter the perceived host specificity of parasitoids, the parasitoid load of host species, and the web-wide role of predators with a cryptic lifestyle. As the northernmost network of host-parasitoid interactions quantified, our data point exerts high leverage on global comparisons of food web structure. However, how we view its structure will depend on what information we use: compared with variation among networks quantified at other sites, the properties of our web vary as much or much more depending on the techniques used to reconstruct it. We thus urge ecologists to combine multiple pieces of evidence in assessing the structure of interaction webs, and suggest that current perceptions of interaction structure may be strongly affected by the methods used to construct them.

One fly to rule them all-muscid flies are the key pollinators in the Arctic.

Global change is causing drastic changes in the pollinator communities of the Arctic. While arctic flowers are visited by a wide range of insects, flies in family Muscidae have been proposed as a pollinator group of particular importance. To understand the functional outcome of current changes in pollinator community composition, we examined the role of muscids in the pollination of a key plant species, the mountain avens (Dryas). We monitored the seed set of Dryas across 15 sites at Zackenberg, northeast Greenland, and used sticky flower mimics and DNA barcoding to describe the flower-visiting community at each site. To evaluate the consequences of shifts in pollinator phenology under climate change, we compared the flower visitors between the early and the late season. Our approach revealed a diverse community of insects visiting Dryas, including two-thirds of all insect species known from the area. Even against this diverse background, the abundance of muscid flies emerged as a key predictor for seed set in Dryas, whereas overall insect abundance and species richness had little or no effect. With muscid flies as the main drivers of the pollinating function in the High Arctic, a recently observed decline in their abundances offers cause for concern.