DNA at the whim of the water: environmental DNA as a course-based undergraduate research experience.

Course-based undergraduate research experiences (CUREs) increase student access to high impact research experiences. CUREs engage students in the scientific process by learning how to pose scientific questions, develop hypotheses, and generate data to test them. Environmental DNA (eDNA) is a growing field of research that is gaining accessibility through decreasing laboratory costs, which can make a foundation for multiple, engaging CUREs. This manuscript describes three case studies that used eDNA in an upper year undergraduate course. The first focusses on a systematic literature review of eDNA metadata reporting. The second describes the biomonitoring of brook trout in southern Ontario using eDNA. The third involves eDNA metabarcoding for freshwater fish detection in southern Ontario. Undergraduates were involved in the development and execution of experiments, scientific communication, the peer review process, and fundraising. Through this manuscript, we show the novel application of eDNA CUREs and provide a roadmap for other instructors interested in implementing similar projects. Interviews with seven students from these courses indicate the benefits experienced from taking these courses. We argue that the use of eDNA in CUREs should be expanded in undergraduate biology programs due to the benefit to students and the increasing accessibility of this technology.

Mitochondrial genome sequencing, mapping, and assembly benchmarking for Culicoides species (Diptera: Ceratopogonidae).

BACKGROUND: Mitochondrial genomes are the most sequenced genomes after bacterial and fungal genomic DNA. However, little information on mitogenomes is available for multiple metazoan taxa, such as Culicoides, a globally distributed, megadiverse genus containing 1,347 species. AIM:  Generating novel mitogenomic information from single Culicoides sonorensis and C. biguttatus specimens, comparing available mitogenome mapping and de novo assembly tools, and identifying the best performing strategy and tools for Culicoides species. RESULTS: We present two novel and fully annotated mitochondrial haplotypes for two Culicoides species, C. sonorensis and C. biguttatus. We also annotated or re-annotated the only available reference mitogenome for C. sonorensis and C. arakawae. All species present a high similarity in mitogenome organization. The general gene arrangement for all Culicoides species was identical to the ancestral insect mitochondrial genome. Only short spacers were found in C. sonorensis (up to 30 bp), contrary to C. biguttatus (up to 114 bp). The mitochondrial genes ATP8, NAD2, NAD6, and LSU rRNA exhibited the highest nucleotide diversity and pairwise interspecific p genetic distance, suggesting that these genes might be suitable and complementary molecular barcodes for Culicoides identification in addition to the commonly utilized COI gene. We observed performance differences between the compared mitogenome generation strategies. The mapping strategy outperformed the de novo assembly strategy, but mapping results were partially biased in the absence of species-specific reference mitogenome. Among the utilized tools, BWA performed best for C. sonorensis while SPAdes, MEGAHIT, and MitoZ were among the best for C. biguttatus. The best-performing mitogenome annotator was MITOS2. Additionally, we were able to recover exogenous mitochondrial DNA from Bos taurus (biting midges host) from a C. biguttatus blood meal sample. CONCLUSIONS: Two novel annotated mitogenome haplotypes for C. sonorensis and C. biguttatus using High-Throughput Sequencing are presented. Current results are useful as the baseline for mitogenome reconstruction of the remaining Culicoides species from single specimens to HTS and genome annotation. Mapping to a species-specific reference mitogenome generated better results for Culicoides mitochondrial genome reconstruction than de novo assembly, while de novo assembly resulted better in the absence of a closely related reference mitogenome. These results have direct implications for molecular-based identification of these vectors of human and zoonotic diseases, setting the basis for using the whole mitochondrial genome as a marker in Culicoides identification.

The Multiple States of Environmental DNA and What Is Known about Their Persistence in Aquatic Environments.

