Skip to main content
  • Submit
  • About
    • Editorial Board
    • PNAS Staff
    • FAQ
    • Accessibility Statement
    • Rights and Permissions
    • Site Map
  • Contact
  • Journal Club
  • Subscribe
    • Subscription Rates
    • Subscriptions FAQ
    • Open Access
    • Recommend PNAS to Your Librarian
  • Log in
  • My Cart

Main menu

  • Home
  • Articles
    • Current
    • Special Feature Articles - Most Recent
    • Special Features
    • Colloquia
    • Collected Articles
    • PNAS Classics
    • List of Issues
  • Front Matter
  • News
    • For the Press
    • This Week In PNAS
    • PNAS in the News
  • Podcasts
  • Authors
    • Information for Authors
    • Editorial and Journal Policies
    • Submission Procedures
    • Fees and Licenses
  • Submit
  • About
    • Editorial Board
    • PNAS Staff
    • FAQ
    • Accessibility Statement
    • Rights and Permissions
    • Site Map
  • Contact
  • Journal Club
  • Subscribe
    • Subscription Rates
    • Subscriptions FAQ
    • Open Access
    • Recommend PNAS to Your Librarian

User menu

  • Log in
  • My Cart

Search

  • Advanced search
Home
Home

Advanced Search

  • Home
  • Articles
    • Current
    • Special Feature Articles - Most Recent
    • Special Features
    • Colloquia
    • Collected Articles
    • PNAS Classics
    • List of Issues
  • Front Matter
  • News
    • For the Press
    • This Week In PNAS
    • PNAS in the News
  • Podcasts
  • Authors
    • Information for Authors
    • Editorial and Journal Policies
    • Submission Procedures
    • Fees and Licenses

New Research In

Physical Sciences

Featured Portals

  • Physics
  • Chemistry
  • Sustainability Science

Articles by Topic

  • Applied Mathematics
  • Applied Physical Sciences
  • Astronomy
  • Computer Sciences
  • Earth, Atmospheric, and Planetary Sciences
  • Engineering
  • Environmental Sciences
  • Mathematics
  • Statistics

Social Sciences

Featured Portals

  • Anthropology
  • Sustainability Science

Articles by Topic

  • Economic Sciences
  • Environmental Sciences
  • Political Sciences
  • Psychological and Cognitive Sciences
  • Social Sciences

Biological Sciences

Featured Portals

  • Sustainability Science

Articles by Topic

  • Agricultural Sciences
  • Anthropology
  • Applied Biological Sciences
  • Biochemistry
  • Biophysics and Computational Biology
  • Cell Biology
  • Developmental Biology
  • Ecology
  • Environmental Sciences
  • Evolution
  • Genetics
  • Immunology and Inflammation
  • Medical Sciences
  • Microbiology
  • Neuroscience
  • Pharmacology
  • Physiology
  • Plant Biology
  • Population Biology
  • Psychological and Cognitive Sciences
  • Sustainability Science
  • Systems Biology
Research Article

Estimating species distribution and abundance in river networks using environmental DNA

Luca Carraro, Hanna Hartikainen, Jukka Jokela, Enrico Bertuzzo, and Andrea Rinaldo
PNAS November 13, 2018 115 (46) 11724-11729; first published October 29, 2018; https://doi.org/10.1073/pnas.1813843115
Luca Carraro
aLaboratory of Ecohydrology, École Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, Switzerland;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • For correspondence: andrea.rinaldo@epfl.ch luca.carraro@epfl.ch
Hanna Hartikainen
bAquatic Ecology Group, Swiss Federal Institute of Aquatic Science and Technology, CH-8600 Dübendorf, Switzerland;
cInstitute of Integrative Biology, Eidgenössiche Technische Hochschule Zürich, CH-8092 Zurich, Switzerland;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Jukka Jokela
bAquatic Ecology Group, Swiss Federal Institute of Aquatic Science and Technology, CH-8600 Dübendorf, Switzerland;
cInstitute of Integrative Biology, Eidgenössiche Technische Hochschule Zürich, CH-8092 Zurich, Switzerland;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Enrico Bertuzzo
dDipartimento di Scienze Ambientali, Informatica e Statistica, Università Ca’ Foscari Venezia, IT-30170 Venice, Italy;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Andrea Rinaldo
aLaboratory of Ecohydrology, École Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, Switzerland;
eDipartimento di Ingegneria Civile, Edile ed Ambientale, Università di Padova, IT-35131 Padua, Italy
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • For correspondence: andrea.rinaldo@epfl.ch luca.carraro@epfl.ch
  1. Contributed by Andrea Rinaldo, September 26, 2018 (sent for review August 15, 2018; reviewed by Gentile Francesco Ficetola and Christopher Jerde)

  • Article
  • Figures & SI
  • Info & Metrics
  • PDF
Loading

Significance

