Nucleic Acid Photolysis by UV254 and the Impact of Virus

Aug 14, 2018 - Nucleic Acid Photolysis by UV254 and the Impact of Virus Encapsidation. Zhong Qiao , Yinyin ... Environmental Science & Technology. Gar...
2 downloads 0 Views 682KB Size
Subscriber access provided by University of Winnipeg Library

Environmental Processes

Nucleic acid photolysis by UV254 and the impact of virus encapsidation Zhong Qiao, Yinyin Ye, Pin Hsuan Chang, Devibaghya Thirunarayanan, and Krista R. Wigginton Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.8b02308 • Publication Date (Web): 14 Aug 2018 Downloaded from http://pubs.acs.org on August 14, 2018

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 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 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.

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

Nucleic acid photolysis by UV254 and the impact of

2

virus encapsidation

3

4

Zhong Qiao,† Yinyin Ye,† Pin Hsuan Chang,† Devibaghya Thirunarayanan ,† and Krista R.

5

Wigginton†* †

6 7

Department of Civil and Environmental Engineering,

University of Michigan, Ann Arbor, Michigan 48109, USA

8 9 10 11 12 13 14 15 16 17

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

18

198 words in Abstract + 5174 words in Text + 900 words in Figures + 300 words in Table + 45

19

words in Acknowledgements = Total 6617

Word Counts

20

1 ACS Paragon Plus Environment

Environmental Science & Technology

21

Abstract:

22

Determining the influence of higher order structure on UVC photolysis will help inform

23

predictions of nucleic acid fate and microorganism inactivation. We measured the direct UV254

24

photolysis kinetics of four model viral genomes composed of single-stranded and double-

25

stranded RNA (ssRNA and dsRNA, respectively), as well as single-stranded and double-stranded

26

DNA (ssDNA and dsDNA, respectively), in ultrapure water, in phosphate buffered saline (PBS),

27

and encapsidated in their native virus particles. The photolysis rate constants of naked nucleic

28

acids measured by qPCR (RT-qPCR for RNA) and normalized by the number of bases measured

29

in a particular sequence exhibited the following trend: ssDNA > ssRNA ≈ dsDNA > dsRNA. In

30

PBS, naked ssRNA bases reacted, on average, 24× faster than the dsRNA bases, whereas naked

31

ssDNA bases reacted 4.3× faster than dsDNA bases. Endogenous indirect photolysis involving

32

1

33

measured rate constants with rate constants reported in the literature shows a general agreement

34

amongst the nucleic acid UV254 direct photolysis kinetics. Our results underscore the high

35

resistance of dsRNA to UVC photolysis and demonstrate the role that nucleic acid structure and

36

solution chemistry play in photoreactivity.

O2 and •OH was ruled out as a major contributing factor in the reactions. A comparison of our

UV254

37

MS2 (ssRNA)

Φ6 (dsRNA)

ΦX174 (ssDNA)

T3 (dsDNA)

2 ACS Paragon Plus Environment

Page 2 of 27

Page 3 of 27

Environmental Science & Technology

38

Introduction

39

Nucleic acids are ubiquitous biopolymers that carry the code for all domains of life, yet several

40

types of nucleic acids are undesirable in the aquatic environment due to their potential impacts

41

on human and ecological health. Disinfection processes often target nucleic acids of problematic

42

microorganisms, but extracellular nucleic acids also pose a potential risk. Antimicrobial

43

resistance genes (ARGs) released from municipal wastewaters and animal farming activities, for

44

example, can transfer resistance to other microorganisms.1 Double-stranded RNA (dsRNA)

45

pesticides in genetically modified crops have the potential to be released into surface waters.2

46

Furthermore, many infectious environmental viruses consist of nucleic acids protected by a

47

simple protein coat.

48 49

The nitrogenous bases in DNA and RNA molecules strongly absorb UV radiation, leading to a

50

number of photochemical reactions. This high photoreactivity of nucleic acids is the main reason

51

why UVC is commonly employed for disinfection in water, air, or on surfaces. The effectiveness

52

of UV at inactivating microorganisms, along with the low levels of disinfection byproducts that

53

form relative to chemical disinfectants, has led to a sustained growth of the UV disinfection

54

industry. By 2019, the UV disinfection market is expected to reach nearly $1 billion, or 1/3 of

55

the total disinfection market.3 Due to the widespread application of UV disinfection, the reaction

56

kinetics that drive the inactivation of microorganisms and nucleic acids should be well-defined.

57 58

The pyrimidine bases in nucleic acids (thymine (T), cytosine (C), and uracil (U)) are more

59

photoreactive with UVC than purine bases (adenine (A), guanine (G)), and the reactions that take

60

place between neighboring thymine bases have been the most extensively studied to-date. The

3 ACS Paragon Plus Environment

Environmental Science & Technology

61

TT dimer products from neighboring pyrimidine bases include the cyclobutane pyrimidine dimer

62

(CPD), pyrimidine pyrimidone photoadducts (i.e., 6−4 products), and their Dewar valence

63

isomers.4 Of these, the CPD products have the highest quantum yields.5,6 In addition to TT

64

dimers, TC, CT, CC, and UU dimer products have also been characterized.5,6 Hydration, protein-

65

nucleic acid linkages, covalent crosslinks between complimentary strands, and backbone breaks

66

also take place in nucleic acids exposed to UVC.7 Pyrimidine photohydrate reactions are the

67

most prevalent of these, and form readily on uracil and cytosine bases, but not on thymine

68

bases.8 Consequently, hydration reactions are more important in RNA than in DNA.

69 70

The physical orientation of bases and the structures that surround them at the time of irradiation

71

impact the base photoreactivity. Dimer formation, for example, requires neighboring pyrimidine

72

bases to have well-aligned double bonds and the correct intrabond separation at the time of UV

73

exposure.9,10 Several studies have characterized the role that the flanking bases of neighboring

74

pyrimidines have on dimer quantum yields,11-13 which directly influence photoreactivity. Purine

75

neighboring bases, in general, and guanine bases at the 5’ end, in particular, reduce dimerization

76

quantum yields.12 Photohydrate formation is not as impacted by the chemistry of the neighboring

77

bases, but the reactions are impacted by access to water molecules.14 The presence of bases in

78

either a single-strand or double-strand can also affect quantum yields; in general, the increased

79

order of double-stranded nucleic acids decreases pyrimidine quantum yields.10,15 An early study

80

reported that viral ssRNA in solution is more susceptible to UVC inactivation than the

81

encapsidated ssRNA genome,16 but to our knowledge, the work has not been extended to other

82

viruses or genome types.

4 ACS Paragon Plus Environment

Page 4 of 27

Page 5 of 27

Environmental Science & Technology

83

Predicting the kinetics of nucleic acid UVC photolysis is complicated by the fact that many

