Environmental DNA enables detection of terrestrial mammals from forest pond water

Terrestrial animals must have frequent contact with water to survive, implying that environmental DNA (eDNA) originating from those animals should be detectable from places containing water in terrestrial ecosystems. Aiming to detect the presence of terrestrial mammals using forest water samples, we applied a set of universal PCR primers (MiMammal, a modified version of fish universal primers) for metabarcoding mammalian eDNA. The versatility of MiMammal primers was tested in silico and by amplifying DNAs extracted from tissues. The results suggested that MiMammal primers are capable of amplifying and distinguishing a diverse group of mammalian species. In addition, analyses of water samples from zoo cages of mammals with known species composition suggested that MiMammal primers could successfully detect mammalian species from water samples in the field. Then, we performed an experiment to detect mammals from natural ecosystems by collecting five 500‐ml water samples from ponds in two cool‐temperate forests in Hokkaido, northern Japan. MiMammal amplicon libraries were constructed using eDNA extracted from water samples, and sequences generated by Illumina MiSeq were subjected to data processing and taxonomic assignment. We thereby detected multiple species of mammals common to the sampling areas, including deer (Cervus nippon), mouse (Mus musculus), vole (Myodes rufocanus), raccoon (Procyon lotor), rat (Rattus norvegicus) and shrew (Sorex unguiculatus). Many previous applications of the eDNA metabarcoding approach have been limited to aquatic/semiaquatic systems, but the results presented here show that the approach is also promising even for forest mammal biodiversity surveys.     

Investigating the Potential Use of Environmental DNA (eDNA) for Genetic Monitoring of Marine Mammals

The exploitation of non-invasive samples has been widely used in genetic monitoring of terrestrial species. In aquatic ecosystems, non-invasive samples such as feces, shed hair or skin, are less accessible. However, the use of environmental DNA (eDNA) has recently been shown to be an effective tool for genetic monitoring of species presence in freshwater ecosystems. Detecting species in the marine environment using eDNA potentially offers a greater challenge due to the greater dilution, amount of mixing and salinity compared with most freshwater ecosystems. To determine the potential use of eDNA for genetic monitoring we used specific primers that amplify short mitochondrial DNA sequences to detect the presence of a marine mammal, the harbor porpoise, Phocoena phocoena, in a controlled environment and in natural marine locations. The reliability of the genetic detections was investigated by comparing with detections of harbor porpoise echolocation clicks by static acoustic monitoring devices. While we were able to consistently genetically detect the target species under controlled conditions, the results from natural locations were less consistent and detection by eDNA was less successful than acoustic detections. However, at one site we detected long-finned pilot whale, Globicephala melas, a species rarely sighted in the Baltic. Therefore, with optimization aimed towards processing larger volumes of seawater this method has the potential to compliment current visual and acoustic methods of species detection of marine mammals.     

Environmental DNA metabarcoding of wild flowers reveals diverse communities of terrestrial arthropods

Terrestrial arthropods comprise the most species‐rich communities on Earth, and grassland flowers provide resources for hundreds of thousands of arthropod species. Diverse grassland ecosystems worldwide are threatened by various types of environmental change, which has led to decline in arthropod diversity. At the same time, monitoring grassland arthropod diversity is time‐consuming and strictly dependent on declining taxonomic expertise. Environmental DNA (eDNA) metabarcoding of complex samples has demonstrated that information on species compositions can be efficiently and non‐invasively obtained. Here, we test the potential of wild flowers as a novel source of arthropod eDNA. We performed eDNA metabarcoding of flowers from several different plant species using two sets of generic primers, targeting the mitochondrial genes 16S rRNA and COI. Our results show that terrestrial arthropod species leave traces of DNA on the flowers that they interact with. We obtained eDNA from at least 135 arthropod species in 67 families and 14 orders, together representing diverse ecological groups including pollinators, parasitoids, gall inducers, predators, and phytophagous species. Arthropod communities clustered together according to plant species. Our data also indicate that this experiment was not exhaustive, and that an even higher arthropod richness could be obtained using this eDNA approach. Overall, our results demonstrate that it is possible to obtain information on diverse communities of insects and other terrestrial arthropods from eDNA metabarcoding of wild flowers. This novel source of eDNA represents a vast potential for addressing fundamental research questions in ecology, obtaining data on cryptic and unknown species of plant‐associated arthropods, as well as applied research on pest management or conservation of endangered species such as wild pollinators.     