Organisms leave traces of DNA in their environment [environmental DNA (eDNA)], such as cells in mucus or feces. When extracted from water or soil, eDNA can be used to track the presence of a target species or the composition of entire communities. In rivers, eDNA dynamics are modulated by transport and decay. Here, we use hydrologically based models to reconstruct the upstream distribution and abundance of target species throughout a river network from eDNA measurements. We validate our method by estimating the catchment-wide biomass distribution of a sessile invertebrate and its parasite, causing disease in salmonids. This work unlocks the power of eDNA for monitoring biodiversity across broad geographies in a way hitherto unfeasible with traditional survey approaches.

Abstract

All organisms leave traces of DNA in their environment. This environmental DNA (eDNA) is often used to track occurrence patterns of target species. Applications are especially promising in rivers, where eDNA can integrate information about populations upstream. The dispersion of eDNA in rivers is modulated by complex processes of transport and decay through the dendritic river network, and we currently lack a method to extract quantitative information about the location and density of populations contributing to the eDNA signal. Here, we present a general framework to reconstruct the upstream distribution and abundance of a target species across a river network, based on observed eDNA concentrations and hydro-geomorphological features of the network. The model captures well the catchment-wide spatial biomass distribution of two target species: a sessile invertebrate (the bryozoan Fredericella sultana) and its parasite (the myxozoan Tetracapsuloides bryosalmonae). Our method is designed to easily integrate general biological and hydrological data and to enable spatially explicit estimates of the distribution of sessile and mobile species in fluvial ecosystems based on eDNA sampling.

  • species distribution model
  • ecohydrology
  • proliferative kidney disease

Environmental DNA (eDNA), present as loose fragments, as shed cells (1, 2), or in microscopic organisms (3, 4), can be extracted from matrices such as water or soil and used to track the presence of target species or the composition of entire communities (5, 6). Approaches using eDNA for qualitative species detection have proved their value in management and conservation, improving the measurement of biodiversity in a replicable and consistent manner (7) and facilitating the detection of rare, invasive, or parasitic species (8⇓⇓⇓⇓–13).

Environmental DNA in river water carries a record of the species present upstream, but the interpretation of this signal is a complex issue (14). Once released to the environment, eDNA undergoes selective decay. Nucleic acids incur progressive damage [e.g., due to biological activity, temperature, or pH (15, 16)] during hydrological advection, retention, and resuspension (17, 18). These processes result in alterations that affect eDNA detection in environmental samples. The magnitude of the decay is highly dependent on the nature of the flow regime and the substrate type (19). Furthermore, eDNA has polydisperse properties due to its origin from diverse organic sources (e.g., spores, cells, tissues, feces), which complicates the evaluation of decay rates (20). The eDNA sampled at any point within a dendritic network of sources is the outcome of diffuse eDNA release from points upstream, modified by decay processes during transport that are governed by network connectivity, in which each path to the observation point may be described by different hydro-morphological conditions. As a result, while it is straightforward to link a positive PCR test with the presence of the target species at some (unknown) distance upstream, quantification of species densities and the location of populations is currently impossible because, besides a number of potentially confounding factors affecting eDNA shedding [e.g., animal behavior, movement, physiology, and size (21, 22)], it requires consideration of the effects of the dynamics of eDNA transport along river branches and the deconvolution of the hierarchical aggregation of the various network branches. Here, we establish a generally applicable framework to interpret quantitative eDNA point measurements in rivers and relate them to the spatial distribution of the DNA sources, jointly with estimates of the density distribution of the target species throughout the river basin.

The proposed framework stems from fundamental mass balance relationships. It is intended for use in river networks discretized into “nodes,” i.e., river stretches of suitable length within which hydrological conditions, as well as the target species density and hence its eDNA production, can be considered homogeneous (23). Within such nodes, the basis for the spatially explicit model contrasting measured eDNA concentrations is given byCj^=1Qj¯∑i∈γ(j)AS,i⁡exp−Lijvij¯τpi,[1]where Cj^ [NL−3] is the eDNA concentration at node j predicted by the model; Qj¯ [L3 T−1] is a characteristic (e.g., median) water discharge at node j; γ(j) indexes the set of nodes upstream of j connected by the river network; AS,i [L2] is the source area (or the spatial extent of the habitat) of the target species pertaining to node j; Lij [L] is the along-stream length of the path connecting i to j; τ [T] is the inverse of the decay rate (i.e., a characteristic decay time) subsuming the damage rates to genetic material experienced during hydrologic transport (18, 19), possibly constant for all eDNA fragments irrespective of the different hydrological and environmental conditions along all paths to the sampling site (24); vij¯ [LT−1] is the average flow velocity along the path connecting i to j; and pi [NL−2 T−1] is the eDNA production at i per unit habitat area and unit time. We assume that pi is proportional to the target species density at node i (12, 24). In the case of sessile species dwelling in riverine habitats, AS,i can be considered equal to the riverbed area of the river stretch i; for terrestrial species, such as certain diatoms that colonize unchanneled areas and have been used to track surface runoff (25), AS,i=AL,i, where AL,i is the directly contributing area (26) to node i (that is, Aj=∑i∈γ(j)AL,i). Eq. 1 is based on the hypothesis that eDNA undergoes first-order exponential decay along the downstream path from the source i to the measurement node j (17). As an alternative characterization of decay, we may introduce the parameterization λ=vij¯τ, where λ is a decay length that is assumed to be constant irrespective of hydrological regimes and heterogeneities in morphological conditions across the watershed (23). This alternative formulation, allowing one to avoid the effective calculation of vij¯, is justified by the observation that water velocities in catchments for a large range of flow regimes generally show modest longitudinal gradients, as high velocities potentially prompted by steeper slopes in the upper reaches are limited by increased flow resistance (27). In previous studies, estimated decay lengths range from the order of magnitude of a few meters for experimental flumes (18) to that of kilometers observed in real catchments (19). Such scale dependence seems to indicate a role for the correlation scale of heterogeneous features, relevant to eDNA transport, likely to obey scaling arguments in river networks (28). Suffice here to note that all estimates of decay times (converted to travel lengths as noted above) from field measurements (19, 23) and mesocosm experiments (17, 24) point to values much greater than the mean distance between significant confluences in real catchments, at least in runoff-producing areas (26). This implies that a plurality of sources could be contributing to detectable eDNA concentrations at each measurement site, thus prompting the need to resort to approaches like Eq. 1 that take into account the structure of the network.