84

quantum yield values available in the literature were collected with DNA oligomers with

85

repeated base sequences or under varied experimental conditions (e.g., UV wavelength, pH,

86

ionic strength, etc.). Furthermore, RNA photochemistry is comparatively much less represented

87

in the literature than DNA photochemistry, making it difficult to draw conclusions about the

88

overall relative reactivity of single stranded and double stranded RNA and DNA molecules.

89

RNA photochemistry is of particular interest for disinfection purposes because viruses can

90

contain genomes with ssRNA and dsRNA, in addition to ssDNA and dsDNA. In viruses, the

91

dsDNA viruses are typically considered more resistant to UVC radiation than viruses with other

92

genomes types.17 Interestingly, this is not due to relatively lower photoreactivity of dsDNA, but

93

due to the ability of host cells to repair dsDNA.18

94

To address the role that structure plays in nucleic acid photochemistry, we investigated the

95

reactions that take place in single-stranded and double-stranded DNA and RNA viral genomes

96

with UV254. As model systems, we compared photolysis rates of four viral genomes within and

97

outside of virus particles to assess the impact of virus particles incorporation on nucleic acid

98

reactivity. We employed quantitative polymerase chain reaction (qPCR) techniques to target

99

specific regions of the genomes, as these methods are commonly used for detecting nucleic acids

100

in water and tracking microorganism fate through water treatment processes. Finally, we

101

compared our measured reaction rate constants to those predicted with quantum yield values in

102

the literature.

103 104

EXPERIMENTAL METHODS

105

Virus preparation 5 ACS Paragon Plus Environment

Environmental Science & Technology

106

Four model viruses were selected to represent four types of nucleic acid genomes, namely MS2

107

(ss RNA), φ6 (ds RNA), φX174 (ssDNA), and T3 (dsDNA). The characteristics of these viruses

108

and their genomes are provided in Table 1.

109 110

Escherichia viruses MS2 and T3, and their corresponding E. coli hosts were purchased from

111

American Type Culture Collection (ATCC 15597 and 11303, respectively). MS2 and T3 stock

112

were propagated, enriched, and enumerated based on previously published methods.19 Enriched

113

virus stock solutions were purified using a Fast Protein Liquid Chromatography system (Econo,

114

Bio-Rad) equipped with a HiPrep Sephacryl S-400 HR column (GE). The purified stocks were

115

filter sterilized with 0.22 µm polyethersulfone (PES) membrane filters (Millipore). The final

116

MS2 and T3 virus stocks (~1011 PFU/mL) were stored in phosphate buffered saline (PBS; 5 mM

117

NaH2PO4, 10mM NaCl, pH 7.5) at 4 oC.

118 119

Pseudomonas virus φ6 and its host Pseudomonas syringae pv. phaseolicola were kindly

120

provided by Dr. Linsey Mar at Virginia Tech. φ6 was added to Pseudomonas syringae grown to

121

OD640 of 0.1 in Luria-Bertani (LB) medium (5 g L-1 NaCl) at 26 °C with a multiplicity of

122

infection (MOI) of 2, and the virus and host mixture was incubated for 7 to 9 hours. Due to its

123

lipid envelope, φ6 was concentrated and purified differently than the other three bacteriophages.

124

Cell lysates were filtered through 0.22 µm PES membranes, and were concentrated with a lab-

125

scale tangential flow filtration system (Millipore) with a 30 kDa cellulous filter. The φ6

126

concentrate was purified by centrifuging at 65,700 × g, 4 °C in a 10-40% (w/v) step sucrose

127

gradient for 1.5 hours, then by centrifuging at 65,700 × g, 4 °C in a 40-60% (w/v) linear sucrose

128

gradient for 15 hours overnight. The φ6 virus band was collected with a needle and was

129

transferred to PBS with a 100 kDa centrifugal ultrafilter (Millipore). The final φ6 virus stocks 6 ACS Paragon Plus Environment

Page 6 of 27

Page 7 of 27

Environmental Science & Technology

130

(~1012 PFU/mL) were filter sterilized through 0.22 µm PES membranes, and stored in PBS at -80

131

°C.

132 133

Escherichia virus φX174 and its corresponding bacterial host E. coli ATCC 13706 were kindly

134

provided by Dr. Charles Gerba at the University of Arizona. φX174 virus was propagated by the

135

agar overlay technique. Specifically, soft tryptic soy agar (TSA; 0.7% w/v agar) was mixed with

136

φX174 virus and host bacteria, was overlaid on hard TSA (1.5% w/v agar), and was incubated at

137

37 °C overnight. The soft agar layer was collected and diluted in PBS to release φX174 from the

138

agar. The agar was removed by centrifuging at 3,000 × g for 10 minutes. The supernatant

139

containing φX174 was treated with a chloroform extraction, filtered through a 0.22 µm PES

140

membrane, and stored in PBS at 4 °C.

141 142

UV254 Photolysis

143

UV photolysis of viral genomes were studied when the genomes were within the virus capsids

144

(i.e., encapsidated) and when they were extracted from the virus capsids (i.e., naked nucleic acids

145

in ultrapure nuclease-free water and in nuclease-free PBS). The prepared PBS was heated to 90

146

°C for 30 minutes to eliminate the activity of nucleases. For each UV treatment, the nucleic acids

147

were extracted from the virus stocks immediately before the UV254 treatment with Maxwell 16

148

Viral Total Nucleic Acid Purification Kits (Promega) according to the manufacturer’s

149

instructions. The concentrations of the extracted virus nucleic acids were determined with a

150

Qubit 2.0 (Life Technologies) in ng/mL and converted to genome copies (gc)/µL using the

151

genome’s molecular mass and Avogadro’s constant, and then diluted to a final concentration of

152

~106 gc/µL with ultrapure nuclease free water (ThermoFisher) or nuclease free PBS. The

153

encapsidated nucleic acids were prepared by directly diluting the virus stocks to a similar 7 ACS Paragon Plus Environment

Environmental Science & Technology

154

genome concentration ~106 gc/µL with PBS. At this concentration, the absorbance of our

155

experimental solutions at 254 nm was below 0.01; it was therefore not necessary to correct the

156

UV254 radiation doses for shielding.

157 158

Samples of the virus in PBS and the extracted nucleic acids in either ultrapure water or PBS (50

159

µL) were added to the wells of 96-well plate with flat bottoms (Costar). The plate was placed

160

approximately 25 cm below four 15 W germicidal UV lamps (model G15T8, Philips) inside a

161

collimated beam unit. The UV irradiance at 254 nm (0.17 mW/cm2) was determined by chemical

162

actinometry measurements.20 The virus and genome solutions were irradiated up to doses of 408

163

mJ/cm2 at room temperature. Aliquots of the experimental solutions were collected from the

