Degradation of Extracellular Antibiotic Resistance Genes with UV254

May 5, 2017 - Disinfected wastewater effluent contains a complex mixture of biomolecules including DNA. If intact genes conveying antibiotic resistanc...
0 downloads 0 Views 2MB Size
Subscriber access provided by UniSA Library

Article 254

Degradation of extracellular antibiotic resistance genes with UV treatment Pin Hsuan Chang, Brianna Juhrend, Terese M Olson, Carl F. Marrs, and Krista R. Wigginton Environ. Sci. Technol., Just Accepted Manuscript • Publication Date (Web): 05 May 2017 Downloaded from http://pubs.acs.org on May 7, 2017

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 27

Environmental Science & Technology

1

Degradation of extracellular antibiotic resistance

2

genes with UV254 treatment

3

Pin Hsuan Chang,† Brianna Juhrend,† Terese M. Olson, † Carl F. Marrs, ‡ and Krista R.

4

Wigginton†* †

5

University of Michigan, Ann Arbor, Michigan 48109, United States

6 7 8

Department of Civil and Environmental Engineering,



Epidemiology Department, University of Michigan, 1415 Washington Heights, Ann Arbor, Michigan 48109, United States

9 10 11 12 13 14 15 16 17 18 19

_____________________________________ *corresponding author: [email protected] Tel. (734) 763-9661 Fax. (734) 764-4292

Word Counts 243 words in Abstract + 4,326 words in Text + 1,500 words in Figures + 600 words in Tables = 6,669

1 ACS Paragon Plus Environment

Environmental Science & Technology

20

Abstract: Disinfected wastewater effluent contains a complex mixture of biomolecules including

21

DNA. If intact genes conveying antibiotic resistance survive the disinfection process,

22

environmental bacteria may take them up. We treated plasmid pWH1266, which contains

23

ampicillin resistance gene blaTEM-1 and tetracycline resistance gene tetA, with UV254 doses up to

24

430 mJ/cm2 and studied the ability of those genes to be acquired by Acinetobacter baylyi. The

25

plasmids required approximately 20-25 mJ/cm2 per log10 loss of transformation efficiency. We

26

monitored plasmid DNA degradation using gel electrophoresis and qPCR with both short

27

amplicons (~200 bps, representative of ARG amplicon lengths commonly used for

28

environmental monitoring) and long amplicons (800-1200 bps, designed to cover the entire

29

resistance genes). The rate of transformability loss due to UV254 treatment was approximately

30

20× and 2× larger than the rate of gene degradation measured with the short and long amplicons

31

qPCR, respectively. When extrapolated to account for the length of the entire pWH1266 plasmid,

32

the qPCR rate constants were 2–7× larger than the rate constants measured with transformation

33

assays. Gel electrophoresis results confirmed that DNA cleavage was not a major inactivating

34

mechanism. Overall, our results demonstrate that qPCR conservatively measures the potential for

35

a gene to be transformed by environmental bacteria following UV254 treatment.

36 37 38 39

2 ACS Paragon Plus Environment

Page 2 of 27

Page 3 of 27

Environmental Science & Technology

40 41 42 43

Introduction The proliferation of bacterial resistance to antibiotics results in more than two million

44

illnesses and 23,000 deaths each year in the U.S. alone.1 To combat the continued emergence and

45

spread of resistant bacteria, the World Health Organization (WHO) announced a global action

46

plan in 2015 that urges international participants to take action in controlling and monitoring the

47

spread of all forms of antimicrobial resistance.2

48

Although links between the release of antibacterial resistance genes into the environment and

49

the emergence of new resistant pathogens in the clinic are not fully developed, resistant bacteria

50

can transfer resistance genes to human pathogenic species, compromising antibiotic treatment

51

effectiveness.3,4 Experts therefore recommend minimizing the release of antibacterial resistance

52

to the environment. Municipal wastewater and livestock wastes have been identified as

53

significant sources of antibiotic resistant bacteria and their associated antibiotic resistance genes

54

(ARGs) in the environment,5,6 with wastewater treatment plant (WWTP) effluent leading to

55

increased levels of resistance genes and resistant bacteria downstream of effluent discharges.7–9

56

Antibiotic resistance is passed by either vertical gene transmission, where genetic information

57

is inherited from parent cells, or horizontal gene transfer (HGT), where a bacterium without

58

resistance acquires the necessary genes from mobile genetic elements. HGT can occur by

3 ACS Paragon Plus Environment

Environmental Science & Technology

59

conjugation (in which DNA is passed from a donor cell to an acceptor cell through direct cell-

60

cell contact), transduction (where bacteriophage introduce ARGs into microbial cells), or

61

transformation (where competent microbes take up free DNA from their surroundings).10,11

62

Wastewater carries different forms of ARGs, including DNA carried within bacteria and viruses,

63

and extracellular DNA.12 Each form has the potential to transfer resistance by the different

64

mechanisms described above. For either vertical gene transmission or conjugation to occur, the

65

bacterium carrying the gene must be viable so that it can pass the gene onto its daughter cells or

66

to another organism. Likewise, for transduction to occur, the virus containing the genes must be

67

infective.13 Transformation, however, does not require a viable or infective donor

68

microorganism; in this case, competent bacteria in the environment may pick up a gene present

69

in extracellular DNA. Consequently, wastewater treatment processes that kill bacteria containing

70

ARGs do not necessarily eliminate the potential for ARG transfer downstream. Indeed, most

71

wastewater treatment and animal waste processes are not designed to destroy nucleic acids in the

72

waste.

73

