parole chiave: DNA, ambienti marini
Our results show evidence for a strong spatial effect of eDNA degradation in the natural marine environment, with eDNA degrading around 1.6 times faster in the terrestrially influenced inshore environment than the ocean-influenced offshore environment. We found that eDNA also degraded around 1.1 times faster in winter than in summer, although this difference was not statistically significant. These results placed in the context of our review of eDNA decay rates in the literature (Table 4), appear to contradict the notion that eDNA degrades faster in marine environments than freshwater 18,20,29. In fact, degradation rates appear to be slower in many cases, with only marine studies or the freshwater studies at low temperature or using non-natural water sources, having a half-life of greater than 30 h (Table 4). The fastest rates in freshwater systems assessed so far are of acidic stream environments (<1.2 h of half-life5), while the fastest marine decay rate was 6.9 h, from anchovy eDNA in Californian inshore waters at 22 °C20. Most marine eDNA decay rates appear, however, to have been estimated at between 10 and 50 h, and with the lowest rates corresponding to the coldest water temperatures: 63 h at 4 °C19 and 71 h at −1 °C15. Rates above 71 h were from freshwater studies using sanitised or purified water from non-natural sources (Table 4).
Compared with freshwater, marine systems are generally characterised by higher salinity and ionic content, typically higher pH, and more stable temperatures. These are factors which have been shown to promote DNA preservation, and tend to correspond to the lowest observed degradation rates5,14,23,25,30,31. Our artificial spatial gradient varied from an offshore treatment with high pH and salinity to an inshore treatment with a lower salinity and a slightly lower pH. This was designed to capture the abiotic heterogeneity that could be expected across the Western English Channel region over the period of a year, a magnitude of variation that will apply to other coastal temperate locations. We found salinity to be a better predictor of eDNA decay than pH, and with salinity varying more between locations and pH varying more over seasons (Table 3), this agrees with the finding that the spatial signal was stronger than the temporal signal, and is reflected in the correlation matrix of the combined predictor-covariate model.
The lack of a statistically significant difference over season may be due to the relatively low degree of variation in pH and temperature. Seawater pH measured in this experiment was between around 8 and 8.6, which may not have any direct impact on DNA hydrolysis, and likewise, temperature ranges in this temperate marine system (10–15 °C) were narrower than those typically studied in terrestrial systems (e.g. 5–35 °C12).
As well as abiotic factors engaging in DNA degradation via oxidisation and hydrolysis by depurination, biotic factors are also likely to play a major role in eDNA persistence dynamics via extracellular DNases produced by heterotrophic microbes4,30. While we found support for faster degradation rates in our inshore and mixed treatments (Fig. 3), this difference did not appear to be proportional to the quantity of inshore water used in the treatment—the two-third offshore treatment tended to be closer to the 100% inshore treatment than the 100% offshore treatment—suggesting that biotic rather than abiotic factors are of stronger influence. Salinity itself may not be therefore entirely responsible for the difference in decay rate, rather that it is associated with particular abundances or communities of microbes. Gilbert et al. 32 reported that microbial community structure in the Western English Channel was highly dynamic seasonally. Free DNA is thought to represent an important organic phosphorus source in marine systems 29, and seasonal phosphate limitation has been identified as a key driver of eDNA turnover rates over abiotic factors such as temperature, pH and salinity 26.
Therefore, the lack of seasonal difference in eDNA decay that we report may also be explained by organic phosphorous or carbon concentrations 14,26.
Taken together with the survey of rates from the literature, this implies that abiotic and biotic factors are co-implicated in eDNA degradation. Assessing the covariance and contribution among these parameters is an area that needs to be addressed, along with more sophisticated analyses of microbial communities incorporating a greater degree of spatial replication. A number of systematic biases were identified as being potentially problematic for our inferences. PCR inhibition in the samples from the inshore site could explain the faster degradation rates from that location. However, we assessed amplification efficiency of the qPCR in a serial dilution experiment, and these were near the expected 100% across treatments and season. Values well above 100% would indicate inhibition. Other studies have also indicated low instances of PCR inhibition when using kits with dedicated inhibitor removal steps such as the PowerWater kit that we used 33,34. Although not significantly different, we found that degradation rates were overall around 1.2 times faster in the crab assay than the shanny assay. This is most likely explained by differences in fragment length between the two assays (153 vs. 132 bp), with longer fragments being shown to decay at a faster rate than shorter fragments 22. It was also noted that despite using a similar mass of crabs and shannys to create the eDNA, initial measured concentrations were roughly an order of magnitude lower in the crab assay (Fig. 1), perhaps indicating that the exoskeleton of the crustaceans, as well as their behaviour and breeding condition at particular times of the year may limit eDNA output 35.
The treatment of qPCR non-amplifications in low-template analyses is an important source of error at the analytical stage. Due to the proportion of non-amplifications at the 192-h sample (0.56 in winter and 0.68 in summer)— i.e. outside of the experimental limit of quantification—and the influence of this time point in estimating the regression slopes, our eDNA decay model was sensitive to how these missing data were treated. Excluding them, or fixing them to the limit-of-detection value resulted in the effects of season and assay becoming statistically significant. However, treating the non-amplifications this way is problematic as these missing data are not randomly distributed; the missing values will tend to be from samples of lower concentrations, and therefore the remaining positive values will then become overestimated 36. Our conservative approach was to follow Ellison et al. 37 and fix their value, although we used the lowest detectable concentration of the assay (13.7 copies/L) rather than fixing the values at zero. Unfortunately, fixing values in this way is also problematic, creating a potential underestimate of concentration, and may interfere with the assumptions of linear regression. A better future strategy may be to avoid estimating decay rates from low copy-number time series, or to impute the missing data 36.
Related to the issue of missing data is that of starting concentrations. Despite normalising each time sample as the proportion of the t = 0 starting concentration, we included in our model the initial value and found it to be a statistically significant predictor associated with faster degradation rates. The summer experiment and the crab assay had lower starting concentrations than the winter experiment and the shanny assay respectively (Fig. 1), but although the average crab-assay decay rate was faster than shanny, the average winter rates were faster than that of summer. Therefore, while they may not have influenced the results overall, a low starting concentration of eDNA resulted in the lower resolution of the summer crab experiment in particular, as qPCR quantification is increasingly stochastic and unreliable at low-template concentrations 38.
In terms of implications for marine ecology, how do eDNA half-lives or decay rate constants relate to detectability of a given organism?
As suggested by Sassoubre et al.20, reporting the duration of time until the detection limit is reached is misleading, as this value will depend upon the starting concentration of eDNA and the sensitivity of the assay; most studies use eDNA starting concentrations far higher than typical natural concentrations in order to generate reliable decay curves with less noise. Our negative biological controls provide an insight into natural concentrations. Sutton Harbour (our inshore treatment) is well populated with common shore crabs, and as expected, we recovered this species at approximate concentrations of 263 copies/L (winter) and 270 copies/L (summer). As the detection rate of the crab assay was 37% at 83 copies/L, and the eDNA half-life inshore was around 24 h, it is estimated that the chance of detection with three PCR replicates would be below the threshold after just two half-life periods (~48 h). However, we did detect eDNA in at least one qPCR replicate from this control at all time points up to 192 h (winter) and 48 h (summer), indicating that eDNA detectability will be difficult to predict at very low concentrations. Quantitative PCR is known to be more sensitive than standard PCR combined with metabarcoding39. Thomsen et al.18 estimated similar values of natural eDNA to ours (535 copies/L for flounder, 120 copies/L for stickleback), and a similar detection limit (63 copies/L). However, it must be noted that we did not consider the loss of DNA in the extraction process, which can be considerable with commercial kits that incorporate steps to remove PCR inhibitors34,40, or any loss of eDNA at the filtration stage, and therefore, real values are likely to be higher and comparisons among studies using different methodologies may be questionable.
Ultimately, how eDNA persists and moves through an environment can have important repercussions for making meaningful ecological inferences, and it is important to document and understand the patterns and processes involved41,42. The combined issues of degradation, transportation and dilution of eDNA are of particular importance in the marine environment, due to the effects of tides and large water volumes 9,18,43. Fortunately, eDNA metabarcoding studies of marine systems have reported a strong local eDNA signal, either closely matching lists of expected fauna 18,44,45 or reporting an expected turnover in diversity over short spatial or temporal scales9,43,46. Most evidence therefore points to eDNA surveys offering a contemporaneous representation of a community, even over the variation encountered on a daily tide9. However, there are cases where non-resident freshwater species have been detected in marine eDNA studies47, and while this source of error can easily be discarded as clearly a riverine input, currents transporting possibly co-occurring marine species eDNA may cause a less obvious source of systematic bias. These biases may become more serious when eDNA is used in applications beyond determining occurrence, for example to monitor the spread of marine invasive species48 or correlating with animal biomass estimates49. By incorporating eDNA degradation rates in different types of water body with oceanographic modelling of tidal currents, it will be possible to build well-informed predictive probability maps of organismal distribution44,48,50. Until these are available, to our knowledge, we show for the first time that it is reasonable to assume large variation in eDNA persistence according to local factors such as salinity gradients over relatively short local scales corresponding to marine environmental stability.
Study species were the shanny (Teleostei: Blenniidae: Lipophrys pholis) and the common shore crab (Decapoda: Portunidae: Carcinus maenas). These species were chosen because they are abundant hardy organisms amenable to transport and experimental manipulation. Reference specimens of shanny (eight individuals) and shore crab (four individuals) were obtained from the Gann estuary, Pembrokeshire, Wales (51.715, −5.173). Using standard molecular methods, we obtained DNA barcodes (COI; 5′ mitochondrial cytochrome c oxidase I gene) for both species using the FishF1/R1 primer set51. Additional sequence data for crabs (149 individuals) were obtained from GenBank; no GenBank COI sequences were available for shanny. Primers and hydrolysis probes were designed using Primer3 v1.1.452,53 under default settings adjusted to aim for an amplicon length between 50 and 170 bp. The resulting 12 candidate primer pairs were tested in silico for general specificity against a dataset of sequences from species present in the United Kingdom. To generate a list of fishes and Malacostraca recorded from the United Kingdom, we searched the Global Biodiversity Information Facility (https://www.gbif.org) using the rgbif v0.9.9 package for R54. COI sequences for these species were then retrieved from GenBank and annotated using rentrez v1.2.155 and traits v0.3.0.931056. Each candidate primer pair was tested in an in silico PCR using MFEprimer v2.057 using liberal settings (k = 5). The final primers were then chosen based on a combination of amplicon length, specificity and melting temperatures, and are reported in Supplementary Table 1. The reporter dye for the shanny assay was FAM, and for the crab assay HEX; both were quenched using BHQ.
The experiment was repeated twice, first in winter (water collected on 17 February, 2017) and once in late summer (water collected on 26 September, 2017). All treatments were set up in a dedicated temperature-controlled aquarium room held at temperatures consistent with natural seawater temperatures at that time of the year (10 °C, winter experiment; 15 °C, summer experiment). Animals were collected 2 days before the start of each experiment (also from the Gann estuary, Pembrokeshire) and placed in a separate and aerated holding tanks for each species (shanny, 50 L of synthetic seawater; crabs, 25 L). Approximately 300 g of animal mass per species were collected (winter, 24 shannys at 343 g of total weight, 18 crabs at 288 g; summer, 25 shannys at 316 g, 9 crabs at 304 g). All animals were euthanised after the experiment was completed, and were formalin fixed and 70% alcohol preserved as voucher specimens for a reference dataset. All experiments were carried out in accordance with the University of Bristol ethical approval (UIN reference UB/16/012).
A total of 24 aquariums at the University of Bristol Animal Services Unit were each filled with 9 L of experimental water. The tanks were initially mixed but not aerated and were maintained under 12 h of light/dark LED room lighting. Five experimental water treatments were carried out as follows: 100% offshore sea surface water—from herein referred to as ‘offshore’—collected from Western Channel Observatory station E1 ~40 km from Plymouth, Devon, UK (50.033, −4.367; Supplementary Fig. 2); inshore urban water—from herein referred to as ‘inshore’—collected from Sutton Harbour, Plymouth Sound, a site located between the estuaries of the rivers Plym and Tamar (50.370, −4.133; Supplementary Fig. 2); a two-thirds/one-third mixture of offshore to inshore water; a one-third/two-thirds mixture of offshore to inshore water; and synthetic seawater made using a proprietary aquarium salt mix. Each of the five treatments had four biological replicates (=20 tanks), plus four no-treatment controls (2× synthetic seawater, 1× offshore and 1× inshore.
After turning off aeration and allowing detritus to settle for an hour, 500 mL of eDNA-rich surface water from both the shanny and crab stock tanks was then added to each experimental tank at the start of the experiment. At each subsequent time point, eDNA was filtered from 600 mL of experimental tank water with a peristaltic pump and Sterivex 0.22-μm PES filters (Millipore part no. SVGP01050)58. Measurements were taken at six intervals from the same tanks (0, 12, 24, 48, 96 and 192 h), resulting in 144 filtered water samples (24×6). After being cleared of water, filters were frozen immediately at −20 °C. DNA was subsequently extracted from the Sterivex filters using the PowerWater DNA isolation kit (MoBio/Qiagen part no. 14900-100-NF) following manufacturers’ instructions, but with 50 μL of final elution volume. Extractions were carried out in a dedicated pre-PCR extraction room regularly decontaminated with 10% bleach and UV sterilisation. Environmental covariates were also measured from each tank with a Hach HQ40D multimeter, and included salinity (conductivity in mS/cm), pH, and temperature at source. As a proxy for biological activity, we also recorded total background double-stranded DNA (dsDNA) concentration from 2 L of source water with a Qubit 3 fluorometer (ThermoFisher) assay (filtered and extracted in the same way as the experimental treatments).
Quantitative PCR reactions were conducted as per the manufacturer’s instructions, in multiplex, on a PCRmax Eco48 machine in 48-well plates of 5 μL per reaction, with ROX normalisation. Each reaction volume comprised 2.5 μL of mastermix (qPCRBIO Lo-Rox Probe mix; part no. PB20.21-05); 0.5 μL of shanny–crab primer-probe mix (optimised reaction concentration for shanny assay: 600 nM each primer, 200 nM probe; crab assay: 600 nM each primer, 300 nM probe); 1 μL of water and 1 μL of eDNA template. The cycling parameters comprised 3 min at 95 °C polymerase activation followed by 42 cycles of denaturation at 95 °C for 5 s and combined extension/annealing at 60 °C for 30 s.
Each plate of 48 reactions comprised: eight extracted water samples of the experimental tanks, with three technical replicates per sample (8 × 3 = 24 reactions); a six-fold standard-curve serial dilution of 1–1 million copies/μL, in triplicate (=21 reactions); and three no-template controls (=three reactions). To allow low-copy-number templates, an increased opportunity to amplify, PCRs were repeated a further three times for each sample when there was no amplification in any of the three initial technical replicates (excluding negative controls). The standard curve stock solutions were generated by PCR-amplifying and purifying tissue extractions of genomic DNA in a standard PCR using the primers in Supplementary Table 1, and were subsequently diluted and quantified using a Qubit assay, with the number of copies estimated at a standard dsDNA molar mass of 650 g 28. We tested for PCR inhibitors by performing triplicate qPCRs on three serial dilutions of the 0 h replicates from three treatments (synthetic, inshore, offshore) over both seasons (total 24 samples). If inhibitors were co-extracted, the cycle threshold (Ct) values at each tenfold dilution point would deviate from the expected increase of 3.3 PCR cycles, and therefore the expected efficiency values of 90–110% 28.
Cycle threshold values and target DNA concentrations were calculated on the Eco48 machine software using the default settings, and converted from copies per reaction (=copies/μL given a 1-μL template volume) to copies/L of initial sample water (given a 600 mL filtration volume and a 50 μL elution volume). All amplifications were checked manually in the log plot view and any amplifications that crossed the baseline threshold, but that did not represent a clean, obviously exponentially increasing reaction, were excluded. The final eDNA concentrations for each sample were averaged over the technical replicates, with non-amplifications included as an arbitrarily low but non-zero value of 13.7 copies per litre of sample water (Ct = 38.5; the lowest concentration that the assay reliably detected).
Statistical analyses were conducted in R v3.5.159. Decay of eDNA was modelled using a linear mixed-effects model as implemented in the lme function of nlme v3.1-13760. The response variable was natural loge transformed eDNA concentration normalised as a proportion of starting concentration, i.e. the value at time t = x divided by the value at time t = 0. We specified time, treatment, season, assay, and the natural log of eDNA starting concentration as predictor variables (our fixed effects), while the individual tank used in each biological replicate was treated as a random effect. To minimise heteroscedasticity—i.e. the increasing variance of regression residuals over time—we excluded the normalised zero-hour (t = 0) data, which had no variance. The synthetic water control was also excluded from the main model—this was a control for reference rather than to investigate its biological effect—and decay rates for this subset were calculated separately (following the same procedure as outlined below).
We determined the optimal model to fit our data according to the procedure of Zuur et al.61. We started with a full model containing all fixed effects and their interactions, and determined the optimal variance weighting for different treatment–season–assay combinations by AIC comparison (given by the form weights = varIdent (form = 1|treatment*season)). We then determined the optimal random structure for the full model with this variance weighting by AIC comparison (given by the form random = 1 + time|tank). Finally, we determined the optimal fixed effects structure using the ‘drop1’ approach and specifying method = ’ML’ until all terms in the model were statistically significant. We switched to method = ’REML’ and performed model validation to ensure that the model residuals were approximately normal and homogeneously distributed (see Supplementary Fig. 3). The fixed effects structure and output for the full model and the optimal model are also presented in Supplementary Note 1.
The first-order decay-rate constant k for each treatment–season–assay combination was calculated from the estimated marginal mean of regression slopes using the emtrends function of emmeans v1.2.362. To test the importance of predictor variables on the degradation rate, pairwise post hoc Tukey tests were carried out on the marginal mean regression slopes, again using emmeans. To explore the environmental covariates we constructed a simple lme model with time, assay, pH, conductivity, natural log transformed starting concentration, and background DNA concentration as fixed effects, and tank as a random effect. For this model, we excluded the treatment and season predictors—which were deliberately chosen for their heterogeneity—as we assumed these to be correlated with the environmental covariates. We additionally included them in a combined model to estimate the degree of correlation between the predictors and covariates.LINK TO PART I LINK TO THE ORIGINAL ARTICLE
The code generated during and/or analysed during the current study is available in the Figshare repository 65
The datasets generated during and/or analysed during the current study are available in the Figshare repository 65. New sequence data generated here were deposited in the GenBank nucleotide archive under the accessions MH931374:MH931388.
1. Bohmann, K. et al. Environmental DNA for wildlife biology and biodiversity monitoring. Trends Ecol. Evol. 29, 358–367 (2014).
2. Rees, H. C., Maddison, B. C., Middleditch, D. J., Patmore, J. R. M. & Gough, K. C. The detection of aquatic animal species using environmental DNA – a review of eDNA as a survey tool in ecology. J. Appl. Ecol. 51, 1450–1459 (2014).
3. Deiner, K. et al. Environmental DNA metabarcoding: Transforming how we survey animal and plant communities. Mol. Ecol. 26, 5872–5895 (2017).
4. Barnes, M. A. & Turner, C. R. The ecology of environmental DNA and implications for conservation genetics. Conserv. Genet. 17, 1–17 (2016).
5. Seymour, M. et al. Acidity promotes degradation of multi-species eDNA in lotic mesocosms. Commun. Biol. 1, 4 (2018).
6. Bakker, J. et al. Environmental DNA reveals tropical shark diversity in contrasting levels of anthropogenic impact. Sci. Rep. 7, 16886 (2017).
7. Kelly, R. P. et al. Genetic signatures of ecological diversity along an urbanization gradient. PeerJ 4, e2444 (2016).
8. Deiner, K. & Altermatt, F. Transport distance of invertebrate environmental DNA in a natural river. PLoS One 9, e88786 (2014).
9. Kelly, R., Gallego, R. & Jacobs-Palmer, E. The effect of tides on nearshore environmental DNA. PeerJ 6, e4521 (2018).
10. Dejean, T. et al. Persistence of environmental DNA in freshwater ecosystems. PLoS One 6, e23398 (2011).
11. Barnes, M. A. et al. Environmental conditions influence eDNA persistence in aquatic systems. Environ. Sci. Technol. 48, 1819–1827 (2014).
12. Strickler, K. M., Fremier, A. K. & Goldberg, C. S. Quantifying effects of UV-B, temperature, and pH on eDNA degradation in aquatic microcosms. Biol. Conserv. 183, 85–92 (2015).
13. Lance, R. et al. Experimental observations on the decay of environmental DNA from bighead and silver carps. Manag. Biol. Invasions 8, 343–359 (2017).
14. Eichmiller, J. J., Best, S. E. & Sorensen, P. W. Effects of temperature and trophic state on degradation of environmental DNA in lake water. Environ. Sci. Technol. 50, 1859–1867 (2016).
15. Cowart, D. A., Murphy, K. R. & Cheng, C.-H. C. Metagenomic sequencing of environmental DNA reveals marine faunal assemblages from the West Antarctic Peninsula. Mar. Genom. 37, 148–160 (2018).
16. Minamoto, T. et al. Environmental DNA reflects spatial and temporal jellyfish distribution. PLoS One 12, e0173073 (2017).
17. Sigsgaard, E. E. et al. Population characteristics of a large whale shark aggregation inferred from seawater environmental DNA. Nat. Ecol. & Evol. 1, 0004 (2016).
18. Thomsen, P. F. et al. Detection of a diverse marine fish fauna using environmental DNA from seawater samples. PLoS One 7, e41732 (2012).
19. Weltz, K. et al. Application of environmental DNA to detect an endangered marine skate species in the wild. PLoS One 12, e0178124 (2017).
20. Sassoubre, L. M., Yamahara, K. M., Gardner, L. D., Block, B. A. & Boehm, A. B. Quantification of environmental DNA (eDNA) shedding and decay rates for three marine fish. Environ. Sci. Technol. 50, 10456–10464 (2016).
21. Andruszkiewicz, E. A., Sassoubre, L. M. & Boehm, A. B. Persistence of marine fish environmental DNA and the influence of sunlight. PLoS One 12, e0185043 (2017).
22. Jo, T. et al. Rapid degradation of longer DNA fragments enables the improved estimation of distribution and biomass using environmental DNA. Mol. Ecol. Resour. 38, 3218–3221 (2017).
23. Okabe, S. & Shimazu, Y. Persistence of host-specific Bacteroides-Prevotella 16S rRNA genetic markers in environmental waters: Effects of temperature and salinity. Appl. Microbiol. Biotechnol. 76, 935–944 (2007).
24. Bae, S. & Wuertz, S. Rapid decay of host-specific fecal Bacteroidales cells in seawater as measured by quantitative PCR with propidium monoazide. Water Res. 43, 4850–4859 (2009).
25. Schulz, C. J. & Childers, G. W. Fecal Bacteroidales diversity and decay in response to variations in temperature and salinity. Appl. Environ. Microbiol. 77, 2563–2572 (2011).
26. Salter, I. Seasonal variability in the persistence of dissolved environmental DNA (eDNA) in a marine system: the role of microbial nutrient limitation. PLoS One 13, e0192409 (2018).
27. Smyth, T. J. et al. A broad spatio-temporal view of the Western English Channel observatory. J. Plankton Res. 32, 585–601 (2010).
28. Agersnap, S. et al. Monitoring of noble, signal and narrow-clawed crayfish using environmental DNA from freshwater samples. PLoS One 12, e0179261 (2017).
29. Dell’Anno, A. & Corinaldesi, C. Degradation and turnover of extracellular DNA in marine sediments: ecological and methodological considerations. Appl. Environ. Microbiol. 70, 4384–4386 (2004).
30. Torti, A., Lever, M. A. & Jørgensen, B. B. Origin, dynamics, and implications of extracellular DNA pools in marine sediments. Mar. Genom. 24, 185–196 (2015).
31. Tsuji, S., Ushio, M., Sakurai, S., Minamoto, T. & Yamanaka, H. Water temperature-dependent degradation of environmental DNA and its relation to bacterial abundance. PLoS One 12, e0176608 (2017).
32. Gilbert, J. A. et al. Defining seasonal marine microbial community dynamics. ISME J. 6, 298–308 (2012).
33. Cox, A. M. & Goodwin, K. D. Sample preparation methods for quantitative detection of DNA by molecular assays and marine biosensors. Mar. Pollut. Bull. 73, 47–56 (2013).
34. Eichmiller, J. J., Miller, L. M. & Sorensen, P. W. Optimizing techniques to capture and extract environmental DNA for detection and quantification of fish. Mol. Ecol. Resour. 16, 56–68 (2016).
35. Dunn, N., Priestley, V., Herraiz, A., Arnold, R. & Savolainen, V. Behavior and season affect crayfish detection and density inference using environmental DNA. Ecol. Evol. 7, 7777–7785 (2017).
36. McCall, M. N., McMurray, H. R., Land, H. & Almudevar, A. On non-detects in qPCR data. Bioinformatics 30, 2310–2316 (2014).
37. Ellison, S. L. R., English, C. A., Burns, M. J. & Keer, J. T. Routes to improving the reliability of low level DNA analysis using real-time PCR. BMC Biotechnol. 6, 1–11 (2006).
38. Hunter, M. E. et al. Detection limits of quantitative and digital PCR assays and their influence in presenceâ€“absence surveys of environmental DNA. Mol. Ecol. Resour. 17, 221–229 (2017).
39. Harper, L. R. et al. Needle in a haystack? A comparison of eDNA metabarcoding and targeted qPCR for detection of great crested newt (Triturus cristatus). Ecol. Evol. 8, 6330–6341 (2018).
40. Deiner, K., Walser, J. C., Mächler, E. & Altermatt, F. Choice of capture and extraction methods affect detection of freshwater biodiversity from environmental DNA. Biol. Conserv. 183, 53–63 (2015).
41. Thomsen, P. F. & Willerslev, E. Environmental DNA – an emerging tool in conservation for monitoring past and present biodiversity. Biol. Conserv. 183, 4–18 (2015).
42. Goldberg, C. S. et al. Critical considerations for the application of environmental DNA methods to detect aquatic species. Methods Ecol. Evol. 7, 1299–1307 (2016).
43. Port, J. A. et al. Assessing vertebrate biodiversity in a kelp forest ecosystem using environmental DNA. Mol. Ecol. 25, 527–541 (2016).
44. Thomsen, P. F. et al. Environmental DNA from seawater samples correlate with trawl catches of subarctic, deepwater fishes. PLoS One 11, e0165252 (2016).
45. Stoeckle, M. Y., Soboleva, L. & Charlop-Powers, Z. Aquatic environmental DNA detects seasonal fish abundance and habitat preference in an urban estuary. PLoS One 12, e0175186 (2017).
46. O’Donnell, J. L. et al. Spatial distribution of environmental DNA in a nearshore marine habitat. PeerJ 5, e3044 (2017).
47. Yamamoto, S. et al. Environmental DNA metabarcoding reveals local fish communities in a species-rich coastal sea. Sci. Rep. 7, 40368 (2017).
48. Richardson, M. F., Sherman, C. D., Lee, R. S., Bott, N. J. & Hirst, A. J. Multiple dispersal vectors drive range expansion in an invasive marine species. Mol. Ecol. 25, 5001–5014 (2016).
49. Yamamoto, S. et al. Environmental DNA as a’snapshot’ of fish distribution: A case study of Japanese jack mackerel in Maizuru Bay, Sea of Japan. PLoS One 11, e0149786 (2016).
50. Sansom, B. J. & Sassoubre, L. M. Environmental DNA (eDNA) shedding and decay rates to model freshwater mussel eDNA transport in a river. Environ. Sci. Technol. 51, 14244–14253 (2017).
51. Ward, R. D., Zemlak, T. S., Innes, B. H., Last, P. R. & Hebert, P. D. N. DNA barcoding Australia’s fish species. Philos. Trans. R. Soc. Lond. B Biol. Sci. 360, 1847–1857 (2005).
52. Rozen, S. & Skaletsky, H. Primer3 on the WWW for General Users and for Biologist Programmers. In Bioinformatics Methods and Protocols, 365–386 (Humana Press, New Jersey, 1999).
53. Untergasser, A. et al. Primer3-new capabilities and interfaces. Nucleic Acids Res. 40, e115 (2012).
54. Chamberlain, S. rgbif: Interface to the Global’Biodiversity’ Information Facility API (2017). URL https://CRAN.R-project.org/package=rgbif. R package version 0.9.9.
55. Winter, D. J. rentrez: an R package for the NCBI eUtils API. R. J. 9, 520–526 (2017).
56. Chamberlain, S., Foster, Z., Bartomeus, I., LeBauer, D. & Harris, D. traits: species trait data from around the Web (2017). URL https://cran.r-project.org/package=traits. R package version 0.3.0.9310.
57. Qu, W. et al. MFEprimer-2.0: A fast thermodynamics-based program for checking PCR primer specificity. Nucleic Acids Res. 40, 205–208 (2012).
58. Spens, J. et al. Comparison of capture and storage methods for aqueous macrobial eDNA using an optimized extraction protocol: advantage of enclosed filter. Methods Ecol. Evol. 8, 635–645 (2017).
59. R Core Team. R: A Language and Environment for Statistical Computing (2017). URL https://www.r-project.org/.
60. Pinheiro, J., Bates, D., DebRoy, S., Sarkar, D. & R Core Team. nlme: linear and nonlinear mixed effects models (2017).
61. Zuur, A. F., Ieno, E. N., Walker, N. J., Saveliev, A. A. & Smith, G. M. Mixed Effects Models and Extensions in Ecology with R. (Springer, New York, 2009).
62. Lenth, R. emmeans: Estimated Marginal Means, aka Least-Squares Means (2018). URL https://CRAN.R-project.org/package=emmeans. R package version 1.2.3.
63. Maruyama, A., Nakamura, K., Yamanaka, H., Kondoh, M. & Minamoto, T. The release rate of environmental DNA from juvenile and adult fish. PLoS One 9, e114639 (2014).
64. Pilliod, D. S., Goldberg, C. S., Arkle, R. S. & Waits, L. P. Factors influencing detection of eDNA from a stream-dwelling amphibian. Mol. Ecol. Resour. 14, 109–116 (2014).
65. Collins, R. A. edna-persistence. doi:10.6084/m9.figshare.7111376.v1 (2018)
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.