164

samples periodically and stored at 4 °C in the dark. Dark control samples were stored at room

165

temperature in the dark for the duration of the longest UV exposure to capture background

166

nucleic acid decay. The experiments were repeated at least four times for every genome type.

167

Following UV treatment, the encapsidated viral nucleic acids were extracted from the virus

168

samples with Maxwell 16 Viral Total Nucleic Acid Purification Kits. The UV-treated samples

169

and dark controls were analyzed by qPCR or RT-qPCR immediately after the completion of the

170

UV experiments.

171 172

qPCR Assay for UV-treated viral genomes

173

Two sets of PCR primers were designed to target two regions of approximately 500 bases or base

174

pairs on each viral genome (Table 1, Table S1). The entire genome of ssRNA of MS2 and

175

dsRNA of φ6 were directly extracted from the purified virus stocks, and were used as standards;

176

gBlock standards containing two target regions (Integrated DNA Technologies) were purchased

177

for φX174 ssDNA and T3 dsDNA. Both experimental samples and standards were quantified in 8 ACS Paragon Plus Environment

Page 8 of 27

Page 9 of 27

Environmental Science & Technology

178

parallel using qPCR or RT-qPCR assays in a RealPlex2 Mastercycler system (Eppendorf). The

179

reaction rate constants of the target regions were calculated based on a first-order reaction model: ln

𝐶 = 𝑘!"# ∙ 𝐷!"!"# 𝐶!

180

Where 𝐶! is the initial concentration of viral genome segment; 𝐶 is the concentration of the viral

181

genome segment after UV254 treatment; 𝑘!"# is the first-order rate constant; 𝐷!"!"# is the UV254

182

dose (mJ/cm2).

183 184

RT-qPCR assay for MS2 ssRNA

185

The 20 µL one-step RT-qPCR reactions consisted of 10 µL of 2× GoTaq qPCR Master Mix

186

(Promega), 0.4 µL of 50× GoScript RT Mix, 0.6 µL of 10 µM forward primer, 0.6 µL of 10 µM

187

reverse primer, 6.4 µL of nuclease-free water, and 2 µL of RNA sample. The following

188

thermocycling conditions were used: 15 min at 40 °C, 10 min at 95 °C, 45 cycles of 95 °C for 15

189

s, 55 °C for 30 s, and 72 °C for 45 s, followed by a melting curve analysis from 68 to 95 °C for 5

190

minutes.

191 192

RT-qPCR assays for φ6 dsRNA

193

For φ6 genome samples, dsRNA samples were mixed with 10 µM forward primer and 10 µM

194

reverse primer at a volume ratio of 10:1.5:1.5. The sample-primer mixture was heated at 99 °C

195

for 5 minutes, and immediately chilled on ice before the RT-qPCR assays. The 20 µL RT-qPCR

196

reactions consisted of 10 µL of 2× GoTaq qPCR Master Mix, 0.4 µL of 50× GoScript RT Mix,

197

5.2 µL of the pre-treated sample-primer mixture, 4 µL of 5 M Betaine (Sigma-Aldrich), and 0.4

198

µL of nuclease-free water. The following thermocycling conditions were used: 15 min at 40 °C,

9 ACS Paragon Plus Environment

Environmental Science & Technology

199

10 min at 95 °C, 40 cycles of 95 °C for 15 s, 59 °C for 30 s, and 72 °C for 45 s, followed by a

200

melting curve analysis from 60 to 95 °C for 10 minutes.

201 202

qPCR assays for φX174 ssDNA and T3 dsDNA

203

For both φX174 ssDNA and T3 dsDNA, each qPCR reaction was run in 10 µL total volume

204

consisting of 5 µL of EvaGreen qPCR Master Mix (Biotium), 0.25 µL of 25 mg/mL bovine

205

serum albumin (BSA) solution (Sigma-Aldrich), 0.4 µL of 10 µM forward primer, 0.4 µL of 10

206

µM reverse primer, 2.95 µL of nuclease-free water, and 1 µL of DNA samples. The following

207

thermocycling conditions were used for φX174 ssDNA: 5 min at 95 °C, 40 cycles of 95 °C for

208

20 s, 55 °C for 20 s, and 72 °C for 20 s, followed by a melting curve analysis from 68 to 95 °C

209

for 5 minutes. The following thermocycling conditions were used for T3 dsDNA: 2 min at 95 °C,

210

35 cycles at 95 °C for 5 s, 60 °C for 5 s, and 72 °C for 25 s, followed by a melting curve analysis

211

from 55 to 95 °C for 5 minutes.

212 213

Statistical analysis

214

Linear regressions and statistical analyses were conducted using Prism 7 (GraphPad Software,

215

Inc. La Jolla, CA) on pooled data from replicate experiments of each tested region (n ≥ 4).

216

Analysis of Covariance (ANCOVA) tests were conducted to determine whether there were

217

significant differences between the rate constants of two groups (i.e. naked genome vs

218

encapsidated genome, ssDNA vs dsDNA). When p was less than 0.05, we concluded that the rate

219

constants of two groups were significantly different.

220 221

RESULTS AND DISCUSSION

10 ACS Paragon Plus Environment

Page 10 of 27

Page 11 of 27

Environmental Science & Technology

222

Photolysis of naked nucleic acids in ultrapure water. The 500 base (b) regions of the naked

223

ssRNA and ssDNA and the 500 base pair (bp) regions of the naked dsRNA and dsDNA in

224

ultrapure water generally reacted according to first order kinetics, although tailing was observed

225

in region B of MS2 ssRNA and regions A and B of T3 dsDNA (Figure S1). Tailing has been

226

repeatedly observed when nucleic acids are exposed to UV25421-24 and a number of explanations

227

have been suggested, including the presence of aggregated particles (for viruses) and the

228

presence of UV-resistant nucleic acid sequences. We propose that the commonly observed

229

tailing effect is due to pyrimidine dimer reactions reaching photostationary state, as described in

230

early UV photolysis literature.25,26 The reason tailing was observed in region B of MS2 after the

231

same UV254 dose as region A may be due to the greater number of neighboring uracil bases in the

232

sequence (32 in region B, 19 in region A, Table S2). Further research will be necessary to

233

characterize the extent of dimer reversion in different nucleic acids with varied sequences and

234

structures.

235 236

For each region in ultrapure water, we pooled data from replicate experiments (n ≥ 7) and

237

conducted linear regressions on the data that exhibited first order kinetics (Figure 1; Table S3).

238

Control samples stored in the dark did not react over the timeframe of the experiments,

239

suggesting that the background decay was negligible (Figures S1). The first order rate constants

240

of the two measured regions in naked dsRNA were not statistically different (p = 0.15), nor were

241

the two regions of naked dsDNA (p = 0.32; Table S4). The two regions analyzed in ssRNA had