To track antibacterial resistance through wastewater treatment, both culture-based and

74

molecular-based methods are employed.14 Culture methods, like broth dilution, disk diffusion,

75

and selective plating, detect the presence of culturable antibiotic resistance bacteria.15,16 These

76

methods do not detect ARGs in viruses, extracellular DNA, or nonculturable bacteria.

77

Molecular-based quantitative PCR (qPCR), on the other hand, detects pieces of the ARGs,

78

regardless of whether they are in a viable cell.17 However, qPCR amplicons rarely encompass the

79

entire ARG, not to mention regions outside of the gene that are necessary for gene transfer

80

pathways.18 Furthermore, qPCR may not detect certain damaged ARGs that can be repaired by

81

the bacteria that pick up the DNA. Finally, enzymes involved with the ARG transfer pathways in

4 ACS Paragon Plus Environment

Page 4 of 27

Page 5 of 27

Environmental Science & Technology

82

the cells may have different fidelities to damaged DNA than polymerase enzymes used in

83

qPCR.19

84

The relationship between concentrations of ARGs detected by qPCR and the concentration of

85

ARGs present in extracellular DNA capable of transformation by bacteria in the environment is

86

unknown. To address this, we compared the kinetics of ARG inactivation by ultraviolet (UV254)

87

disinfection using qPCR and transformation assays. UV254 is increasingly used to treat

88

wastewater effluent and is known to target nucleic acids. We studied plasmid transformation in

89

Acinetobacter species as a proxy for transformation of extracellular DNA in the environment.

90

Acinetobacter species are ubiquitous in soil and aquatic environments, undergo spontaneous

91

gene transformation at high frequency,20,21 and can be opportunistic human pathogens.22

92 93

Experimental Methods

94

General Experimental Approach

95

Solutions of purified plasmids that contain antibiotic resistance genes (ARGs) were exposed to

96

different doses of UV254 (Figure 1). The treated plasmids were then subjected to transformation

97

assays to establish the reaction kinetics for the loss of ARG transformation ability. The same

98

samples were assayed by qPCR to track the reaction kinetics of the plasmid ARG regions. Short

99

amplicon qPCR was conducted to represent qPCR assays that are commonly conducted to track

100

ARGs in water treatment and the environment. Long amplicon qPCR was conducted to track the

101

modification of the entire resistance genes. The same samples were analyzed by gel

102

electrophoresis to detect plasmid nicks and breaks following UV treatment. Ultimately, the

103

kinetics of the ARG transformation loss and qPCR amplicon reactions were compared to assess

104

the effectiveness of qPCR at tracking the destruction of ARG transformation ability.

5 ACS Paragon Plus Environment

Environmental Science & Technology

105 106

Model transformation system

107

We adopted a model transformation system to mimic environmental transformation events

108

downstream of WWTP effluents. The system (kindly provided by Dr. Chuanwu Xi, University of

109

Michigan), was comprised of the A. baylyi strain AC811, a derivative of Acinetobacter sp. strain

110

BD413,20,23 and plasmid pWH1266 which is a construct of the E. coli plasmid pBR322 (4.4

111

kbps) and the cryptic plasmid pWH1277 (4.5 kbps) isolated from A. lwoffi (Figure 2).24 The

112

pBR322 portion of the plasmid carries the tetA tetracycline resistance gene and the ampicillin

113

resistance gene blaTEM-1 that encodes TEM-1 beta-lactamase (Figure 2).24–26 The pWH1266 was

114

propagated in E. coli strain TOP10.

115 116

Plasmid Extraction

117

Host E. coli TOP10 cells carrying the plasmid were incubated in 5 mL LB broth (Sigma-

118

Aldrich, St. Louis, MO) containing 20 µg/L tetracycline with shaking (180 rpm) at 37 °C until

119

they reached stationary growth phase (~ 16 hr). The cells were then treated with QIAprep Spin

120

Miniprep Kits (Qiagen, Valencia, CA)15 per the manufacturer’s instruction. Several extractions

121

were conducted concurrently to prepare the plasmid stocks. The resulting purified plasmid

122

solutions contained 20-60 ng/µL plasmid DNA, as measured by a NanoDrop ND-1000 (SI

123

Figure 1; ThermoFisher Scientific, Waltham, MA).

124 125 126

Transformation assays Transformation assays were conducted to quantify the ability A. baylyi to acquire antibiotic

127

resistance from the plasmid.15 Transformation efficiency of the plasmid/A. baylyi system was

128

determined as the ratio of the antibiotic resistant bacteria (ARB) colonies detected on selective 6 ACS Paragon Plus Environment

Page 6 of 27

Page 7 of 27

Environmental Science & Technology

129

plates to the total colonies detected on nonselective plates, as expressed by the following

130

equation:  

131

    =

132

where CFU stands for colony forming units.

133

Experiments were conducted to establish the optimal transformation efficiencies, including

134

plasmid concentration, the point in the A. baylyi growth curve the cells were inoculated with

135

plasmids, and the time to incubate the plasmids with the bacteria. Details of these experiments

136

are reported in the SI. Once optimized, the transformation assays were always conducted in the

137

same manner. In brief, aliquots of A. baylyi strain AC811 stored at -80 °C were thawed and

138

streaked on LB agar plates, and incubated overnight at 30 °C. Three colonies of A. baylyi were

139

inoculated in 5 mL LB broth and incubated for 18 hours at 30 °C and 140 rpm in order to reach

140

stationary growth phase (SI Figure S2). Subsequently, 0.5 mL of the cell suspension was diluted

141