Comparing diversity levels in environmental samples: DNA sequence capture and metabarcoding approaches using 18S and COI genes

Environmental DNA studies targeting multiple taxa using metabarcoding provide remarkable insights into levels of species diversity in any habitat. The main drawbacks are the presence of primer bias and difficulty in identifying rare species. We tested a DNA sequence‐capture method in parallel with the metabarcoding approach to reveal possible advantages of one method over the other. Both approaches were performed using the same eDNA samples and the same 18S and COI regions, followed by high throughput sequencing. Metabarcoded eDNA libraries were PCR amplified with one primer pair from 18S and COI genes. DNA sequence‐capture libraries were enriched with 3,639 baits targeting the same gene regions. We tested amplicon sequence variants (ASVs) and operational taxonomic units (OTUs) in silico approaches for both markers and methods, using for this purpose the metabarcoding data set. ASVs methods uncovered more species for the COI gene, whereas the opposite occurred for the 18S gene, suggesting that clustering reads into OTUs could bias diversity richness especially using 18S with relaxed thresholds. Additionally, metabarcoding and DNA sequence‐capture recovered 80%–90% of the control sample species. DNA sequence‐capture was 8x more expensive, nonetheless it identified 1.5x more species for COI and 13x more genera for 18S than metabarcoding. Both approaches offer reliable results, sharing ca. 40% species and 72% families and retrieve more taxa when nuclear and mitochondrial markers are combined. eDNA metabarcoding is quite well established and low‐cost, whereas DNA‐sequence capture for biodiversity assessment is still in its infancy, is more time‐consuming but provides more taxonomic assignments.     

Targeted multiplex next‐generation sequencing: advances in techniques of mitochondrial and nuclear DNA sequencing for population genomics

Next‐generation sequencing (NGS) is emerging as an efficient and cost‐effective tool in population genomic analyses of nonmodel organisms, allowing simultaneous resequencing of many regions of multi‐genomic DNA from multiplexed samples. Here, we detail our synthesis of protocols for targeted resequencing of mitochondrial and nuclear loci by generating indexed genomic libraries for multiplexing up to 100 individuals in a single sequencing pool, and then enriching the pooled library using custom DNA capture arrays. Our use of DNA sequence from one species to capture and enrich the sequencing libraries of another species (i.e. cross‐species DNA capture) indicates that efficient enrichment occurs when sequences are up to about 12% divergent, allowing us to take advantage of genomic information in one species to sequence orthologous regions in related species. In addition to a complete mitochondrial genome on each array, we have included between 43 and 118 nuclear loci for low‐coverage sequencing of between 18 kb and 87 kb of DNA sequence per individual for single nucleotide polymorphisms discovery from 50 to 100 individuals in a single sequencing lane. Using this method, we have generated a total of over 500 whole mitochondrial genomes from seven cetacean species and green sea turtles. The greater variation detected in mitogenomes relative to short mtDNA sequences is helping to resolve genetic structure ranging from geographic to species‐level differences. These NGS and analysis techniques have allowed for simultaneous population genomic studies of mtDNA and nDNA with greater genomic coverage and phylogeographic resolution than has previously been possible in marine mammals and turtles.     

An empirical evaluation of the African elephant as a focal species for connectivity planning in East Africa

