Quantification of Environmental DNA (eDNA) Shedding and Decay

Aug 31, 2016 - We thank Ian Rowbotham, Nick Mendoza, and Alex Norton for their help setting up the experiments and caring for the fish, Oliver Fringer...
0 downloads 0 Views 599KB Size
Subscriber access provided by Northern Illinois University

Article

Quantification of environmental DNA (eDNA) shedding and decay rates for three marine fish Lauren M Sassoubre, Kevan M Yamahara, Luke D. Gardner, Barbara A. Block, and Alexandria B. Boehm Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.6b03114 • Publication Date (Web): 31 Aug 2016 Downloaded from http://pubs.acs.org on August 31, 2016

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

Environmental Science & Technology is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 23

Environmental Science & Technology

ACS Paragon Plus Environment

Environmental Science & Technology

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29

Page 2 of 23

Quantification of environmental DNA (eDNA) shedding and decay rates for three marine fish

Lauren M. Sassoubrea†*, Kevan M. Yamaharab, Luke D. Gardnerc, Barbara A. Blockc, Alexandria B. Boehma a Department of Civil and Environmental Engineering, Stanford University, Stanford, CA 94305 USA † Present Address: Department of Civil, Structural, and Environmental Engineering, The State University of New York at Buffalo, Buffalo, New York 14260 USA b Monterey Bay Aquarium Research Institute, Moss Landing, CA 95039 USA c Hopkins Marine Station, Department of Biology, Stanford University, Pacific Grove, CA 93950 USA KEYWORDS: environmental DNA, eDNA, shedding, decay, DNA persistence, marine waters, fish, qPCR * Corresponding author: Lauren Sassoubre University at Buffalo Civil, Structural and Environmental Engineering Tel: (716) 645-1810 Email: [email protected]

ACS Paragon Plus Environment

1

Page 3 of 23

30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48

Environmental Science & Technology

ABSTRACT Analysis of environmental DNA (eDNA) to identify macroorganisms and biodiversity has the potential to significantly augment spatial and temporal biological monitoring in aquatic ecosystems. Current monitoring methods relying on the physical identification of organisms can be time consuming, expensive and invasive. Measuring eDNA shed from organisms provides detailed information on the presence and abundance of communities of organisms. However, little is known about eDNA shedding and decay in aquatic environments. In the present study, we designed novel Taqman® qPCR assays for three ecologically and economically important marine fish – Engraulis mordax (Northern Anchovy), Sardinops sagax (Pacific Sardine), and Scomber japonicas (Pacific Chub Mackerel). We subsequently measured fish eDNA shedding and decay rates in seawater mesocosms. eDNA shedding rates ranged from 165 to 3368 pg DNA per hour per gram of biomass. First order decay rate constants ranged from 0.055 to 0.070 per hour. We also examined the size fractionation of eDNA and concluded eDNA is both intra- and extracellular. Finally, we derived a simple mass-balance model to estimate fish abundance from eDNA concentration. The mesocosm derived shedding and decay rates inform the interpretation of eDNA concentrations measured in environmental samples and future use of eDNA as a monitoring tool.

ACS Paragon Plus Environment

2

Environmental Science & Technology

49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94

Page 4 of 23