While not all of the above assumptions are equally valid in the general case (see SI Appendix for a digression on caveats on model assumptions), the framework represents a flexible general theory of spatially explicit eDNA source tracking. Some of these assumptions may be easily relaxed. For instance, the purely convective treatment of the decay of genetic material outlined above may be the subject of more refined formulations of the transport problem. In particular, travel times from source to measurement site can be made explicitly dependent on the hydrodynamic and geomorphological dispersion induced by the hierarchical nature of the network (27).

To derive estimates of biomass distribution, we coupled the general source area model described above with a species distribution model. This combined approach provides a versatile method to capture the influence of the catchment-wide ecological, hydro-morphological, or geological drivers promoting eDNA production (and reflecting species density) within the defined river stretches. Local eDNA production can be expressed by means of the exponential linkpi=p0⁡exp[βTX(i)],[2]where p0 is a baseline eDNA production value constant in space, X(i) is a vector of covariates evaluated at site i, and β is a vector of parameters requiring calibration that identify the effect of such covariates on eDNA production and thereby on the distribution of the target species. Covariates included in the vector X(i) will depend on the particular system and data available. Their extent can be either local (e.g., pertaining only to the area directly contributing to a given stretch) or nonlocal, i.e., related to the whole catchment area upstream of the stretch.

We tested the framework with joint field measurements of eDNA concentrations of the myxozoan parasite Tetracapsuloides bryosalmonae and its primary host, the freshwater bryozoan Fredericella sultana, across various locations within the Wigger watershed (Switzerland) (Fig. 1). T. bryosalmonae is the causative agent of proliferative kidney disease (PKD), a high-mortality disease affecting salmonid fish populations. PKD is recognized as one of the leading causes of declines in brown trout populations in Europe, it affects diverse salmonid populations in North America, and it is a major aquaculture disease (29, 30). The water samples used for eDNA detection of both F. sultana and T. bryosalmonae were collected at roughly monthly intervals at 15 sites during 12 mo (one 500-mL sample per sampling occasion and site) (Fig. 1C). T. bryosalmonae eDNA is likely to be largely derived from spores shed into the environment. Parasite spores, released into water by infected bryozoans, infect brown trout through skin and gills and proliferate in the kidney. To complete the life cycle, spores infective to bryozoans are excreted in the urine of infected fish. These two types of spores are genetically indistinguishable yet differentiated in terms of function, which poses further challenges for modeling. The T. bryosalmonae eDNA concentration may thus be a product of the genomic contents of the two types of spore originating from very different transport sources, i.e., from an immobile source (bryozoans) coupled with a mobile source (fish). In this particular case, a comparative analysis of field-measured eDNA for both F. sultana (sessile source of eDNA) and T. bryosalmonae (eDNA that could jointly originate from sessile and mobile hosts) will prove particularly instructive as a demonstration of the potential of the proposed framework.

Fig. 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 1.

Environmental DNA sampling of F. sultana and T. bryosalmonae. (A) Map and color-coded digital elevation map of the study region showing the extracted river network and location of the eDNA sampling sites. (B) Location of the study region within Switzerland. (C) Measured eDNA concentrations of F. sultana and T. bryosalmonae at the 15 sampling sites during the period May 2014–May 2015 (LOQ, limit of quantification).

