Abstract
Environmental DNA analysis has emerged as a key component of biodiversity and environmental monitoring. However, the state and fate of eDNA in natural environments is still poorly understood for many ecological systems. Here we assess the state and fate of eDNA derived from the water flea, Daphnia magna, using a full factorial mesocosm experiment. We measured the quantity and degradation of eDNA over a two month period across a range of filters differing in pore size (0, 0.2, 1 and 10âµm), which spans the range of eDNA source material including subcellular, cellular and tissue. We also used two primer sets targeting mitochondrial (COI) and nuclear (18S) genomic regions. Our findings demonstrated that eDNA was most prevalent in the effluent water, but also reliably detected on the 0.2 μm filter, suggesting subcellular material is the predominate state of eDNA. Temporal eDNA quantity dynamics followed an exponential decay function over the course of 6-17 days, demonstrating a predictable decline in eDNA concentration. Nuclear eDNA was more abundant than mitochondrial eDNA, which may be a result of greater primer affinity, or indicate greater availability of nuclear eDNA gene targets in the environment. In contrast to two previous size-sorting experiments, which utilizing fish eDNA, our findings suggest that the state of invertebrate eDNA is much smaller than previously suspected. Overall, our data suggest that the detection of eDNA greatly depends on our knowledge of the state and fate of eDNA, which differ among species, and likely across environmental conditions.
Similar content being viewed by others
Introduction
There is no other source of biodiversity information like environmental DNA (eDNA), which can detect whole communities from any environment. Consequently, eDNA analysis has been named as a âgame changerâ for biodiversity sampling1,2. Environmental DNA refers to macrobial DNA that is extracted from an environmental sample (e.g. water, soil or air), without targeting a particular organism3. Combined with genetic amplification and sequencing, eDNA analysis enables a wide range of research questions, across a wide range of disciplines; including molecular ecology, palaeontology, conservation/invasion biology, ecology and environmental sciences4,5. The non-invasive means of detecting species via eDNA provides a reliable biomonitoring approach that avoids disturbing the ecology of the target species3,6,7. In case of aqueous environments, eDNA of macro-organisms offers a simple and sensitive standardized means of sample collection and identification, often with better detection of diverse fauna, compared to traditional direct sampling methods8,9. Aquatic eDNA has been utilized for monitoring ecological communities10,11, assessing invasion dynamics12, conservation monitoring13 and assessing localized extinction14. Furthermore, eDNA can be collected from any type of aquatic environment, and has been used to assess species living in lakes, rivers, groundwater and marine environments11,15,16,17. Though eDNA stands as an important tool to complement traditional methods, there are still some shortcomings in our understanding of the very nature of eDNA that can impede practical eDNA study designs18,19.
The ability to relate eDNA information to its original source, either a target species or whole community diversity, is limited by our understanding of the temporal, physical, and chemical factors that influence eDNA detectability in natural environments. The detection of diverse aquatic macrofauna using eDNA requires the knowledge of origin (e.g. physiological source), state (e.g. intra- or subcellular particles), fate (e.g. suspension time) and transport dynamics of eDNA20. The state of eDNA is an important domain that directly affects the ability of current sampling methods to detect eDNA. Our current understanding is that eDNA state predominates as necromass (dead biomass), of various concentrations, in water, soil and sediments21,22. In aqueous environments, eDNA necromass can remain in many forms such as intracellular DNA, subcellular DNA, soluble DNA, non-soluble DNA, DNA within dead but structurally intact cells, organically or inorganically complex DNA and DNA absorbed by sediment minerals23. Structurally, eDNA is found in the environment in both intracellular (e.g. tissue or cells) and subcellular (e.g. mitochondria, ribosomes or free floating nucleotide strands) forms with intracellular eDNA transforming to subcellular eDNA over time as cells degrade24.
What remains unclear is what eDNA state (intracellular or subcellular) dominates an environment when eDNA sampling takes place. Present aqueous eDNA sampling methods primarily utilize a range of water filtering approaches, with limited consideration for the particle size (i.e. state) of the targeted eDNA of interest. The range for nominal filtering for eDNA studies varies widely across taxonomic groups, ranging from <0.2 μm to â¥180 μm25,26,27,28,29,30. However, eDNA particle sizes likely differ among species and bodily sources of origin, such as eDNA originating from whole body decomposition versus eDNA from defecated material. It is likely that different ranges of eDNA particle sizes exist per species, and across seasons within species20. Therefore, samplers looking to maximize the capture of the most abundant particle sizes across species or environments generally opt for small filter sizes23. However, smaller filters also capture additional non-target material, which can dilute or inhibit the intended eDNA target31. Knowledge of the size and distribution of the optimal eDNA particle size for sampling and genetic analyses are key variables to consider when designing an eDNA based study. Additionally, studies should consider the genetic region used to identify eDNA particles, as the genetic state of eDNA in the natural environment is composed of nuclear and mitochondrial DNA, which differ in their structure, abundance and availability32,33,34.
Whether eDNA exists primarily as mitochondrial (mtDNA) or nuclear (cDNA) DNA is unclear. Presently, many eDNA studies utilize genetic tools to target mitochondrial genes due, in part, to existing sequence databases being predominately based on the COI mitochondrial gene3. Additionally, mtDNA is expected to occur in higher density compared to cDNA since each cell has 2â10 mitochondria per cell, compared to 2 copies of cDNA inhabiting the nucleaus35. However, tandemly repeated ribosomal cDNA-based genes, including the small subunit (18S) and the large subunit (28S), occur at similar or higher copies per cell compared to mtDNA34,36. From forensic-based research there are known benefits to utilizing short nuclear ribosomal markers, due to higher PCR amplification success of degraded ribosomal cDNA33,37. Recent eDNA studies have also suggested that ribosomal gene markers are optimal for differentiating closely related species38, and for increased detectability of fish species from eDNA samples, compared to mtDNA markers34. However, mtDNA has a more stable cellular structure compared to ribosomal cDNA, rendering mtDNA likely better for detecting and assessing the fate of eDNA over longer time spans33.
The fate (i.e. persistence) of eDNA in aqueous environments varies between hours and a month, depending on the system being sampled and the underlying abiotic and transport dynamics associated with the sampled ecosystem39. Fate dynamics of eDNA are perhaps the most studied, yet one of the least understood aspects of aqueous eDNA research due to the complex dynamics involved in lentic40,41, lotic42,43 and marine environments44. In general, once released into an environment, eDNA starts to decay or leave the system, approximately following a first order exponential decay function, due to abiotic conditions (temperature, oxygen, pH)43, transport factors (flow rate, turbidity)15,45, substrate absorption (substrate and biofilm)46 and possibly biotic interactions (microbes and extracellular enzymes)20. Additionally, variations in physical properties of DNA, such as length, conformation, and structure, alter the susceptibility of eDNA decay to abiotic and biotic factors20. Less understood is whether various states of eDNA including particle size or nuclear versus mitochondrial sources degrade similarly under common environmental conditions.
Only by knowing the state of macrobial eDNA in the environment will we be able to select the most efficient eDNA sampling method. In this study we conducted a mesocosm experiment (Fig. 1) to assess the size distribution (state) and temporal dynamics (fate) of aquatic macrobial eDNA particles originating from mitochondrial and nuclear genetic fragments. Using eDNA originating from Daphnia magna we addressed three main objectives: (i) quantify the eDNA particle size distribution of eDNA using filter size sorting; (ii) compare eDNA quantifications between mitochondrial and nuclear derived material via quantitative PCR (qPCR); (iii) calculate the decay rate across detectable eDNA particle size classes for mitochondrial and nuclear derived eDNA. Collectively, this study enhances the understanding of the state of macrobial eDNA in aquatic environments to assists researchers in adapting eDNA sampling and analytical methods.
Results
Mean starting (day 0) eDNA concentrations, quantified via qPCR, were 75.61â±â31.67 copy numbers (copies) for the 18S water effluent (size 0), 98.77â±â31.30 copies for 18S 0.2 μm filter, 51.58â±â22.63 copies for COI water effluent, and 42.19â±â0.30.25 copies for COI 0.2 μm filter. Total DNA concentrations, quantified via Qbit fluorometer, were lowest in the water effluent 0.544â±â0.504âng/ul compared to 3.02â±â3.28 for the 0.2 μm filters, 4.15â±â6.92 in the 1 μm filers and 2.66â±â4.32 in the 10 μm filters (Fig. 2). The negative controls samples showed zero amplification across all filter sizes, replicates and time points. Successful amplification (across all replicates) from water effluent and 0.2âµm samples were recorded for both mitochondrial (COI) and nuclear (18s) markers for days 0, 1 and 3. Environmental DNA decay occurred primarily from day 0 to day 6 (one week) for COI and 18S markers for all water effluent and 0.2âµm filter samples, with inconsistent amplification occurring for day 17 and no amplification occurring for day 31 (Fig. 3). Amplification generally failed for samples collected on the 1âµm and 10âµm filters, for both COI and 18S markers, apart from isolated amplification occurring up to day 6 for the 1âµm and 10âµm filters (Fig. 2).
From our repeated measures ANOVA, which included an AR1 correlation structure to account for the temporal autocorrelation for each tank across time47, we found significant effects of filter size (p-valueâ<â0.001) and marker (p-valueâ=â0.022) on eDNA copy number quantification, and a non-significant filter x marker interaction (p-value 0.124) (Table 1). The model employed had the lowest AIC score of all alternative models, indicating a better model fit over alternative models, including a linear regression model that excluded the AR1 correlation structure (Table 2). Decay rates differed primarily among filter types within the 18S amplified samples, with a decay rate of â0.360 per day for the effluent water samples versus â0.544 per day for the 18S 0.2âµm filter samples and â0.628 for COI effluent water samples versus â0.501 per day for the COI 0.2âµm filter samples (Table 3, Fig. 3).
Discussion
Our findings show that eDNA fate followed an exponential decay function. We found that the eDNA particle size for D. magna is predominantly between 0 and 0.2âµm, suggesting eDNA state is predominately subcellular in nature, meaning cell organelles or free-floating DNA strands. Furthermore, nuclear derived eDNA quantifications were greater than mitochondrial derived eDNA quantifications, suggesting nuclear eDNA was more readily available for capture in our environmental samples.
Degradation of eDNA in our study occurred primarily over a week, with highly variable eDNA detection persisting for up to a month. Under natural environmental conditions eDNA is detectable up to a month in lentic systems40, around 24âhours in the lotic system19, and ranging from hours43 to months under experimental settings48. Additionally, the persistence of eDNA across different species has been experimentally shown to vary from hours to months28,31,49. Hence, the rate of eDNA decay observed was comparable to previous studies, particularly for eutrophic environments44 where biotic interactions are suspected to promote eDNA degradation rates20. We did observe algal growth within the mesocosms overtime, likely facilitated by the constant 20â°C climate, which suggested that biotic activity could have played a role in degradation20. Additionally, the pH of the water source, though not measured here, is categorized as moderately soft (low mineral), which infers neutral to slightly acidic pH (~6.5â7). The low pH, which has been found to promote eDNA degradation43, may also explain the eDNA decay observed here, which is similar to rates of eDNA decay observed at neutral and slightly acid pH conditions in other lentic experimental studies50,51. However, further testing is required to verify the exact mechanisms of eDNA degradation in our particular system. We also observed unequal degradation rates for different size eDNA particles, specifically for the 0.2âµm and 0âµm (water effluent) particle size fractions, whereby decay was greater over time for the 0.2âµm size fraction compared to the 0âµm size. As suggested by Turner et al. 2014, the breakdown of sub-cellular structures (0.2âµm) into free nucleotide strands (0âµm) is likely prolonging the decay of the 0âµm size class, thereby allowing the extended observation of the smaller size class through the extension of natural degradation processes28.
Total copy numbers from the 18S derived nuclear eDNA samples were ten times greater than the COI derived mtDNA eDNA, indicating nuclear eDNA was more readably amplified from our eDNA samples. The temporal decay dynamics and size fractioning (i.e. state) were similar for the 18S and COI derived eDNA, indicating that nuclear eDNA was generally more abundant, but that eDNA detectability was not affected by nuclear or mitochondrial origin. In general, technical performance differences among PCR-based genetic markers/primers is greatly associated with GC content and melting temperature, however these aspects should be tested in-silico during the primer design stage and do not reflect differences observed in this study. Specifically, the primers utilized in this study have similar melting points. The GC content is lower in the COI marker despite the 18S marker having greater performance, when higher GC content is expected to reduce PCR performance due to the higher GC bond strength31. The underlying differences between the COI and 18S markers observed in this study are possibly due to nuclear eDNA being more sensitive to eDNA detection compared to mitochondrial eDNA34,36. Forensic based molecular biology suggest that nuclear DNA may be more readily available for PCR amplification, while also degrading more rapidly compared to mitochondrial DNA33. The increased degradation rate and increased availability of nuclear DNA, are possibly due to the chromatin structure of nuclear DNA being more susceptible to exonucleases and protease compared to the more protective and less accessible circular structure of mitochondrial DNA33,52.
Animal mitochondria range from 0.2 to 1.2âµm in diameter, with length between 1â8 µm25. Our findings suggest that given the lack of eDNA detection from the 1 and 10âµm filters, that the primary state of accessible eDNA originating from D. magna was subcellular. At the start of the experiment, the most abundant eDNA particle size was on the 0.2âµm filters, for both nuclear and mitochondrial derived eDNA. In contrast, the water effluent samples (size 0âµm) generally retained the highest amount of eDNA over time, which can be associated with free-floating nucleotide strands. The high initial quantity of 0.2âµm eDNA combined with the longer persistence of smaller fragments exhibits a distinct temporal dynamic of eDNA particle state and fate. Aqueous macrobial eDNA typically originates from excreted urine, faeces, cellular by-products, and epidermal tissues. Once released from its bodily source eDNA initially exist as large/bulk particles (>1000âµm)28. However, large particles readily break apart due to physical (e.g. water currents, abrasion) and environmental factors, such as temperature, oxygen, pH, salinity, light, and decomposition via microbial and extracellular enzymatic activity24. Potentially, once released and exposed to the environment, decaying cellular matter rapidly degrades due to hydrolysis and apoptosis destroying cell walls53. Subsequently, apoptotic released mitochondria, ribosomes and other subcellular material accumulate, due to endonucleases resistance of intercellular structures, thereby prolong the persistance of short eDNA fragments54,55. In aqueous environments the double membrane of the mitochondria protects the mitochondrial cell and inner mitochondrial nucleoid from rapid lysis, whereas nuclear DNA decays more rapidly due to its chromatin structure exposing it to exonuclease activity56.
Daphnia magna eDNA was primarily captured on the 0.2âµm filter and from the water effluent, whereas a paucity of D. magna eDNA was capture on the larger 1âµm and 10âµm filters. The only two previous size sorting experiments found eDNA predominantly on filters of larger pore size (1â1.2âµm)26,28. The overarching conclusion from previous studies, was that fish sourced eDNA was of a cellular nature and likely comprised of loosely aggregated cells and tissue, which are captured on filters with pore sizes greater than 1 µm26,28. Why then the differences between the D. magna and fish eDNA capture studies? Both of the previous experiments utilized fish sourced eDNA, specifically Common carp (Cyprinus carpio) and Brook trout (Salvelinus fontinalis) at high eDNA concentrations (i.e. fish biomass to water volume), which may have included sampling of fish tissue or eDNA material released as a result of animal stress. Conversely, our study did not include D. magna in the experimental tanks and relied on a diluted source of eDNA, thereby avoiding the likelihood of filter clogging from non-targeted fragments. Two working hypotheses can then be postulated to address the subtle, yet important differences between the present and previous studies. 1) The nature of eDNA particles from fish and macroinvertebrate taxa differ in their size distributions. Whole mitochondrial cells range between 0.2âµm and 1.2âµm in diameter with length between 1â8âµm, and would be captured primarily by the 0.2âµm. Whereas, eukaryotic cellular structures (>0.8âµm) are larger than the 0.2âµm filter sizes and would accumulate on 1 or 10âµm filters25. Fisheries derived eDNA may feature larger agglomerations of cellular material, thereby explaining the preferential captured efficiently on 1âµm filters26,28. Conversely, invertebrate eDNA, originating from smaller-bodied individuals compared to fish, may consist primarily of fragmented epithelial cells or mitochondrial DNA. 2) Alternatively, the nature of the eDNA in the two previous and our present studies are of a similar size distribution, however the volume of eDNA and organic particulates was greater in previous studies that utilized eDNA from larger biomass (fish) and water sourced from natural ponds and streams. Such that the suspended particulate matter and increased eDNA concentration in Turner et al. (2014) and Wilcox et al. (2015) may have clogged the larger filter, which allowed subcellular eDNA to be captured on the 1âµm filter. Further experiments, exploring different levels of sediment loading, or particulate organic matter could distinguish between the two hypotheses, providing definitive insights into the genuine state of macrobial eDNA.
There is a general need to understand the reliability of eDNA-based detection methods3,31,57. Environmental DNA fate and state, as presented here, encompasses two core components required to enhance future eDNA-based studies58. Knowledge of the state of eDNA in natural environments, including particle size distribution and the predominantly available gene, are important for efficient capture of eDNA material and for modelling eDNA dynamics28. Furthermore, if eDNA sampling design fails to account for the correct eDNA state there is increased potential for false negatives, which can compromise research efforts. Our study suggests that sampling frequency, filter size selection, and decay dynamics play a significant role in the detection of eDNA particles, though such variables are likely to be influenced by the background environment or the study species. Specifically, we show that efficient capture of invertebrate eDNA from water samples requires smaller pore sizes than previously suggested. The larger implications of our findings to wider eDNA research being that eDNA particle state differs among species or environmental conditions, which need to be considered for effective eDNA study design. Overall, our findings yield much needed insights into the state, fate and dynamics of eDNA, with clear implications and guidance for future eDNA research aiming to improve species detection and quantification.
Methods
Experimental setup and sampling
Daphnia magna were cultivated for two weeks in mesocosms (initial density ~200 individuals/L) at Bangor University, UK. Animals originated from a single clone provided by Birmingham University. Environmental DNA-rich water from D. magna was collected by sieving individuals from the water using a 45 μm sieve into 3 sterilized 10L tanks containing 10L unchlorinated tap water (the same water used to culture D. magna), to serve as experimental replicates. We did not collect water chemistry for the water, which is categorized as moderately soft (low mineral), which infers a below neutral pH (~6.5â7). A fourth tank filled with 20L unchlorinated tap water served as the experimental control, which serves to detect contamination during the course of the experiment. Replicate and control tanks were placed in a walk-in climate controlled room set to 20â°C throughout the experiment. The climate room was located in a separate room and building from where the D. magna was cultivated and sieving of D. magna did not occur in the climate-controlled room to avoid potential contamination and to ensure no live animals were included in the experiment. Lids were placed on the tanks and only removed when sampling occurred for each tank. Water samples (500âml) were taken at 0, 1, 3, 6, 17 and 31 days after the start of the experiment. For each 500âml sample, we sequentially filtered the water through 10, 1, and 0.2âµm filters and retained 50âml of the final filtered effluent water (size 0) for eDNA extraction. If filters clogged during filtration, additional filters were used such that 500âml of water was passed through each filter size. Filters and water effluent samples were immediately placed in â20â°C for subsequent DNA extraction (Fig. 1). Filters were stored in 15âml conical tubes while water effluent samples were stored in 50âml conical tubes. No additional preservatives were added to the samples prior to storages at â20â°C. Extractions were performed for all samples after the experiment was completed.
DNA extraction and qPCR analysis
DNA was extracted from both the filter and retained water samples. For DNA extracted from filters, we used a modified Qiagen blood and tissue extraction protocol59,60. DNA from water samples were first concentrated by adding 3M sodium acetate (10% per water sample volume) and ethanol (2:1 ethanol to sample volume) then storing at â20â°C for 2âhours. Concentrated samples were then centrifugation at maximum speed for 30âmin to form a pellet and the effluent removed61. DNA was then extracted from the pellet using the standard Qiagen blood and tissue kit following the manufacturerâs protocol.
Quantification of extracted eDNA was performed using COI and 18S species-specific targeted qPCR assays (Table 4), developed by Primer Design Ltd (Southampton, UK). As the primers were purchased by a company, the exact design methods were not disclosed, however we can provide a brief overview of the general methods employed to generate the primer sets. In short, reference sequences for the targeted gene regions are queried for potential amplicons between 50â150âbp (e.g. using NCBI primer blast). Potential amplicon sequences are then assessed in-silico for species specificity by querying a sequence database with each potential amplicon. Once suitable amplicons are found the respective primers and probes are tested against template DNA originating from the species of interest (here D. magna) to verify amplification. It is important that any primer set be assessed for specificity within a given study. For the purposed of this experiment, we confirmed that the primers were not successful in amplifying algae originating from biofilm samples in the climate control or cultivation rooms. Furthermore, the control tanks, which utilized the same source water as the treatment tanks, served as direct comparison against non-targeted amplification. Triplicate reactions were performed for each sample and primer combination. Each 20âμL reaction contained 1âμL primer/probe mix (300ânM), 10âμL (2X) PrecisionPLUS Mastermix (Primer Design Ltd.), 2âμL DNA extract, and 7âμL DNAse-free water. Reactions were run on a QuantStudio Flex 6 Real-Time PCR System (Applied Biosystems, USA) with the following protocol: 2âmin at 95â°C, followed by 40 cycles of 10âs at 95â°C and 60âs at 60â°C. For the standard curve, 363 and 341âbp synthesized gBlock fragments (Integrated DNA Technologies Inc.) were designed for the COI and 18S primer sets, respectively, that encapsulated the targeted amplified regions of the aforementioned primer pairs. For each qPCR plate, a five-fold dilution of the appropriate standard curve dilution was added, including 1, 10, 100, 1000, 10000 copies. We also quantified the total DNA collected for each sample using a Qubit fluorimeter (Invitrogen, Waltham, USA).
In addition, we tested qPCR inhibition by spiking control samples from time points 1, 3 and 31(days) for each filter size (water effluent, 0.2âµm, 1âµm, 10âµm) as well as all samples at time day 31, which did not amplify for any samples. We added 2 ul of our 10000 size standard to 18 ul of each tested sample, resulting in a one-step dilution, which we then compared to the CT value of the 1000 size standard. Inhibition in the control samples spiked size standard was negligible (<3% divergence)31 with percent divergence being on average 0.633% (SDâ=â0.56%) among control tanks, 0.973% (SDâ=â0.415%) among time 31 tanks and 0.962% (SDâ=â0.435%) across each filter size. We also checked COI inhibition by spiking control tank samples, as above, but were limited due to low remaining extraction volume, but still found negligible inhibition with percent divergence for across control tanks being 1.25% (SDâ=â0.117) and across filter types being 1.69% (SDâ=â1.25).
Statistical analysis
For statistical analyses, we did not include eDNA quantities where no amplification occurred across the group (filter sizes 1âµm and10 µm and days 17 and 31) to avoid zero inflation of the data and potential type I errors. All statistical analyses employed R, version 3.5.162. To assess the effect of filter size and gene marker choice (i.e. the response variables) on DNA quantification (i.e. the explanatory variable), we used a repeated measures ANOVA via the gls function in the nlme package, with a AR1 correlation structure to account for the temporal autocorrelation for each tank across time (i.e. the repeated measure)47. Generalized least squares (gls) allows within group residual structuring to directly model variance-covariance structure for the response by extending the ordinary least squares model used in linear regression47. The gls residuals structuring is particularly useful for dealing with temporal or spatial autocorrelation, which can lead to type I errors when unaccounted for in linear regression models63. The AR1 implements an autoregressive correlation structure for a sequential set of integers in a supplied vector, which makes it suitable for structuring residual spreads over times64. We assessed the fitness of the gls model by comparing the fitted model to the null and parsimonious versions of the model as well as alternative fitted linear model using AIC47. Additionally, we assessed temporal dynamics of the water effluent and 0.2âµm filter data over time using a first order decay model for each set of filter and gene marker data combination65.
Data Availability
Data associated with the study are available on Figshare; https://doi.org/10.6084/m9.figshare.9699143.
Change history
04 August 2021
A Correction to this paper has been published: https://doi.org/10.1038/s41598-021-95567-6
References
Lawson Handley, L. How will the âmolecular revolutionâ contribute to biological recording? Biol. J. Linn. Soc. 115, 750â766 (2015).
Seymour, M. Rapid progression and future of environmental DNA research. Commun. Biol. 2, 80 (2019).
Deiner, K. et al. Environmental DNA metabarcoding: Transforming how we survey animal and plant communities. Mol. Ecol., https://doi.org/10.1111/mec.14350 (2017).
Foote, A. D. et al. Investigating the Potential Use of Environmental DNA (eDNA) for Genetic Monitoring of Marine Mammals. PLoS One 7, e41781 (2012).
Rees, H. C., Maddison, B. C., Middleditch, D. J., Patmore, J. R. M. & Gough, K. C. REVIEW: 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).
Baird, D. J. & Hajibabaei, M. Biomonitoring 2.0: a new paradigm in ecosystem assessment made possible by next-generation DNA sequencing. Mol. Ecol. 21, 2039â2044 (2012).
Kelly, R. P. et al. Harnessing DNA to improve environmental management. Science (80-.). 344, 1455â1456 (2014).
Lodge, D. M. et al. Conservation in a cup of water: estimating biodiversity and population abundance from environmental DNA. Mol. Ecol. 21, 2555â2558 (2012).
McDevitt, A. D. et al. Environmental DNA metabarcoding as an effective and rapid tool for fish monitoring in canals. J. Fish Biol. 0 (2019).
Hänfling, B. et al. Environmental DNA metabarcoding of lake fish communities reflects long-term data from established survey methods. Mol. Ecol. 25, 3101â3119 (2016).
Bista, I. et al. Annual time-series analysis of aqueous eDNA reveals ecologically relevant dynamics of lake ecosystem biodiversity. Nat. Commun. 8, 14087 (2017).
Jerde, C. L. et al. Detection of Asian carp DNA as part of a Great Lakes basin-wide surveillance program. Can. J. Fish. Aquat. Sci. 70, 522â526 (2013).
Biggs, J. et al. Using eDNA to develop a national citizen science-based monitoring programme for the great crested newt (Triturus cristatus). Biol. Conserv. 183, 19â28 (2015).
Weltz, K. et al. Application of environmental DNA to detect an endangered marine skate species in the wild. PLoS One 12, e0178124 (2017).
Deiner, K. & Altermatt, F. Transport Distance of Invertebrate Environmental DNA in a Natural River. PLoS One 9, e88786 (2014).
Tsuji, S. et al. Real-time multiplex PCR for simultaneous detection of multiple species from environmental DNA: an application on two Japanese medaka species. Sci. Rep. 8, 9138 (2018).
Yamamoto, S. et al. Environmental DNA metabarcoding reveals local fish communities in a species-rich coastal sea. Sci. Rep. 7, 40368 (2017).
Pawlowski, J. et al. The future of biotic indices in the ecogenomic era: Integrating (e)DNA metabarcoding in biological assessment of aquatic ecosystems. Sci. Total Environ. 637â638, 1295â1310 (2018).
Deiner, K., Fronhofer, E. A., Mächler, E., Walser, J.-C. & Altermatt, F. Environmental DNA reveals that rivers are conveyer belts of biodiversity information. Nat. Commun. 7, 12544 (2016).
Barnes, M. A. & Turner, C. R. The ecology of environmental DNA and implications for conservation genetics. Conserv. Genet. 17, 1â17 (2016).
DellâAnno, A. & Danovaro, R. Extracellular DNA Plays a Key Role in Deep-Sea Ecosystem Functioning. Science (80-.). 309, 2179 LP â 2179 (2005).
Pietramellara, G. et al. Extracellular DNA in soil and sediment: fate and ecological relevance. Biol. Fertil. Soils 45, 219â235 (2009).
Torti, A., Lever, M. A. & Jørgensen, B. B. Origin, dynamics, and implications of extracellular DNA pools in marine sediments. Mar. Genomics 24, 185â196 (2015).
Levy-Booth, D. J. et al. Cycling of extracellular DNA in the soil environment. Soil Biol. Biochem. 39, 2977â2991 (2007).
Flindt, R. Amazing numbers in biology. (Springer Science & Business Media, 2006).
Wilcox, T. M., McKelvey, K. S., Young, M. K., Lowe, W. H. & Schwartz, M. K. Environmental DNA particle size distribution from Brook Trout (Salvelinus fontinalis). Conserv. Genet. Resour. 7, 639â641 (2015).
Matsui, K., Ishii, N. & Honjo, M. Use of the SYBR Green I fluorescent dye and a centrifugal filter device for rapid determination of dissolved DNA concentration in fresh water. Aquat. Microb. Ecol. 36, 99â105 (2004).
Turner, C. R. et al. Particle size distribution and optimal capture of aqueous macrobial eDNA. Methods Ecol. Evol. 5, 676â684 (2014).
Piggott, M. P. Evaluating the effects of laboratory protocols on eDNA detection probability for an endangered freshwater fish. Ecol. Evol. 6, 2739â2750 (2016).
Muha, T. P., Robinson, C. V., Garcia de Leaniz, C. & Consuegra, S. An optimised eDNA protocol for detecting fish in lentic and lotic freshwaters using a small water volume. PLoS One 14, e0219218 (2019).
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).
Stewart, K. A. Understanding the effects of biotic and abiotic factors on sources of aquatic environmental DNA. Biodivers. Conserv. 28, 983â1001 (2019).
Foran, D. R. Relative Degradation of Nuclear and Mitochondrial DNA: An Experimental Approach*. J. Forensic Sci. 51, 766â770 (2006).
Minamoto, T. et al. Nuclear internal transcribed spacer-1 as a sensitive genetic marker for environmental DNA studies in common carp Cyprinus carpio. Mol. Ecol. Resour. 17, 324â333 (2017).
Taberlet, P., Coissac, E., Hajibabaei, M. & Rieseberg, L. H. Environmental DNA. Mol. Ecol. 21, 1789â1793 (2012).
Dysthe, J. C., Franklin, T. W., McKelvey, K. S., Young, M. K. & Schwartz, M. K. An improved environmental DNA assay for bull trout (Salvelinus confluentus) based on the ribosomal internal transcribed spacer I. PLoS One 13, e0206851 (2018).
Holland, M. M. & Parsons, T. J. Mitochondrial DNA Sequence Analysis - Validation and Use for Forensic Casework. Forensic Sci. Rev. 11, 21â50 (1999).
Toju, H., Tanabe, A. S., Yamamoto, S. & Sato, H. High-Coverage ITS Primers for the DNA-Based Identification of Ascomycetes and Basidiomycetes in Environmental Samples. PLoS One 7, e40863 (2012).
Buxton, A. S., Groombridge, J. J. & Griffiths, R. A. Is the detection of aquatic environmental DNA influenced by substrate type? PLoS One 12, e0183371 (2017).
Dejean, T. et al. Persistence of Environmental DNA in Freshwater Ecosystems. PLoS One 6, e23398 (2011).
Goldberg, C. S., Sepulveda, A., Ray, A., Baumgardt, J. & Waits, L. P. Environmental DNA as a new method for early detection of New Zealand mudsnails (Potamopyrgus antipodarum). Freshw. Sci. 32, 792â800 (2013).
Wilcox, T. M. et al. Understanding environmental DNA detection probabilities: A case study using a stream-dwelling char Salvelinus fontinalis. Biol. Conserv. 194, 209â216 (2016).
Seymour, M. et al. Acidity promotes degradation of multi-species environmental DNA in lotic mesocosms. Commun. Biol. 1, 4 (2018).
Collins, R. A. et al. Persistence of environmental DNA in marine systems. Commun. Biol. 1, 185 (2018).
Levi, T. et al. Environmental DNA for the enumeration and management of Pacific salmon. Mol. Ecol. Resour. 19, 597â608 (2019).
Jerde, C. L. et al. Influence of Stream Bottom Substrate on Retention and Transport of Vertebrate Environmental DNA. Environ. Sci. Technol. 50, 8770â8779 (2016).
Pinheiro, J. C. & Bates, D. M. Mixed-effects models in S and S-plus. (Springer, 2000).
Goldberg, C. S., Strickler, K. M. & Fremier, A. K. Degradation and dispersion limit environmental DNA detection of rare amphibians in wetlands: Increasing efficacy of sampling designs. Sci. Total Environ. 633, 695â703 (2018).
Jane, S. F. et al. Distance, flow and PCR inhibition: eDNA dynamics in two headwater streams. Mol. Ecol. Resour. 15, 216â227 (2015).
Goldberg, C. S., Strickler, K. M. & Pilliod, D. S. Moving environmental DNA methods from concept to practice for monitoring aquatic macroorganisms. Biol. Conserv. 183, 1â3 (2015).
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).
Alaeddini, R., Walsh, S. J. & Abbas, A. Forensic implications of genetic analyses from degraded DNAâA review. Forensic Sci. Int. Genet. 4, 148â157 (2010).
Matsui, K. & Honjo, M. Estimation of the fate of dissolved DNA in thermally stratified lake water from the stability of exogenous plasmid DNA. Aquat. Microb. Ecol. 26, 95â102 (2001).
Murgia, M. et al. Mitochondrial DNA is not fragmented during apoptosis. J. Biol. Chem. 267, 10939â10941 (1992).
Tepper, C. G. & Studzinski, G. P. Resistance of mitochondrial DNA to degradation characterizes the apoptotic but not the necrotic mode of human leukemia cell death. J. Cell. Biochem. 52, 352â361 (1993).
Rickwood, D. & Chambers, J. A. A. Evidence for protected regions of DNA in the mitochondrial nucleoid of Saccharomyces cerevisiae. FEMS Microbiol. Lett. 12, 187â190 (1981).
Evans, N. T. & Lamberti, G. A. Freshwater fisheries assessment using environmental DNA: A primer on the method, its potential, and shortcomings as a conservation tool. Fish. Res. 197, 60â66 (2018).
Sutherland, W. J. et al. Horizon scan of global conservation issues for 2011. Trends Ecol. Evol. 26, 10â16 (2011).
Kirshtein, J. D., Anderson, C. W., Wood, J. S., Longcore, J. E. & Voytek, M. A. Quantitative PCR detection of Batrachochytrium dendrobatidis DNA from sediments and water. Dis. Aquat. Organ. 77, 11â15 (2007).
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).
Ficetola, G. F., Miaud, C., Pompanon, F. & Taberlet, P. Species detection using environmental DNA from water samples. Biol. Lett. 4, 423â425 (2008).
R Core Team. R: A Language and Environment for Statistical Computing (2018).
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, 2009).
Legendre, P. & Legendre, L. Numerical ecology. (Elsevier Science BV, 2012).
Crawley, M. J. The R book. (John Wiley & Sons, 2012).
Acknowledgements
Funding from Bangor University and NERC Highlight Topic grant NE/N006216/1. Thank you to Mike Hayle for assistance with climate control and experimental setup, Simon Webster for D. magna cultivation advise, Luisa Orsini for constructive feedback, John Colbourne for initial provision of the D. magna clones, and two anonymous reviewers for their valued input.
Author information
Authors and Affiliations
Contributions
M.S., S.C. and G.W. designed the experiment. G.W. conducted the experiment; R.M., G.W. and M.S. performed the laboratory analyses. M.S. and R.M. performed the statistical analyses. M.S., R.M., G.R.C. and S.C. drafted the manuscript. M.S. coordinated the editing and revisions. All authors comment on the final draft.
Corresponding author
Ethics declarations
Competing Interests
The authors declare no competing interests.
Additional information
Publisherâs note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The original online version of this Article was revised: This original version of this Article contained errors. Modifications have been made to the legends of Figure 2 and 3, and the headings in Table 4. Full information regarding the corrections made can be found in the correction for this Article.
Rights and permissions
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/.
About this article
Cite this article
Moushomi, R., Wilgar, G., Carvalho, G. et al. Environmental DNA size sorting and degradation experiment indicates the state of Daphnia magna mitochondrial and nuclear eDNA is subcellular. Sci Rep 9, 12500 (2019). https://doi.org/10.1038/s41598-019-48984-7
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-019-48984-7
This article is cited by
-
Detection, movement and persistence of invertebrate eDNA in groundwater
Scientific Reports (2024)
-
Larger particle size distribution of environmental RNA compared to environmental DNA: a case study targeting the mitochondrial cytochrome b gene in zebrafish (Danio rerio) using experimental aquariums
The Science of Nature (2024)
-
Simple extraction and analysis of environmental DNA using glass fibers in suspension form
Limnology (2023)
-
Persistence and degradation dynamics of eDNA affected by environmental factors in aquatic ecosystems
Hydrobiologia (2022)
-
Environmental DNA Methods for Ecological Monitoring and Biodiversity Assessment in Estuaries
Estuaries and Coasts (2022)