INTRODUCTION Marine biodiversity is threatened by climate change, pollution, human population growth, overfishing, invasive species and habitat loss/alteration 1, 2. Monitoring and quantifying the impacts of these threats, however, is challenging due to a lack of biological data describing the abundance, spatial and temporal distribution of organisms in the marine realm 1. A number of techniques are currently employed to monitor marine organisms including underwater visual counts by divers or video, trawls or netting, tagging, and traditional fishing 3. These techniques rely on the physical identification and counting of organisms. They can therefore be: time-consuming; expensive when taxonomic experts and shiptime are required; invasive or destructive when netting; trapping or electrofishing methods are used; dependent on chance encounters with potentially rare organisms and limited to waters accessible by ships or divers. To address limitations associated with traditional biological monitoring, researchers have proposed the use of molecular techniques to analyze DNA extracted from environmental samples (termed environmental DNA or eDNA) to census the presence of fish and other aquatic organisms 4-11. Organisms shed DNA into their environments in the form of sloughed tissue or cells, waste products, gametes, saliva or other secretions. eDNA can be isolated from water samples and used in qPCR assays targeting species-specific DNA sequences or in next generation sequencing (NGS) that identifies communities of organisms 7. The concentration of eDNA in a water sample is controlled by a number of environmental processes. As shown in Figure 1, the eDNA concentration depends on its sources (e.g., shedding from organisms) and sinks (e.g., decay) in a given parcel of water and the advection and dispersion of that water parcel in the aquatic system. Thus, the first step to infer organism presence and population densities from eDNA concentrations is obtaining information on eDNA shedding and decay rates. This information can then inform the use of eDNA for biological monitoring, biodiversity assessments and ultimately policy decisions. The eDNA shedding rates of aquatic macroorganisms depend on a number of factors: (1) the type of organism or species 5, 12, (2) organism size13, (3) number of organisms 5, 13-15, (4) life stage 5, 12, 16, (5) skin/scale properties (6) the stress an organism is under 17 and (7) water temperature 15. Few studies have directly measured how much eDNA an organism sheds into a water body over time 17-19 and these studies focused on freshwater amphibian larvae, salamanders and carp. A microcosm study with amphibian larvae found that eDNA concentrations increased over time and were higher in microcosms with more larval species, however, shedding rates differed between larval species 5. eDNA decay is influenced by (1) DNA characteristics (e.g.,, whether it is extracellular or cellular) 20, (2) abiotic factors (e.g.,, sunlight, temperature, pH, salinity, flow rate or residence time) 21-23, and (3) biotic factors (e.g.,, extracellular enzymes and microorganisms) 22, 24. The majority of research on eDNA decay has focused on freshwater fish and amphibians. These studies reported eDNA concentrations falling below detection limits in 4.2 to 54 days 5, 14, 17, 21-23, 25, 26. Thus far, eDNA reportedly decays faster in marine waters than freshwater 27, 28, however, very few eDNA decay studies have been carried out in marine waters. Studies conducted with marine organisms

ACS Paragon Plus Environment

3

Page 5 of 23

95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140

Environmental Science & Technology

in experimental tanks report eDNA decaying below detection limits in 0.9 days for European flounder (Platichthys flesus) 27, 6.7 days for three-spined stickleback (Gasterosteus aculeatus) 27, and a decay rate constant of 0.104 ± 0.047 per hour for bluegill sunfish (Lepomis macrochirus) 16. The implementation of eDNA in biological monitoring and management practices will likely require eDNA abundance to be representative of organismal abundance. A few studies have been published revealing a positive relationship between biomass or species abundance and eDNA concentration or sequence abundance 5, 13, 15, 19, 24, 29-31. However, there are limited studies in the marine environment and none that present an empirical equation for estimating fish abundance from eDNA concentrations that incorporate the fate and transport of eDNA in the marine environment. The goal of the present study is to determine eDNA shedding and decay rates for three ecologically and economically important marine fish in the temperate waters of the eastern Pacific Ocean 32-34. We developed a novel set of sensitive and specific Taqman qPCR assays for Northern Anchovy (Engraulis mordax), Pacific Sardine (Sardinops sagax), and Pacific Chub Mackerel (Scomber japonicas). We then conducted seawater mesocosm experiments to quantify fish eDNA shedding and decay rates. We also investigated eDNA size fractionation to gain insight into the source of eDNA (e.g., intracellular or extracellular). Finally, we utilized the eDNA shedding rates and decay rate constants for Northern Anchovy in a simplified mass balance model that estimates fish abundance from eDNA concentrations. The results of this research are critical to inform the use of eDNA as a biological monitoring tool for fisheries management and marine conservation efforts.

MATERIALS AND METHODS Genus-specific assay design and qPCR optimization. Genus-specific primers and Taqman® probes were designed to amplify Northern Anchovy (Engraulis mordax), Pacific Sardine (Sardinops sagax), and Pacific Chub Mackerel (Scomber japonicas). Primers and probes were designed using PrimerBlast 35 (see SI). The specificity of potential primers/probes sequences was assessed in silco using PrimerBlast and using an alignment of closely related sequences in NCBI. Primer/probe sequences showing specificity in silico were obtained from Integrated DNA Technologies (IDT, Coralville, Iowa) and tested for cross-reactivity using genomic DNA (1-4 ng) extracted from the tissue samples of organisms found in the same marine habitats as anchovies, sardines and mackerel. Genomic DNA was extracted using the Qiagen DNeasy Blood and Tissue extraction kit according to manufacturers instructions (Qiagen, Valencia, CA). Genomic DNA from the following fish and one cephalopod was tested against each primer set: sardines (Sardinops), anchovies (Engraulis), mackerel (Scomber), herring (Clupea), rockfish (Sebastes), tuna (Thunnus) and squid (Doryteuthis). Primer/probe sets were considered specific if no amplification (within 40 cycles) was observed for any of the non-target genomic DNA. The final primer and probe sequences and concentrations, DNA targets, and primer annealing temperatures are shown in Table 1.