We implemented and calibrated the model as described in Materials and Methods. Note that we explicitly accounted for the possibility that samples with low eDNA concentration may be interpreted as zeros (Fig. 1C) by introducing the nondetection probability ϕj, a monotonically decreasing function of Cj^. The covariates included in vector X(i) were local elevation, contributing area, and the fractions of contributing area covered by moraine, peat, or superficial water (e.g., lakes, ponds, or wetlands) upstream of site i. These five covariates were chosen as representative of morphological and geological features of the catchment (23). Overall, our model proved to be considerably good at reproducing the observed data, with regard to both observed eDNA concentrations (Figs. 2 and 3A) and frequency of nondetection of the target eDNA (compared with ϕj distributions in SI Appendix, Fig. S1). The median eDNA concentration of the different sampling sites (Fig. 3A) is well reproduced by the model (Pearson’s correlation coefficient r=0.84 for F. sultana, r=0.75 for T. bryosalmonae). According to two-sample Kolmogorov–Smirnov tests (31) (Fig. 2), the null hypothesis that observed positive eDNA concentrations are drawn from the predicted distributions cannot be rejected at the 5% confidence level for 33% (F. sultana) and 71% (T. bryosalmonae) of the sampling sites. At the 1% confidence level, these percentages rise to 73% and 100%, respectively.

Fig. 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 2.

Comparison between the observed cumulative frequencies of measured eDNA concentrations at the 15 sampling sites for F. sultana (Fs; purple) and T. bryosalmonae (Tb; green) and the cumulative distribution function obtained by the model. α values (color coded to match the solid lines) indicate the confidence level at which, according to a two-sample Kolmogorov–Smirnov test, the null hypothesis that the two samples (modeled and observed) come from the same distribution cannot be rejected. Higher values of α indicate a better fit. Tested values for α were 0.05, 0.01, 0.005, 0.001. N/A, not applicable (i.e., at site 15 no positive values of eDNA concentration for Tb were detected). Note that some of the cases where the two distributions are not different at α=0.05 are characterized by a limited number of positive eDNA detections.

Fig. 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 3.

Correlation between eDNA concentrations, predicted species density, and covariates. (A) Correlation between observed and modeled eDNA concentrations for F. sultana (Fs; purple) and T. bryosalmonae (Tb; green). Each circle corresponds to a sampling site. x and y values, respectively, are taken from the medians of the observed and modeled cumulative distribution functions displayed in Fig. 2. (B) Correlation between medians of predicted F. sultana and T. bryosalmonae eDNA production. Each circle represents one river stretch. The trend line is displayed in red. Pearson’s correlation coefficients r are also reported. (C) Correlations between medians of predicted eDNA production for the two target species and the covariates representing the contributing catchment area and fraction of the upstream catchment covered by moraine. Correlations with other covariates are reported in SI Appendix, Fig. S2B; maps of covariate values are presented in SI Appendix, Fig. S3. (D) Posterior distributions for the decay time τ and values of the β coefficients associated with the covariates reported on the axis labels. The posterior distributions for the other covariates are shown in SI Appendix, Fig. S2A.

Predicted maps of eDNA production for F. sultana (Fig. 4A) identify the southeastern portion of the watershed as a hotspot for bryozoans. This is mainly due to the positive correlation between the presence of moraines and the production of bryozoan eDNA (Fig. 3 C and D). This correlation was uncovered in previous work (23), although it is not clear yet what the underlying causes of this positive association are. The upstream bryozoan reservoir explains the eDNA concentration patterns observed at the downstream sites (Fig. 1C), because the estimated values of the decay time (Fig. 3D) allow for the detection of eDNA material at distances comparable with the maximum length from source to outlet of the river system. The predicted distribution of T. bryosalmonae (Fig. 4B) mirrors that of F. sultana in headwaters and upper sites, while higher values of production are estimated toward the outlet. This is shown by a positive shift in the posterior distribution of the parameter expressing the effect of contributing area (Fig. 3 C and D). The correlation between predicted densities of F. sultana and T. bryosalmonae is strong (Fig. 3B) and suggests that bryozoans release disproportionately more spores compared with fish hosts. Thus, we suggest that a full description of the spatial distribution of PKD-infected fish (23) might be unnecessary to understand the bulk of the distribution of the parasite sources when the eDNA signal is dominated by locally abundant colonies of infected bryozoans like in the case at hand. We argue that the much stronger correlation observed between predicted densities of T. bryosalmonae and contributing area in comparison with that for F. sultana (Fig. 3C) posits that the density of overtly infected (i.e., spore producing) bryozoans tends to increase along downstream directions. In fact, the positive correlation between PKD prevalence and total contributing area in river networks has been shown to be a byproduct of network connectivity (32). The fact that T. bryosalmonae is mostly shed by bryozoans rather than fish seems plausible in this specific host–parasite system, as parasite maturation in fish kidney tubuli is observed relatively rarely, compared with the prolific spore production within large parasite sacs inside the bryozoan host (33). We note further that the spatial match of the bryozoan and fish populations is unlikely to drive this relationship as the biomass of fish is expected to be higher in deeper, more downstream sections (34). Finally, estimated median values of the decay time (Fig. 3D) were 4.0 h for T. bryosalmonae (with a 25–75% range of the posterior distribution of 2.7–7.0 h) and 6.9 h (25–75% range: 5.0–11.1 h) for F. sultana, corresponding to decay lengths of 14 km (25–75% range: 10–25 km) and 25 km (25–75% range: 18–40 km), respectively (obtained by assuming an average flow velocity of 1 ms−1), in agreement with previous findings (19).

Fig. 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 4.