10× in LB broth and incubated at 30 °C and 140 rpm for another 8 hours to reach late

142

exponential phase. In a sterile glass tube, 480 µL of the bacteria was combined with 20 µL

143

untreated or UV-treated plasmid solution (10 ng/µL), resulting in a final ARG concentration of

144

0.4 ng/µL. The mixture was incubated at 30 °C and 140 rpm for 24 hours. Serial dilutions of the

145

suspension were plated onto selective media plates (10 µg/mL tetracycline or 100 µg/mL

146

ampicillin) for transformant counts27,28 and onto LB plates for total cell counts.

  

(1)

147 148

Gel Electrophoresis

149

Potential plasmid breakage after UV254 treatment and enzyme digestion was assessed with DNA

150

gel electrophoresis. The pWH1266 plasmids were incubated with restriction enzymes BamHI

151

and PvuI (Thermo Fisher Scientific, Waltham, MA) at 37 °C for 1 hour to cut the plasmids inside

7 ACS Paragon Plus Environment

Environmental Science & Technology

152

and outside of the tetA gene, respectively. Gel electrophoresis of the plasmid DNA before and

153

after UV254 treatment and restriction enzyme treatment was conducted on 0.5% agarose gels at 3

154

V/cm for 60 minutes. The bands were visualized by SYBR Safe DNA gel staining (Life

155

Technologies, Carlsbad, CA) alongside the GeneRuler 1kb DNA ladders (ThermoFisher

156

Scientific, Waltham, MA).

157 158

qPCR measurements

159

The tetA gene (1191 bps) and blaTEM-1 gene (861 bps) on the pWH1266 plasmid were quantified

160

with both short amplicon and long amplicon qPCR (Table 1). The short amplicons covered a

161

fraction of the resistance genes (~200 bps), whereas the long amplicons were designed to cover

162

the entire tetA and blaTEM-1 gene sequences. Amplicon and primer sequences were identified

163

using the plasmid pBR322 sequence, which is the portion of plasmid pWH1266 that carries both

164

resistance genes. Sequences were acquired from the NCBI GenBank database (J01749.1). All

165

primers were designed with NCBI Primer-BLAST tool.

166 167

The qPCR measurements were performed on a Mastercycler RealPlex 2 (Eppendorf, Hamburg,

168

Germany) using a Fast EvaGreen qPCR Master Mix (Biotium). Standard curves were conducted

169

in triplicate, with 10-fold dilutions covering 5 orders of magnitude. Each 10 µL qPCR reaction

170

contained 5 µL of 2X Master Mix, 0.5 µL of forward and reverse primers at 100 µM, 0.1 µL of

171

DNA template, and 4.4 µL sterile DNase free water. The temperature profile included one cycle

172

at 95 °C for 2 minutes, 40 cycles at 95 °C for 5 seconds, the annealing temperature (TA, Table 1)

173

for 15 seconds, 72 °C for either 15 seconds (short amplicons) or 60 seconds (long amplicons),

174

and then a melt curve to verify specificity. qPCR assays amplifying short and long amplicons

8 ACS Paragon Plus Environment

Page 8 of 27

Page 9 of 27

Environmental Science & Technology

175

were successfully developed for both tetA and blaTEM-1 genes. The amplification efficiency was

176

0.90 ± 0.05 for tetA short amplicon (216 bps), 0.82 ± 0.08 for blaTEM-1 short amplicon (209 bps),

177

0.71 ± 0.05 for tetA long amplicon (1200 bps), and 0.74 ± 0.04 for blaTEM-1 long amplicon (861

178

bps). The average R2 value across all qPCR assays was 0.995 ± 0.003 (mean ± SD).

179 180 181

UV disinfection experiments The UV disinfection experiments were conducted in a collimated beam reactor (SI Figure S3)

182

housing four lamps that emit 254 nm germicidal UV254 (Philips TUV G15T8, Amsterdam,

183

Netherlands). The UV irradiance was 0.18 ± 0.01 mW/cm2 (mean ± SD) as determined with

184

chemical actinometry using potassium iodide. 29 The chemical actinometry and UV experiments

185

were conducted in the same 96-well plates and the UV intensity was similar across all plate

186

wells. For UV experiments, the plasmid stocks were diluted to 10 ng/µL in 50 µL DNase free

187

water in 96-well plates. At this concentration and sample depth, the transmittance of UV254 was

188

approximately 93%, thus shielding corrections were not deemed necessary. The samples were

189

exposed to UV254 for 1, 3, 5, 10, 20 and 40 minutes, corresponding to UV254 doses of 11, 32, 54,

190

108, 215, and 430 mJ/cm2. At each sample point, the remaining ARG amplicon concentration

191

was quantified with qPCR and the remaining transformation efficiency was measured with

192

transformation assays. Inactivation kinetics were established for the transformation and qPCR

193

data. Kinetics were expressed as first-order (ln C/C0 versus UV254 dose).

194 195

Statistical analyses

196

Statistical analyses were conducted in Microsoft Excel 2016 and StatPlus (AnalystSoft Inc.,

197

Walnut, CA). Data from replicate experiments were pooled for linear regression analyses. The

9 ACS Paragon Plus Environment

Environmental Science & Technology

198

first order reaction rate constants from different types of experiments were compared with

199

multiple linear regression analyses using StatPlus (AnalystSoft Inc., Walnut, CA). The null

200

hypothesis in all of the experimental comparisons (e.g., transformation assay vs. long amplicon

201

qPCR) was that the first order rate constants were not significantly different. The P values were

202

computed at a confidence level of 95%.

203 204 205