242

statistically different kinetics (p = 0.001), as did the two regions in ssDNA (p < 0.0001; Table

243

S4). This was likely due to the different number of neighboring pyrimidine bases in the two

244

regions of ssRNA and the two regions of ssDNA (Table S2).

11 ACS Paragon Plus Environment

Environmental Science & Technology

245 246

φ6 and T3 are dsRNA and dsDNA viruses, respectively, and MS2 and φX174 are ssRNA and

247

ssDNA viruses, respectively; thus, the analyzed 500 bp regions in φ6 and T3 have twice as many

248

bases as the 500 b regions of MS2 and φX174. It was therefore not possible to directly compare

249

the measured rate constants of the single-stranded regions with the double-stranded regions. By

250

comparing the ssRNA regions with the ssDNA regions, and the dsRNA regions with the dsDNA

251

regions, it was evident that the DNA regions reacted faster than the RNA regions. Specifically,

252

the two ssDNA regions reacted ~2× faster than the two ssRNA regions, and the two dsDNA

253

regions reacted ~8× faster than the two dsRNA regions. A study on the UV254 photolysis of short

254

DNA and RNA hairpin structures with HPLC also concluded that DNA is more photoreactive

255

than RNA.27

256 257

Viral nucleic acids can act as sensitizers with UV254, leading to other reactions in virus

258

particles.28 We therefore assessed the role that 1O2 and ·OH played in the UV254 photolysis of

259

naked T3 dsDNA in ultrapure water (see SI). The reaction kinetics of the naked dsDNA in D2O

260

were no different than it was in ultrapure water (Figure S2, p = 0.83), and no difference in

261

reaction kinetics was observed between dsDNA with methanol and dsDNA in ultrapure water

262

(Figure S2, p = 0.95). These results demonstrate that neither 1O2 nor ·OH was contributing

263

significantly to the decay of the naked dsDNA. We therefore concluded that the observed nucleic

264

acid reactions in ultrapure water were due to direct photolysis.

265 266

Photolysis of naked nucleic acids in PBS. For the naked ssRNA, dsRNA, and dsDNA regions,

267

the photolysis rate constants in PBS were substantially lower than in ultrapure water (Figure 1,

12 ACS Paragon Plus Environment

Page 12 of 27

Page 13 of 27

Environmental Science & Technology

268

Table S3). Specifically, the ssRNA rate constants were, on average, 1.7× higher in ultrapure

269

water than PBS, the dsRNA rate constants were, on average, 5.3× higher in ultrapure water, and

270

the dsDNA rate constants were, on average, 2.2× higher in ultrapure water. For ssDNA, only

271

region A had faster kinetics in ultrapure water based on ANCOVA analysis, and the difference in

272

this case was minor (1.1×, p = 0.04). The reason that the photoreactivities of MS2 ssRNA, φ6

273

dsRNA, and T3 dsDNA were much more impacted by solution chemistry than φX174 ssDNA is

274

not clear, but may involve the circularity of the φX174 ssDNA genome. The MS2 and T3

275

genomes are linear, and the φ6 genome is linear and segmented.

276 277

Our PBS buffer consisted of 5 mM NaH2PO4 and 10 mM NaCl (pH 7.5), whereas our

278

suspensions of nucleic acids in ultrapure water solutions were approximately pH 5.8. Pyrimidine

279

dimer quantum yields are generally not impacted by pH6,14 and early work with cytosine and

280

uracil hydrates suggest their quantum yields are relatively consistent above pH 5.29 It is therefore

281

unlikely that our observed differences in ultrapure water and PBS were due to pH. In terms of

282

ionic strength, Douki observed a 20% decrease of TT dimer quantum yields and a 300% increase

283

in TC dimer quantum yields when 20 mM NaCl was added to calf thymus DNA in pure water.6

284

In that study, the DNA dimer quantum yields were consistent when the NaCl concentrations

285

increased above 20 mM. Pyrimidine photohydrate quantum yields are also impacted by ionic

286

strength, as well as buffer type.29 These observed effects of ionic strength and buffer type on

287

quantum yield are likely related to nucleic acid structure. Hydrate formation, for example,

288

requires water molecules to access bases; this occurs more readily when double helix RNA

289

structures are denatured.29 Dimer formation requires that neighboring pyrimidines are favorably

290

aligned at the time of excitation, and this can be inhibited when structures are highly ordered.10

13 ACS Paragon Plus Environment

Environmental Science & Technology

291

Here, the nucleic acids in PBS likely had more ordered structures than the nucleic acids in

292

ultrapure water, thus inhibiting dimer and hydrate formation.

293 294

Predicted photolysis rate constants. We next predicted the theoretical photolysis rate constants,

295

𝑘!"#$ , for the eight nucleic acid regions based on a base composition method. Specifically, we

296

used the sequences in our ssRNA, dsRNA, ssDNA, and dsDNA regions and the RNA and DNA

297

extinction coefficients and quantum yields in the literature (Table S2), and summed the rate

298

constants of the major potential reactive bases in the regions:

299 𝑘!"#$,! =

2.3 𝑏 ∙ 𝜀!"# ∙ Φ 𝑈 !

!

𝑘!"#$ =

𝑘!"#$,! = !

!

2.3 𝑏 ∙ 𝜀!"# ∙ Φ 𝑈

300 301

Here, 𝑘!"#$ (cm2 mJ-1) is the predicted rate constant of the nucleic acid regions, 𝑘!"#$,! (cm2 mJ-

302

1

303

UU, and CC for RNA, and C, TT, CT, TC, and CC for DNA, 𝜀!"# (M-1 cm-1) is the extinction

304

coefficient of the nucleic acid monomer or dimer at 254 nm in M-1 cm-1, Φ (mol eins-1) is the

305

quantum yield of a particular reaction, U is a constant of the number of joules per einstein at 254

306

nm (4.7 × 105 J eins-1), and b is the number of reactive nucleic acid monomers or dimers in the

307

nucleic acid region.30 For the dsRNA and dsDNA regions, the sequences of both strands were

308

included in the 𝑘!"#$ calculation.

) is the predicted rate constant of the reactive nucleic acid monomer or dimer 𝑖, including C, U,

309

14 ACS Paragon Plus Environment

Page 14 of 27

Page 15 of 27

Environmental Science & Technology

310

The resulting predicted rate constants of the eight regions exhibited trends similar to the

311