Increased use of environmental DNA (eDNA) analysis for indirect species detection has spurred the need to understand eDNA persistence in the environment. Understanding the persistence of eDNA is complex because it exists in a mixture of different states (e.g., dissolved, particle adsorbed, intracellular, and intraorganellar), and each state is expected to have a specific decay rate that depends on environmental parameters. Thus, improving knowledge about eDNA conversion rates between states and the reactions that degrade eDNA in different states is needed. Here, we focus on eukaryotic extraorganismal eDNA, outline how water chemistry and suspended mineral particles likely affect conversion among each eDNA state, and indicate how environmental parameters affect persistence of states in the water column. On the basis of deducing these controlling parameters, we synthesized the eDNA literature to assess whether we could already derive a general understanding of eDNA states persisting in the environment. However, we found that these parameters are often not being measured or reported when measured, and in many cases very few experimental data exist from which to draw conclusions. Therefore, further study of how environmental parameters affect eDNA state conversion and eDNA decay in aquatic environments is needed. We recommend analytic controls that can be used during the processing of water to assess potential losses of different eDNA states if all were present in a water sample, and we outline future experimental work that would help determine the dominant eDNA states in water.

Assessment of stream macroinvertebrate communities with eDNA is not congruent with tissue-based metabarcoding.

Freshwater biomonitoring programmes routinely sample aquatic macroinvertebrates. These samples are time-consuming to collect, as well as challenging and costly to identify reliably genus or species. Environmental DNA (eDNA) metabarcoding has emerged as a surrogate to traditional collection techniques and has been used in whole-community approaches across several taxa and ecosystems. However, the usefulness of eDNA-based detection of freshwater macroinvertebrates has not been extensively explored. Few studies have directly compared bulk sample and eDNA metabarcoding at a local scale to assess how effective each method is at characterizing aquatic macroinvertebrate communities. Here, we collected both eDNA and kicknet samples at the same sample transect locations across nine different streams in southern Ontario, Canada. We observed minimal overlap in community composition between these paired samples. Bulk tissue metabarcoding resulted in a greater proportion of sequences belonging to metazoan taxa (over 99%) than eDNA (12%) and had higher OTU richness for macroinvertebrate taxa. We suggest that degenerate primers are not effective for eDNA metabarcoding due to the high degree of nontarget amplification and subsequently low yield of target DNA. While both bulk sample and eDNA metabarcoding had the power to detect differences between stream communities, eDNA did not represent local communities. Bulk tissue metabarcoding thus provides a more accurate representation of local stream macroinvertebrate communities and is the preferred method if smaller-scale spatial resolution is an important factor in data analyses.

Environmental DNA bioassays corroborate field data for detection of overwintering species at risk Blanding's turtles (Emydoidea blandingii).

Environmental DNA (eDNA) is gaining traction in conservation ecology as a powerful tool for detecting species at risk. We developed a quantitative polymerase chain reaction assay to detect a DNA amplicon fragment of the mitochondrial nicotinamide adenine dinucleotide locus of the Blanding's turtle (Emydoidea blandingii) for detecting overwintering individuals. Seventy-eight water samples were collected from 17 wetland sites in Ontario, Canada. We used traditional field data to identify a priori positive and negative control sites. Fifty percent of positive control sites amplified. Detection was related to the number of individuals estimated from field observations in at least one region surveyed. Positive control sites had lower total dissolved solids and electrical conductivity in relation to negative control sites. Shedding rates were within the same order of magnitude for brumating and active turtles. We recommend collecting additional samples at a larger number of locations to maximize detection. Recommended sampling design changes may overshadow the additional effects of water chemistry and low eDNA shedding rates. eDNA offers tremendous potential to practitioners conducting species at risk assessments in environmental consulting by providing a faster, more efficient method of detection compared with traditional surveys.

Environmental DNA detection of endangered and invasive species in Kejimkujik National Park and Historic Site.