Predicted species distributions. (A–D) Maps of predicted production p (A and B) and concentration C^ (C and D) of eDNA for F. sultana (Fs) and T. bryosalmonae (Tb). Estimates were obtained as the medians of the distributions predicted by the model.

Our framework proved capable of interpreting both eDNA data incidentally shed from benthic populations (F. sultana, with likely sources of eDNA from fecal pellets and sloughed cells) and eDNA from spores released into the water (T. bryosalmonae). Although different forms of eDNA may be differently impacted by environmental factors such as temperature and pH, the choice of formulation involving a single parameter expressing the decay time for both species appeared satisfactory for capturing the integrated eDNA transport dynamics at the catchment scale. Further considerations on this aspect are presented in SI Appendix. Another strength of our approach is the possibility of applying adequate parameterizations for p in an explicit manner (as in Eq. 2), to accommodate the nature of the link between the target species density and its biological and environmental filters along hydrologic pathways, such as the environmental conditions or the density of species with which it interacts. Such parameterization provides a simple and versatile means to assimilate field data and to integrate population or species distribution models.

Accurate field validations of the current assumptions are needed to generalize this framework, and this could be achieved by relatively simple experimental designs. For instance, the displacement and decay of genetic material from nonnative known biomasses placed in well-differentiated positions (say, within a catchment where hydrologic and geomorphologic drivers are known) could be key. Subsequent sampling at downstream sites, where eDNA would be contributed by sources at known distances, could then be used to assess the strengths and weaknesses of each assumption underlying the proposed approach.

Tracking the source area and the local biomass density of target species via downstream eDNA measurement is possible, provided that a suitable spatially explicit framework is used to interpret the field data. Key is accounting for the filtering produced by the progressive damage occurring during hydrological transport and harnessing it to recover spatial information on species distributions. The integration of quantitative eDNA measurements, hydro-geomorphological scaling, and ecological models presented here reveals another direction in ecohydrological studies by unlocking the huge potential of remote monitoring using eDNA.

Materials and Methods

eDNA Data Collection.

Stream water samples were collected in 15 locations (Fig. 1) along the river network of the Wigger watershed, Switzerland. For each site, a total of 21 500-mL samples were taken at approximately biweekly (or monthly during December, January, and February) intervals (except site 5, a connected pond which was artificially drained after 12 samples were taken).

Presterilized (10% bleach followed by UV-B treatment) plastic bottles were used to collect water from the river by submerging the bottle with a gloved hand. The samples were transported to the laboratory on ice and filtered using gentle vacuum within the same day onto 5-cm diameter, 0.45-μm pore size individually packaged sterile membrane filters (Merck Millipore). A vacuum pump with a borosilicate glass filtration setup was used and sterilized between each sample in 10% bleach followed by three clean water rinses. Negative controls were created by filtering MilliQ water through a sterile filter at the start and end of each filtration session, as well as once during the filtration (after sample 7). Filter papers were placed in 2-mL bead beating tubes (obtained from the kit described below) and frozen at −80 oC until extraction. Before extraction, filter papers were cut with sterilized scissors to break them up. eDNA was extracted from all filter papers, including controls, using a PowerSoil DNA kit (MO BIO Laboratories) in a dedicated clean laboratory (free of PCR products). The kit includes a bead beating step and a separate inhibitor removal step. The eDNA was eluted in 60 μL of Solution C6 and subsequently preserved at −20 oC. Samples were removed from the freezer only for analysis and remained at room temperature for a maximum of 2 h.

Details on F. sultana and T. bryosalmonae eDNA measurement and characterization of the study area are reported in SI Appendix.

Choice of Covariates and Model Settings.

Details on covariates and model settings are reported in SI Appendix. Chosen covariates were checked for multicollinearity. All variance inflation factors for the five considered covariates were below the rule-of-thumb threshold (35) of 10. The three geological covariates were obtained from the vectorized geological map of Switzerland provided by the Swiss Federal Office of Topography (Swisstopo). Covariates were normalized (i.e., linearly transformed into vectors within the range [−1; 1]). SI Appendix, Fig. S3 shows the values of the covariates plotted on the maps of the river network.

The average velocity vij¯ along the path between nodes i and j was calculated as vij¯=∑k∈Pi→jlk/∑k∈Pi→j(lk/vk), where Pi→j indexes the path connecting i to j, while lk and vk are, respectively, the length and the average water velocity of stretch k. Velocities vk were calculated by assuming nearly uniform flow conditions for each stretch and at all times via Manning’s equation. The following further assumptions were made: Water discharges measured by the Swiss Federal Office for the Environment in Zofingen (corresponding to site 3 in Fig. 1A) were used to calculate discharges across all network stretches based on the assumed proportionality between discharge and contributing area; river cross-sections were assumed to be rectangular with width much larger than depth; river widths wk were estimated for all stretches based on aerial images; Manning’s roughness coefficient was taken as n=0.033 m−1/3 s and deemed representative of the flow resistance in the whole river network. Source areas AS,i were assumed equal to river-bed surfaces liwi.

Calibration Algorithm.