ACS Paragon Plus Environment

4

Environmental Science & Technology

141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186

Page 6 of 23

qPCR standards were constructed using genomic DNA (gDNA) extracted from target fish tissue samples and quantified using a QUBIT fluorometer 2.0 (Life Technologies, Grand Island, NY). Standard curves consisted of 1:10 dilutions of gDNA from 1 ng/µl to 0.01 pg/µl and were run in triplicate alongside samples in each 96 well plate. Standard curves were pooled across plates to calculate concentrations of unknown samples 36. All unknowns were amplified in triplicate 20 µl qPCR reactions and each qPCR plate contained triplicate no template controls (NTCs). Reaction chemistry consisted of 1X Taqman Universal Mastermix II, 0.2 mg/ml bovine serum album (BSA), a Taqman probe and forward and reverse primers in optimized concentrations. Inhibition was assessed using serial dilutions (see SI) 37, 38. Based on inhibition tests, samples from the anchovy only experiment were diluted 1:5, samples from sardine only experiment were diluted 1:10, and samples from the mackerel plus sardine experiment were diluted 1:10 for the sardine assay. Cycle quantification (Ct) thresholds were set at 0.02 for anchovy, 0.01 for sardine and 0.01 for mackerel. Experimental design. Experiments were conducted at the Tuna Research Conservation Center (TRCC) at Hopkins Marine Laboratory of Stanford University, in Pacific Grove, California. We performed four seawater mesocosm experiments where we examined the shedding and decay rates of (1) anchovies only, (2) sardines only, (3) mackerel only and (4) mackerel plus sardines. The mesocosm consisted of an indoor, circular fiberglass tank with a capacity of ~5200 liters and a water depth of 1 meter. Seawater from Monterey Bay, CA was used for all experiments. For experiments examining shedding and decay of anchovy and mackerel eDNA, the tank was filled with TRCC seawater (a mixture of seawater from Monterey Bay and seawater supplied from the Monterey Bay Aquarium which is filtered through sand filters and aerated) (see description in SI). For the sardine only and the mackerel plus sardine experiment, the tank was filled with seawater from Monterey Bay rather than TRCC seawater. TRCC seawater was not used for the experiments with sardines to minimize the potential background sardine eDNA signal that may result due to the fact the TRCC tuna are fed sardines. Once the experimental tank was filled, seawater was circulated in a closed loop system to aerate the water throughout the experiments. For the mackerel and sardine experiments, additional aeration was supplied by bubbling oxygen into the tank due to the high biomass content in the tank. Tank water temperatures ranged from 22. 0 ±1°C during the anchovy experiment (due to a malfunction of the water chilling device) and 18.7±1˚C for the sardine only, mackerel only and mackerel plus sardine experiments. The temperatures were within a range suitable for the fish and are similar to temperatures within Monterey Bay, CA (10-15°C)39. Fish were added to the tank and fasted during the experiments (~3 days) to control for potential eDNA input resulting from feeding. The experiment with only anchovies contained 43 anchovies with an average weight of 20 ± 2 grams per fish (determined by displacement). The experiment with only sardines contained 27 sardines with an average weight of 88 ± 11 grams per fish. The experiment with only mackerel contained 20 mackerel with an average weight of 424 ± 4 grams per fish. The experiment with sardines and mackerel contained 30 sardines and 15 mackerel with average weights of

ACS Paragon Plus Environment

5

Page 7 of 23

187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232

Environmental Science & Technology