Aim Large, charismatic and wide‐ranging animals are often employed as focal species for prioritizing landscape linkages in threatened ecosystems (i.e. ‘connectivity conservation’), but there have been few efforts to assess empirically whether focal species co‐occur with other species of conservation interest within potential linkages. We evaluated whether the African elephant (Loxodonta africana), a world‐recognized flagship species, would serve as an appropriate focal species for other large mammals in a potential linkage between two major protected area complexes. Location A 15,400 km² area between the Ruaha and Selous ecosystems in central Tanzania, East Africa. Methods We used walking transects to assess habitat, human activity and co‐occurrence of elephants and 48 other large mammal species (> 1 kg) at 63 sites using animal sign and direct sightings. We repeated a subset of transects to estimate species detectability using occupancy modelling. We used logistic regression and AIC model selection to characterize patterns of elephant occurrence and assessed correlation of elephant presence with richness of large mammals and subgroups. We considered other possible focal species, compared habitat‐based linear regression models of large mammal richness and used circuit theory to examine potential connectivity spatially. Results Elephants were detected in many locations across the potential linkage. Elephant presence was highly positively correlated with the richness of large mammals, as well as ungulates, carnivores, large carnivores and species > 45 kg in body mass (‘megafauna’). Outside of protected areas, both mammal richness and elephant presence were negatively correlated with human population density and distance from water. Only one other potential focal species was more strongly correlated with species richness than elephants, but detectability was highest for elephants. Main conclusions Although African elephants have dispersal abilities that exceed most other terrestrial mammals, conserving elephant movement corridors may effectively preserve habitat and potential landscape linkages for other large mammal species among Tanzanian reserves.     

Detection of the endangered European weather loach (Misgurnus fossilis) via water and sediment samples: Testing multiple eDNA workflows

The European weather loach (Misgurnus fossilis) is classified as highly endangered in several countries of Central Europe. Populations of M. fossilis are predominantly found in ditches with low water levels and thick sludge layers and are thus hard to detect using conventional fishing methods. Therefore, environmental DNA (eDNA) monitoring appears particularly relevant for this species. In previous studies, M. fossilis was surveyed following eDNA water sampling protocols, which were not optimized for this species. Therefore, we created two full factorial study designs to test six different eDNA workflows for sediment samples and twelve different workflows for water samples. We used qPCR to compare the threshold cycle (C_t) values of the different workflows, which indicate the target DNA amount in the sample, and spectrophotometry to quantify and compare the total DNA amount inside the samples. We analyzed 96 water samples and 48 sediment samples from a pond with a known population of M. fossilis. We tested several method combinations for long‐term sample preservation, DNA capture, and DNA extraction. Additionally, we analyzed the DNA yield of samples from a ditch with a natural M. fossilis population monthly over one year to determine the optimal sampling period. Our results showed that the long‐term water preservation method commonly used for eDNA surveys of M. fossilis did not lead to optimal DNA yields, and we present a valid long‐term sample preservation alternative. A cost‐efficient high salt DNA extraction led to the highest target DNA yields and can be used for sediment and water samples. Furthermore, we were able to show that in a natural habitat of M. fossilis, total and target eDNA were higher between June and September, which implies that this period is favorable for eDNA sampling. Our results will help to improve the reliability of future eDNA surveys of M. fossilis.     

Combined hybridization capture and shotgun sequencing for ancient DNA analysis of extinct wild and domestic dromedary camel

