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