As sampled eDNA concentrations for both F. sultana and T. bryosalmonae did not show a clear temporal pattern (Fig. 1C), measured values were considered as realizations of a random variable whose distribution is constant in time. Measured eDNA concentrations Cj,t above the limit of quantification (LOQ) observed at site j and time t were assumed to be log-normally distributed; i.e., lnCj,t=lnCj^+ϵj,t, where ϵj,t∼N(0,σ2), or alternatively Cj,t=Cj^expϵj,t. Samples with low eDNA concentration may fall below the LOQ or induce sampling errors and be interpreted as indicating absence, owing to the small sampling volumes and the lack of replicates for a sample taken at a given site and time. Here, to account for the number of samples where the target eDNA goes undetected, we make the minimalist assumption that the probability φj of not detecting eDNA from a sample collected at site j is a monotonically decreasing function of the eDNA concentration Cj^ predicted by the model: φj=exp−Cj^/C*, where C* is a concentration scale (species specific but constant in space) that requires calibration. The likelihood of nondetection at j is therefore equal to φj, while the likelihood of a positive observation reads (1−φj)ϕlnCj,t/Cj^/σ, where ϕ(⋅) is the standard normal probability density function. Thus, the overall likelihood readsL(C|β,p0,τ,σ,C*)=∏j=1MφjNj1−φjDj∏t=1DjϕlnCj,t/Cj^σ,where C indicates the full set of eDNA concentrations observed at any time and site; M is the number of sampling sites; Nj and Dj are, respectively, the number of null and positive observations at site j; and t=1,…,Dj spans all positive observations at site j. The sampling of the likelihood L was performed by means of a Metropolis-within-Gibbs algorithm (36). For all free parameters, prior distributions were chosen as flat. Further details are reported in SI Appendix.

Acknowledgments

The authors thank Prof. Carlo Gaetan (Università Cà Foscari, Venice) for advice on the statistical calculations. The authors acknowledge the support provided by the Swiss National Science Foundation through the Sinergia project CRSII3_147649: “Temperature driven emergence of Proliferative Kidney Disease in salmonid fish - role of ecology, evolution and immunology for aquatic diseases in riverine landscapes.”

