Subscriber access provided by McMaster University Library
Article
Rapid molecular detection of invasive species in ballast and harbor water by integrating environmental DNA and Light Transmission Spectroscopy Scott Patrick Egan, Erin Grey, Brett Olds, Jeff Feder, Steve Ruggiero, Carol E. Tanner, and David Lodge Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/es5058659 • Publication Date (Web): 16 Feb 2015 Downloaded from http://pubs.acs.org on February 18, 2015
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 38
Environmental Science & Technology
1
TITLE: Rapid molecular detection of invasive species in ballast and harbor water by integrating
2
environmental DNA and Light Transmission Spectroscopy
3 4
Scott P. Egan1,2,3,*, Erin Grey3,4, Brett Olds3, Jeffery L. Feder2,3, Steven T. Ruggiero3,5, Carol E.
5
Tanner3,5, & David M. Lodge2,3
6 7
1
Department of BioSciences, Rice University, Houston, TX 77005
8
2
Department of Biological Sciences, University of Notre Dame, Notre Dame, Indiana 46556
9
3
Environmental Change Initiative, University of Notre Dame, Notre Dame, Indiana 46556
10
4
College of Science, Governors State University, University Park, IL 60484
11
5
Department of Physics, University of Notre Dame, Notre Dame, Indiana 46556
12 13
*Corresponding author: Scott P. Egan, Department of BioSciences, Rice University, Houston,
14
TX 77005;
[email protected], phone: 615-618-6601, fax: 713-348-5232
15 16
Running title: Aquatic invasive species detection using eDNA
ACS Paragon Plus Environment
Environmental Science & Technology
Page 2 of 38
17
ABSTRACT
18
Invasive species introduced via the ballast water of commercial ships cause enormous
19
environmental and economic damage worldwide. Accurate monitoring for these often
20
microscopic and morphologically indistinguishable species is challenging, but critical for
21
mitigating damages. We apply eDNA sampling, which involves the filtering and subsequent
22
DNA extraction of microscopic bits of tissue suspended in water, to ballast and harbor water
23
sampled during a commercial ship’s 1,400km voyage through the North American Great Lakes.
24
Using a lab-based gel electrophoresis assay and a rapid, field-ready Light Transmission
25
Spectroscopy (LTS) assay, we test for the presence of two invasive species: quagga (Dreissena
26
bugensis) and zebra (D. polymorpha) mussels. Furthermore, we spiked a set of uninfested ballast
27
and harbor samples with zebra mussel tissue to further test each assay’s detection capabilities. In
28
unmanipulated samples, zebra mussel was not detected, while quagga mussel was detected in all
29
samples at a rate of 85% for the gel assay and 100% for the LTS assay. In the spiked
30
experimental samples, both assays detected zebra mussel in 94% of spiked samples and 0% of
31
negative controls. Overall, these results demonstrate that eDNA sampling is effective for
32
monitoring ballast-mediated invasions and that LTS has the potential for rapid, field-based
33
detection.
34
Key Words: conservation, environmental DNA, invasive species, light transmission spectroscopy
ACS Paragon Plus Environment
2
Page 3 of 38
35
Environmental Science & Technology
INTRODUCTION
36 37
Invasive species have had a globally damaging effect on freshwater and marine ecosystems,
38
adversely impacting biodiversity and human commerce.1 In the United States alone, invasive
39
species cause environmental damage and economic losses in excess of US $137 billion
40
annually2,3 and are now recognized as one of the leading global threats to native biodiversity and
41
ecosystem functions.4,5 Within the North American Great Lakes, there are over 180 aquatic
42
invasive species that have significantly altered the ecosystem and are estimated to cost up to
43
$800 million annually.6-10 The majority of these aquatic invasive species (>70%) entered through
44
transoceanic shipping routes, including the ballast water discharge of incoming ships.11-14 A
45
general lack of sensitive technologies for early detection can allow harmful invasions to establish
46
to the point where it is too late to intercede for an increasing number of highly damaging
47
species.15
48 49
The ecological and economic threats of aquatic invasions therefore require improved monitoring
50
of the various anthropogenic pathways, such as ships’ ballast water. Unfortunately, ballast water
51
treatment technologies are still mainly in the research and development phase.16,17 Consequently,
52
there is an urgent need for rapid detection methods to identify harmful invasive species in ships’
53
ballast prior to port entry to provide adequate time for initiating appropriate control actions. Even
54
if ballast water treatment systems become widely adopted, the need will still exist for rapid
55
inspection and enforcement tools, as a monitoring firewall against harmful invasives evading the
56
first line of control.
57
ACS Paragon Plus Environment
3
Environmental Science & Technology
Page 4 of 38
58
Traditional microscope-based methods have been used to examine ballast samples for target
59
organisms, but are expensive and slow.18 By the time results from microscopic examination are
60
available, ships have long since reached their destinations, making it too late to intervene.14 In
61
addition, the most often encountered life stage in ballast tanks, larvae, are often not identifiable
62
to the species level using morphological keys.19 Accordingly, a rapid and sensitive on-board
63
detection platform that can identify target species in large volume ballast samples requiring
64
minimal user expertise would greatly facilitate inspection (voluntary or regulatory), enforcement,
65
and control efforts.
66 67
Species detection also represents a fundamental tool for effective conservation and management.
68
Without accurate, up-to-date knowledge of where rare endangered species are and where
69
invasive propagules reside, cost effective management is impossible. Environmental DNA
70
(eDNA) sampling has become increasingly prevalent in such conservation and management
71
efforts20 and has recently experienced rapid development in their application to monitoring
72
biodiversity.21,22 In a management context, there are at least two major advantages of eDNA
73
methods over traditional methods: (1) increased sensitivity22,23 (i.e., increased probability of
74
detecting a species if it is present) and (2) potentially lower costs, especially as eDNA methods
75
increase in portability and automation.24,25
76 77
The overarching goal of our work is to enable accurate, sensitive, timely, and cost effective
78
detection of target species. In the present study, we merge two new, but proven technologies for
79
the purpose of quickly and accurately detecting target species in ship’s ballast and surrounding
80
harbor waters. First, we use DNA filtering methods to isolate and purify field samples of
ACS Paragon Plus Environment
4
Page 5 of 38
Environmental Science & Technology
81
environmental DNA (eDNA) from a ship’s ballast water and the surrounding harbor area. In this
82
regard, analysis of aquatic samples with traditional methods, such as microscopy and
83
morphology, often poses insurmountable logistical challenges for detecting rare and
84
taxonomically-challenging organisms. However, species surveillance using eDNA exploits the
85
fact that the aqueous environment often keeps microscopic bits of tissue from target invasive
86
species in suspension in the water column, making it easier to sample DNA from even rare
87
organisms that are present but invisible to traditional tools. We define eDNA here as all types of
88
tissue samples collected by filtering water from an aquatic environment. These samples can
89
include, but are not limited to: sloughed tissues or cells, larvae or adult organisms that are
90
microscopic, milt, eggs, extracellular DNA from degraded tissues, scales, and invertebrate
91
exoskeletons. The eDNA method has been successfully developed and applied to possible Asian
92
carp invasions into the Great Lakes23,26 and for other species surveillance questions.22,27-31 The
93
eDNA approach to species surveillance is rapidly developing into a burgeoning field of research
94
in molecular ecology and application in conservation and invasive species management.21,32
95 96
Second, we combine the efficacy of eDNA sampling with a new field-based tool for target DNA
97
detection, Light Transmission Spectroscopy (LTS).33,34 Typically, eDNA detection methods have
98
relied on PCR or qPCR to amplify target DNA present in a sample, where amplification is
99
equivalent to detection.35 However, the LTS platform offers the possibility of field-based species
100
detection because it is portable, rapid, and requires low electrical power to operate.24 LTS
101
determines the size, shape, and number of nanoparticles in suspension by measuring wavelength-
102
dependent light transmittance through a sample containing nanoparticles in suspension. 33,34 LTS
103
is applicable for species detection by measuring nanoparticles that specifically bind to target
ACS Paragon Plus Environment
5
Environmental Science & Technology
Page 6 of 38
104
species DNA and therefore grow in size.24,25,34 Specifically, carboxylated nanoparticles of a
105
known size (typically around 200nm in diameter) are bound with species-specific tags
106
representing ~18-26 base pair oligonucleotides that bind to a unique section of target DNA.
107
When the tagged nanoparticles are in solution with the DNA of a target species, the tagged
108
nanoparticles grow in diameter because the tags bind to PCR amplified DNA fragments. LTS
109
can detect the size shift in particle size within seconds (Supplemental Figure 1). The LTS device
110
is portable and can run off a car’s battery. Thus, LTS technology has promise for uses where
111
speed and portability are paramount (e.g. real-time field surveys of ballast water, law
112
enforcement inspections).
113 114
In the present study, we test the effectiveness of merging eDNA sampling with LTS for invasive
115
species detection using two test species that are currently found in the ballast waters of ships
116
from foreign ports entering the Great Lakes: (1.) Dreissena polymorpha (zebra mussel) and (2.)
117
Dreissena bugensis (quagga mussel). Both invasive species are known to be in the Great Lakes
118
with viable populations. We report that combining eDNA sampling with LTS can rapidly and
119
accurately detect target invasive species in two types of environmental samples where the threat
120
of aquatic invasion is the highest, within the ballast water of a large commercial shipping vessel
121
and in the surrounding harbor where these ship-borne aquatic invasions usually begin. Moreover,
122
these data support the continued development of the LTS portable real-time DNA detection
123
system that could be used for onboard monitoring during a ship’s passage prior to ballast
124
discharge to dramatically reduce the impacts of future invasions.12,14
125 126
MATERIALS and METHODS
ACS Paragon Plus Environment
6
Page 7 of 38
Environmental Science & Technology
127 128
Ballast and harbor sampling and eDNA filtration
129 130
Ballast water was collected directly from the intake and discharge of a 1,000-ft. shipping vessel,
131
the M/V Indiana Harbor, during its five-day, 1,400-km voyage from the Indiana Harbor (41° 39'
132
57.1788" N, 87° 26' 43.9038" W) near Gary, Indiana, on Lake Michigan, to the Duluth-Superior
133
Harbor near Duluth, MN (46° 45' 14.3136" N, 92° 5' 59.1426" W), on Lake Superior during mid-
134
October 2012 (Supplemental Figure 2). In parallel, we also collected water from the ballast
135
intake source, Indiana Harbor, IN and the destination port receiving the ballast water discharge,
136
the Duluth-Superior Harbor, MN. All water samples were collected in a two-liter sterile Nalgene
137
bottle as follows: Ballast intake samples (n = 20) were collected during ballast intake from a
138
direct tap connected to the pumping system on October 17th, 2012. Source harbor samples (n = 5)
139
were collected from the surface of the water directly next to the ship during ballast intake.
140
Receiving harbor surface water samples (n = 3) were collected five days later (October 22nd,
141
2012) in the same manner just prior to ballast discharge. Ballast discharge samples (n = 13) were
142
collected from the direct line to the pumping system from three separate tanks (designated tank
143
#3 = 3 samples, tank #678 = 3 samples, and tank #5 = 5 samples). Each ballast tank is
144
independently sealed from the other and thus, can potentially harbor different contaminants
145
(biotic or abiotic) in each.
146 147
For each of the ballast and harbor samples, a 250mL aliquot was filtered through a 1.2-micron
148
polycarbonate membrane filter (GE Water & Process Technologies) to screen for the aquatic
149
invasive species, quagga mussel (Dreissena bugensis), expected to be present and active in the
ACS Paragon Plus Environment
7
Environmental Science & Technology
150
Great Lakes region at the time of sampling36-38 and zebra mussel (Dreissena polymorpha), also
151
present but not expected to be as active as quagga mussels at this time of year. We will refer to
152
these samples throughout as the ‘standard’ ballast and harbor samples.
Page 8 of 38
153 154
Spiked experimental ballast and harbor samples
155 156
To further test the effectiveness of eDNA sampling, we spiked a second 250mL aliquot of ballast
157
water samples with a small fragment of tissue (range in tissue amount: 66 – 520 µg) from the
158
aquatic invasive species zebra mussel (Dreissena polymorpha). As discussed above, zebra
159
mussels are present in the Great Lakes region, but not present in the water column as
160
temperatures cool during the fall, such as when we sampled in October of 2012.37,38 Thus, the
161
prediction is that the spiked samples should serve as a positive control for the negative result
162
expected for zebra mussel from the field samples. All samples (N = 41) were screened for zebra
163
mussels using eDNA, species-specific mitochondrial DNA (mtDNA) markers (Table 1) and gel
164
electrophoresis prior to the spiking experiment, with 0% detection of zebra mussels. Each tissue
165
sample was removed from a preserved adult using fine-tipped dissecting forceps and added to the
166
eDNA filter to imitate the event that this target invasive species was collected during ballast
167
water screening. We spiked sixteen samples distributed equally across the four sampling
168
categories. In addition, we generated eight control eDNA samples where no target zebra mussel
169
tissue was added. The control eDNA extractions were done in parallel with the spiked samples.
170
We will refer to these samples henceforth as the ‘experimental’ ballast samples. All experimental
171
ballast samples were filtered through a 1.2-micron polycarbonate membrane filter (GE Water &
172
Process Technologies).
ACS Paragon Plus Environment
8
Page 9 of 38
Environmental Science & Technology
173 174
DNA extraction from filter and Polymerase Chain Reaction
175 176
eDNA was isolated from filters following the protocol described in of Jerde et al.23 and Egan et
177
al.24. In brief, the polycarbonate membrane filter from standard and experimental ballast and
178
harbor samples was combined with 700 µL of CTAB Buffer and 20 µL of Protinase K, vortexed
179
for 15 seconds, and incubated at 63ºC for 1hour. After incubation, 700 µL of a 24:1 chloroform :
180
isoamylalcohol solution was added to the solution and vortexed for 5 seconds, centrifuged at
181
14000 RPM for 10 minutes, and the supernatent was transferred to a new tube. Next, 1mL of
182
isopropanol was added, the sample was inverted gently to mix, and incubated at -20ºC for ≥4
183
hours. Samples were removed from -20ºC, centrifuged at 14000 RPM for 20 minutes to form a
184
pellet, the supernatent was poured off, and washed twice with 500µL 70% ethanol, pouring off
185
ethanol and maintaining the pellet each time. Pellets were dried for 10 minutes at 45ºC in a
186
vacuum centrifuge and then resuspended in 50µL TE Buffer overnight at 4ºC in a refrigerator.
187 188
We measured the DNA concentration of each sample extracted using a Qubit 2.0 fluorometer
189
(Life Technologies) with either the dsDNA Broad Range or High Sensitivity assay kits
190
(Invitrogen). Sample DNA concentration was first measured with the Broad Range kit, and
191
samples with undetectable DNA concentration were then measured again with the High
192
Sensitivity kit. Two PCRs were performed on each sample. First, a species-specific primer was
193
used in concert with gel-electrophoresis for a standard detection. Second, a universal primer,
194
which amplifies the same section of the mtDNA in most invertebrates (Folmer et al. 1994), was
195
used in concert with the LTS device to test this new, field-ready device.
ACS Paragon Plus Environment
9
Environmental Science & Technology
Page 10 of 38
196 197
The first PCRs were performed using species-specific primers listed in Table 1 that would
198
amplify a section of the cytochrome c oxidase subunit I (COI) gene in the mtDNA for one of the
199
target species only (quagga or zebra mussel). The PCRs were used to confirm the success of
200
eDNA sampling from ballast and harbor water for species-specific detection of the two taxa. The
201
species-specific mtDNA markers were developed for previous studies where they were tested for
202
specificity, even against their closest related sympatric relative (i.e., quagga versus zebra and
203
vice versa), and some fragments were sequence verified for additional verification (Li et al.
204
2011, Mahon et al. 2011, Mahon et al. 2013, Egan et al. 2013). Species-specific PCR was
205
followed by standard gel electrophoresis detection analysis and included positive and negative
206
controls. The lab-based gel detection offered a baseline for detection to which we could compare
207
results from LTS detection. For each sample, nine replicate PCRs were run.
208 209
The second set of PCRs used universal invertebrate primers (Table 1), which amplify a ~600
210
base pair section of the mtDNA COI gene for both target and non-target invertebrate species.
211
The long term strategy for using universal primers is to develop a single PCR from which it
212
would be possible to use LTS to detect the presence or absence of multiple species of interest
213
using nanoparticles functionalized with different species-specific tags (i.e., single stranded
214
oligonucleotides). The current LTS sequence tag is ~340bp downstream of the 5’ end of the
215
~600bp mtDNA fragment amplified by the universal primers. In effect, the universal primers
216
would amplify all invertebrates present in a sample. Then, separate LTS analyses could be
217
conducted on the single PCR reaction using a different species specific set of functionalized
ACS Paragon Plus Environment
10
Page 11 of 38
Environmental Science & Technology
218
nanoparticles for each LTS trial to determine if a particular target species of invertebrate was
219
present in the sample.24,25,34
220
PCRs with both primer sets were performed as follows: (i) 94ºC for 1 minute, (ii) 94ºC for 30
221
seconds, (iii) 41ºC for 45 seconds, (iv) 72ºC for 1 minute, (v) go to step 2 and repeat 30 times,
222
and (vi) 72ºC for 8 minutes.39 The total reaction volume was 25 µL and included: 1µL DNA
223
template, 2.5µL 10X buffer, 2.5µL MgCl2 (25mM), 0.5µL forward primer (10µM), 0.5µL
224
reverse primer (10µM), 0.5µL dNTPs (2.5mM each; Invitrogen), 0.15µL Taq polymerase
225
(5Prime), and 17.35µL ddH2O.
226 227
Test of PCR inhibition
228 229
We also tested for PCR inhibition in each sample. Abiotic substances, such as humic and fulvic
230
acids, are known PCR inhibitors that are likely to be suspended or dissolved in ballast and port
231
waters, where they are found in soils and sediments and get turned up a lot due to boat traffic and
232
repeated dredging within the ports.40 To test for PCR inhibition, we compared the observed and
233
expected concentration of a positive control DNA sequence (Eurogentec Universal Exogenous
234
qPCR Positive Control) using quantitative real-time PCR on an Eppendorf MasterCycler
235
realplex machine. The master mix for each sample included 3.6 µL H2O, 10.0 µL ABI TaqMan
236
Environmental, 2.0 µL positive control mix (Eurogentec), 0.4 µL positive control DNA
237
(Eurogentec), and 4.0 µL template DNA.
238 239
Nanoparticle preparation and hybridization to eDNA samples for LTS
240
ACS Paragon Plus Environment
11
Environmental Science & Technology
Page 12 of 38
241
To detect target species using LTS,25,34 we first bound polystyrene beads with a carboxylated
242
surface to a species-specific oligonucleotide tag that will only bind to DNA fragment amplified
243
in the target species (see Table 1 for species specific tags). To generate these ‘functionalized’
244
nanoparticles for species detection we followed published protocols.24,25,34 Briefly, we combined
245
10 µL of Polybead Carboxylate Microspheres (diameter = 209 nm; [C] = 5.68x1012 beads/mL;
246
Polysciences, Inc.) with 1 mL of ddH2O, 50 µL of Bead Coupling Buffer (Polysciences, Inc.),
247
and 25 µL of an EDC coupling activating solution [EDC = 1-Ethyl-3-(3-dimethylaminopropyl)
248
carbodiimide HCl at 192 g/mol; solution = 1mg EDC in 5mL of ddH2O]. Bead solution was
249
mixed at 300 rpm on an orbital shaker for 15 minutes. Separately, 10 µL of each species-specific
250
oligonucleotide tag ([C] = 100 ng/µL) was combined with 1 mL of ddH2O and 50 µL of Bead
251
Coupling Buffer (Polysciences, Inc.) and mixed for 1 minute at 300 rpm on an orbital shaker.
252
The bead solution and the tag solution were then combined and mixed constantly at 300 rpm for
253
two to four hours. For downstream LTS measurement, a ‘reference’ solution containing all
254
reagents except the beads and tags was also made, where bead and tag volume was replaced by
255
equal amounts of ddH2O.
256 257
PCR products (double stranded DNA from the amplified mtDNA COI region) from each sample
258
was denatured by heating to 95ºC for 2 minutes, then immediately chilled on ice for 2 minutes.
259
PCR product was combined with functionalized bead solution ([C] = 1.046109 x 109 beads/mL)
260
in a 1:2 ratio at 48ºC for one minute.34 If target DNA was present, it would anneal (or bind) to
261
the complementary species-specific tag and the nanoparticle would grow in diameter. We have
262
previously shown that it does not grow when combined with non-target DNA.24,25,33,34 Each
263
PCR-bead hybridization reaction was then placed on ice until LTS measurement.
ACS Paragon Plus Environment
12
Page 13 of 38
Environmental Science & Technology
264 265
Laser Transmission Spectroscopy (LTS) and nanoparticle measurement
266 267
LTS is used in the context of species detection to measure the size of nanoparticles in
268
suspension. For the purpose of DNA diagnostics, LTS can detect the increase in size of the
269
functionalized nanoparticles if the species specific tag has bound to the target species PCR
270
product causing an increase in the diameter of the bead. LTS is based on measuring wavelength-
271
dependent light transmittance through a sample containing nanoparticles in suspension. Briefly,
272
the transmission of light through a sample cell containing particles plus suspension fluid is
273
recorded along with that of a similar cell containing only the suspension fluid. The fundamental
274
data-acquisition process involves measuring the wavelength-dependent transmission of light
275
(quantified as extinction) through an aqueous suspension of nanoparticles. Here, the pertinent
276
wavelength range is from ~300 to 1000 nm. Given the extinction information, and the known
277
wavelength-dependent properties of the beads, Mie theory 33 can be used to accurately determine
278
the bead diameter. The extinction data are analyzed and inverted by a computer algorithm that
279
outputs the particle size distribution. The sensitivity limit of LTS reported in Li et al.33 for 1025
280
nm polystyrene spheres is ~3580 particles/mL (i.e. 3.5x10-17 molar). The precision, accuracy,
281
sensitivity, and resolution of LTS using NIST traceable polystyrene particles are detailed in Li et
282
al.33, where these properties are quantified for the size range important for DNA detection (~50–
283
1000 nm). Further details of LTS theory and operation are detailed in Li et al.33 and Li et al.34
284 285
RESULTS & DISCUSSION
286
ACS Paragon Plus Environment
13
Environmental Science & Technology
Page 14 of 38
287
Species detection probabilities are typically low in aquatic environments, especially for the
288
ballast of large commercial shipping vessels.11-14 Yet, effective management of high-risk
289
invasive species requires the early detection of an incipient invasion, which increases the
290
feasibility of rapid responses to eradicate the species or at least contain its spread.15,41 Here, we
291
tested whether eDNA sampling can be an effective tool for species detection in heterogeneous
292
ballast and harbor water samples, which contain a variety of aquatic species, their eDNA filtered
293
tissues and cells, and a mixture of biotic and abiotic debris.
294 295
Amounts of eDNA isolated from harbor and ballast samples
296 297
First, the concentration of DNA isolated from ballast and harbor water samples was quantified
298
using a Qubit 2.0 flourometer (Life Technologies). Ballast intake and discharge were slightly
299
less than the harbor samples, which is consistent with ballast representing a subset of available
300
eDNA content in each harbor (mean ± SE: Indiana and Duluth harbors = 5.54 ± 1.32 µg/mL,
301
Ballast intake and discharge = 3.28 ± 0.65 µg/mL; t-test assuming unequal sample sizes: tdf=19.45
302
= 2.16, P = 0.04). However, ballast intake and discharge did not differ significantly (Ballast
303
intake = 3.67 ± 0.84 µg/mL, Ballast discharge = 2.68 ± 1.04 µg/mL; t-test assuming unequal
304
sample sizes: tdf=27.81 = 0.81, P = 0.43), nor did samples from the two harbors (Indiana source
305
harbor = 6.42 ± 0.95 µg/mL, Duluth destination harbor = 4.09 ± 1.17 µg/mL; t-test assuming
306
unequal sample sizes: tdf=2.29 = 1.24, P = 0.32). In addition, the DNA concentrations did not differ
307
between the spiked and control samples in our experiment (DNA [C] ± SE: spiked = 1.68 ± 0.21
308
µg/mL, control = 1.26 ± 0.29 µg/mL; t df=22 = 1.21, P = 0.24). These results are consistent with an
ACS Paragon Plus Environment
14
Page 15 of 38
Environmental Science & Technology
309
expectation from real-world ballast and harbor water samples, where the DNA from many
310
different organisms would be filtered and extracted from a sample.
311 312
Tests for PCR inhibition
313 314
Prior to screening eDNA samples for invasive species with each technology, we tested for PCR
315
inhibition since PCR inhibitors are likely present in ballast and harbor waters.40 There was no
316
evidence for systematic biotic or abiotic inhibition of PCRs in the study. We compared the
317
observed and expected concentration of a positive control DNA sequence added to each sample
318
using quantitative real-time PCR. The greater the difference generated between observed from
319
expected concentrations suggests an increased chance of PCR inhibition. None of the 41 standard
320
samples (ballast or harbor) or the 24 experimentally spiked samples (zebra or controls) showed
321
any evidence of PCR inhibition (all values were less than 1.4; mean ± SE: 0.64 ± 0.01).
322 323
Standard lab-based eDNA detection using gel electrophoresis
324 325
Each standard ballast water sample from the field was screened for quagga and zebra mussel
326
eDNA using species specific primers (Table 1) under the prediction that we should detect quagga
327
mussels at a higher rate than zebra mussels in the Great Lakes during October because they are
328
known to be more cold-tolerant.36-38 Based on gel electrophoresis, quagga mussel was detected in
329
35/41 (85%) of the field samples, representing 17/20 of the ballast intake samples, 5/5 of the
330
source harbor samples, 12/13 of the ballast discharge samples and 1/3 of the destination harbor
331
samples (Table 2). The six samples where quagga was not detected may represent: (1) a natural
ACS Paragon Plus Environment
15
Environmental Science & Technology
Page 16 of 38
332
failure to capture quagga mussel in the samples, (2) a failure to isolate quagga DNA in samples
333
in which it was present, but perhaps in low frequency, (3) PCR failure, or (4) a failure to detect
334
DNA by gel electrophoresis in PCRs in which quagga amplification was low. The overall DNA
335
concentration in six negative quagga samples did not, however, differ from the overall mean
336
(failed quagga detection samples: 2.36 ± 1.12 µg/mL; one-sample t-test: tdf=5 = -0.88, P = 0.42).
337
Thus, it was not the case that the negative quagga results for these six samples was due to little or
338
no DNA captured and extracted from the filters.
339 340
Zebra mussel eDNA was not detected in any standard field sample using gel electrophoresis
341
(0/41), as expected based on their natural history. However, zebra mussels were detected in
342
15/16 of the spiked environmental samples compared to 0/8 control samples (Table 2). The one
343
failed zebra detection had undetectable amounts of DNA present based on the Qubit analysis,
344
suggesting that zebra mussel tissue was not captured on the filter in this eDNA sample or that the
345
DNA extraction itself failed.
346 347
Experimental field-based eDNA detection using Light Transmission Spectroscopy
348 349
Prior to LTS species detection tests, nanoparticles were attached with species-specific
350
oligonucleotide tags and measured to be 230 (± 0.4) nm in diameter (Fig. 1A). This LTS size
351
distribution is larger than the original diameter of 209-nm for the carboxylated beads and
352
represents the binding of the oligonucleotide tags to the beads.24,34 Thus, 230 nm is the baseline
353
diameter of beads used to test for size shifts associated with positive detection of target quagga
354
and zebra mussel DNA.
ACS Paragon Plus Environment
16
Page 17 of 38
Environmental Science & Technology
355 356
LTS successfully detected quagga mussel eDNA in all 41 standard field samples, including the
357
six samples where quagga mussel was not detected using gel electrophoresis. The LTS result
358
suggests that the failure to detect quagga in the six cases by gel electrophoresis was due to a
359
reduced sensitivity of electrophoresis versus LTS to detect target DNA at low concentration. We
360
present LTS measurements from four representative positive quagga mussel detections from the
361
source harbor, ballast intake, ballast discharge, and receiving harbor in Figs. 1B-E, respectively.
362
Positive detection is indicated by two peaks in the size distribution of nanoparticles. The peak
363
having smaller diameter size represents functionalized nanoparticles that did not hybridize with
364
target DNA in the reaction, while the peak of larger-size diameter represents instances when the
365
oligonucleotide tag on the nanobead hybridized with target species DNA. Note, that in Figures
366
1B, 1C, and 1E, the estimated diameters of the nanoparticles measured by LTS matched the
367
predicted sizes of functionalized, non-hybridized nanoparticles (~230 nm) determined in control
368
experiments, with positive detection causing a size shift of from 140 to 160 nm for a subset of
369
nanoparticles that bound target species DNA to mean diameters of 370 to 390 nm (see Fig. 1F
370
for summary statistics for all 41 positive quagga mussel results). In comparison, in Figure 1D,
371
the two peaks were also present indicating a positive detection. However, both peaks were
372
shifted larger than the mean diameters expected for non-hybridizing (230 nm) and hybrizing
373
(~380 nm) nanoparticles. Variation from anticipated peak shifts was observed in 5 of the 57
374
positive results (8.8%) in the current study. The presence of the two peaks represented a positive
375
detection in these 5 cases; since we have never observed a double peak in any negative control in
376
this or any previous study.24,25 However, it is possible that interactions between bound and non-
377
bound nanoparticles may occasionally occur in the hybridization reaction to generate the greater
ACS Paragon Plus Environment
17
Environmental Science & Technology
Page 18 of 38
378
than expected shifting of peaks in positive detections. Further work is needed to test this
379
hypothesis. The average size shift (± standard deviation) in the diameter of nanoparticles relative
380
to the baseline beads attached with quagga-specific tags for all 41 standard samples (ballast and
381
harbor) was 139 ± 4.18 nm (Figure 1.F).
382 383
In the experimentally manipulated ballast and harbor water samples, LTS successfully detected
384
zebra mussel eDNA in 15/16 of the spiked samples (94%). The one sample where zebra mussel
385
was not detected by gel electrophoresis coincided with the sample that zebra mussel was not
386
detected by LTS. Thus, either DNA was not successfully isolated from this sample or the PCR
387
failed, which is to be expected with real-world sampling of eDNA. As was the case for gel
388
electrophoresis, all eight control samples were also negative for zebra mussel in LTS tests.
389
Figure 2 shows in panel (A) the LTS measurement of one control sample where no target zebra
390
mussel tissue was added, (B) the positive LTS result for one individual sample spiked with zebra
391
mussel tissue, and (C) the average size shift of 141 ± 6.33 nm standard deviation observed in the
392
diameter of nanoparticles relative to the baseline beads for the 15 cases of positive LTS
393
detection.
394 395
eDNA as an effective tool for species detection in ballast and habor waters
396 397
In unmanipulated field samples, eDNA species detection matched expectations from the natural
398
history of each species.36-38 Zebra mussel, which are known to be less cold-tolerant that quagga
399
mussels, was not detected, while quagga mussel was detected in most samples at a rate of 85%
400
for the gel assay and 100% for the LTS assay. In the spiked experimental samples, both assays
ACS Paragon Plus Environment
18
Page 19 of 38
Environmental Science & Technology
401
detected zebra mussel in 94% of spiked samples and 0% of negative control samples. Overall,
402
these results demonstrate that eDNA sampling is effective for monitoring ballast-mediated
403
invasions and that LTS is a sensitive species-specific assay with the potential for rapid, field-
404
based detection. Moreover, using our eDNA sampling and species detection methods, we found
405
clear evidence of naturally occurring quagga mussel DNA in all four sampling localities (ballast
406
and harbor water from the Indiana and Duluth harbors). Using traditional methods, it is common
407
to find the dreissenid veligers (i.e., planktonic larval forms) in the ballast tanks of large ships,
408
however they are nearly impossible to determine morphologically to species based on standard
409
microscopy.42 Thus, in addition to detection, eDNA can offer greater taxonomic clarity.
410 411
One previous study has applied species-specific PCR to ballast water surveillance43 in which
412
standard PCR was used to screen for starfish larvae from bulk plankton samples, not eDNA.
413
Deagle et al.43 showed it worked with no false positives detected and could detect 1
414
larvae/200mg plankton. In general, the eDNA approach has been applied successfully for species
415
detection to multiple aquatic organisms (fish, amphibians, and invertebrates) worldwide22,23,26,29-
416
31
417
application in conservation and invasive species management.21,32 Our results suggest that eDNA
418
species detection from the ballast water of large transoceanic ships is a natural extension of the
419
eDNA approach.
and is rapidly developing into a burgeoning field of research in molecular ecology and
420 421
Next steps: Moving into rapid field-based detection with LTS
422 423
We also demonstrate that the field-ready LTS is a very sensitive technique for rapid target
ACS Paragon Plus Environment
19
Environmental Science & Technology
Page 20 of 38
424
species detection from ballast and harbor waters. Previous research regarding LTS as a DNA
425
detection tool has shown that it can work effectively in species detection using pure lab samples
426
with one species present,34 when in the presence of the DNA of other closely related species,25,33
427
at low concentrations, 25 and for real-world samples, where there are a variety of species present
428
and the target DNA fragment is amplified with universal invertebrate primers and species
429
specific LTS tagged nano beads.24 In fact, detection capabilities of this LTS system suggest
430
potential detection to the picomolar (10−12) range.25 The LTS device fits into a large briefcase
431
and could run on a car battery in the field, thus LTS offers the opportunity to become a rapid
432
DNA detection method for real-world applications.33,34 Our ultimate objective is to apply the
433
technology for rapid point of sampling detection in nature. To this end, LTS technology is
434
portable, fast, and relatively inexpensive; we have recently developed a beta-model LTS device
435
that together with a field-based PCR unit fit in a “carry-on” sized suitcase, it can run for a day
436
from a rechargeable 12-volt battery, and following PCR amplification can measure a sample in
437
less than 10 seconds. The next major phase of research is two fold: (1.) integrate already
438
developed field DNA extraction techniques 44,45 and field PCR technologies 46,47 to be fully field-
439
ready and (2.) conduct LTS testing on site and eventually to transfer the technology to personnel
440
on the front line of invasive species detection in the field. In this regard, one extremely valuable
441
application of LTS would be to begin rapidly screening the ballast tanks of transoceanic ships
442
coming into the U.S. for invasive species. An urgent need exists for rapid detection methods like
443
LTS to identify harmful invasive species in ballast prior to port entry to provide adequate time
444
for implementing appropriate management actions.
445
ACS Paragon Plus Environment
20
Page 21 of 38
Environmental Science & Technology
446
The general purpose of our work is to enable improved environmental protection, which requires
447
timely data to determine if potentially invasive species are contained in particular ships entering
448
the Great Lakes and other coastal regions of the world. Only then can the risks be understood at
449
the spatial and temporal scales needed to make appropriate decisions to reduce the potential for
450
harm. We feel that the further development and implementation of the technology tested here
451
could provide important information to guide decisions in the public and private sector, including
452
those by government agencies, the shipping industry, insurers, and other stakeholders, as it will
453
become possible to rapidly test for the presence of target organisms in specific ships. For
454
example, the tools described here could enable regulators to ensure compliance in actionable
455
time frames, and could be used by ship operators to document compliance with various enacted
456
and proposed ballast water concentration standards.48 Finally, ballast and harbor water
457
monitoring is only one of many potential applications for the tested technology, including
458
detection of invasive species in ports or other environments, threatened or endangered species,
459
species of public health concern (e.g., pathogens, parasites, or indicator species), and species of
460
concern for biosecurity.
461 462
Acknowledgements
463 464
For providing funding for this research, we would like to thank the U.S. Environmental
465
Protection Agency’s Great Lakes Restoration Initiative, and two interdisciplinary research
466
initiatives funded by the University of Notre Dame, the Environmental Change Initiative and the
467
Advanced Diagnostics and Therapeutics Initiative. We would also like to thank the American
468
Steamship Company and the captain and crew of the M/V Indiana Harbor for access to ballast
ACS Paragon Plus Environment
21
Environmental Science & Technology
469
Page 22 of 38
and harbor water.
470 471
ACS Paragon Plus Environment
22
Page 23 of 38
472
Environmental Science & Technology
REFERENCES
473 474
1. Pejchar, L.; Mooney, H. Invasive species, ecosystem services and human well-being. Trends
475
Ecol. Evol. 2009, 24, 497–504.
476 477
2. Pimentel, D.; Lach, L.; Zuniga, R.; Morrison, D. Environmental and economic costs of
478
nonindigenous species in the United States. Bioscience 1999, 50, 53–65.
479 480
3. Pimentel, D.; Zuniga, R.; Morrison, D. Update on the environmental and economic costs
481
associated with alien-invasive species in the United States. Ecol. Econ. 2005, 52, 273–288.
482 483
4. Wilcove, D.S.; Rothstein, D.; Dubow, J.; Phillips, A.; Losos, E. Quantifying threats to
484
imperiled species in the United States. BioScience 1998, 48, 607–615.
485 486
5. Sala, O.E.; Chapin, F.S.; Armesto, J.J.; Berlow, E.; Bloomfield, J.; et al. Global biodiversity
487
scenarios for the year 2100. Science 2000, 287, 1770–1774.
488 489
6. Mills, E.L.; Leach, J.H.; Carlton, J.T.; Secor, C.L. Exotic species in the Great Lakes: a history
490
of biotic crises and anthropogenic introductions. J. Great Lakes Res. 1993, 19, 1–54.
491 492
7. Ricciardi, A.; MacIsaac, H.J. Recent mass invasion of the North American great lakes by
493
Ponto-Caspian species. Trends Ecol. Evol. 2000, 15, 62–65.
494
ACS Paragon Plus Environment
23
Environmental Science & Technology
495
8. Ricciardi, A. Patterns of invasion of the Laurentian Great Lakes in relation to changes in
496
vector activity. Diversity and Distributions 2006, 12, 425–433.
Page 24 of 38
497 498
9. Keller, R.P.; Lodge, D.M.; Shogren, J.F.; Lewis, M. (2009) Bioeconomics of Invasive Species:
499
Integrating Ecology, Economics, Policy, and Management; Oxford University Press, Oxford,
500
U.K., 2009.
501 502
10. Rothlisberger, J.D.; Finnoff, D.C.; Cooke, R.M.; Lodge, D.M. Ship-borne nonindigenous
503
species diminish Great Lakes ecosystem services. Ecosystems 2012, 15, 472–476.
504 505
11. Duggan, I.C.; Van Overdijk, C.D.A.; Bailey, S.A.; Jenkins, P.T.; Limen, H. et al.
506
Invertebrates associated with residual ballast water and sediments of cargo-carrying ships
507
entering the Great Lakes. Can. J. Fish. Aquat. Sci. 2012, 62, 2463–2474.
508 509
12. Drake, J.M.; Lodge, D.M. Global hotspots of biological invasions: evaluating options for
510
ballast-water management. Proc. Roy. Soc. B 2004, 271, 575–580.
511 512
13. Drake, J.M.; Lodge, D.M. Hull fouling is a risk factor for intercontinental species exchange
513
in aquatic ecosystems. Aquatic Invasions 2007, 2, 121–131.
514 515
14. Drake, J.M.; Lodge, D.M. Rates of species introductions in the Great Lakes via ships’ ballast
516
water and sediments. Can. J. Fish. Aquat. Sci. 2007, 64, 530–538.
517
ACS Paragon Plus Environment
24
Page 25 of 38
Environmental Science & Technology
518
15. Lodge, D.M.; Williams, S.; MacIsaac, H.J.; et al. Biological invasions: recommendations for
519
U.S. policy and management. Ecol. Appl. 2006, 16, 2035–2054.
520 521
16. Tsolaki, E.; Diamadopoulos, E. Technologies for ballast water treatment: a review. J. Chem.
522
Technol. Biot. 2010, 85, 19–32.
523 524
17. Albert, R.J.; Lishman, J.M.; Saxena, J.R. Ballast water regulations and the move toward
525
concentration-based numeric discharge limits. Ecol. Appl. 2013, 23, 289–300.
526 527
18. Lawrence, D.J.; Cordell, J.R. Relative contributions of domestic and foreign sourced ballast
528
water to propagule pressure in Puget Sound, Washington, USA. Biol. Cons. 2010, 143, 700–709.
529 530
19. Bott, N.J.; Ophel-Keller, K.M.; Sierp, M.T.; Herdina; Rowling, K.P. et al. Toward routine,
531
DNA based detection methods for marine pests. Biotechnol. Adv. 2010, 28, 706–714.
532 533
20. Schwartz, M.K.; Luikart, G.; Waples, R.S. Genetic monitoring as a promising tool for
534
conservation and management. Trends Ecol. Evol. 2007, 22, 25–33.
535 536
21. Lodge, D.M.; Turner, C.R.; Jerde, C.L. et al. Conservation in a cup of water: estimating
537
biodiversity and population abundance from environmental DNA. Mol. Ecol. 2012, 21, 2555–
538
2558.
539
ACS Paragon Plus Environment
25
Environmental Science & Technology
540
22. Thomsen, P.F.; Kielgast, J.; Iversen, L.L. et al. Monitoring endangered freshwater
541
biodiversity using environmental DNA. Mol. Ecol. 2012, 21, 2565–2573.
Page 26 of 38
542 543
23. Jerde, C.L.; Mahon, A.R.; Chadderton, W.L.; Lodge, D.M. “Sight-unseen” detection of rare
544
aquatic species using environmental DNA. Conserv. Lett. 2011, 4, 150–157.
545 546
24. Egan, S.P.; Barnes, M.A.; Hwang, C.-T.; Mahon, A.R.; Feder, J.L.; Ruggiero, S.T.; Tanner,
547
C.E.; Lodge, D.M. Rapid invasive species detection by combining environmental DNA with
548
Light Transmission Spectroscopy. Conserv. Lett. 2013, 6, 402–409.
549 550
25. Mahon, A.R.; Barnes, M.A.; Li, F.; et al. DNA-based species detection capabilities using
551
Laser Transmission Spectroscopy. J. Roy. Soc. Int. 2013, 10, 20120637.
552 553
26. Jerde CL, Chadderton WL, Mahon AR, Renshaw MA, Corush J, et al. Detection of Asian
554
carp DNA as part of a Great Lakes basin-wide surveillance program. Can. J. Fish. Aquat. Sci.
555
2013, 70, 522–526.
556 557
27. Amos, W.; Whitehead, H.; Ferrari, M.J.; Glockner-Ferrari, D.A.; Payne, R.; Gordon, J.
558
Restrictable DNA from sloughed cetacean skin: its potential for use in population analysis. Mar.
559
Mam. Sci. 1992, 8, 275–283.
560 561
28. Ficetola, G.F.; Miaud, C.; Pompanon, F.; Taberlet, P. Species detection using environmental
562
DNA from water samples. Biol. Lett. 2008, 4, 423–425.
ACS Paragon Plus Environment
26
Page 27 of 38
Environmental Science & Technology
563 564
29. Dejean, T.; Valentini, A.; Duparc, A. et al. Persistence of environmental DNA in freshwater
565
ecosystems. PLoS ONE 2011, 6, e23398.
566 567
30. Goldberg, C.S.; Pilliod, D.S.; Arkle, R.S.; Waits, L.P. Molecular detection of vertebrates in
568
stream water: a demonstration using Rocky Mountain Tailed Frogs and Idaho Giant
569
Salamanders. PLoS ONE 2011, 6, e22746.
570 571
31. Minamoto, T.; Yamanaka, H.; Takahara, T.; Honjo, M.N.; Kawabata, Z. Surveillance of fish
572
species composition using environmental DNA. Limnology 2012, 13, 193-197.
573 574
32. Taberlet, P.; Coissac, E.; Hajibabaei, M.; Rieseberg, L.H. Environmental DNA. Mol. Ecol.
575
2012, 21, 1789–1793.
576 577
33. Li, F.; Schafer, R.; Hwang, C.-T.; Tanner, C.E.; Ruggiero, S.T. High-precision sizing of
578
nanoparticles by laser transmission spectroscopy. Appl. Opt. 2010, 49, 6602–6611.
579 580
34. Li, F.; Mahon, A.R.; Barnes, M.A.; Feder, J.; Lodge, D.M. et al. Quantitative and rapid DNA
581
detection by Laser Transmission Spectroscopy. PLoS ONE 2010, 6, e29224.
582 583
35. Darling, J.A.; Mahon, A.R. From molecules to management: Adopting DNA-based methods
584
for monitoring biological invasions in aquatic environments. Environ. Res. 2011, 111, 978–988.
585
ACS Paragon Plus Environment
27
Environmental Science & Technology
Page 28 of 38
586
36. MacIsaac, H.G. Comparative growth and survival of Dreissena polymorpha and Dreissena
587
bugensis, exotic mollusks introduced to the Great Lakes. J. Great Lakes Res. 1994, 20, 783–790.
588 589
37. Mills, E.L.; Rosenberg, G.; Spidle, A.P.; Ludyanskiy, M.; Pligin, Y.; May, B. A review of
590
the biology and ecology of the quagga mussel (Dreissena bugensis), a second species of
591
freshwater Dreissenid introduced to North America. Am. Zool. 1996, 36, 271–286.
592 593
38. Claxton,W.T.; Mackie, G.L. Seasonal and depth variations in gametogenesis and spawning
594
of Dreissena polymorpha and Dreissena bugensis in eastern Lake Erie. Can. J. Zool. 1998, 76,
595
2010–2019.
596 597
39. Folmer, O.; Black, M.; Hoeh, W.; Lutz, R.; Vrijenhoek, R. DNA primers for amplification of
598
mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. Mol. Mar.
599
Biol. Biotech. 1994, 3, 294–299.
600 601
40. Hartman, L.J.; Coyne, S.R.; Norwood, D.A. Development of a novel internal positive control
602
for Taqman based assays. Mol. Cell. Probe 2005, 19, 51–59.
603 604
41. Hulme, P.E. Beyond control: wider implications for the management of biological invasions.
605
J. Appl. Ecol. 2006, 43, 835–847.
606 607
42. Domm, S.; McCauley, R.W.; Kott, E. (1993) Physiological and taxonomic separation of two
608
dreissenid mussels in the Laurentian Great Lakes. Can. J. Fish. Aquat. Sci. 1993, 50, 2294–2298.
ACS Paragon Plus Environment
28
Page 29 of 38
Environmental Science & Technology
609 610
43. Deagle, B.E.; Bax, N.; Hewitt, C.L.; Patil, J.G. Development and evaluation of a PCR-based
611
test for detection of Asterias (Echinodermata:Asteroidea) larvae in Australian plankton samples
612
from ballast water. Mar. Freshwater Res. 2003, 54, 709–719.
613 614
44. Tomlinson, J.A.; Dickinson, M.J.; Boonham, N. Rapid detection of Phytophthora ramorum
615
and P. kernoviae by two-minute DNA extraction followed by isothermal amplification and
616
amplicon detection by generic lateral flow device. Phytopathology 2010, 100,143–149.
617 618
45. Kasajima, I.; Ide, Y.; Ohkama-Ohtsu, N.; Hayashi, H.; Yoneyama, T.; Fujiwara, T. A
619
protocol for rapid DNA extraction from Arabidopsis thaliana for PCR analysis. Plant Molecular
620
Biology Reporter 2004, 22, 49-52.
621 622
46. Chang, H.-F. G.; Tsai, Y.-L.; Tsai, C.-F.; Lin, C.-K.; Lee, P.-Y.; Teng, P.-H.; Su, C.; Jeng,
623
C.-C. A thermally baffled device for highly stabilized convective PCR. Biotechnology Journal
624
2012, 7, 1-5.
625 626
47. Miller, S.M.; Tourlousse, D.M.; Stedtfeld, R.D.; Baushke, S.M.; Herzog, A.M.; Wick, L.M.;
627
Rouillard, J.M.; Gulari, E.; Tiedje, J.M.; Hashsham, S.A. An in situ synthesized Virulence and
628
Marker Gene (VMG) biochip for the detection of waterborne pathogens. Applied and
629
Environmental Microbiology 2008, 74, 2200-2209.
630
ACS Paragon Plus Environment
29
Environmental Science & Technology
Page 30 of 38
631
48. IMO. International Convention for the Control and Management of Ships’ Ballast Water and
632
Sediments. http://www.imo.org/, 2004.
633 634
49. Mahon, A.R.; Barnes, M.A.; Senapati, S.; Feder, J.L.; Darling, J.A. et al. Molecular detection
635
of invasive species in heterogeneous mixtures using a microfluidic carbon nanotube platform.
636
PLoS ONE 2011, 6, e17280.
ACS Paragon Plus Environment
30
Page 31 of 38
637
Environmental Science & Technology
Data Accessibility
638 639
All data presented in this paper will be publicly available on the Dryad Digital Repository.
ACS Paragon Plus Environment
31
Environmental Science & Technology
640
Page 32 of 38
TABLE 1. List of species-specific and universal PCR primers and species-specific tags for LTS used in the present study,
641
Oligonucleotide
Sequence
Ref.
Species-specific PCR primers Dreissena bugensis (quagga mussel) forward
5’-CCTTATTATCTGTTCGGCGTTTAG-3’
49
Dreissena bugensis (quagga mussel) reverse
5’-ACTTGTAACACCAATAGAAGTAC-3’
49
Dreissena polymorpha (zebra mussel) forward
5’-GGGATTCGGAAATTGATTGGTAC-3’
49
Dreissena polymorpha (zebra mussel) reverse
5’-GAATCTGGTCACACCAATAGATGTGC-3’
49
Universal invertebrate PCR primers Universal invertebrate forward (LCO-1490)
5’-GGTCAACAAATCATAAAGATATTG-3’
39
Universal invertebrate reverse (HCO-2198)
5’-TAAACTTCAGGGTGACCAAAAAATCA-3’
39
Species-specific tag for LTS Dreissena polymorpha (zebra mussel)
5’-GAATCTGGTCACACCAATAGATGTGC-3’
49
Dreissena bugensis (quagga mussel)
5’-ACAAGTTGGGGGTGGTTTAGGCGGGAGT-3’
49
642 643
ACS Paragon Plus Environment
32
Page 33 of 38
Environmental Science & Technology
644
TABLE 2. Results of quagga mussel (Dreissena bugensis) screens of the standard ballast and harbor water samples for the standard
645
gel and LTS assay methods. There were no positive hits for zebra mussel (D. polymorpha) using either of the two methods. PCR
646
inhibition was also not observed in any of the samples. (Gel assay 1: % of three separate triplicate PCR reactions, where at least one
647
positive reaction occurred out of the three; Gel assay 2: % of all nine total PCR reactions that was positive; LTS assay: % of samples
648
with a clear double peak demonstrating a positive species detection) eDNA sample
Standard
Experimental
Type
DNA Concentration
(number of samples)
(mean µg/mL ± sd)
Gel assay 1
Gel assay 2
LTS assay
Source harbor (5)
6.43 ± 1.01
100% (5/5)
89% (40/45)
100% (5/5)
Ballast intake (20)
3.86 ± 4.90
85% (17/20)
63% (113/180)
100% (20/20)
Receiving harbor (3)
4.01 ± 3.14
33% (1/3)
17% (3/18)
100% (3/3)
Ballast discharge (13)
2.68 ± 2.01
92% (12/13)
54% (63/117)
100% (13/13)
Spiked (16)
1.68 ± 0.79
94% (15/16)
94% (45/48)
94% (15/16)
Control (8)
1.26 ± 1.09
0% (0/8)
0% (0/24)
0% (0/8)
ACS Paragon Plus Environment
33
Environmental Science & Technology
649
Page 34 of 38
FIGURE LEGEND
650 651
Figure 1. Naturally infested ballast and harbor water samples
652 653
Environmental DNA-based quagga mussel detection using LTS from naturally infested ballast
654
water samples from the M/V Indiana Harbor and surrounding source and receiving harbor
655
waters. Density distribution for one LTS measurement of (A.) the control blank, where no eDNA
656
was added; four positive quagga mussel detections from the (B.) source harbor, (C.) ballast
657
intake, (D.) ballast discharge, (E.) receiving harbor; and (F.) the average shift (± standard
658
deviation) in the diameter of nanoparticles relative to the baseline beads attached with quagga-
659
specific tags for all 41 standard samples (ballast and harbor). Dashed line in pane (F.) marks zero
660
peak shift for nanoparticle and represents no detection of target.
661 662
Figure 2. Experimentally spiked ballast and harbor water samples
663 664
Environmental DNA-based zebra mussel detection using LTS from experimentally spiked ballast
665
water samples collected from the M/V Indiana Harbor and surrounding source and receiving
666
harbor waters. (A.) Density distribution of one LTS measurement of a control sample, where no
667
zebra mussel tissue was added, (B.) density distribution of one LTS measurement where zebra
668
mussel tissue was added, and (C.) the average shift (± standard deviation) in the diameter of
669
nanoparticles for ballast samples where zebra mussel was added (+) relative to the samples were
670
zebra mussel was not added (–). Dashed line in pane (C.) marks zero peak shift for nanoparticle
671
and represents no detection of target.
ACS Paragon Plus Environment
34
3x107
2x107
1x107
0x100 0
200
400
600
Quagga mussel positive Source harbor
B.
1x108
5x107
0x100 0
200
E.
1x107
Quagga mussel positive Ballast discharge
8x106 6x106 4x106 2x106 0x100 0
200
400
Diameter (nm)
600
Quagga mussel positive Ballast intake
C.
6x107
4x107
2x107
0x100 0
200
Diameter (nm)
600
4x107
Density distribution (1/mL*nm)
Density distribution (1/mL*nm)
Diameter (nm)
400
Density distribution (1/mL*nm)
Quagga mussel negative Control blank
8x107
Environmental Science & Technology
Quagga mussel positive Receiving harbor
E.
3x107
2x107
.
1x107
0x100 0
ACS Paragon Plus Environment
200
400
Diameter (nm)
400
600
Diameter (nm)
Mean shift in diameter (nm)
A.
2x108
Density distribution (1/mL*nm)
Density distribution (1/mL*nm)
7 4x10 Page 35 of 38
140
F.
E
100
60
20 E -20
600
Quagga eDNA
Control
A.
Zebra mussel negative Control
6x106
2x107
Environmental Technology Zebra Science mussel&positive
Spiked sample 1x107
4x106
8x106
2x106
4x106
0x100
B.
0
200
400
0x100 600 0
Mean shift in diameter (nm)
8x106
140
Page 36 of 38
E
100
60
20 E -20
200 400Environment 600 ACS Paragon Plus
C.
_
+ Zebra mussel
Page 37 of 38
Environmental Science & Technology
Taking an environmental DNA sample. 363x241mm (72 x 72 DPI)
ACS Paragon Plus Environment
Environmental Science & Technology
Laser transmission spectroscopy read-out and briefcase size device. 404x253mm (240 x 240 DPI)
ACS Paragon Plus Environment
Page 38 of 38