experimental rate constants (𝑘!"# ) of the naked nucleic acid regions (Figure 1). For both the

312

experimental and predicted rate constants, the two ssDNA regions had the fastest kinetics, on

313

average, followed by the dsDNA and ssRNA regions, and the dsRNA regions were the least

314

reactive. The trends observed between the measured rate constants of the two regions in ssDNA

315

(region B > region A) and the two regions in ssRNA (regions B > region A) were also predicted

316

based on the different numbers of reactive monomers and dimers in their sequences (Table S2).

317 318

For each region, the predicted rate constant was larger than the measured values in ultrapure

319

water and in PBS. The predictions were closest for ssDNA, with the predicted rate constants

320

within a factor of 2 of the measured rate constants in both ultrapure water and PBS. The

321

predictions for dsRNA were the furthest from the measured values, with up to a 24-fold

322

difference between the predicted rate constant and the measured rate constant. The differences

323

were more pronounced when ssRNA, dsRNA, and dsDNA were in PBS than when in ultrapure

324

water (Figure 1). The discrepancies between predicted rate constants and measured rate constants

325

have several possible explanations. First, many of the quantum yield values available in the

326

literature were measured on short oligomers with repeated bases (e.g. UUUUUUU), and these

327

likely react differently than longer sequences with mixed bases. Quantum yield values available

328

in the literature were measured under a range of temperatures, solution pH, UVC wavelengths, or

329

ionic strengths that differed from our experimental conditions; these parameters can affect

330

photolysis quantum yields.31,32

331

15 ACS Paragon Plus Environment

Environmental Science & Technology

332

Quantum yields can also be impacted by the technique by which they are measured. In many of

333

the early reports, quantum yields were measured using the nucleic acid base chromosphere

334

decay, and this can suffer from interference by the formed products.14 The molecular methods

335

used here may not detect all of the reactions that take place in the RNA and DNA. Updated

336

quantum yield measurements under a variety of experimental conditions and a better

337

characterization of PCR method tolerance to nucleic acid photolysis products would help address

338

some of the observed discrepancies.

339 340

We present each 𝑘!"#$ value with the contributing reactions (𝑘!"#$,! ) identified by color (Figure

341

1; Table S2); this highlights the relative importance of hydration and dimerization reactions in

342

RNA versus DNA, as predicted with quantum yields in the literature. In both the ssRNA and

343

dsRNA regions, most of 𝑘!"#$ (~70%) is due to hydration reactions, with the rest from

344

dimerization reactions. In DNA regions, the opposite is true, with most of 𝑘!"#$ (~85%)

345

resulting from dimerization reactions. Another interesting observation from the photolysis

346

predictions is the relative contributions of TT dimers in the ssDNA and dsDNA reactions. TT

347

dimers are commonly assumed to be the most prevalent UV photoproducts in DNA, but here, the

348

sum of the predicted C hydrate contribution and TC dimer contribution is greater than the TT

349

dimer contribution. We should note that C hydrates can revert back to cytosines under some

350

conditions, and that was not included in our predictions.29 The CT, TC, and CC quantum yields

351

used in 𝑘!"#$ were based on a study in which DNA photoproducts were quantified following the

352

exposure of calf thymus DNA to UV254.5 We used the relative yields of TC, CT, or CC dimer

353

products to TT dimer products in that study, along with the well-established TT cyclobutane

354

dimer quantum yield for UV254 to estimate TC, CT, and CC quantum yields (Table S2). We were 16 ACS Paragon Plus Environment

Page 16 of 27

Page 17 of 27

Environmental Science & Technology

355

unable to locate quantum yield values for UC and CU reactions in RNA; the predicted RNA rate

356

constants would have been even higher if those reactions had been included.

357 358

Single stranded versus double stranded photoreactivity. As mentioned above, the measured

359

rate constants of the single stranded and double stranded nucleic acid regions were not directly

360

comparable due to the regions containing different numbers of bases. We therefore normalized

361

the rate constants by the number of bases in the analyzed regions (Figure 2) to observe the

362

impact of double-strand structures on RNA and DNA reactivity. In ultrapure water, RNA bases

363

reacted, on average, 7.7× faster in the single-stranded regions than in the double-stranded

364

regions; in buffer, the difference in kinetics increased to 24×. For DNA in ultrapure water, the

365

single-stranded DNA regions reacted, on average, 2.1× faster than the double-stranded genome,

366

whereas in buffer, that difference increased to 4.3×. In the case of DNA, some of the observed

367

differences were likely due to the higher proportion of adjacent TT bases in the two ssDNA

368

regions (36 and 48 TT pairs in 500 bases; Table S2) compared to the two dsDNA regions (46 and

369

47 TT pairs in 1000 bases; Table S2). The ss and ds RNA sequences, however, had similar

370

proportions of reactive bases, thus the lower photoreactivity of dsRNA compared to ssRNA was

371

likely due to the impact of the double helix.

372 373

Together, these results demonstrate that the double helix has a greater impact on RNA

374

photoreactivity than on DNA photoreactivity, and that the impact of the double helix was

375

enhanced in buffer solution compared to ultrapure water. An early study of RNA exposed to

376

UV280 found that hydrates formed 10× faster in polyU oligomers in ssRNA than in polyU:polyA

377

in dsRNA, while dimers formed 5× faster in ssRNA compared to dsRNA.15 By comparison, in

17 ACS Paragon Plus Environment

Environmental Science & Technology

378

DNA, reported TT dimer quantum yields differ by less than 2× in ssDNA and dsDNA.10,33 Our

379

results underscore the high resistance of dsRNA to UV254 photolysis compared to other nucleic

380

acids, particularly in buffered solutions. The resistance of dsRNA rotavirus to UVC inactivation

381

has been noted previously.34

382 383

Photolysis of encapsidated nucleic acids. To explore how the incorporation of nucleic acids in

384

viral particles impacts the rate of direct photolysis by UV254, we conducted experiments with the

385

same nucleic acids analyzed above, but encapsidated in virus particles and suspended in PBS

386

(Figure 1). We hypothesized that that genomes compressed inside a protein capsid would be less

387

photoreactive, due to possible decreased availability of H2O molecules for hydrate formation, the

388

restriction of motions necessary for dimer formations, and the potential for energy transfer to

389

virus proteins. In fact, for three of the four viruses tested, namely MS2 (ssRNA), φ6 (dsRNA),

390

and φX174 (ssDNA), encapsidation did not impact the reaction kinetics of the nucleic acids

391

(Table S4; p > 0.05). We observed slightly faster reaction kinetics (~1.2×) for both T3 dsDNA

392

regions when encapsidated, although the ANCOVA analysis p values were barely below our

393

0.05 cutoff for statistical significance (Table S4). These results of little or no effect from

394

encapsidation contradict an early report on the UV photolysis of tobacco mosaic virus ssRNA;

395

that research found that the ssRNA genome was 6x more sensitive to UV when outside of the

396

virus capsid than when inside the capsid.16

397 398

Literature review of nucleic acid photolysis kinetics. A number of previous studies have

399

measured the UV254 photodegradation of viral nucleic acids and extracellular DNA with qPCR or