The performance of hybridization capture combined with next‐generation sequencing (NGS) has seen limited investigation with samples from hot and arid regions until now. We applied hybridization capture and shotgun sequencing to recover DNA sequences from bone specimens of ancient‐domestic dromedary (Camelus dromedarius) and its extinct ancestor, the wild dromedary from Jordan, Syria, Turkey and the Arabian Peninsula, respectively. Our results show that hybridization capture increased the percentage of mitochondrial DNA (mtDNA) recovery by an average 187‐fold and in some cases yielded virtually complete mitochondrial (mt) genomes at multifold coverage in a single capture experiment. Furthermore, we tested the effect of hybridization temperature and time by using a touchdown approach on a limited number of samples. We observed no significant difference in the number of unique dromedary mtDNA reads retrieved with the standard capture compared to the touchdown method. In total, we obtained 14 partial mitochondrial genomes from ancient‐domestic dromedaries with 17–95% length coverage and 1.27–47.1‐fold read depths for the covered regions. Using whole‐genome shotgun sequencing, we successfully recovered endogenous dromedary nuclear DNA (nuDNA) from domestic and wild dromedary specimens with 1–1.06‐fold read depths for covered regions. Our results highlight that despite recent methodological advances, obtaining ancient DNA (aDNA) from specimens recovered from hot, arid environments is still problematic. Hybridization protocols require specific optimization, and samples at the limit of DNA preservation need multiple replications of DNA extraction and hybridization capture as has been shown previously for Middle Pleistocene specimens.     

Detection of invasive freshwater fish species using environmental DNA survey

Invasive species are one of the most significant problem in freshwater ecosystems. Most common non-native freshwater species in Turkish freshwater fish fauna are Prussian Carp (Carassius gibelio), North African Catfish (Clarias gariepinus), Nile Tilapia (Oreochromis niloticus) and Topmouth Gudgeon (Pseudorasbora parva).Recent studies showed that environmental DNA could be used to detect target species inhabiting the ecosystem with higher precision and less effort compared to traditional field surveys. In this study, eDNA approach was used to investigate non-native freshwater fish species from fifteen different locations of Upper Sakarya Basin. eDNA was successfully extracted from the water samples of locations where the species were visually observed. Mean amplification rate of eDNA was calculated as 77.03%.This study is the first environmental DNA study used in detection of four of the most common invasive freshwater fish species. Results clearly indicating that eDNA surveys could be used as an important molecular tool to monitor invasive fish species in freshwater ecosystems.     

Quantification of mesocosm fish and amphibian species diversity via environmental DNA metabarcoding

Freshwater fauna are particularly sensitive to environmental change and disturbance. Management agencies frequently use fish and amphibian biodiversity as indicators of ecosystem health and a way to prioritize and assess management strategies. Traditional aquatic bioassessment that relies on capture of organisms via nets, traps and electrofishing gear typically has low detection probabilities for rare species and can injure individuals of protected species. Our objective was to determine whether environmental DNA (eDNA) sampling and metabarcoding analysis can be used to accurately measure species diversity in aquatic assemblages with differing structures. We manipulated the density and relative abundance of eight fish and one amphibian species in replicated 206‐L mesocosms. Environmental DNA was filtered from water samples, and six mitochondrial gene fragments were Illumina‐sequenced to measure species diversity in each mesocosm. Metabarcoding detected all nine species in all treatment replicates. Additionally, we found a modest, but positive relationship between species abundance and sequencing read abundance. Our results illustrate the potential for eDNA sampling and metabarcoding approaches to improve quantification of aquatic species diversity in natural environments and point the way towards using eDNA metabarcoding as an index of macrofaunal species abundance.     

The Use of Genus-Specific Amplicon Pyrosequencing to Assess Phytophthora Species Diversity Using eDNA from Soil and Water in Northern Spain

Phytophthora is one of the most important and aggressive plant pathogenic genera in agriculture and forestry. Early detection and identification of its pathways of infection and spread are of high importance to minimize the threat they pose to natural ecosystems. eDNA was extracted from soil and water from forests and plantations in the north of Spain. Phytophthora-specific primers were adapted for use in high-throughput Sequencing (HTS). Primers were tested in a control reaction containing eight Phytophthora species and applied to water and soil eDNA samples from northern Spain. Different score coverage threshold values were tested for optimal Phytophthora species separation in a custom-curated database and in the control reaction. Clustering at 99% was the optimal criteria to separate most of the Phytophthora species. Multiple Molecular Operational Taxonomic Units (MOTUs) corresponding to 36 distinct Phytophthora species were amplified in the environmental samples. Pyrosequencing of amplicons from soil samples revealed low Phytophthora diversity (13 species) in comparison with the 35 species detected in water samples. Thirteen of the MOTUs detected in rivers and streams showed no close match to sequences in international sequence databases, revealing that eDNA pyrosequencing is a useful strategy to assess Phytophthora species diversity in natural ecosystems.     

