Creating, curating and evaluating a mitogenomic reference database to improve regional species identification using environmental DNA

Environmental DNA: The next chapter

Molecular tools are an indispensable part of ecology and biodiversity sciences and implemented across all biomes. About a decade ago, the use and implementation of environmental DNA (eDNA) to detect biodiversity signals extracted from environmental samples opened new avenues of research. Initial eDNA research focused on understanding population dynamics of target species. Its scope thereafter broadened, uncovering previously unrecorded biodiversity via metabarcoding in both well-studied and understudied ecosystems across all taxonomic groups. The application of eDNA rapidly became an established part of biodiversity research, and a research field by its own. Here, we revisit key expectations made in a land-mark special issue on eDNA in Molecular Ecology in 2012 to frame the development in six key areas: (1) sample collection, (2) primer development, (3) biomonitoring, (4) quantification, (5) behaviour of DNA in the environment and (6) reference database development. We pinpoint the success of eDNA, yet also discuss shortfalls and expectations not met, highlighting areas of research priority and identify the unexpected developments. In parallel, our retrospective couples a screening of the peer-reviewed literature with a survey of eDNA users including academics, end-users and commercial providers, in which we address the priority areas to focus research efforts to advance the field of eDNA. With the rapid and ever-increasing pace of new technical advances, the future of eDNA looks bright, yet successful applications and best practices must become more interdisciplinary to reach its full potential. Our retrospect gives the tools and expectations towards concretely moving the field forward.

Using low volume eDNA methods to sample pelagic marine animal assemblages

Environmental DNA (eDNA) is an increasingly useful method for detecting pelagic animals in the ocean but typically requires large water volumes to sample diverse assemblages. Ship-based pelagic sampling programs that could implement eDNA methods generally have restrictive water budgets. Studies that quantify how eDNA methods perform on low water volumes in the ocean are limited, especially in deep-sea habitats with low animal biomass and poorly described species assemblages. Using 12S rRNA and COI gene primers, we quantified assemblages comprised of micronekton, coastal forage fishes, and zooplankton from low volume eDNA seawater samples (n = 436, 380-1800 mL) collected at depths of 0-2200 m in the southern California Current. We compared diversity in eDNA samples to concurrently collected pelagic trawl samples (n = 27), detecting a higher diversity of vertebrate and invertebrate groups in the eDNA samples. Differences in assemblage composition could be explained by variability in size-selectivity among methods and DNA primer suitability across taxonomic groups. The number of reads and amplicon sequences variants (ASVs) did not vary substantially among shallow (<200 m) and deep samples (>600 m), but the proportion of invertebrate ASVs that could be assigned a species-level identification decreased with sampling depth. Using hierarchical clustering, we resolved horizontal and vertical variability in marine animal assemblages from samples characterized by a relatively low diversity of ecologically important species. Low volume eDNA samples will quantify greater taxonomic diversity as reference libraries, especially for deep-dwelling invertebrate species, continue to expand.

Population genomics of the invasive Northern Giant Hornet Vespa mandarinia in North America and across its native range

The northern giant hornet Vespa mandarinia (NGH) is a voracious predator of other insect species, including honey bees. NGH's native range spans subtropical and temperate regions across much of east and southeast Asia and, in 2019, exotic populations of the species were discovered in North America. Despite this broad range and invasive potential, investigation of the population genomic structure of NGH across its native and introduced ranges has thus far been limited to a small number of mitochondrial samples. Here, we present analyses of genomic data from NGH individuals collected across the species' native range and from exotic individuals collected in North America. We provide the first survey of whole-genome population variation for any hornet species, covering this species' native and invasive ranges, and in doing so confirm likely origins in Japan and South Korea for the two introductions. We additionally show that, while this introduced population exhibited strongly elevated levels of inbreeding, these signatures of inbreeding are also present in some long-standing native populations, which may indicate that inbreeding depression alone is insufficient to prevent the persistence of NGH populations. As well as highlighting the importance of ongoing monitoring and eradication efforts to limit the spread of this species outside of its natural range, our data will serve as a foundational database for future genomic studies into introduced hornet populations.

rCRUX: A Rapid and Versatile Tool for Generating Metabarcoding Reference libraries in R

The sequencing revolution requires accurate taxonomic classification of DNA sequences. Key to making accurate taxonomic assignments are curated, comprehensive reference barcode databases. However, the generation and curation of such databases has remained challenging given the large and continuously growing volumes of both DNA sequence data and novel reference barcode targets. Monitoring and research applications require a greater diversity of specialized gene regions and targeted taxa then are currently curated by professional staff. Thus there is a growing need for an easy to implement computational tool that can generate comprehensive metabarcoding reference libraries for any bespoke locus. We address this need by reimagining CRUX from the Anacapa Toolkit and present the rCRUX package in R which, like it's predecessor, relies on sequence homology and PCR primer compatibility instead of keyword-searches to avoid limitations of user-defined metadata. The typical workflow involves searching for plausible seed amplicons (get_seeds_local() or get_seeds_remote()) by simulating in silico PCR to acquire a set of sequences analogous to PCR products containing a user-defined set of primer sequences. Next, these seeds are used to iteratively blast search seed sequences against a local copy of the National Center for Biotechnology Information (NCBI) formatted nt database using a taxonomic-rank based stratified random sampling approach ( blast_seeds() ). This results in a comprehensive set of sequence matches. This database is dereplicated and cleaned (derep_and_clean_db()) by identifying identical reference sequences and collapsing the taxonomic path to the lowest taxonomic agreement across all matching reads. This results in a curated, comprehensive database of primer-specific reference barcode sequences from NCBI. Databases can then be compared (compare_db()) to determine read and taxonomic overlap. We demonstrate that rCRUX provides more comprehensive reference databases for the MiFish Universal Teleost 12S, Taberlet trnl, fungal ITS, and Leray CO1 loci than CRABS, MetaCurator, RESCRIPt, and ecoPCR reference databases. We then further demonstrate the utility of rCRUX by generating 24 reference databases for 20 metabarcoding loci, many of which lack dedicated reference database curation efforts. The rCRUX package provides a simple to use tool for the generation of curated, comprehensive reference databases for user-defined loci, facilitating accurate and effective taxonomic classification of metabarcoding and DNA sequence efforts broadly.