400

RT-qPCR.24,34-43 It is normally difficult to compare these results because the regions they target

18 ACS Paragon Plus Environment

Page 18 of 27

Page 19 of 27

Environmental Science & Technology

401

have a range of sizes. We therefore normalized the literature rate constants by the number of

402

bases in the analyzed regions and compared the normalized rate constants from different studies

403

and from different genome types (Table S5). The nucleic acid UV254 photolysis rate constants in

404

the literature generally agree with one another and with our data (Figure 3). The compiled data

405

confirms that the ssDNA encapsidated genomes are the most reactive with UV254 and that the

406

encapsidated dsRNA genomes are the least reactive with UV254. The average rate constants per

407

base (mean ± 95% C.I.) for encapsidated ssRNA, dsRNA, ssDNA, and dsDNA were 3.1 × 10-5

408

± 1.0 × 10-5, 4.3 × 10-6 ± 5.8 × 10-6, 1.8 × 10-4 ± 1.9 × 10-5, and 1.7 × 10-5 ± 8.3 × 10-6 cm2 mJ-1

409

base-1, respectively. The average rate constants per base for naked dsDNA was 4.8 × 10-5 ± 2.2 ×

410

10-5 cm2 mJ-1 base-1. The scatter observed for each genome type is likely due to variations in the

411

analyzed region sequences (e.g., number of TT sequences) and differences in experimental

412

conditions between laboratories (e.g., buffer, temperature, etc.).

413 414

Environmental Implications. Our research demonstrates the relatively high resistance of

415

dsRNA compared to the other genome types. The dsRNA rate constants measured here were

416

~80× lower than the ssDNA rate constants, which exhibited the highest rate constants on a per

417

base basis. The dsRNA rate constants from this work will help with future efforts to disinfect

418

dsRNA viruses and help predict the environmental fate of dsRNA plant incorporated protectants,

419

which are increasingly applied to combat agricultural pests. Our results suggest that dsRNA

420

viruses and dsRNA plant-incorporated protectants likely persist after other forms of nucleic acids

421

have degraded from photochemical reactions.

422

19 ACS Paragon Plus Environment

Environmental Science & Technology

423

Another important conclusion of this work is the large impact that higher order structure and

424

solution chemistry plays in nucleic acid photochemistry, particularly for RNA. This complicates

425

efforts to predict the UV254 inactivation rate constants of virus genomes using published nucleic

426

acid quantum yields in combination with the genome size and sequence. The ssDNA of φX174

427

was the only genome unaffected by the solution chemistry. Future research will be needed to

428

determine if this is the case for all ssDNA viruses or for all circular genomes.

429 430

Finally, it is worth discussing the implications of employing qPCR and RT-qPCR methods to

431

quantify intact RNA and DNA, respectively. These two approaches employ different enzymes to

432

detect damage. Namely, qPCR uses polymerase to copy intact DNA regions and RT-qPCR uses

433

reverse transcriptase to convert intact RNA regions to complimentary DNA. The relative

434

tolerances of polymerase and reverse transcriptase to photolysis products is not currently known,

435

but our earlier work demonstrated that the photolysis kinetics of short RNA oligomers were

436

faster when measured with a mass spectrometry method than when measured with RT-qPCR.44

437

We selected qPCR and RT-qPCR methods for this research, as opposed to HPLC or LC-MS

438

methods, in order to obtain region-specific rate constants in the genomes; this information is lost

439

when genomes are digested and individual bases are quantified by HPLC or LC-MS.

440

Furthermore, the hydrate products that form with UV254 can undergo facile dehydration, a

441

reaction that is acid-catalyzed.8 Consequently, acid digestion steps and acidic mobile phases in

442

LC separations will affect the measured UV254 reaction rates and this is avoided with qPCR and

443

RT-qPCR. Future research should characterize the specific reactions measured by the

444

polymerases used in qPCR and the reverse transcriptases used in RT-qPCR, as well as the

445

photoproducts that impact the biological activity of RNA and DNA.

20 ACS Paragon Plus Environment

Page 20 of 27

Page 21 of 27

Environmental Science & Technology

446 447

SUPPORTING INFORMATION

448

The Supporting Information is available free of charge via the internet at http://pubs.acs.org.

449

Additional materials and methods on experiments probing the role of reactive oxygen

450

species; figures of UV254 photolysis kinetics of encapsidated and naked viral nucleic

451

acids; UV254 photolysis kinetics of naked dsDNA genomes in ultrapure water, D2O and

452

methanol; table that lists the viruses, primer sets, and the target genome regions; table

453

that describes the prediction of the UV photolysis rate constants; table with reaction rate

454

constants from linear regressions; p values obtained from ANCOVA tests; UV254

455

photolysis rate constants of encapsidated viral nucleic acids and plasmids reported in

456

previous literature.

457 458

AUTHOR INFORMATION

459

Corresponding author

460

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

461 462

ACKNOWLEDGMENT

463

This research was supported by National Science Foundation (NSF) BRIGE Award 1329576 and

464

NSF PIRE Award 1545756. We thank Professor Rose Cory from the University of Michigan for

465

her technical assistance and Kathryn Langenfeld from the University of Michigan for her

466

assistance with virus preparations.

467 468

REFERENCES:

21 ACS Paragon Plus Environment

Environmental Science & Technology

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 496 497 498 499 500 501 502 503 504 505 506 507 508 509 510 511 512 513 514

(1)

(2) (3) (4) (5) (6) (7) (8) (9) (10) (11)

(12) (13) (14) (15) (16) (17) (18)

