Rapid Molecular Detection of Invasive Species in ... - ACS Publications

Feb 16, 2015 - requires low electrical power to operate.24 LTS determines the size, shape, and number of ..... Dryad Digital Repository. □ REFERENCE...
0 downloads 0 Views 742KB Size
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