The effects of spatial and temporal replicate sampling on eDNA metabarcoding

Background The heterogeneous nature of environmental DNA (eDNA) and its effects on species detection and community composition estimates has been highlighted in several studies in the past decades. Mostly in the context of spatial distribution over large areas, in fewer occasions looking at spatial distribution within a single body of water. Temporal variation of eDNA, similarly, has mostly been studied as seasonality, observing changes over large periods of time, and often only for small groups of organisms such as fish and amphibians. Methods We analyzed and compared small-scale spatial and temporal variation by sampling eDNA from two small, isolated dune lakes for 20 consecutive weeks. Metabarcoding was performed on the samples using generic COI primers. Molecular operational taxonomic unit (MOTUs) were used to assess dissimilarities between spatial and temporal replicates. Results Our results show large differences between samples taken within one lake at one point in time, but also expose the large differences between temporal replicates, even those taken only 1 week apart. Furthermore, between-site dissimilarities showed a linear correlation with time frame, indicating that between-site differences will be inflated when samples are taken over a period of time. We also assessed the effects of PCR replicates and processing strategies on general patterns of dissimilarity between samples. While more inclusive PCR replicate strategies lead to higher richness estimations, dissimilarity patterns between samples did not significantly change. Conclusions We conclude that the dissimilarity of temporal replicates at a 1 week interval is comparable to that of spatial replicate samples. It increases, however, for larger time intervals, which suggests that population turnover effects can be stronger than community heterogeneity. Spatial replicates alone may not be enough for optimal recovery of taxonomic diversity, and cross-comparisons of different locations are susceptible to inflated dissimilarities when performed over larger time intervals. Many of the observed MOTUs could be classified as either phyto- or zooplankton, two groups that have gained traction in recent years as potential novel bio-indicator species. Our results, however, indicate that these groups might be susceptible to large community shifts in relatively short periods of time, highlighting the need to take temporal variations into consideration when assessing their usability as water quality indicators.