Results and Discussion After optimization, the transformation efficiencies of both tetracycline resistance from tetA

206

and ampicillin resistance from blaTEM-1 were consistently between 10-4 and 10-5. In other words,

207

one out of 10,000 or 100,000 A. baylyi cells attained resistance when exposed to the plasmids.

208 209 210

Loss of transformation efficiency during UV treatment. Transformation efficiencies of both tetracycline and ampicillin resistance decreased with

211

increasing UV254 doses (Figure 3). At 54 mJ/cm2, the transformation efficiencies dropped to

212

approximately 10-7 or lower, and beyond this dose, they were below the assay detection limit

213

(90% reductions in concentration

283

(Figure S6). For this reason, rate constants for the qPCR experiments were determined using

284

only data points where C/C0 > 0.1. Deviations from first order kinetics measured with qPCR

285

have previously been reported with UV treatment.17,39 Reaction kinetics measured with qPCR

286

did not change with replicate experiments on different days. For this reason, qPCR data was

287

pooled to find the first order rate constants.

288

For both genes, the LA first-order rate constants were larger than the SA first-order rate

289

constants (Figure 5). This was expected due to the larger number of targets in the LA sequence

13 ACS Paragon Plus Environment

Environmental Science & Technology

290

compared to the SA sequence. Both SAs and LAs had reaction kinetics that were slower than the

291

inactivation kinetics of antibiotic resistance transformation. For example, at a dose of 108

292

mJ/cm2, the tetA SA had decreased about 0.25 log10 units, whereas the LA had decreased about

293

1.5 log10 units; extrapolating the inactivation kinetics of tetracycline resistance transformation to

294

this dose results in a 5-log10 loss in the gene transformation efficiency. These results highlight

295

that qPCR of the entire resistance gene will underestimate the loss in a gene’s ability to be

296

transformed. Ultimately, the results demonstrate that using qPCR results to track ARGs through

297

a UV treatment process will provide conservative conclusions on the ability of DNA to transfer

298

resistance. For the case of our two ARGs, this is true even when the entire ARG is covered by

299

the selected amplicon.

300

The slower LA-qPCR reaction kinetics compared to that of the antibiotic resistance

301

transformation could be explained by the fact that modifications outside the tetA and blaTEM-1

302

genes are critical for ARG transformation and expression. In an earlier study of the shuttle

303

plasmid pWH1266, deletion and insertion analyses were conducted to identify specific plasmid

304

regions that were necessary for the transformation of ampicillin and tetracycline resistance in A.

305

calcoaceticus BD413.23 They identified a 1350-bp region in the cryptic pWH1277 portion of the

306

pW1266 plasmid construct that was critical for transformation of blaTEM-1 and tetA in A.

307

calcoaceticus BD413. The study did not attempt to identify elements on the pBR322 region of

308

the pW1266 plasmid construct that were critical for transformation in A. calcoaceticus BD413.

309

The rate of transformation efficiency loss with UV treatment was the same for ampicillin and

310

tetracycline resistance genes on the same plasmid, despite the fact that the blaTEM-1 and tetA

311

genes degraded at different rates, based on LA-qPCR results (Figures 5, Table S2). This suggests

312

that the transformation inactivation mechanism may be the same for both resistance types. 14 ACS Paragon Plus Environment

Page 14 of 27

Page 15 of 27

Environmental Science & Technology

313

To test if modifications anywhere in the plasmid (and detectable by qPCR) would obstruct

314

antibiotic resistance transformation, we extrapolated the damage measured by qPCR in a specific

315

amplicon to the theoretical damage across the entire plasmid based on the following relationship:

 C ′ C plasmid size  log  = log × C0 amplicon size  C0 

316

 "

(2)



317

where  ! is the extrapolated value of normalized ARG copies, and  is the

318

normalized values of detected amplicon copies. This approach has been applied to track the

319

infectivity of RNA viruses with RT-qPCR following UV254 treatment.18 The extrapolation

320

approach assumes that modifications in the plasmid are randomly distributed, or that damage in

321

the measured amplicon is representative of damage anywhere in the plasmid. By covering

322

approximately 25% of the plasmid sequence with our two long amplicons, we aimed to negate

323

issues raised when different regions have different susceptibilities to UV254.18 Indeed, our results

324

suggested that the blaTEM-1 gene reacted faster with UV than the tetA according to LA-qPCR

325

(Figure 3), despite the fact that the tetA gene was 339 bps larger than the blaTEM-1 gene.

326

We applied the extrapolation approach on qPCR data from all four of the measured amplicons

327

(tetA SA, tetA LA, blaTEM-1 SA, blaTEM-1, LA) to assess if any of the extrapolated qPCR rate

328

constants were in the general range of the transformation efficiency rate constants. For the

329

extrapolations, each qPCR data point was treated with equation 2 and linear regressions were

330

conducted on the extrapolated data, resulting in four new predicted rate constants (tetA SA-

331

extrap, tetA LA-extrap, blaTEM-1 SA-extrap, blaTEM-1 LA-extrap) and associated standard errors

332

(Figures 5, Table S1). For both genes, extrapolated LA-qPCR rate constants were larger than the

333

transformation assay rate constants (Figure 5,Table S2). Although the rate constants of the





15 ACS Paragon Plus Environment

Environmental Science & Technology

334

extrapolated SA-qPCR data were not statistically different than the transformation rate constants

335

for both genes (Table S2, p > 0.05), t-tests on replicate data collected at the largest tested dose

336

(54 mJ/cm2) did suggest that the extrapolated SA-qPCR overestimated the loss of gene

337

transformation (p-values = 3.4×10-3 for tetA comparison and 5.0×10-3 for blaTEM-1 comparison).

338

In summary, extrapolating any of the four amplicon first order rate constants measured by qPCR

339

to the entire plasmid overestimated the loss in a gene’s ability to be transformed. The

340

overestimation of the extrapolated qPCR results in predicting ARG transformation loss could be

341

due to certain regions of the plasmid not being critical for the effective transformation of the

342

antibiotic resistance genes.23 Another possible explanation is that the 25% of the plasmid

343

measured by qPCR is much more sensitive to UV254 than the rest of the plasmid, and thus our

344

extrapolated reaction rate constants overestimate the reaction rate constant of the entire plasmid.

345

It is also possible that the A. baylyi cells repair the UV-damaged DNA upon uptake; if this is the

346

case, then the damage we detect with qPCR does not necessarily inactivate the plasmid’s ability

347

to be transformed in A. baylyi. Bacteriophage DNA damaged by UV254 can be repaired when it

348

enters its host cell40 and the same may occur with plasmid DNA that enters the bacterial cell. A

349

number of potential DNA repair mechanisms have been identified in Acinetobacter species.41,42

350

When qPCR is applied to predict the inactivation of ARGs, we assume that the polymerase

351

enzymes in the qPCR technique (Taq polymerase in this study) detect the same DNA

352

modifications that inhibit the cellular mechanisms responsible for plasmid transformation and

353

ARG expression. Research on the chemistry of DNA photolysis with UV254 suggests that

354

cyclobutane-pyrimidine dimers (CPDs) and 6-4 photoproduct (6-4PP) are the major DNA

355

products.43 Both modifications occur between two adjacent pyrimidines, with thymine-thymine

356

(T-T) CPDs being the most prevalent product.43,44 The ability of polymerase enzymes in qPCR to 16 ACS Paragon Plus Environment

Page 16 of 27

Page 17 of 27

Environmental Science & Technology

357

read over DNA adducts depends on the type of modification, with T-T dimers having the biggest

358

impact on amplification efficiency.45 It is also dependent on the polymerase used for

359

amplification; Taq polymerases, which were used here, do not read over pyrimidine dimers.46

360

The impacts of specific photoproducts on DNA transformation mechanisms have not been

361

identified.

362

To summarize the qPCR and transformation results, the ARG transformation inactivation rate

363

constants did not correspond to the ARG rate constants measured by LA-qPCR. In this case, the

364

LA-qPCR results underestimated ARG transformation loss. The transformation inactivation rate

365

constants also did not correspond to the qPCR reaction rate constants when the results were

366

extrapolated to the size of entire plasmid. In this case, the extrapolated qPCR results

367

overestimated ARG transformation loss.

368

As mentioned above, the blaTEM-1 gene first order reaction rate measured with qPCR were faster

369

than the first order reaction rate of the tetA gene, despite the fact that the tetA gene is larger than

370

the blaTEM-1 gene. Previously, McKinney and Pruden studied reactions in ARGs due to UV254

371

radiation using qPCR and demonstrated that the disappearance of ARG amplicons correlated

372

strongly with the number of adjacent T-T bases in the amplicon targets (r = -0.93). 17 In their

373

analysis, neighboring thymidine bases on one of the two DNA strands were enumerated. We

374

applied the same approach, but summed up adjacent pyrimidine bases on both DNA strands due

375

to the fact that qPCR measures both strands (Table 2). Interestingly, the blaTEM-1 gene contains

376

~1.5× the number of adjacent T-T bases than tetA, but is only 72% the length of tetA. The higher

377

T-T base content of blaTEM-1 explains the LA-qPCR results; the blaTEM-1 gene reacted ~1.2 faster

378

than tetA (Table 2; p < 0.05). Our analysis demonstrates that the T-T base content of DNA is

379

arguably more important than DNA size when predicting a genome’s reactivity with UV254. 17 ACS Paragon Plus Environment

Environmental Science & Technology

Page 18 of 27

380

Environmental Implications

381

In this study, qPCR underestimated the loss of the DNA’s ability to be transformed into the

382

Acinetobacter species, even when the entire resistance gene was amplified. In other words, qPCR

383

detection of an ARG represents a conservative assessment of the gene’s ability to be transformed

384

when it is released into the environment. When the qPCR results were extrapolated to the length

385

of the entire plasmid, the results overestimated the loss of transformation efficiency. Future

386

research should expand to other transformation models to see if the inactivation kinetics differs

387

amongst bacterial species or types of extracellular DNA carrying the genes. When qPCR is

388

employed to track the loss of ARGs capable of being transformed into competent bacterial cells,

389

amplicons should therefore cover the entire gene sequence when possible. Although this will

390

result in conservative estimates of the ability of the gene to be transformed, it will be far less

391

conservative than when only a fragment of the gene is monitored.

392 393

Supporting Information Available. Details on how plasmid concentrations were measured

394

including figure of QubitTM and NanodropTM measurement comparison; details of the

395

transformation assay optimization experiments including figure of A. baylyi growth curve, figure

396

of plasmid concentration vs. transformation efficiency, and figure of transformation efficiency

397

vs. incubation time; schematic of the UV disinfection apparatus; figure of ARG reaction kinetics

398

measured with qPCR; tables of measured first-order rate constants and corresponding P-values

399

from multiple linear regressions. This information is available free of charge via the Internet at

400

http://pubs.acs.org.

401

Author Information

402

Corresponding author

18 ACS Paragon Plus Environment

Page 19 of 27

403

Environmental Science & Technology

*Phone: +1- (734)-763-9661; e-mail: [email protected]

404 405

Acknowledgments

406

Funding for this research was provided through the University of Michigan MCubed program.

407

We would also like to acknowledge the laboratory of Dr. Chuanwu Xi’s for assistance with

408

transformation assays. The findings in this study do not represent the views of the sponsors.

409 410 411

References

412 413 414 415 416 417 418 419 420 421 422 423 424 425 426 427 428 429 430 431 432 433 434 435 436 437 438 439 440 441 442 443 444

(1) (2) (3) (4) (5)

(6) (7)

(8) (9)

(10) (11) (12) (13)

(14) (15)

CDC. Antibiotic resistance threats in the United Sates, 2013; 2013. World Health Organization. Global action plan on antimicrobial resistance. 2015, Geneva: World Health Organization. Threlfall, E. J.; Ward, L. R.; Frost, J. A.; Willshaw, G. A. The emergence and spread of antibiotic resistance in food-borne bacteria. Int. J. Food Microbiol. 2000, 62 (1–2), 1–5. Barton, M. D. Impact of antibiotic use in the swine industry. Curr. Opin. Microbiol. 2014, 19, 9– 15. Munir, M.; Wong, K.; Xagoraraki, I. Release of antibiotic resistant bacteria and genes in the effluent and biosolids of five wastewater utilities in Michigan. Water Res. 2011, 45 (2), 681– 693. Heuer, H.; Schmitt, H.; Smalla, K. Antibiotic resistance gene spread due to manure application on agricultural fields. Curr. Opin. Microbiol. 2011, 14 (3), 236–243. Pei, R.; Kim, S. C.; Carlson, K. H.; Pruden, A. Effect of River Landscape on the sediment concentrations of antibiotics and corresponding antibiotic resistance genes (ARG). Water Res. 2006, 40 (12), 2427–2435. Czekalski, N.; Gascón Díez, E.; Bürgmann, H. Wastewater as a point source of antibioticresistance genes in the sediment of a freshwater lake. ISME J. 2014, 8 (7), 1381–1390. Amos, G. C. A.; Gozzard, E.; Carter, C. E.; Mead, A.; Bowes, M. J.; Hawkey, P. M.; Zhang, L.; Singer, A. C.; Gaze, W. H.; Wellington, E. M. H. Validated predictive modelling of the environmental resistome. ISME J. 2015, 9 (6), 1467–1476. Thomas, C. M.; Nielsen, K. M. Mechanisms of, and barriers to, horizontal gene transfer between bacteria. Nat. Rev. Microbiol. 2005, 3 (September), 711–721. Levy, S. B.; Marshall, B. Antibacterial resistance worldwide: causes, challenges and responses. Nat.Med. 2004, 10 (1078–8956 (Print)), S122–S129. Colomer-Lluch, M.; Jofre, J.; Muniesa, M. Antibiotic resistance genes in the bacteriophage DNA fraction of environmental samples. PLoS One 2011, 6 (3), e17549. Vettori, C.; Gallori, E.; Stotzky, G. Clay minerals protect bacteriophage PBS1 of Bacillus subtilis against inactivation and loss of transducing ability by UV radiation. Can. J. Microbiol. 2000, 46 (8), 770–773. Bouki, C.; Venieri, D.; Diamadopoulos, E. Detection and fate of antibiotic resistant bacteria in wastewater treatment plants: a review. Ecotoxicol. Environ. Saf. 2013, 91, 1–9. Modali, L. A. The Use of Laboratory and Agent-based Models to Evaluate the Role of Natural Transformation in Biofilms in the Formation and Spread of Antibiotic Resistant Bacteria in 19 ACS Paragon Plus Environment

Environmental Science & Technology

445 446 447 448 449 450 451 452 453 454 455 456 457 458 459 460 461 462 463 464 465 466 467 468 469 470 471 472 473 474 475 476 477 478 479 480 481 482 483 484 485 486 487 488 489 490 491 492 493 494 495

(16) (17)

(18) (19)

(20)

(21)

(22)

(23) (24)

(25) (26) (27)

(28)

(29) (30) (31) (32)

(33) (34)

Page 20 of 27

Water Systems. 2014. Palmen, R.; Hellingwerf, K. J. Uptake and processing of DNA by Acinetobacter calcoaceticus a review. Gene 1997, 192 (1), 179–190. McKinney, C. W.; Pruden, A. Ultraviolet disinfection of antibiotic resistant bacteria and their antibiotic resistance genes in water and wastewater. Environ. Sci. Technol. 2012, 46 (24), 13393–13400. Pecson, B. M.; Ackermann, M.; Kohn, T. Framework for using quantitative PCR as a nonculture based method to estimate virus infectivity. Environ. Sci. Technol. 2011, 45, 2257–2263. Palmen, R.; Buijsman, P.; Hellingwerf, K. J. Physiological regulation of competence induction for natural transformation in Acinetobacter calcoaceticus. Arch. Microbiol. 1994, 162 (5), 344– 351. Vaneechoutte, M.; Young, D. M.; Ornston, L. N.; Baere, T. De; Nemec, A.; Reijden, T. Van Der; Carr, E.; Tjernberg, I.; Dijkshoorn, L. Naturally transformable Acinetobacter sp. Strain ADP1 belongs to the newly described species Acinetobacter baylyi. Appl. Environ. Microbiol. 2006, 72 (1), 932–936. Fondi, M.; Bacci, G.; Brilli, M.; Papaleo, M. C.; Mengoni, A.; Vaneechoutte, M.; Dijkshoorn, L.; Fani, R. Exploring the evolutionary dynamics of plasmids: the Acinetobacter panplasmidome. Bmc Evol. Biol. 2010, 10(1), 59. Chen, T. L.; Siu, L. K.; Lee, Y. T.; Chen, C. P.; Huang, L. Y.; Wu, R. C. C.; Cho, W. L.; Fung, C. P. Acinetobacter baylyi as a pathogen for opportunistic infection. J. Clin. Microbiol. 2008, 46 (9), 2938–2944. Wu, J.; Xi, C. Evaluation of different methods for extracting extracellular DNA from the biofilm matrix. Appl. Environ. Microbiol. 2009, 75 (16), 5390–5395. Hunger, M.; Schmucker, R.; Kishan, V.; Hillen, W. Analysis and nucleotide sequence of an origin of DNA replication in Acinetobacter calcoaceticus and its use for Escherichia coli shuttle plasmids. Gene 1990, 87 (1), 45–51. Cooksey, R. C.; Clark, N. C.; Thornsberry, C. A gene probe for TEM type β-lactamases. Antimicrob. Agents Chemother. 1985, 28 (1), 154–156. McNicholas, P.; Chopra, I.; Rothstein, D. M. Genetic analysis of the tetA(C) gene on plasmid pBR322. J. Bacteriol. 1992, 174 (24), 7926–7933. Palmen, R.; Vosman, B.; Kok, R.; van der Zee, J. R.; Hellingwerf, K. J. Characterization of transformation-deficient mutants of Acinetobacter calcoaceticus. Mol. Microbiol. 1992, 6 (13), 1747–1754. Palmen, R.; Vosman, B.; Buijsman, P.; Breek, C. K.; Hellingwerf, K. J. Physiological characterization of natural transformation in Acinetobacter calcoaceticus. J. Gen. Microbiol. 1993, 139 (2), 295–305. Rahn, R. O. Potassium Iodide as a Chemical Actinometer for 254 nm Radiation: Use of Iodate as an Electron Scavenger. Photochem. Photobiol. 1997, 66 (6), 885–885. Thurston-Enriquez, J. a; Haas, C. N.; Riley, K.; Gerba, C. P.; Jacangelo, J. Inactivation of Feline Calicivirus and Adenovirus Type 40 by UV Radiation. 2003, 69 (1), 577–582. Gerba, C. P.; Gramos, D. M.; Nwachuku, N. Comparative inactivation of enteroviruses and adenovirus 2 by UV light. Appl. Environ. Microbiol. 2002, 68 (10), 5167–5169. Calgua, B.; Carratalà, A.; Guerrero-Latorre, L.; de Abreu Corrêa, A.; Kohn, T.; Sommer, R.; Girones, R. UVC innactivation of dsDNA and ssRNA viruses in water: UV fluences and a qPCR-based approach to evaluate decay on viral infectivity. Food Environ. Virol. 2014, 260– 268. Görner, H. Photochemistry of DNA and related biomolecules: Quantum yields and consequences of photoionization. Journal of Photochemistry and Photobiology, B: Biology. 1994, pp 117–139. Dietrich, G. J.; Szpyrka, A.; Wojtczak, M.; Dobosz, S.; Goryczko, K.; Zakowski; Ciereszko, A. Effects of UV irradiation and hydrogen peroxide on DNA fragmentation, motility and fertilizing ability of rainbow trout (Oncorhynchus mykiss) spermatozoa. Theriogenology 2005, 64 (8), 20 ACS Paragon Plus Environment

Page 21 of 27

496 497 498 499 500 501 502 503 504 505 506 507 508 509 510 511 512 513 514 515 516 517 518 519 520 521 522 523 524 525 526 527 528 529 530 531

Environmental Science & Technology

(35)

(36) (37)

(38) (39)

(40)

(41)

(42) (43) (44) (45)

(46)

1809–1822. Peak, J. G.; Peak, M. J. Ultraviolet light induces double-strand breaks in DNA of cultured human P3 cells as measured by neutral filter elution. Photochem. Photobiol. 1990, 52 (2), 387– 393. Opel, K. L.; Chung, D.; McCord, B. R. A study of PCR inhibition mechanisms using real time PCR. J. Forensic Sci. 2010, 55 (1), 25–33. Jung, D.; Cho, Y.; Meyer, J. N.; Giulio, R. T. Di. The long amplicon quantitative PCR for DNA damage assay as a sensitive method of assessing DNA damage in the environmental model, Atlantic killfish (Fundulus heteroclitus). 2010, No. 919, 1–13. Jia, H.; Guo, Y.; Zhao, W.; Wang, K. Long-range PCR in next-generation sequencing: comparison of six enzymes and evaluation on the MiSeq sequencer. Sci Rep 2014, 4, 5737. Zhuang, Y.; Ren, H.; Geng, J.; Zhang, Y.; Zhang, Y.; Ding, L.; Xu, K. Inactivation of antibiotic resistance genes in municipal wastewater by chlorination, ultraviolet, and ozonation disinfection. Environ. Sci. Pollut. Res. 2014, 7037–7044. Kadavy, D. R.; Shaffer, J. J.; Lott, S. E.; Thomas, A.; Bolton, C. E.; Gallimore, W. H.; Eugene, L.; Nickerson, K. W.; Kokjohn, T. A.; Wolf, T. A.; et al. Influence of infected cell growth state on bacteriophage reactivation levels. Society 2000, 66 (12), 5206–5212. Albarracín, V. H.; Pathak, G. P.; Douki, T.; Cadet, J.; Borsarelli, C. D.; Gärtner, W.; Farias, M. E. Extremophilic Acinetobacter strains from high-altitude lakes in argentinean puna: remarkable UV-B resistance and efficient DNA damage repair. Orig. Life Evol. Biosph. 2012, 42 (2), 201– 221. Norton, M. D.; Spilkia, A. J.; Godoy, V. G. Antibiotic resistance acquired through a DNA damage-inducible response in Acinetobacter baumannii. J. Bacteriol. 2013, 195 (6), 1335–1345. Sinha, R. P.; Häder, D.-P. UV-induced DNA damage and repair: a review. Photochem. Photobiol. Sci. 2002, 1 (4), 225–236. Brash, D. E.; Haseltine, W. A. UV-induced mutation hotspots occur at DNA damage hotspots. Nature 1982, 298 (5870), 189–192. Sikorsky, J. A.; Primerano, D. A.; Fenger, T. W.; Denvir, J. Effect of DNA damage on PCR amplification efficiency with the relative threshold cycle method. Biochem. Biophys. Res. Commun. 2004, 323 (3), 823–830. McDonald, J. P.; Hall, A.; Gasparutto, D.; Cadet, J.; Ballantyne, J.; Woodgate, R. Novel thermostable Y-family polymerases: Applications for the PCR amplification of damaged or ancient DNAs. Nucleic Acids Res. 2006, 34 (4), 1102–1111.

21 ACS Paragon Plus Environment

Environmental Science & Technology

532

Page 22 of 27

Table 1. Primers used in this study. Primer

Sequence (5’- 3’)

Primer Length

Amplicon Size

Annealing Temperature

Melting Temperature (℃ ℃)

(℃ ℃) tetA Long FW

CGTGTATGAAATCTAACAATGCGCT

25

tetA Long RV

CCATTCAGGTCGAGGTGGC

19

tetA Short FW

GACTATCGTCGCCGCACTTA

20

tetA Short RV

ATAATGGCCTGCTTCTCGCC

20

blaTEM-1 Long FW

TTACCAATGCTTAATCAGTGAGGC

24

blaTEM-1 Long RV

ATGAGTATTCAACATTTCCGTGTCG

25

blaTEM-1 Short FW

AATAAACCAGCCAGCCGGAA

20

blaTEM-1 Short RV

TTGATCGTTGGGAACCGGAG

20

1200

51.9

56.9 60.6

216

53.9

58.9 59.8

861

51.4

56.6 56.4

209

53.8

58.8 59.2

533 Table 2. Number of neighboring pyrimidine bases in amplicons. Dimer count*

Dimer count/amplicon length (%)

Amplicon Name (size bps)

TT

TC

CT

CC

total

TT

TC

CT

CC

Total

blaTEM-1 short (209)

30

21

35

41

127

14.4

10.0

16.7

19.6

60.8

tetA short (216)

18

32

26

26

102

8.3

14.8

12.0

12.0

47.2

blaTEM-1 long (861)

123

106

105

97

431

14.3

12.3

12.2

11.3

50.1

tetA long (1200)

80

149

132

199

560

6.7

12.4

11.0

16.6

46.7

* The TTTT sequence was counted as 3 neighboring thymidines (T1T2T3T)

22 ACS Paragon Plus Environment

Page 23 of 27

Environmental Science & Technology

Figure 1. Schematic depicting experimental approach. Purified plasmids containing ARGs are treated with UV254. The treated plasmids are subjected to transformation assays in Acinetobacter bayli to measure the extent that UV254 treatment has reduced the transformation efficiency of the ARGs on the plasmids. The treated plasmid samples are also analyzed with short amplicon- and long amplicon-qPCR to measure the reaction kinetics of ARGs on the plasmids. The qPCR and transformation reaction kinetics are compared. Gel electrophoresis is applied to assess the extent of plasmid DNA breaks following UV254 treatment.

23 ACS Paragon Plus Environment

Environmental Science & Technology

Page 24 of 27

Figure 2. Plasmid 1266, which is a construct of E. coli plasmid pBR322 and A. calcoaceticus cryptic plasmid pWH1277 (grey portion). 534 535 536 537 538 539 540 541

24 ACS Paragon Plus Environment

Page 25 of 27

Environmental Science & Technology

Figure 3. tetA and blaTEM-1 degradation versus UV dose measured with qPCR short amplicon (SA), long amplicon (LA), and transformation assays. The error bars indicate one standard deviation from the mean (n≥3). For visibility, data has been pooled. 542

25 ACS Paragon Plus Environment

Environmental Science & Technology

Figure 4. DNA electrophoresis gel of pHW1266 plasmids treated with 0 (control), 11, 32, 54, 108 mJ/cm2 UV or with restriction enzymes BamHI and PvuI. The transformation efficiencies of the samples are labeled at the bottom of the gel. 543

26 ACS Paragon Plus Environment

Page 26 of 27

Page 27 of 27

Environmental Science & Technology

544 1.0

1.0

tetA

0.8

0.4

k (cm2/mJ)

k (cm2/mJ)

0.8 0.6

blaTEM-1

0.6 0.4 0.2

0.2

0.0

0.0

Figure 5. Rate constants corresponding to SA-qPCR, LA-qPCR, extrapolated qPCR, and transformation assays. The values of k are included in the SI Table S1. The error bars represent 95% confidence intervals. qPCR rate constants were analyzed over the first 90% decrease in copy number concentrations. The rate constant values and statistical differences between rate constants are presented in Tables S1 and S2. 545 546

27 ACS Paragon Plus Environment