The use of environmental DNA (eDNA) allows the early detection of aquatic species at low densities (e.g., elusive and invasive species), which otherwise could be challenging to monitor using conventional techniques. Here, we assess the ability of eDNA sampling to detect the presence or absence of one species at risk (Blanding's turtle) and two invasive species (chain pickerel and smallmouth bass) in Kejimkujik National Park and National Historic Site, Nova Scotia, where the aquatic system is highly acidic and rich in organic compounds. Five replicates of 1 L water samples were taken per sampling site. Water filtration and eDNA extractions were performed on-site, while qPCR reactions were performed in the laboratory using species-specific assays. Samples were treated with an inhibition removal kit and analyzed pre- and post-inhibition removal. Despite the low pH and PCR inhibitors in water samples, our results showed positive eDNA detections in almost all expected positive sites (except in one site for Blanding's turtle). Detections of the target species were also observed at sites where their presence was previously unknown. Our study supports the advantage of eDNA to monitor species at low densities, revealing new distributions or recently invaded areas. We also demonstrate how eDNA can directly instruct management strategies in Kejimkujik.

Rapid cooling via dry ice preserves the genetic and morphological integrity of fish embryos.

Larval fishes provide a valuable metric for assessing and monitoring species, populations, and ecosystem trends and condition. However, taxonomic resolution for this life stage is inherently problematic because of their individual sizes, limited morphological characteristics and high tissue degradation rates. There is little research on methods that rapidly preserve larval tissues for later morphological and molecular identification. The goal of this study was to test methods of rapidly killing fish embryos that maintain both morphological and molecular integrity. Rapid cooling with dry ice successfully maintained morphological and molecular integrity and may offer a simple and cost-effective approach for larval fish identification.

Assessing off-target cytotoxicity of the field lampricide 3-trifluoromethyl-4-nitrophenol using novel lake sturgeon cell lines.

Lampricides are currently being applied to streams and rivers to control the population of sea lamprey, an invasive species, in the Great Lakes. The most commonly used lampricide agent used in the field is 3-trifluoromethyl-4-nitrophenol (TFM), which targets larval sea lamprey in lamprey-infested rivers and streams. The specificity of TFM is due to the relative inability of sea lamprey to detoxify the agent relative to non-target fishes. There is increasing concern, however, about non-target effects on fishes, particularly threatened populations of juvenile lake sturgeon (LS; Acipenser fulvescens). There is therefore a need to develop models to better define lake sturgeon's response to TFM. Here we report the establishment of five LS cell lines derived from the liver, gill, skin and intestinal tract of juvenile LS and some of their cellular characteristics. All LS cell lines grew well at 25 °C in Leibovitz's (L)- 15 medium supplemented with 10% FBS. All cell lines demonstrated high senescence-associated ß-galactosidase activity and varying levels of Periodic acid Schiff-positive polysaccharides, indicating substantial production of glycoproteins and mucosubstances by the cells. Comparative toxicity of TFM in the five LS cell lines was assessed by two fluorescent cell viability dyes, Alamar Blue and CFDA-AM, in conditions with and without serum and at 24 or 72 h exposure. Deduced EC50 values were compared between the cell lines and to the reported in vivo LC50s. Tissues sensitive to the effects of TFM in vivo correlated with cell lines from the same tissues being most sensitive to TFM in vitro. EC50 values for the LSliver-e cells was significantly lower than the EC50 for the rainbow trout (RBT) liver cells RTL-W1, reaffirming the in vivo observation that LS was generally more TFM-sensitive than rainbow trout. Our data suggests that whole-fish sensitivity of LS to TFM is likely attributable to sensitivity at the cellular level. Thus, LS cell lines, as well as those of RBT, can be used to screen and evaluate the toxicity of the next generation of lampricides on non-target fish such as lake sturgeon.

DNA barcodes identify medically important tick species in Canada.