Vikesland, P. J.; Pruden, A.; Alvarez, P. J. J.; Aga, D.; Bürgmann, H.; Li, X.-D.; Manaia, C. M.; Nambi, I.; Wigginton, K.; Zhang, T.; Zhu, Y.-G. Toward a comprehensive strategy to mitigate dissemination of environmental sources of antibiotic resistance. Environ. Sci. Technol. 2017, 51 (22), 13061–13069. Parker, K. M.; Sander, M. Environmental fate of insecticidal plant-incorporated protectants from genetically modified crops: knowledge gaps and research opportunities. Environ. Sci. Technol. 2017, 51 (21), 12049–12057. Disinfection market set for billion dollar growth up to 2019. waterworld.com. 2018. Douki, T. The variety of UV-induced pyrimidine dimeric photoproducts in DNA as shown by chromatographic quantification methods. Photochem. Photobiol. Sci. 2013, 12 (8), 1286–1302. Douki, T.; Cadet, J. Individual determination of the yield of the main UV-induced dimeric pyrimidine photoproducts in DNA suggests a high mutagenicity of CC photolesions. Biochemistry 2001, 40 (8), 2495–2501. Douki, T. Low ionic strength reduces cytosine photoreactivity in UVC-irradiated isolated DNA. Photochem. Photobiol. Sci. 2006, 5 (11), 1045–1051. Görner, H. New trends in photobiology: Photochemistry of DNA and related biomolecules: Quantum yields and consequences of photoionization. J. Photoch. Photobio. B. 1994, 26 (2), 117–139. Lemaire, D. E.; Ruzsicska, B. P. Quantum yields and secondary photoreactions of the photoproducts of dTpdT, dTpdC and dTpdU. Photochem. Photobiol. 1993, 57 (5), 755– 757. McCullagh, M.; Hariharan, M.; Lewis, F. D.; Markovitsi, D.; Douki, T.; Schatz, G. C. Conformational control of TT dimerization in DNA conjugates. A molecular dynamics study. J. Phys. Chem. B. 2010, 114 (15), 5215–5221. Schreier, W. J.; Schrader, T. E.; Koller, F. O.; Gilch, P.; Crespo-Hernández, C. E.; Swaminathan, V. N.; Carell, T.; Zinth, W.; Kohler, B. Thymine dimerization in DNA is an ultrafast photoreaction. Science 2007, 315 (5812), 625–629. Pan, Z.; Hariharan, M.; Arkin, J. D.; Jalilov, A. S.; McCullagh, M.; Schatz, G. C.; Lewis, F. D. Electron donor–acceptor interactions with flanking purines influence the efficiency of thymine photodimerization. J. Am. Chem. Soc. 2011, 133 (51), 20793– 20798. Pan, Z.; McCullagh, M.; Schatz, G. C.; Lewis, F. D. Conformational control of thymine photodimerization in purine-containing trinucleotides. J. Phys. Chem. Lett. 2011, 2 (12), 1432–1438. Law, Y. K.; Forties, R. A.; Liu, X.; Poirier, M. G.; Kohler, B. Sequence-dependent thymine dimer formation and photoreversal rates in double-stranded DNA. Photochem. Photobiol. Sci. 2013, 12 (8), 1431–1439. Fisher, G. J.; Johns, H. E. Pyrimidine photodimers; 1976; pp 225–294. Pearson, M.; Johns, H. E. Suppression of hydrate and dimer formation in ultravioletirradiated poly (A + U) relative to poly U. J. Mol. Biol. 1966, 20 (2), 215–229. McLaren, A. D.; Takahashi, W. N. Inactivation of infectious nucleic acid from tobacco mosaic virus by ultraviolet light (2537 A). Radiat. Res. 1957, 6 (5), 532. 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. Weitzman, M. D.; Carson, C. T.; Schwartz, R. A.; Lilley, C. E. Interactions of viruses 22 ACS Paragon Plus Environment

Page 22 of 27

Page 23 of 27

515 516 517 518 519 520 521 522 523 524 525 526 527 528 529 530 531 532 533 534 535 536 537 538 539 540 541 542 543 544 545 546 547 548 549 550 551 552 553 554 555 556 557 558 559 560

Environmental Science & Technology

(19)

(20) (21) (22) (23) (24) (25) (26) (27) (28) (29) (30) (31) (32) (33) (34) (35)

with the cellular DNA repair machinery. DNA Repair 2004, 3 (8-9), 1165–1173. Pecson, B. M.; Martin, L. V.; Kohn, T. Quantitative PCR for determining the infectivity of bacteriophage MS2 upon inactivation by heat, UV-B radiation, and singlet oxygen: advantages and limitations of an enzymatic treatment to reduce false-positive results. Appl. Environ. Microbiol. 2009, 75 (17), 5544–5554. 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. Beck, S. E.; Rodríguez, R. A.; Hawkins, M. A.; Hargy, T. M.; Larason, T. C.; Linden, K. G. Comparison of UV-induced inactivation and RNA damage in MS2 phage across the germicidal UV spectrum. Appl. Environ. Microbiol. 2016, 82 (5), 1468–1474. Mattle, M. J.; Kohn, T. Inactivation and tailing during UV254 disinfection of viruses: contributions of viral aggregation, light shielding within viral aggregates, and recombination. Environ. Sci. Technol. 2012, 120830135926005. McKinney, C. W.; Pruden, A. Ultraviolet disinfection of antibiotic resistant bacteria and their antibiotic resistance genes in water and wastewater. J. Am. Chem. Soc. 2012, 46 (24), 13393–13400. Chang, P. H.; Juhrend, B.; Olson, T. M.; Marrs, C. F.; Wigginton, K. R. Degradation of extracellular antibiotic resistance genes with UV254 treatment. Environ. Sci. Technol. 2017, 51 (11), 6185–6192. Helleiner, C. W.; Pearson, M. L.; Johns, H. E. The ultraviolet photochemistry of deoxyuridylyl (3′→5′) deoxyuridine. Proc. Natl. Acad. Sci. U.S.A. 1963, 50 (4), 761– 767. Johns, H. E.; Pearson, M. L.; LeBlanc, J. C.; Helleiner, C. W. The ultraviolet photochemistry of thymidylyl-(3′→5′)-thymidine. J. Mol. Biol. 1964, 9 (2), 503–IN1. Kundu, L. M.; Linne, U.; Marahiel, M.; Carell, T. RNA is more UV resistant than DNA: The formation of UV‐induced DNA lesions is strongly sequence and conformation dependent. Chem. Eur. J. 2004, 10 (22), 5697–5705. Wigginton, K. R.; Menin, L.; Sigstam, T.; Gannon, G.; Cascella, M.; Hamidane, H. B.; Tsybin, Y. O.; Waridel, P.; Kohn, T. UV radiation induces genome‐mediated, site‐ specific cleavage in viral proteins. ChemBioChem 2012, 13 (6), 837–845. Fisher, G. J.; Johns, H. E. Pyrimidine Photohydrates. In Photochemistry and Photobiology of Nucleic Acids; Elsevier, 1976; pp 169–224. Schwarzenbach, R. P.; Gschwend, P. M.; Imboden, D. M. Environmental organic Chemistry 2 nd Edition. Brown, I. H.; Johns, H. E. Photochemistry of uracil. Intersystem crossing and dimerization in aqueous solution. Photochem. Photobiol. 1968, 8 (4), 273–286. Evans, N. A.; McLaren, A. D. Quantum yields for inactivation of tobacco mosaic virus nucleic acid by ultraviolet radiation (254 nm.). J. Gen. Virol. 1969, 4 (1), 55–63. Tarmer, Z.; Wierzchowski, K. L.; Shugar, D. Influence of polynucleotide secondary structure on thymine photodimerization. Acta Biochim Pol 1696, 16 (1), 83–107. Li, D.; Gu, A. Z.; He, M.; Shi, H.-C.; Yang, W. UV inactivation and resistance of rotavirus evaluated by integrated cell culture and real-time RT-PCR assay. Water Res. 2009, 43 (13), 3261–3269. Ho, J.; Seidel, M.; Niessner, R.; Eggers, J.; Tiehm, A. Long amplicon (LA)-qPCR for the discrimination of infectious and noninfectious phix174 bacteriophages after UV inactivation. Water Res. 2016, 103, 141–148. 23 ACS Paragon Plus Environment