Rapid degradation of longer DNA fragments enables the improved estimation of distribution and biomass using environmental DNA

The advent of environmental DNA (eDNA) analysis methods has enabled rapid and wide‐range ecological monitoring in aquatic ecosystems, but there is a dearth of information on eDNA degradation. The results of previous studies suggest that the decay rate of eDNA varies depending on the length of DNA fragments. To examine this hypothesis, we compared temporal change in copy number of long eDNA fragments (719 bp) with that of short eDNA fragments (127 bp). First, we isolated rearing water from a target fish species, Japanese Jack Mackerel (Trachurus japonicus), and then quantified the copy number of the long and short eDNA fragments in 1 L water samples after isolating the water from the fish. Long DNA fragments showed a higher decay rate than short fragments. Next, we measured the eDNA copy numbers of long and short DNA fragments using field samples, and compared them with fish biomass as measured by echo intensity. Although a previous study suggested that short eDNA fragments could be overestimated because of nontarget eDNA from a nearby fish market and carcasses, the eDNA concentrations of long fragments were correlated with echo intensity. This suggests that the concentration of longer eDNA fragments reflects fish biomass more accurately than the previous study by removing the effects of the fish market and carcasses. The length‐related differences in eDNA have a substantial potential to improve estimation of species biomass.     

Optimal extraction methods for the simultaneous analysis of DNA from diverse organisms and sample types

Using environmental DNA (eDNA) to assess the distribution of micro‐ and macroorganisms is becoming increasingly popular. However, the comparability and reliability of these studies is not well understood as we lack evidence on how different DNA extraction methods affect the detection of different organisms, and how this varies among sample types. Our aim was to quantify biases associated with six DNA extraction methods and identify one which is optimal for eDNA research targeting multiple organisms and sample types. We assessed each methods’ ability to simultaneously extract bacterial, fungal, plant, animal and fish DNA from soil, leaf litter, stream water, stream sediment, stream biofilm and kick‐net samples, as well as from mock communities. Method choice affected alpha‐diversity for several combinations of taxon and sample type, with the majority of the differences occurring in the bacterial communities. While a single method performed optimally for the extraction of DNA from bacterial, fungal and plant mock communities, different methods performed best for invertebrate and fish mock communities. The consistency of methods, as measured by the similarity of community compositions resulting from replicate extractions, varied and was lowest for the animal communities. Collectively, these data provide the first comprehensive assessment of the biases associated with DNA extraction for both different sample types and taxa types, allowing us to identify DNeasy PowerSoil as a universal DNA extraction method. The adoption of standardized approaches for eDNA extraction will ensure that results can be more reliably compared, and biases quantified, thereby advancing eDNA as an ecological research tool.     

Impact of enrichment conditions on cross‐species capture of fresh and degraded DNA