Medically important ticks (Acari: Ixodidae) are often difficult to identify morphologically. A standardized, molecular approach using a 658 base pair DNA barcode sequence (from the 5' region of the mitochondrial cytochrome c oxidase subunit I gene) was evaluated for its effectiveness in discriminating ticks in North America, with an emphasis on Canadian ticks. DNA barcodes were generated for 96 of 154 specimens representing 26 ixodid species. A genetic cluster analysis was performed on the barcode sequences, which separated specimens into haplogroups closely corresponding with morphologically identified species. The tree topology was further supported by a BIN analysis. COI sequences generated were found to have a mean maximum intraspecific divergence of 1.59% and a mean nearest neighbour divergence of 12.8%, indicating a significant "barcode gap". This study also revealed possible cryptic diversity among specimens morphologically identified as Ixodes soricis and Ixodes texanus. A PCR-based test for Borrelia burgdorferi determined that 18.1% of Lyme-competent ticks in this study were positive. This study is also the first to record a B. burgdorferi-positive exoskeleton. In conclusion, DNA barcoding is a powerful tool that clinicians can use to determine the identification of tick specimens which can help them to suggest whether an attached tick is a potential health risk.

DNA-based identification of invasive alien species in relation to Canadian federal policy and law, and the basis of rapid-response management.

Managing invasive alien species in Canada requires reliable taxonomic identification as the basis of rapid-response management. This can be challenging, especially when organisms are small and lack morphological diagnostic features. DNA-based techniques, such as DNA barcoding, offer a reliable, rapid, and inexpensive toolkit for taxonomic identification of individual or bulk samples, forensic remains, and even environmental DNA. Well suited for this requirement, they could be more broadly deployed and incorporated into the operating policy and practices of Canadian federal departments and should be authorized under these agencies' articles of law. These include Fisheries and Oceans Canada, Canadian Food Inspection Agency, Transport Canada, Environment Canada, Parks Canada, and Health Canada. These efforts should be harmonized with the appropriate provisions of provincial jurisdictions, for example, the Ontario Invasive Species Act. This approach necessitates that a network of accredited, certified laboratories exists, and that updated DNA reference libraries are readily accessible. Harmonizing this approach is vital among Canadian federal agencies, and between the federal and provincial levels of government. Canadian policy and law must also be harmonized with that of the USA when detecting, and responding to, invasive species in contiguous lands and waters. Creating capacity in legislation for use of DNA-based identifications brings the authority to fund, train, deploy, and certify staff, and to refine further developments in this molecular technology.

Genetic calibration of species diversity among North America's freshwater fishes.

Freshwater ecosystems are being heavily exploited and degraded by human activities all over the world, including in North America, where fishes and fisheries are strongly affected. Despite centuries of taxonomic inquiry, problems inherent to species identification continue to hamper the conservation of North American freshwater fishes. Indeed, nearly 10% of species diversity is thought to remain undescribed. To provide an independent calibration of taxonomic uncertainty and to establish a more accessible molecular identification key for its application, we generated a standard reference library of mtDNA sequences (DNA barcodes) derived from expert-identified museum specimens for 752 North American freshwater fish species. This study demonstrates that 90% of known species can be delineated using barcodes. Moreover, it reveals numerous genetic discontinuities indicative of independently evolving lineages within described species, which points to the presence of morphologically cryptic diversity. From the 752 species analyzed, our survey flagged 138 named species that represent as many as 347 candidate species, which suggests a 28% increase in species diversity. In contrast, several species of parasitic and nonparasitic lampreys lack such discontinuity and may represent alternative life history strategies within single species. Therefore, it appears that the current North American freshwater fish taxonomy at the species level significantly conceals diversity in some groups, although artificially creating diversity in others. In addition to providing an easily accessible digital identification system, this study identifies 151 fish species for which taxonomic revision is required.

Demographic trends in mixed Bemisia tabaci (Hemiptera: Aleyrodidae) cryptic species populations in commercial poinsettia under biological control- and insecticide-based management.

Bemisia tabaci (Gennadius) is an economically important pest of agricultural and ornamental plants worldwide and is now widely recognized as a cryptic species complex. In North America, B. tabaci is a particularly important pest of greenhouse poinsettia. In poinsettia production, two cryptic species from the B. tabaci complex, Mediterranean and Middle East Minor 1, often infest crops simultaneously. Differences in pesticide susceptibility between these two cryptic species have the potential to influence growers' management decisions, including the use of biological control or insecticides, and the choice of insecticide active ingredient. However, the demographic behavior of mixed-species infestations in commercial greenhouses has yet to be investigated. We conducted a survey of B. tabaci populations in commercial greenhouses in Ontario, Canada, and provide evidence that under biological control-based management, Middle East Minor 1 can displace Mediterranean, whereas under insecticide-based management Mediterranean populations will persist. Furthermore, we comment on implications of this behavior on the management of B. tabaci, and comment on methods used to identify B. tabaci cryptic species.

Delimiting Species with Single-Locus DNA Sequences.

DNA sequences are increasingly used for large-scale biodiversity inventories. Because these genetic data avoid the time-consuming initial sorting of specimens based on their phenotypic attributes, they have been recently incorporated into taxonomic workflows for overlooked and diverse taxa. Major statistical developments have accompanied this new practice, and several models have been proposed to delimit species with single-locus DNA sequences. However, proposed approaches to date make different assumptions regarding taxon lineage history, leading to strong discordance whenever comparisons are made among methods. Distance-based methods, such as Automatic Barcode Gap Discovery (ABGD) and Assemble Species by Automatic Partitioning (ASAP), rely on the detection of a barcode gap (i.e., the lack of overlap in the distributions of intraspecific and interspecific genetic distances) and the associated threshold in genetic distances. Network-based methods, as exemplified by the REfined Single Linkage (RESL) algorithm for the generation of Barcode Index Numbers (BINs), use connectivity statistics to hierarchically cluster-related haplotypes into molecular operational taxonomic units (MOTUs) which serve as species proxies. Tree-based methods, including Poisson Tree Processes (PTP) and the General Mixed Yule Coalescent (GMYC), fit statistical models to phylogenetic trees by maximum likelihood or Bayesian frameworks.Multiple webservers and stand-alone versions of these methods are now available, complicating decision-making regarding the most appropriate approach to use for a given taxon of interest. For instance, tree-based methods require an initial phylogenetic reconstruction, and multiple options are now available for this purpose such as RAxML and BEAST. Across all examined species delimitation methods, judicious parameter setting is paramount, as different model parameterizations can lead to differing conclusions. The objective of this chapter is to guide users step-by-step through all the procedures involved for each of these methods, while aggregating all necessary information required to conduct these analyses. The "Materials" section details how to prepare and format input files, including options to align sequences and conduct tree reconstruction with Maximum Likelihood and Bayesian inference. The Methods section presents the procedure and options available to conduct species delimitation analyses, including distance-, network-, and tree-based models. Finally, limits and future developments are discussed in the Notes section. Most importantly, species delimitation methods discussed herein are categorized based on five indicators: reliability, availability, scalability, understandability, and usability, all of which are fundamental properties needed for any approach to gain unanimous adoption within the DNA barcoding community moving forward.

A Measure of the DNA Barcode Gap for Applied and Basic Research.

DNA barcoding has largely established itself as a mainstay for rapid molecular taxonomic identification in both academic and applied research. The use of DNA barcoding as a molecular identification method depends on a "DNA barcode gap"-the separation between the maximum within-species difference and the minimum between-species difference. Previous work indicates the presence of a gap hinges on sampling effort for focal taxa and their close relatives. Furthermore, both theory and empirical work indicate a gap may not occur for related pairs of biological species. Here, we present a novel evaluation approach in the form of an easily calculated set of nonparametric metrics to quantify the extent of proportional overlap in inter- and intraspecific distributions of pairwise differences among target species and their conspecifics. The metrics are based on a simple count of the number of overlapping records for a species falling within the bounds of maximum intraspecific distance and minimum interspecific distance. Our approach takes advantage of the asymmetric directionality inherent in pairwise genetic distance distributions, which has not been previously done in the DNA barcoding literature. We apply the metrics to the predatory diving beetle genus Agabus as a case study because this group poses significant identification challenges due to its morphological uniformity despite both relative sampling ease and well-established taxonomy. Results herein show that target species and their nearest neighbor species were found to be tightly clustered and therefore difficult to distinguish. Such findings demonstrate that DNA barcoding can fail to fully resolve species in certain cases. Moving forward, we suggest the implementation of the proposed metrics be integrated into a common framework to be reported in any study that uses DNA barcoding for identification. In so doing, the importance of the DNA barcode gap and its components for the success of DNA-based identification using DNA barcodes can be better appreciated.

Fecal DNA metabarcoding helps characterize the Canada jay's diet and confirms its reliance on stored food for winter survival and breeding.

Accurately determining the diet of wild animals can be challenging if food items are small, visible only briefly, or rendered visually unidentifiable in the digestive system. In some food caching species, an additional challenge is determining whether consumed diet items have been previously stored or are fresh. The Canada jay (Perisoreus canadensis) is a generalist resident of North American boreal and subalpine forests with anatomical and behavioural adaptations allowing it to make thousands of arboreal food caches in summer and fall that are presumably responsible for its high winter survival and late winter/early spring breeding. We used DNA fecal metabarcoding to obtain novel information on nestling diets and compiled a dataset of 662 published and unpublished direct observations or stomach contents identifications of natural foods consumed by Canada jays throughout the year. We then used detailed natural history information to make informed decisions on whether each item identified to species in the diets of winter adults and nestlings was best characterized as 'likely cached', 'likely fresh' (i.e., was available as a non-cached item when it appeared in a jay's feces or stomach), or 'either possible'. Of the 87 food items consumed by adults in the winter, 39% were classified as 'likely cached' and 6% were deemed to be 'likely fresh'. For nestlings, 29% of 125 food items identified to species were 'likely cached' and 38% were 'likely fresh'. Our results support both the indispensability of cached food for Canada jay winter survival and previous suggestions that cached food is important for late winter/early spring breeding. Our work highlights the value of combining metabarcoding, stomach contents analysis, and direct observations to determine the cached vs. non-cached origins of consumed food items and the identity of food caches, some of which could be especially vulnerable to degradation through climate change.

Glacial cycles as an allopatric speciation pump in north-eastern American freshwater fishes.

Allopatric speciation may be the principal mechanism generating new species. Yet, it remains difficult to judge the generality of this process because few studies have provided evidence that geographic isolation has triggered the development of reproductive isolation over multiple species of a regional fauna. Here, we first combine results from new empirical data sets (7 taxa) and published literature (9 taxa) to show that the eastern Great Lakes drainage represents a multispecies suture zone for glacial lineages of freshwater fishes with variable levels of genetic divergence. Second, we performed amplified fragment length polymorphism analyses among four pairs of lineages. Results indicate that lineages with relatively deep levels of mtDNA 5' COI (barcode) sequence divergence (>2%) developed strong reproductive barriers, while lineages with lower levels of divergence show weaker reproductive isolation when found in sympatry. This suggests that a threshold of 2% sequence divergence at mtDNA could be used as a first step to flag cryptic species in North American freshwater fishes. By describing different levels of divergence and reproductive isolation in different co-occurring fishes, we offer strong evidence that allopatric speciation has contributed significantly to the diversification of north-eastern American freshwater fishes and confirm that Pleistocene glacial cycles can be viewed as a 'speciation pump' that played a predominant role in generating biodiversity.

Can DNA barcoding accurately discriminate megadiverse Neotropical freshwater fish fauna?

BACKGROUND: The megadiverse Neotropical freshwater ichthyofauna is the richest in the world with approximately 6,000 recognized species. Interestingly, they are distributed among only 17 orders, and almost 80% of them belong to only three orders: Characiformes, Siluriformes and Perciformes. Moreover, evidence based on molecular data has shown that most of the diversification of the Neotropical ichthyofauna occurred recently. These characteristics make the taxonomy and identification of this fauna a great challenge, even when using molecular approaches. In this context, the present study aimed to test the effectiveness of the barcoding methodology (COI gene) to identify the mega diverse freshwater fish fauna from the Neotropical region. For this purpose, 254 species of fishes were analyzed from the Upper Parana River basin, an area representative of the larger Neotropical region. RESULTS: Of the 254 species analyzed, 252 were correctly identified by their barcode sequences (99.2%). The main K2P intra- and inter-specific genetic divergence values (0.3% and 6.8%, respectively) were relatively low compared with similar values reported in the literature, reflecting the higher number of closely related species belonging to a few higher taxa and their recent radiation. Moreover, for 84 pairs of species that showed low levels of genetic divergence (<2%), application of a complementary character-based nucleotide diagnostic approach proved useful in discriminating them. Additionally, 14 species displayed high intra-specific genetic divergence (>2%), pointing to at least 23 strong candidates for new species. CONCLUSIONS: Our study is the first to examine a large number of freshwater fish species from the Neotropical area, including a large number of closely related species. The results confirmed the efficacy of the barcoding methodology to identify a recently radiated, megadiverse fauna, discriminating 99.2% of the analyzed species. The power of the barcode sequences to identify species, even with low interspecific divergence, gives us an idea of the distribution of inter-specific genetic divergence in these megadiverse fauna. The results also revealed hidden genetic divergences suggestive of reproductive isolation and putative cryptic speciation in some species (23 candidates for new species). Finally, our study constituted an important contribution to the international Barcoding of Life (iBOL.org) project, providing barcode sequences for use in identification of these species by experts and non-experts, and allowing them to be available for use in other applications.

Hidden diversity: DNA metabarcoding reveals hyper-diverse benthic invertebrate communities.

BACKGROUND: Freshwater ecosystems, such as streams, are facing increasing pressures from agricultural land use and recent literature stresses the importance of robust biomonitoring to detect trends in insect decline globally. Aquatic insects and other macroinvertebrates are often used as indicators of ecological condition in freshwater biomonitoring programs; however, these diverse groups can present challenges to morphological identification and coarse-level taxonomic resolution can mask patterns in community composition. Here, we incorporate molecular identification (DNA metabarcoding) into a stream biomonitoring sampling design to explore the diversity and variability of aquatic macroinvertebrate communities at small spatial scales. While individual stream reaches can be very heterogenous, most community ecology studies focus on larger, landscape-level patterns of community composition. A high degree of community variability at the local scale has important implications for both biomonitoring and ecological research, and the incorporation of DNA metabarcoding into local biodiversity assessments will inform future sampling protocols. RESULTS: We sampled twenty streams in southern Ontario, Canada, for aquatic macroinvertebrates across multiple time points and assessed local community variability by comparing field replicates taken ten meters apart within the same stream. Using bulk-tissue DNA metabarcoding, we revealed that aquatic macroinvertebrate communities are highly diverse at small spatial scales with unprecedented levels of local taxonomic turnover. We detected over 1600 Operational Taxonomic Units (OTUs) from 149 families, and a single insect family, the Chironomidae, contained over one third of the total number of OTUs detected in our study. Benthic communities were largely comprised of rare taxa detected only once per stream despite multiple biological replicates (24-94% rare taxa per site). In addition to numerous rare taxa, our species pool estimates indicated that there was a large proportion of taxa that remained undetected by our sampling regime (14-94% per site). Our sites were located across a gradient of agricultural activity, and while we predicted that increased land use would homogenize benthic communities, this was not supported as within-stream dissimilarity was unrelated to land use. Within-stream dissimilarity estimates were consistently high for all levels of taxonomic resolution (invertebrate families, invertebrate OTUs, chironomid OTUs), indicating stream communities are very dissimilar at small spatial scales.

Unveiling invasive insect threats to plant biodiversity: Leveraging eDNA metabarcoding and saturated salt trap solutions for biosurveillance.

The negative global impacts of invasive alien species (IAS) on biodiversity are second only to habitat loss. eDNA metabarcoding allows for a faster and more comprehensive evaluation of community species composition, with a higher taxonomic resolution and less taxonomic expertise required than traditional morphological-based biosurveillance. These advantages have positioned eDNA metabarcoding as the standard method for molecular-based detection of invasive alien species, where fast and accurate detectability allows prompt responses to mitigate their adverse effects. Here, eDNA metabarcoding is used for biosurveillance of invasive alien species regulated by Canada in high-risk areas with four main objectives: i) validate the effectiveness of eDNA metabarcoding of salt trap solutions as a molecular technique for IAS detection, ii) compare detection from DNA extracts obtained from filter quarters versus whole filters, iii) benchmark two different bioinformatic pipelines (MetaWorks and mBRAVE), and iv) compare canopy and ground level trapping. eDNA from up to five IAS (Agrilus planipennis, Daktulosphaira vitifoliae, Lymantria dispar, Popillia japonica, and Trichoferus campestris) were successfully detected across years from 2017 to 2022 in southern Ontario, Canada, with successful morphological validation for all except Lymantria dispar and Trichoferus campestris. Analysis of filter quarters in contrast to whole filters was demonstrated to be insufficient for effective IAS detection in each sample. All IAS were detected in only one filter quarter, suggesting a patchy eDNA distribution on the filter. The MetaWorks and mBRAVE bioinformatics pipelines proved effective in identifying IAS, with MetaWorks yielding a higher success rate when comparing molecular and morphological identifications. Ground-level and canopy-level sampling showed differential IAS recovery rates based on the molecular detection, which also varied per collection year, with all found IAS detected at the canopy level in 2022 while only one (Lymantria dispar) in 2020. The present study ratifies the efficacy and importance of eDNA-based detection in a regulatory context and the utility of adding eDNA metabarcoding of saturated salt trap solutions, a critical tool for IAS detection.