Environmental Science & Technology

561 562 563 564 565 566 567 568 569 570 571 572 573 574 575 576 577 578 579 580 581 582 583 584 585 586 587 588 589 590 591

(36)

(37)

(38) (39) (40) (41) (42) (43) (44)

Beck, S. E.; Rodríguez, R. A.; Linden, K. G.; Hargy, T. M.; Larason, T. C.; Wright, H. B. Wavelength dependent UV inactivation and DNA damage of adenovirus as measured by cell culture infectivity and long range quantitative PCR. Environ. Sci. Technol. 2013, 48 (1), 591–598. Calgua, B.; Carratalà, A.; Guerrero-Latorre, L.; de Abreu Corrêa, A.; Kohn, T.; Sommer, R.; Girones, R. UVC inactivation of dsDNA and ssRNA viruses in water: UV fluences and a qPCR-based approach to evaluate decay on viral infectivity. Food Environ. Virol. 2014, 6 (4), 260–268. Simonet, J.; Gantzer, C. Inactivation of poliovirus 1 and F-specific RNA phages and degradation of their genomes by UV irradiation at 254 nanometers. Appl. Environ. Microbiol. 2006, 72 (12), 7671–7677. Dancho, B. A.; Chen, H.; Kingsley, D. H. Discrimination between infectious and noninfectious human norovirus using porcine gastric mucin. Int. J. Food Microbiol. 2012, 155 (3), 222–226. Yoon, Y.; Chung, H. J.; Wen Di, D. Y.; Dodd, M. C.; Hur, H.-G.; Lee, Y. Inactivation efficiency of plasmid-encoded antibiotic resistance genes during water treatment with chlorine, UV, and UV/H2O2. Water Res. 2017, 123, 783–793. Wigginton, K. R.; Pecson, B. M.; Sigstam, T.; Bosshard, F.; Kohn, T. Virus inactivation mechanisms: impact of disinfectants on virus function and structural integrity. Environ. Sci. Technol. 2012, 46 (21), 12069–12078. Ye, Y.; Chang, P. H.; Hartert, J.; Wigginton, K. R. Reactivity of enveloped virus genome, proteins, and lipids with free chlorine and UV254. Environ. Sci. Technol. 2018, 52 (14), 7698–7708. Vazquez-Bravo, B.; Gonçalves, K.; Shisler, J. L.; Mariñas, B. J. Adenovirus replication cycle disruption from exposure to polychromatic ultraviolet irradiation. Environ. Sci. Technol. 2018, 52 (6), 3652–3659. Qiao, Z.; Wigginton, K. R. Direct and indirect photochemical reactions in viral RNA measured with RT-qPCR and mass spectrometry. Environ. Sci. Technol. 2016, 50 (24), 13371–13379.

24 ACS Paragon Plus Environment

Page 24 of 27

Page 25 of 27

Environmental Science & Technology

Naked, water

Naked, PBS

Encapsidated, PBS

C hydrates U hydrates CC dimers UU dimers TT dimers TC dimers CT dimers

ns

*

0.08

ns

0.16 ****

****

0.06 0.04 0.02 0.00

592

0.20

ns

0.12 ****

*

*

ns

**** ns

A ssRNA (MS2)

B

**** ns

**** ns

A

B

dsRNA (φ6)

0.08 0.04

A

B

ssDNA (φX174)

A

B

kpred (cm2 mJ-1)

kexp (cm2 mJ-1)

0.10

0.00

dsDNA (T3)

593

Figure 1. Measured UV254 rate constants of virus genome regions A and B in naked forms

594

(white and black bars) and in encapsidated forms (grey bars) and predicted rate constants of the

595

same genome regions (bars with patterns). Error bars represent standard error for replicate

596

experiments (n ≥ 4). Asterisks indicate rate constants are significantly different (one asterisk for

597

0.01 < p < 0.05, four asterisks for p < 0.0001) and ns indicates that rate constants are not

598

significantly different (p > 0.05). The ANCOVA analysis p values are presented in Table S4.

599

25 ACS Paragon Plus Environment

kexp (cm2 mJ-1 base-1)

Environmental Science & Technology

Naked, water

1.8×10-4

Naked, PBS

Page 26 of 27

Encapsidated, PBS

1.2×10-4

6.0×10-5

0.0

A

B

ssRNA (MS2)

600

A

B

A

dsRNA (φ6)

B

ssDNA (φX174)

A

B

dsDNA (T3)

601

Figure 2. Experimental first order rate constants of target viral genome regions normalized by

602

the number of bases in corresponding region. Error bars represent standard error.

k (cm2 mJ-1 base-1)

10-3

10-4

10-5

603

Encapsidated

ds D N A

ds D N A

A N D ss

ds R N A

ss

R

N

A

10-6

Naked

26 ACS Paragon Plus Environment

Page 27 of 27

Environmental Science & Technology

604

Figure 3. Comparison of the first order rate constants measured in this study on two regions of

605

each genome with first order rate constants reported in the literature normalized by the number

606

of bases in the PCR amplicon. Black points represent data from this study conducted in PBS,

607

gray points represent data from this study conducted in ultrapure water, and white points

608

represent data from the literature.

609 610

Table 1. Four model viruses used in this study, including details of the genome regions analyzed. Region Location

Region size (bases or base pairs)

Region A Region B

944-1439 2693-3189

496 b 497 b

Enveloped

Region A Region B

S1141-S1639 L1510-L1993

499 bp 484 bp

Nonenveloped

Region A Region B

571-1074 1717-2209

504 b 493 b

Family/Genus*

Genome type

MS2

Leviviridae/ Levivirus

ssRNA (linear)

3.6 kb

~25

Nonenveloped

φ6

Cystoviridae/ Cystovirus

dsRNA (segmented)

2.9 kbp (S) 4.1 kbp (M) 6.4 kbp (L)

~85

5.4 kb

~25

Virus

φX174

Microviridae/ ssDNA Phix174microvirus (circular)

Genome Particle size Enveloped size (nm) /nonenveloped

Podoviridae/ dsDNA ~50 × 20 Region A 1678-2186 509 bp 38.2 kbp Nonenveloped T7virus (linear) (tailed) Region B 11826-12324 499 bp * Family and genus were based on the International Committee on Taxonomy of Viruses (ICTV), https://talk.ictvonline.org/taxonomy/ T3

611

Regions Analyzed

27 ACS Paragon Plus Environment