By combining high‐throughput sequencing with target enrichment (‘hybridization capture’), researchers are able to obtain molecular data from genomic regions of interest for projects that are otherwise constrained by sample quality (e.g. degraded and contamination‐rich samples) or a lack of a priori sequence information (e.g. studies on nonmodel species). Despite the use of hybridization capture in various fields of research for many years, the impact of enrichment conditions on capture success is not yet thoroughly understood. We evaluated the impact of a key parameter – hybridization temperature – on the capture success of mitochondrial genomes across the carnivoran family Felidae. Capture was carried out for a range of sample types (fresh, archival, ancient) with varying levels of sequence divergence between bait and target (i.e. across a range of species) using pools of individually indexed libraries on Agilent SureSelect^? arrays. Our results suggest that hybridization capture protocols require specific optimization for the sample type that is being investigated. Hybridization temperature affected the proportion of on‐target sequences following capture: for degraded samples, we obtained the best results with a hybridization temperature of 65 °C, while a touchdown approach (65 °C down to 50 °C) yielded the best results for fresh samples. Evaluation of capture performance at a regional scale (sliding window approach) revealed no significant improvement in the recovery of DNA fragments with high sequence divergence from the bait at any of the tested hybridization temperatures, suggesting that hybridization temperature may not be the critical parameter for the enrichment of divergent fragments.     

Moving eDNA surveys onto land: Strategies for active eDNA aggregation to detect invasive forest insects

The use of environmental DNA (eDNA) surveys to monitor terrestrial species has been relatively limited, with successful implementations still confined to sampling DNA from natural or artificial water bodies and soil. Sampling water for eDNA depends on proximity to or availability of water, whereas eDNA from soil is limited in its spatial scale due to the large quantities necessary for processing and difficulty in doing so. These challenges limit the widespread use of eDNA in several systems, such as surveying forests for invasive insects. We developed two new eDNA aggregation approaches that overcome the challenges of above‐ground terrestrial sampling and eliminate the dependency on creating or utilizing pre‐existing water bodies to conduct eDNA sampling. The first, “spray aggregation,” uses spray action to remove eDNA from surface substrates and was developed for shrubs and other understorey vegetation, while the second, “tree rolling,” uses physical transfer via a roller to remove eDNA from the surface of tree trunks and large branches. We tested these approaches by surveying for spotted lanternfly, Lycorma delicatula, a recent invasive pest of northeastern USA that is considered a significant ecological and economic threat to forests and agriculture. We found that our terrestrial eDNA surveys matched visual surveys, but also detected L. delicatula presence ahead of visual surveys, indicating increased sensitivity of terrestrial eDNA surveys over currently used methodology. The terrestrial eDNA approaches we describe can be adapted for use in surveying a variety of forest insects and represent a novel strategy for surveying terrestrial biodiversity.     

Environmental DNA metabarcoding: Transforming how we survey animal and plant communities

The genomic revolution has fundamentally changed how we survey biodiversity on earth. High‐throughput sequencing (“HTS”) platforms now enable the rapid sequencing of DNA from diverse kinds of environmental samples (termed “environmental DNA” or “eDNA”). Coupling HTS with our ability to associate sequences from eDNA with a taxonomic name is called “eDNA metabarcoding” and offers a powerful molecular tool capable of noninvasively surveying species richness from many ecosystems. Here, we review the use of eDNA metabarcoding for surveying animal and plant richness, and the challenges in using eDNA approaches to estimate relative abundance. We highlight eDNA applications in freshwater, marine and terrestrial environments, and in this broad context, we distill what is known about the ability of different eDNA sample types to approximate richness in space and across time. We provide guiding questions for study design and discuss the eDNA metabarcoding workflow with a focus on primers and library preparation methods. We additionally discuss important criteria for consideration of bioinformatic filtering of data sets, with recommendations for increasing transparency. Finally, looking to the future, we discuss emerging applications of eDNA metabarcoding in ecology, conservation, invasion biology, biomonitoring, and how eDNA metabarcoding can empower citizen science and biodiversity education.     

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