VLF: An R package for the analysis of very low frequency variants in DNA sequences.

Here, we introduce VLF, an R package to determine the distribution of very low frequency variants (VLFs) in nucleotide and amino acid sequences for the analysis of errors in DNA sequence records. The package allows users to assess VLFs in aligned and trimmed protein-coding sequences by automatically calculating the frequency of nucleotides or amino acids in each sequence position and outputting those that occur under a user-specified frequency (default of p = 0.001). These results can then be used to explore fundamental population genetic and phylogeographic patterns, mechanisms and processes at the microevolutionary level, such as nucleotide and amino acid sequence conservation. Our package extends earlier work pertaining to an implementation of VLF analysis in Microsoft Excel, which was found to be both computationally slow and error prone. We compare those results to our own herein. Results between the two implementations are found to be highly consistent for a large DNA barcode dataset of bird species. Differences in results are readily explained by both manual human error and inadequate Linnean taxonomy (specifically, species synonymy). Here, VLF is also applied to a subset of avian barcodes to assess the extent of biological artifacts at the species level for Canada goose (Branta canadensis), as well as within a large dataset of DNA barcodes for fishes of forensic and regulatory importance. The novelty of VLF and its benefit over the previous implementation include its high level of automation, speed, scalability and ease-of-use, each desirable characteristics which will be extremely valuable as more sequence data are rapidly accumulated in popular reference databases, such as BOLD and GenBank.