Footnotes

  • ↵1To whom correspondence may be addressed. Email: andrea.rinaldo{at}epfl.ch or luca.carraro{at}epfl.ch.
  • Author contributions: L.C., H.H., J.J., E.B., and A.R. designed research; L.C. and H.H. performed research; L.C., J.J., E.B., and A.R. analyzed data; and L.C., H.H., J.J., E.B., and A.R. wrote the paper.

  • Reviewers: G.F.F., University of Milan; and C.J., University of California, Santa Barbara.

  • The authors declare no conflict of interest.

  • Data deposition: Malacosporean mitochondrial COI sequences have been deposited in the GenBank database, (accession nos. MH986597–MH986599). Matlab codes have been deposited on GitHub (available at https://github.com/lucarraro/edna-species-distribution).

  • This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10.1073/pnas.1813843115/-/DCSupplemental.

  • Copyright © 2018 the Author(s). Published by PNAS.

This open access article is distributed under Creative Commons Attribution-NonCommercial-NoDerivatives License 4.0 (CC BY-NC-ND).

View Abstract

References

  1. ↵
    1. Taberlet P,
    2. Coissac E,
    3. Hajibabaei M,
    4. Rieseberg LH
    (2012) Environmental DNA. Mol Ecol 21:1789–1793.
    OpenUrlCrossRefPubMed
  2. ↵
    1. Thomsen PF,
    2. Willerslev E
    (2015) Environmental DNA - An emerging tool in conservation for monitoring past and present biodiversity. Biol Conserv 183:4–18.
    OpenUrl
  3. ↵
    1. Pace NR
    (1997) A molecular view of microbial diversity and the biosphere. Science 276:734–740.
    OpenUrlAbstract/FREE Full Text
  4. ↵
    1. Bass D,
    2. Stentiford GD,
    3. Littlewood DTJ,
    4. Hartikainen H
    (2015) Diverse applications of environmental DNA methods in parasitology. Trends Parasitol 31:499–513.
    OpenUrlCrossRefPubMed
  5. ↵
    1. Bohmann K, et al.
    (2014) Environmental DNA for wildlife biology and biodiversity monitoring. Trends Ecol Evol 29:358–367.
    OpenUrlCrossRefPubMed
  6. ↵
    1. Kelly RP, et al.
    (2014) Harnessing DNA to improve environmental management. Science 344:1455–1456.
    OpenUrlAbstract/FREE Full Text
  7. ↵
    1. Yoccoz NG
    (2012) The future of environmental DNA in ecology. Mol Ecol 21:2031–2038.
    OpenUrlCrossRefPubMed
  8. ↵
    1. Ficetola GF,
    2. Miaud C,
    3. Pompanon F,
    4. Taberlet P
    (2008) Species detection using environmental DNA from water samples. Biol Lett 4:423–425.
    OpenUrlCrossRefPubMed
  9. ↵
    1. Jerde CL,
    2. Mahon AR,
    3. Chadderton WL,
    4. Lodge DM
    (2011) “Sight-unseen” detection of rare aquatic species using environmental DNA. Conserv Lett 4:150–157.
    OpenUrl
  10. ↵
    1. Dejean T, et al.
    (2012) Improved detection of an alien invasive species through environmental DNA barcoding: The example of the American bullfrog Lithobates catesbeianus. J Appl Ecol 49:953–959.
    OpenUrlCrossRef
  11. ↵
    1. Mächler E,
    2. Deiner K,
    3. Steinmann P,
    4. Altermatt F
    (2014) Utility of environmental DNA for monitoring rare and indicator macroinvertebrate species. Freshwater Sci 33:1174–1183.
    OpenUrl
  12. ↵
    1. Takahara T,
    2. Minamoto T,
    3. Yamanaka H,
    4. Doi H,
    5. Kawabata Z
    (2012) Estimation of fish biomass using environmental DNA. PLoS One 7:e35868.
    OpenUrlCrossRefPubMed
  13. ↵
    1. Huver JR,
    2. Koprivnikar J,
    3. Johnson PTJ,
    4. Whyard S
    (2015) Development and application of an eDNA method to detect and quantify a pathogenic parasite in aquatic ecosystems. Ecol Appl 25:991–1002.
    OpenUrl
  14. ↵
    1. Olds B, et al.
    (2016) Estimating species richness using environmental DNA. Ecol Evol 6:4214–4226.
    OpenUrl
  15. ↵
    1. Barnes M, et al.
    (2014) Environmental conditions influence eDNA persistence in aquatic systems. Environ Sci Technol 48:1819–1827.
    OpenUrlCrossRefPubMed
  16. ↵
    1. Lance RF, et al.
    (2017) Experimental observations on the decay of environmental DNA from bighead and silver carps. Manage Biol Invasions 8:343–359.
    OpenUrl
  17. ↵
    1. Jerde CL, et al.
    (2016) Influence of stream bottom substrate on retention and transport of vertebrate environmental DNA. Environ Sci Technol 50:8770–8779.
    OpenUrl
  18. ↵
    1. Shogren AJ, et al.
    (2017) Controls on eDNA movement in streams: Transport, retention, and resuspension. Sci Rep 7:5065.
    OpenUrl
  19. ↵
    1. Deiner K,
    2. Altermatt F
    (2014) Transport distance of invertebrate environmental DNA in a natural river. PLoS One 9:e88786.
    OpenUrlCrossRefPubMed
  20. ↵
    1. Wilcox TM,
    2. McKelvey KS,
    3. Young MK,
    4. Lowe WH,
    5. Schwartz MK
    (2015) Environmental DNA particle size distribution from brook trout (Salvelinus fontinalis). Conserv Genet Resour 7:639–641.
    OpenUrlCrossRef
  21. ↵
    1. Klymus K,
    2. Richter C,
    3. Chapman D,
    4. Paukert C
    (2015) Quantification of eDNA shedding rates from invasive bighead carp Hypophthalmichthys nobilis and silver carp Hypophthalmichthys molitrix. Biol Conserv 183:77–84.
    OpenUrl
  22. ↵
    1. Bylemans J, et al.
    (2017) An environmental DNA-based method for monitoring spawning activity: A case study, using the endangered Macquarie perch (Macquaria australasica). Methods Ecol Evol 8:646–655.
    OpenUrl
  23. ↵
    1. Carraro L, et al.
    (2017) Integrated field, laboratory, and theoretical study of PKD spread in a Swiss prealpine river. Proc Natl Acad Sci USA 114:11992–11997.
    OpenUrlAbstract/FREE Full Text
  24. ↵
    1. Sansom BJ,
    2. Sassoubre LM
    (2017) Environmental DNA (eDNA) shedding and decay rates to model freshwater mussel eDNA transport in a river. Environ Sci Technol 51:14244–14253.
    OpenUrl
  25. ↵
    1. Pfister L, et al.
    (2009) The rivers are alive: On the potential for diatoms as a tracer of water source and hydrological connectivity. Hydrol Process 23:2841–2845.
    OpenUrl
  26. ↵
    1. Rodriguez-Iturbe I,
    2. Rinaldo A
    (2001) Fractal River Basins. Chance and Self-Organization (Cambridge Univ Press, New York).
  27. ↵
    1. Pilgrim DH
    (1977) Isochrones of travel time and distribution of flood storage from a tracer study on a small watershed. Water Resour Res 13:587–595.
    OpenUrl
  28. ↵
    1. Rinaldo A,
    2. Marani A,
    3. Rigon R
    (1991) Geomorphological dispersion. Water Resour Res 27:513–525.
    OpenUrl
  29. ↵
    1. Burkhardt-Holm P, et al.
    (2005) Where have all the fish gone? Environ Sci Technol 39:441A–447A.
    OpenUrlPubMed
  30. ↵
    1. Okamura B,
    2. Hartikainen H,
    3. Schmidt-Posthaus H,
    4. Wahli T
    (2011) Life cycle complexity, environmental change and the emerging status of salmonid proliferative kidney disease. Freshwater Biol 56:735–753.
    OpenUrlCrossRef
  31. ↵
    1. James F
    (2006) Statistical Methods in Experimental Physics (World Scientific, Singapore).
  32. ↵
    1. Carraro L,
    2. Mari L,
    3. Gatto M,
    4. Rinaldo A,
    5. Bertuzzo E
    (2018) Spread of proliferative kidney disease in fish along stream networks: A spatial metacommunity framework. Freshwater Biol 63:114–127.
    OpenUrl
  33. ↵
    1. Tops S,
    2. Okamura B
    (2003) Infection of bryozoans by Tetracapsuloides bryosalmonae at sites endemic for salmonid proliferative kidney disease. Dis Aquat Org 57:221–226.
    OpenUrlPubMed
  34. ↵
    1. Heggenes J,
    2. Baglinière JL,
    3. Cunjak RA
    (1999) Spatial niche variability for young Atlantic salmon (Salmo salar) and brown trout (S. trutta) in heterogeneous streams. Ecol Freshwater Fish 8:1–21.
    OpenUrl
  35. ↵
    1. Hastie T,
    2. Tibshirani R,
    3. Friedman J
    (2001) The Elements of Statistical Learning, Springer Series in Statistics(Springer, New York).
  36. ↵
    1. Roberts GO,
    2. Rosenthal JS
    (2009) Examples of adaptive MCMC. J Comput Graph Stat 18:349–367.
    OpenUrlCrossRef
PreviousNext
Back to top
Article Alerts
Email Article

Thank you for your interest in spreading the word on PNAS.

NOTE: We only request your email address so that the person you are recommending the page to knows that you wanted them to see it, and that it is not junk mail. We do not capture any email address.

Enter multiple addresses on separate lines or separate them with commas.
Estimating species distribution and abundance in river networks using environmental DNA
(Your Name) has sent you a message from PNAS
(Your Name) thought you would like to see the PNAS web site.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Citation Tools
Estimating species distribution and abundance in river networks using environmental DNA
Luca Carraro, Hanna Hartikainen, Jukka Jokela, Enrico Bertuzzo, Andrea Rinaldo
Proceedings of the National Academy of Sciences Nov 2018, 115 (46) 11724-11729; DOI: 10.1073/pnas.1813843115

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Request Permissions
Share
Estimating species distribution and abundance in river networks using environmental DNA
Luca Carraro, Hanna Hartikainen, Jukka Jokela, Enrico Bertuzzo, Andrea Rinaldo
Proceedings of the National Academy of Sciences Nov 2018, 115 (46) 11724-11729; DOI: 10.1073/pnas.1813843115
Digg logo Reddit logo Twitter logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Mendeley logo Mendeley
Proceedings of the National Academy of Sciences: 115 (46)
Table of Contents

Submit

Sign up for Article Alerts

Article Classifications

  • Physical Sciences
  • Environmental Sciences
  • Biological Sciences
  • Ecology

Jump to section

  • Article
    • Abstract
    • Materials and Methods
    • Acknowledgments
    • Footnotes
    • References
  • Figures & SI
  • Info & Metrics
  • PDF

You May Also be Interested in

Abstract depiction of a guitar and musical note
Science & Culture: At the nexus of music and medicine, some see disease treatments
Although the evidence is still limited, a growing body of research suggests music may have beneficial effects for diseases such as Parkinson’s.
Image credit: Shutterstock/agsandrew.
Scientist looking at an electronic tablet
Opinion: Standardizing gene product nomenclature—a call to action
Biomedical communities and journals need to standardize nomenclature of gene products to enhance accuracy in scientific and public communication.
Image credit: Shutterstock/greenbutterfly.
One red and one yellow modeled protein structures
Journal Club: Study reveals evolutionary origins of fold-switching protein
Shapeshifting designs could have wide-ranging pharmaceutical and biomedical applications in coming years.
Image credit: Acacia Dishman/Medical College of Wisconsin.
White and blue bird
Hazards of ozone pollution to birds
Amanda Rodewald, Ivan Rudik, and Catherine Kling talk about the hazards of ozone pollution to birds.
Listen
Past PodcastsSubscribe
Goats standing in a pin
Transplantation of sperm-producing stem cells
CRISPR-Cas9 gene editing can improve the effectiveness of spermatogonial stem cell transplantation in mice and livestock, a study finds.
Image credit: Jon M. Oatley.

Similar Articles

Site Logo
Powered by HighWire
  • Submit Manuscript
  • Twitter
  • Facebook
  • RSS Feeds
  • Email Alerts

Articles

  • Current Issue
  • Latest Articles
  • Archive

PNAS Portals

  • Anthropology
  • Chemistry
  • Classics
  • Front Matter
  • Physics
  • Sustainability Science
  • Teaching Resources

Information

  • Authors
  • Editorial Board
  • Reviewers
  • Librarians
  • Press
  • Site Map
  • PNAS Updates

Feedback    Privacy/Legal

Copyright © 2021 National Academy of Sciences. Online ISSN 1091-6490