Species detection using eDNA is revolutionizing global capacity to monitor biodiversity. However, the lack of regional, vouchered, genomic sequence information—especially sequence information that includes intraspecific variation—creates a bottleneck for management agencies wanting to harness the complete power of eDNA to monitor taxa and implement eDNA analyses. eDNA studies depend upon regional databases of mitogenomic sequence information to evaluate the effectiveness of such data to detect and identify taxa. We created the Oregon Biodiversity Genome Project to create a database of complete, nearly error-free mitogenomic sequences for all of Oregon's fishes. We have successfully assembled the complete mitogenomes of 313 specimens of freshwater, anadromous and estuarine fishes representing 24 families, 55 genera and 129 species and lineages. Comparative analyses of these sequences illustrate that many regions of the mitogenome are taxonomically informative, that the short (~150 bp) mitochondrial ‘barcode’ regions typically used for eDNA assays do not consistently diagnose for species and that complete single or multiple genes of the mitogenome are preferable for identifying Oregon's fishes. This project provides a blueprint for other researchers to follow as they build regional databases, illustrates the taxonomic value and limits of complete mitogenomic sequences and offers clues as to how current eDNA assays and environmental genomics methods of the future can best leverage this information.     

Responses of foraminifera communities to aquaculture‐derived organic enrichment as revealed by environmental DNA metabarcoding

Current monitoring methods to assess benthic impacts of marine finfish aquaculture are based on complex biological indices and/or geochemistry data. The former requires benthic macrofauna morpho‐taxonomic characterization that is time‐ and cost‐intensive, while the latter provides rapid assessment of the organic enrichment status of sediments but does not directly measure biotic impacts. In this study, sediment samples were collected from seven stations at six salmon farms in British Columbia, Canada, and analyzed for geochemical parameters and by eDNA metabarcoding to investigate linkages between geochemistry and foraminifera. Sediment texture across farm sites ranged from sand to silty loam, while the maximum sediment pore‐water sulphide concentration at each site ranged from 1,000 to 13,000 μM. Foraminifera alpha diversity generally increased with distance from cage edge. Adonis analyses revealed that farm site explained the most variation in foraminifera community, followed by sediment type, enrichment status, and distance from cage edge. Farm‐specific responses were observed in diversity analyses, taxonomic difference analyses, and correlation analyses. Results demonstrated that species diversity and composition of foraminifera characterized by eDNA metabarcoding generated signals consistent with benthic biodiversity being impacted by finfish farming activities. This substantiates the validity of eDNA metabarcoding for augmenting current approaches to benthic impact assessments by providing more cost‐effective and practicable biotic measures than traditional morpho‐taxonomy. To capitalize on this potential, further work is needed to design a new nomogram that combines eDNA metabarcoding data and geochemistry data to enable accurate monitoring of benthic impacts of fish farming in a time‐ and cost‐efficient way.     

Environmental DNA metabarcoding of lake fish communities reflects long‐term data from established survey methods

Organisms continuously release DNA into their environments via shed cells, excreta, gametes and decaying material. Analysis of this ‘environmental DNA’ (eDNA) is revolutionizing biodiversity monitoring. eDNA outperforms many established survey methods for targeted detection of single species, but few studies have investigated how well eDNA reflects whole communities of organisms in natural environments. We investigated whether eDNA can recover accurate qualitative and quantitative information about fish communities in large lakes, by comparison to the most comprehensive long‐term gill‐net data set available in the UK. Seventy‐eight 2L water samples were collected along depth profile transects, gill‐net sites and from the shoreline in three large, deep lakes (Windermere, Bassenthwaite Lake and Derwent Water) in the English Lake District. Water samples were assayed by eDNA metabarcoding of the mitochondrial 12S and cytochrome b regions. Fourteen of the 16 species historically recorded in Windermere were detected using eDNA, compared to four species in the most recent gill‐net survey, demonstrating eDNA is extremely sensitive for detecting species. A key question for biodiversity monitoring is whether eDNA can accurately estimate abundance. To test this, we used the number of sequence reads per species and the proportion of sampling sites in which a species was detected with eDNA (i.e. site occupancy) as proxies for abundance. eDNA abundance data consistently correlated with rank abundance estimates from established surveys. These results demonstrate that eDNA metabarcoding can describe fish communities in large lakes, both qualitatively and quantitatively, and has great potential as a complementary tool to established monitoring methods.