110 ± 29 and 449 ± 90 grams per fish, respectively. The number of fish used in the experiments was based on the number of fish available and the capacity of the tank. Sampling Procedure. All water samples were collected in two 500 ml acid-washed (10% HCl) polypropylene bottles at each time point as biological replicates. Two samples were collected from the tank before fish were added (t=0). The fish were then added to the tank. Water was sampled frequently while the fish were in the tank (3-8 times per day depending on the experiment and results from previous experiments). On the third day of the experiment, the fish were removed from the tanks and the water was sampled 2-3 times per day for 3-4 more days (Table S1). Each day, a filtration blank was created with 50 ml molecular grade water (Sigma Aldrich, St. Louis, MO) to test for contamination during sample filtration. New gloves were worn during each sample collection. Samples were placed on ice and immediately processed at the laboratory. 250-500 ml was filtered in duplicate through 0.2 µm-pore size 47 mm diameter track-etched polycarbonate filters (Nucleopore Track-Etch Membrane, Whatman, GE Healthcare Bio-Sciences, Pittsburgh, PA) in disposable filter funnels (ThermoScientific, Waltham, MA). Filters were stored at -20˚ C until DNA extraction (approximately 1-4 weeks after the experiment). Size fractionation of environmental DNA. The particle size distributions of mackerel and sardine eDNA was determined during the mackerel plus sardine experiment. Fortyone hours after sardine and mackerel were introduced into the tank, triplicate 20 ml and 1 liter seawater samples were collected directly from the tank. Each 1 liter sample was sequentially filtered through a 10 µm, followed by a 1 µm, followed by a 0.2 µm poresize, 47 mm diameter polycarbonate filter (Nucleopore Track-Etch Membrane, Whatman, GE Healthcare Bio-Sciences, Pittsburgh, PA). The full 1 liter was sequentially filtered through the 10 and 1 µm filters. Due to filter clogging, 250 ml of the 1 µm pore size filtrate was filtered through the 0.2 µm filter. All filters were stored at -20° C until extracted. The triplicate 20 ml unfiltered tank seawater samples were prepared for DNA precipitation. The filtrate from the 0.2 µm pore size filter was also prepared for DNA precipitation. Briefly, 20 ml samples were combined with 20 ml of molecular grade isopropanol (Fisher Scientific, Fair Lawn, NJ) and 1.5 ml sodium acetate (3M, pH 5.2, CALBIOCHEM, EMD Biosciences, La Jolla, CA) and saved at -20 °C until extraction 40. A DNA precipitation blank made with 20 ml of molecular grade water was also prepared and processed alongside the samples. Filter and precipitate DNA extraction methods are described below. DNA precipitation and extraction. DNA was precipitated from unfiltered seawater from the tank, filtrate from the 0.2 µm pore size filter and the DNA precipitation blank. The filtrate was centrifuged at 6000 x g for 45 minutes at room temperature 39. The supernatant was carefully removed from the pellet. DNA was extracted from pellets and filters using the Qiagen DNeasy Blood and Tissue extraction kit (Qiagen, Valencia, CA) with modifications described in SI. A DNA extraction blank was extracted with each set of samples to test for contamination in the reagents. DNA was eluted in 100 µl warmed Buffer AE (10 mM Tris-Cl and 0.5 mM EDTA, Qiagen, Valencia, CA) for the filters and 50 µl warmed Buffer AE for the DNA precipitations.

ACS Paragon Plus Environment

6

Environmental Science & Technology

Page 8 of 23

233 234 235 236 237 238 239 240 241 242 243 244 245 246

Inhibition. Samples from each experiment were investigated for inhibition using serial dilutions 37, 38. Briefly, samples were diluted in two series, five-fold (1:5, 1:25) and tenfold (1:10, 1:100) to find the minimal dilution needed to dilute out inhibitors but not lower the target concentration below the detection limit. A five-fold dilution is expected to result in a Ct change of 2.32 cycles (log2(5)=2.32) and a ten-fold dilution is expected to result in a Ct change of 3.32 cycles (log2(10)=3.32), assuming 100% efficiency. We considered the sample inhibited if the difference between the sample and the dilution was 0.5 cycles less than the expected change. We also added bovine serum album (BSA) to each QPCR to help with inhibition.

247

V

248 249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270 271 272 273 274 275 276 277

Data analysis eDNA shedding. The tank with fish was modeled as a completely mixed batch reactor for each experiment:

dC = S − kCV dt

(Eqn. 1)

where V is the volume of the tank with units of liters, C is the concentration of eDNA with units of pg/mL, t is the time since the start of the experiment with units of hours, S is the eDNA shedding rate with units of pg/hr, and k is the first order decay rate constant with units of per hour (see decay rate constant calculation in eDNA decay section below). Eqn. 1 assumes that the tank is well mixed and that decay is first order. Given that the fish were constantly swimming in the tank, we believe this is a reasonable assumption. Steady-state was reached after 17-25 hours depending on the experiment (Table S1) and lasted until the fish were removed from the system. At steady state, dC/dt = 0 and S = kCV. The error associated with the shedding rate was determined by propagating errors associated with k, C, and V. A z-test was used to compare shedding rates between fish species and experiments and a significant difference was determined by p