Induction of Microbial Oxidative Stress as a New Strategy to

Jul 29, 2019 - Center for Urban Waters, Tacoma, Washington 98421, United States ..... fitting the data using Solver in MS Excel. ..... regulatory syst...
2 downloads 0 Views 1MB Size
Subscriber access provided by BUFFALO STATE

Environmental Processes

Induction of microbial oxidative stress as a new strategy to enhance the enzymatic degradation of organic micropollutants in synthetic wastewater Amrita Bains, Octavio Perez-Garcia, Gavin Lear, David R. Greenwood, Simon Swift, Martin J Middleditch, Edward P. Kolodziej, and Naresh Singhal Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.9b02219 • Publication Date (Web): 29 Jul 2019 Downloaded from pubs.acs.org on July 31, 2019

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 26

Environmental Science & Technology

1

Induction of microbial oxidative stress as a new strategy to enhance the enzymatic

2

degradation of organic micropollutants in synthetic wastewater

3 4

Amrita Bains †, Octavio Perez-Garcia†, Gavin Lear‡, David Greenwood‡, Simon Swift§, Martin

5

Middleditch‡, Edward P. Kolodziej┴♯ꜞ, Naresh Singhal*†

6 7

†Department

8

‡School

9

§Faculty

of Civil and Environmental Engineering, University of Auckland, New Zealand

of Biological Sciences, University of Auckland, New Zealand of Medical and Health Sciences, University of Auckland, New Zealand

10

┴Division

11

United States

12

♯Department

13

United States

14

ꜞCenter

of Sciences and Mathematics, University of Washington-Tacoma, Washington, of Civil & Environmental Engineering, University of Washington, Washington,

for Urban Waters, Tacoma, Washington, United States.

15 16

ABSTRACT:

17

Organic micropollutants (OMPs) are pervasive anthropogenic contaminants of receiving

18

waters where they can induce various adverse effects to aquatic life. Their ubiquitous

19

environmental occurrence is primarily attributed to discharge from wastewater treatment plants

20

due to incomplete removal by common biological wastewater treatment processes. Here, we

21

assess a new strategy for promoting the degradation of six representative OMPs (i.e.

22

sulfamethoxazole, carbamazepine, tylosin, atrazine, naproxen and ibuprofen) by intentionally

23

stimulating the production of microbial oxidoreductases to counter oxidative stress caused by

24

oxygen perturbations. Mixed microbial cultures from a dairy farm wastewater were subjected

25

to cyclic perturbations of dissolved oxygen (DO). A distance-based redundancy analysis was

26

used to show that DO perturbations correlate with the abundance of Pseudomonadaceae and

27

Rhodocyclaceae families, activities of peroxidases and cytochromes, and the degradation of

28

OMPs. DO perturbation of 0.25 and 0.5 cycles/hr led to most abundance of Pseudomonadaceae

29

and Rhodocyclaceae families, showed higher activity of peroxidase and cytochrome, and gave

30

largest removal of OMPs (removal of 92±3% for sulfamethoxazole, 84±3% for naproxen,

31

82±3% for ibuprofen, 66±2% for carbamazepine, 57±15% for tylosin, and 88±1% for atrazine).

32

1 ACS Paragon Plus Environment

Environmental Science & Technology

33

INTRODUCTION

34

The discharge of trace levels of pharmaceuticals, industrial chemicals, pesticides and personal

35

care products into aquatic environments is of global environmental concern.1 The resistance of

36

such organic micropollutants (OMPs) to conventional wastewater treatment and their

37

environmental persistence risks adverse ecological impacts to biota via direct toxicity and

38

endocrine and metabolic disruption2, DNA fragmentation and immune deficiency.3–5 For

39

example, exposure to the potent synthetic estrogen 17α-ethinylestradiol and the antidepressant

40

fluoxetine even at concentrations below 5 ng/L and 28 ng/L, respectively, causes reproductive

41

endocrine malfunction (e.g., production of vitellogenin mRNA and protein and disruption of

42

gonad development) in fathead minnow (Pimephales promelas) males.6,7 OMP exposure

43

effects may be additive; for example, exposure of female zebrafish (Danio rerio) to

44

pharmaceutical mixtures of acetaminophen, carbamazepine, gemfibrozil and venlafaxine

45

significantly reduces fecundity.8 Such concerns about the occurrence of OMPs have triggered

46

attempts to enhance their biodegradation during municipal wastewater treatment.9 Although

47

conventional activated sludge processes are not specifically optimised to eliminate OMPs 10,

48

exploiting the diverse enzymatic potential of microbial consortia has been identified as a

49

critical pathway for enhancing OMP removal.11

50

In the period since the first oxidative conditions arose on Earth due to early photosynthesis,

51

microbes have evolved to survive exposure to a variety of harmful oxidative stresses.12 For

52

instance, oxidative stress generated by high concentrations of intracellular reactive oxygen

53

species (ROS), including superoxide (O2-), hydrogen peroxide (H2O2) and hydroxyl radicals

54

(OH.), induces oxidative enzyme gene expression 13 and consequent synthesis of antioxidative

55

enzymes such as oxidoreductases (peroxidases and cytochromes) to protect against such

56

oxidative stress.14 From an engineered wastewater treatment perspective, intentional

57

manipulations of DO concentrations may be able to induce the production of such enzymes,

58

some of which are also capable of degrading pollutants like polycyclic aromatic hydrocarbons

59

and organophosphorus contaminants.15–17 Therefore, DO perturbations might pre-dispose cells,

60

through transcriptional responses of their oxidoreductases to protect against iterative oxidative

61

stress, to enhance OMP degradation. This proposed mechanism implies that simple oxygen

62

control strategies within wastewater treatment bioreactors may provide a cost-effective

63

solution for enhanced OMP removal.

64

In this study, we investigate the impact of dynamic variation in DO concentration on OMP

65

biotransformation by representative microbial consortia. We hypothesise that varying DO 2 ACS Paragon Plus Environment

Page 2 of 26

Page 3 of 26

Environmental Science & Technology

66

concentrations alter the intracellular metabolic pathways to disrupt the prevailing equilibrium

67

in electron flow, leading to enhanced ROS formation and thereby upregulating ROS

68

scavenging and defence enzymes. We evaluate OMP removal under three oxygenation patterns

69

– constant aeration, continuous cyclic aeration, and intermittent cyclic aeration for both high

70

and low (~10 and ~4 mg-DO/L, respectively) aerobic conditions. Our objective is to assess

71

whether cyclic patterns enhance OMP degradation by inducing oxidative stress and forcing

72

microbial oxidoreductase production (superoxide dismutase (SOD), catalases, peroxidases

73

etc.).

74

MATERIALS AND METHODS

75 76

Reactor operation and sampling Fed-batch experiments were performed using bioreactors inoculated with a mixed microbial

77

culture (sludge from a dairy wastewater settling pond, Alfriston, New Zealand). This inoculum

78

represents a complex source of microbial communities potentially exposed to antibiotics and

79

different carbon sources typical of intensive animal agricultural operations. The sludge was

80

diluted to 3 g/L volatile suspended solids (VSS) with distilled water. A 10x concentrated stock

81

synthetic wastewater with the following composition was used: 4000 mg-COD/L acetate, 1000

82

mg-COD/L methanol, and 600 mg-N/L ammonium along with 1 mg/L of each OMP

83

(sulfamethoxazole, carbamazepine, tylosin, atrazine, naproxen and ibuprofen) (detailed in

84

supplementary information). Each 1 L bioreactor was continuously fed with this synthetic

85

wastewater at 0.034 mL/min over 48 hours using a low flow syringe pump (KDS Legato 210-

86

788210 Syringe Pumps, KD Scientific, USA). Following other studies2,18, a higher than

87

environmentally realistic nominal OMP concentration of 100 µg/L was used in the bioreactors

88

to more easily monitor their concentration changes during the study and account for the

89

relatively limited selectivity of the LC/MS detector.

90

Oxidative stress was imposed on cultures by iteratively exposing them to ON-OFF cycles of

91

an oxygen (20%) enriched air supply during the 48h cultivation period. (Table 1, Figure 1). Air

92

delivery to reactors was controlled by a flow meter in combination with Millenium 3 CD20

93

logic controller (Crouzet, New Zealand) connected to solenoid valve operating on compressed

94

air lines. DO concentrations in reactors were continuously measured using Visiferm DO ARC

95

425 probes (Hamilton Bonaduz AG, Switzerland) connected to a PC using Hamilton Device

96

Manager 1.0.0 software (Hamilton Bonaduz AG, Switzerland). Additionally, constant (non-

97

perturbed) high (10 mg-DO/L) and low (4 mg-DO/L) DO conditions and autoclaved (dead)

98

biomass were evaluated as controls under constant aeration at 10 mg-DO/L. 3 ACS Paragon Plus Environment

Environmental Science & Technology

99

Samples were collected at 1, 24 and 48 hours, with sampling at the end of any perturbation

100

cycles. The samples collected at 48 hours were analysed for enzyme activity, proteins,

101

microbial speciation, and OMP concentrations, while the samples collected at 1 and 24 hours

102

were analysed for enzyme activity (Table S1). Each analysis was performed in duplicate (two

103

biological and two analytical replicates). Additional experimental details are described in the

104

supplementary information.

105 106

Enzymes activity assays

107

Oxidoreductase activity in crude biomass extracts was determined spectrophotometrically by

108

measuring the degradation (oxidation and hydrolysis) of various chromogenic substrates used

109

as surrogate xenobiotics (Table S2). The oxidoreductase (lignin peroxidase, horseradish

110

peroxidase, laccase and cytochrome P450) and hydrolase (-glucosaminidase, -glucosidase)

111

activities in bioreactors were detected spectrophotometrically by the oxidation of their

112

respective chromogenic substrates (Methylene Blue [MB], Azure B [AB], 3,4-Dihydroxy-L-

113

phenylalanine [L-DOPA], 2,2’-Azino-bis-(3-ethylbenzothiazoline-6-sulfonic acid)[ABTS],

114

Sudan Orange G [SO], 4-nitrophenyl-dodecanoate [pNP-12], Indole and 4-Aminoantipyrine

115

[4-AAP], 4-nitrophenyl N-acetyl-β-D-glucosaminide [pNP-A] and 4-nitrophenyl- β-D-

116

glucopyranoside [pNP-G]

117

selected as they are synthesised during intracellular oxidative stress (generated by cyclic

118

oxygen shocks here) and their ability to degrade OMPs22. Oxidoreductases catalyse OMP

119

biotransformation by both oxidation and reduction reactions22–24. Except for Cyt P450 activity,

120

estimated using three substrates, and hydrolases β-glu and β-glcNAc, estimated using one

121

substrate, all enzyme activities were estimated using two chromogenic substrates. Standard

122

normal distribution (N(0,1) was used to obtain the autoscaled enzyme activities (Figure 2a). In

123

brief, 20 ml aliquots of microbial cultures samples were centrifuged at 20,000 x g for 15 min.

124

at 4˚C in 50 ml falcon tubes. The pellets from each tube were resuspended in either 10 ml

125

acetate buffer or 10 ml phosphate buffer (detailed in supporting information). Each well of a

126

96 well microplate was filled with 100 µl aliquots of the respective buffers, 100 µl chromogenic

127

dye and 100 µl resuspended biomass. Dyes and samples resuspended in double the amount of

128

buffer served as controls. The microplates were mixed and incubated at 30˚C for 1 hour 10 µl

129

of 0.3% H2O2 at 30% as primary oxidant was added to start the reactions for lignin peroxidase

130

(AB, MB) and horseradish peroxidase (ABTS, L-DOPA) enzymes25,26 and 10 µl of 1M NaOH

131

was added to quench the para-nitrophenol dye reactions19 (i.e. -glucosaminidase (pNP-A), -

19–21

(Table S2). These oxidoreductases and hydrolases were

4 ACS Paragon Plus Environment

Page 4 of 26

Page 5 of 26

Environmental Science & Technology

132

glucosidase (pNP-G) and cytochrome P450 (pNP-12)) only. Changes in absorbance caused by

133

chromogenic reactions were quantified using a Victor X3 Multimode Plate Reader

134

(PerkinElmer, USA) at different wavelengths (refer to supplementary information).

135

Protein Analysis

136

Proteomic analysis was performed to screen for oxidative stress associated enzymes as well as

137

non-ROS enzymes. Protein extraction was performed by cell lysis followed by precipitation of

138

proteins and separation of low and high stringency protein fractions by 1D SDS PAGE. Briefly,

139

30 ml of the DO treated sample was washed with 50 ml of 0.9% sodium chloride and spun

140

down at 20,000g (20 min at -4°C). The pellet was rewashed in 40 ml Tris-HCl (pH 7) and

141

centrifuged at 20,000g (20 min at -4°C) to get a another pellet, which was resuspended in

142

sample buffer, pulse-vortexed and placed on ice for 2 h with regular mixing at 15 min intervals.

143

Samples were sonicated for 15 s (six rounds on ice) and centrifuged at 20,000g (3 min at 4°C).

144

The supernatant was mixed with trichloroacetic acid (100% w/v) and centrifuged at 23,400g

145

(30 min at -4°C) after which the pellet was washed with 5 ml cold acetone twice. The final

146

pellet was heat dried and re-suspended in 400 µl low stringency buffer (see supplementary

147

information) to generate a low stringency fraction (LSF) in the supernatant after,

148

centrifugation, then extracting the pellet in 400 µl high stringency buffer to generate a soluble

149

high stringency fraction (HSF) after recentrifugation.

150

Proteins in both fractions were quantified spectrophotometrically by fluorescence. Crude

151

separation of LSF and HSF proteins was achieved using 1D SDS PAGE. Individual bands were

152

excised with a sharp razor blade and placed into low-binding, siliconized microcentrifuge tubes

153

for destaining reduction, alkylation and trypsinolysis. Finally, 10 µL of the generated tryptic

154

peptides of the microbial proteins were injected onto a SCIEX 6600 tripleTOF mass

155

spectrometer. Protein identification was carried out by comparing the obtained peptide

156

sequences against those of the UniProt database. Additional details can be found in the

157

Supplementary Information.

158

Organic micropollutant analysis

159

The individual 1L bioreactors were continuously fed with OMPs mixture (sulfamethoxazole,

160

carbamazepine, tylosin, atrazine, naproxen and ibuprofen) at a feed rate of 0.034 mL/min to

161

the bioreactors for 48 h. Residual OMPs were extracted using solid phase extraction (SPE) on

162

OASIS HLB cartridges (Waters, Milford, MA, USA) following the described method.27

163

Briefly, 200 mL of aqueous bioreactor culture was collected at the end of the cultivation cycle 5 ACS Paragon Plus Environment

Environmental Science & Technology

Page 6 of 26

164

and centrifuged to remove suspended particles (20,400 g, 20 min, -4°C). The supernatant was

165

passed through the cartridges (pre-conditioned with 5 mL each of tert-methyl butyl ether,

166

methanol (Merck; laboratory grade purity) and deionised water) and washed with water. The

167

OMPs were eluted with methanol (10 mL) under vacuum.

168

OMP analysis was performed by liquid chromatography coupled with mass spectrometry (LC-

169

MS) using a Shimadzu 2020 Series LC-MS (Shimadzu, Japan) with an Agilent ZORBAX

170

Eclipse Plus C18 column (2.1 mm × 100 mm, particle size 1.8 µm, Agilent Technologies,

171

Germany) following EPA Method 1694.28 A binary gradient system of mobile phase A, 0.1%

172

formic acid in deionised water, and mobile phase B, 0.1% formic acid in acetonitrile were used

173

to separate analytes in ESI+ mode, while 5 mM ammonium acetate, pH-5.5 (mobile phase A)

174

and methanol (mobile phase B) were used for analysis in ESI- mode. The solvent gradient for

175

ESI+ mode was 10% B to 60% B in 24 min and then maintained at 100% B to 30 min. For

176

ESI- mode, the gradient started at 40% B and was increased linearly to 100% B over 13.5 min,

177

and then maintained at 40% B to 17 min. The flow rate was 0.2 ml/min. in both modes and the

178

injection volume was set to 2 µl and 10 µl for +/- polarity modes, respectively. Limits of

179

detection and quantification were determined using signal/noise ratios of 3 and 10,

180

respectively. OMP removal efficiencies were calculated as [(OMP0 – OMP48)/OMP0] *100,

181

where OMP0 and OMP48 are the OMP concentrations at the start of the study (0 hour) and after

182

48 hours. The quality assurance and quality control were checked within each measurement

183

series by spike recovery experiments (spiking at 100 µg/L, n ≥ 3) and repeated injections of

184

matrix recovery samples. Recovery of the matrix spikes ranged from 54 to 100% and was

185

consistent between analytes. Replicate samples showed a maximum range of relative standard

186

deviation (RSD%) within 15%. OMPs were never detected in distilled water blanks.

187

Estimating the removal efficiency and degradation rate constant

188

The OMP removal efficiencies (%) for fed-batch mode involving a continuous addition of

189

OMPs to bioreactors, were calculated as follows:

190

 48hr    Qin Cin dt  VCresidual   x100 % Removal =  0 48hr   Q C dt   in in 0  

191

Where, Qin is the flow rate and Cin is the OMP concentration in the feed from syringe pump, V

192

is the volume of synthetic media, and Cresidual is the OMP concentration after 48 hours of reactor

(Eq. 1)

6 ACS Paragon Plus Environment

Page 7 of 26

Environmental Science & Technology

193

operation. The first-order degradation rate constants were estimated for the fed batch mode

194

using the following equation, obtained by applying the mass balance principle.

195

Ct 

196

where, Ct is the OMP concentration in reactor after time t and k is the first-order degradation

197

rate constant (/hr).

198

The estimates for first-order rate constants were obtained by fitting the data using Solver in MS

199

Excel. These values were then used to predict the decay that would have occurred in batch

200

mode, i.e., for a spike addition of OMPs to bioreactors, using the following equation.

201

Qin Cin 1  e  k*t kV





(Eq. 2)

C t  C0 e  k*t

(Eq. 3)

202

where, C0 is the initial OMP concentration (100 µg/L).

203

Microbial DNA isolation and bacterial species identification

204

A PowerSoil DNA isolation kit (MoBio, Carlsbad, USA) was used for the isolation of bacterial

205

total genomic DNA extracted from samples (1 mL) following the manufacturer’s protocol.

206

Bacterial community composition was characterised by amplifying and sequencing bacterial

207

16S ribosomal RNA (rRNA) gene fragments with the universal 16S PCR forward primer (5'-

208

TCGTCGGCAGCGTCAGATGTGTATAAGAGACAGCCTACGGGNGGCWGCAG-3')

209

and

210

GTCTCGTGGGCTCGGAGATGTGTATAAGAGACAGGACTACHVGGGTATCTAA

211

TCC-3') following standard protocols29; nucleotide bases in bold are Illumina overhang adapter

212

sequences for high throughput sequencing. Amplified PCR products were purified using an

213

AMPure XP beads kit (Beckman Coulter Inc.) and concentrations recorded using a Qubit®

214

dsDNA HS Assay Kit (Life Technologies) and sequenced using an Illumina MiSeq machine

215

(New Zealand Genomics Ltd., Auckland, New Zealand) using 2-by-300 bp chemistry. Prior to

216

DNA sequencing, the sequencing provider attached a unique combination of Nextera XT dual

217

indices (Illumina Inc., USA) to the DNA from each sample to allow multiplex sequencing. The

218

resulting paired-end read DNA sequence data were merged and quality filtered using the

219

USEARCH sequence analysis tool.30 Data were dereplicated so that only one copy of each

220

sequence was reported. Sequence data were then checked for chimeric sequences and clustered

221

into groups of operational taxonomic units based on a sequence identity threshold equal to or

222

greater than 97% (thereafter referred to as 97% OTUs) using the clustering pipeline UPARSE

223

30

16S

PCR

reverse

primer

(5'-

in QIIME v.1.6.0, as described.31 Prokaryote phylotypes were classified to their 7 ACS Paragon Plus Environment

Environmental Science & Technology

224

corresponding taxonomy by implementing the RDP classifier routine 32 in QIIME v. 1.6.0 33 to

225

interrogate the Greengenes 13˙8 database.34 All sequences of chloroplast and mitochondrial

226

DNA were removed. Finally, DNA sequence data were rarefied to a depth of 5,600 randomly

227

selected reads per sample and two samples per treatment to achieve a standard number of

228

sequencing reads across all samples.

229

Statistical analyses

230

Enzyme activities were plotted using the heat map function within the R package ‘gplots’.

231

Similarities in bacterial community composition, enzyme activity and OMP remaining in

232

cultures were investigated using multivariate statistical analyses. Bray-Curtis distance matrices

233

of relative bacterial abundance, enzyme synthesis and OMP residual concentrations were

234

calculated and significant differences across all sample groups assessed using permutational

235

multivariate analysis of variance (PERMANOVA). The RELATE function within the

236

PRIMER package was used to calculate the Spearman rank correlation between data matrices

237

constructed from the bacterial community, enzyme activity and OMP residual concentrations.

238

All multivariate statistical analyses were performed in PRIMER 6 (Plymouth Marine

239

Laboratory). Significant differences in OMP removal and residual concentrations were

240

calculated by GraphPad Prism 7 using one-way ANOVA and post-hoc Tukey tests.

241 242

RESULTS

243

Influence of oxygen perturbations on microbial enzymatic activity

244

The cyclic DO exposure to representative mixed microbial communities induces the activation

245

of different carbon metabolic pathways and oxidative stress. The shift in the intracellular

246

metabolism and ROS generation results in the synthesis of oxidoreductases and other enzymes

247

essential for metabolic regulation. To test this hypothesis, five frequencies (2, 1, 0.5, 0.25 and

248

0.16 cycles per hour) and constant controls (0 cycles/hr) were tested under two aeration regimes

249

(high-DO, up to 10 mg/L DO and low-DO, up to 4 mg/L DO) that generated cyclic patterns

250

(both continuous and intermittent) of DO concentrations in bioreactors (Figure 1, Table 1). The

251

activities of oxidoreductases (lignin peroxidase, horseradish peroxidase, laccase and

252

cytochrome P450) and hydrolases (-glucosidase and -glucosaminidase) in bioreactors were

253

determined spectrophotometrically.

254

Cultures perturbed with 0.25 and 0.5 cycles/hr at both low and high-DO regimes showed

255

significantly higher (ANOVA, p < 0.05) cytochrome P450 (Indole, 4-AAP) and horseradish

256

peroxidase (L-DOPA, ABTS) compared to constant DO non-perturbed controls. The 8 ACS Paragon Plus Environment

Page 8 of 26

Page 9 of 26

Environmental Science & Technology

257

horseradish peroxidase (L-DOPA, ABTS) activity remained significantly higher in 0.16

258

cycles/hr DO perturbed samples but decreased in cultures perturbed with 1 or 2 cycle/hr at

259

high-DO aeration where -glucosidase (pNP-G) activity was high. This indicated that the DO

260

transitions activate the enzymes associated with regulation of carbon metabolism and the

261

oxidative stress caused by ROS formation. Lignin peroxidase (AB, MB) activity was evident

262

in all DO conditions (perturbed and non-perturbed). By contrast, under constant non-perturbed

263

oxygen concentrations, the high-DO aerobic cultures showed greater (ANOVA, p < 0.05)

264

cytochrome P450 (pNP-12), -glucosaminidase (pNP-A), -glucosidase (pNP-G) activities

265

while the low-DO aerobic samples exhibited greater laccase (SO) and horseradish peroxidase

266

(ABTS) activity (Figure 2a).

267

The influence of varying DO conditions on protein abundance in biomass samples was

268

illustrated by proteomic analysis. Proteins were resolved for perturbed (0.25 and 0.5 cycles/hr

269

at low-DO regime) and non-perturbed (constant high and low-DO aerobic regime; Figure 2b).

270

Overall, 436 targeted proteins (>1.3 unused score and 95% confidence interval) were identified,

271

of which 245 proteins were detected in perturbed (0.25 and 0.5 cycle/hr low-DO aerobic)

272

cultures and 191 proteins in non-perturbed (constant high and low-DO aerobic) cultures.

273

Relatively higher abundance of expressed cytochrome c, cbb3 cytochrome oxidase,

274

flavocytochrome c, NADH-quinone oxidoreductase and putative oxidoreductase were

275

observed in cultures perturbed with 0.25 and 0.5 cycles/hr at low-DO aerobic regime. This

276

observation demonstrates that exposing microbes to transitioning aerobic-microaerobic-anoxic

277

conditions will induce the activation of different cellular central metabolic pathways leading

278

to differential enzyme synthesis35. The non-perturbed constant high and low-DO aerobic

279

cultures showed higher abundances for superoxide dismutase (SOD) and putative

280

oxidoreductases, synthesised by the intracellular ROS formation (Figure 2b). SOD, as the

281

primary ROS scavenging enzyme indicates the first line of defence against oxidative stress

282

caused by ROS formation14,36. Table S3 and S4 provides the autoscaled enzyme activities and

283

proteins relative abundance at 48 hours under different DO perturbed and non-perturbed

284

conditions, respectively. Figure S1 shows the change in the enzyme activities in DO perturbed

285

and non-perturbed cultures at 1 and 24 hours.

286

Micropollutant biotransformation by mixed cultures

287

Minimal (3±3%) OMP removal was observed in autoclaved controls (Figure S1) and is

288

attributed to abiotic factors such as hydrophobic partitioning and sorption to solids. In

289

biologically active reactors, we reasoned that synthesised oxidoreductases might enhance OMP 9 ACS Paragon Plus Environment

Environmental Science & Technology

290

biotransformation and increase removals. To further validate this hypothesis, we observed that

291

OMP removal (at 100 µg/L initial nominal concentrations) was significantly (ANOVA test at

292

p

307

0.05 by post-hoc Tukey test) for sulfamethoxazole, carbamazepine and tylosin between DO

308

perturbed and non-perturbed constant high and low-DO aerobic cultures. Tables S5, S6 and S7

309

present OMP removal efficiencies, standard deviation using pair of two replicates and quality

310

control measures, respectively. The estimated first-order rate constants for the observed OMP

311

removals are presented in Table S8 and the corresponding half-lives of OMP in DO perturbed,

312

non-perturbed, autoclaved biomass and wastewater treatment systems are presented in Table

313

S9. Figure S2 shows predictions for OMP removals in reactors operated in batch mode. Figure

314

S3 represents the OMP residual concentrations after 48 hours in the bioreactors operated in

315

fed-batch mode.

316

Influence of oxidative stress on microbial speciation

317

Next, we assessed the effect of oxidative stress on the microbial community composition and

318

whether bacterial species presence correlated with oxidoreductases synthesis and OMP

319

biotransformation. A total of 527 operational taxonomic units (97% OTUs) were identified in

320

biomass samples, of which 513 represented Bacteria and 13 were Archaea. The occurrence of

321

bacteria belonging to the phylum Proteobacteria, which are linked to the biotransformation of

322

recalcitrant contaminants in wastewaters, 37, predominated and were correlated to ecological 10 ACS Paragon Plus Environment

Page 10 of 26

Page 11 of 26

Environmental Science & Technology

323

variables like process operations (changing oxygen conditions), influent wastewater

324

characteristics or geographical location.38 Comamonadaceae (40%) was the dominant family

325

in cultures perturbed with 0.25 cycles/hr in high-DO regime while a higher relative abundance

326

of Methylophilaceae (25%) was found in 1 and 2 cycles/hr treated cultures at high-DO range

327

(Figure S2). Cultures exposed to 0.16 cycles/hr in low-DO regime showed a higher relative

328

abundance of Rhodocyclaceae (15%). Flavobacteriaceae (10-12%) and Pseudomonadaceae

329

(14%) were dominant in cultures perturbed with 0.25 and 0.5 cycles/hr at low-DO regime. The

330

non-perturbed constant high and low-DO regime subjected cultures resulted in higher

331

abundances (30-45%) for Comamonadaceae (Figure S4). A dbRDA constrained ordination

332

based on Bray-Curtis similarity was conducted to determine the extent to which bacterial

333

community variations at family level are explained by the DO conditions (Figure 4). A

334

relationship was confirmed between dependent variables (bacterial families, enzyme activity

335

and residual OMP concentration) and changes in predictor variables (perturbed and non-

336

perturbed DO conditions). This indicates that DO modulations can cause physiological and

337

biochemical changes, greatly affecting bacterial community composition, enzymatic activities

338

and OMP biotransformations. The DO treatments (perturbed and non-perturbed)

339

[permutational multivariate ANOVA (PERMANOVA): F = 3.26; df = 10, 11; p < 0.001] had

340

significant effects on the bacterial community composition (Figure 4). The variation in

341

bacterial

342

Sphingobacteriaceae, Methylophilaceae, Xanthomonadaceae, Campylobacteriacea and

343

Aeromonadaceae), expressed enzyme activities (laccase, lignin, horseradish peroxidase and

344

cytochrome P450) and OMPs (sulfamethoxazole, carbamazepine, naproxen and ibuprofen)

345

removal were assessed in both perturbed and constant non-perturbed DO conditions (Figure 4).

346

These data indicate that microbial species that from oxidative stress showed more enzyme

347

activities and OMP biotransformation. Variation in the composition of bacterial communities

348

was positively correlated with microbial enzyme activity (relate analysis, R = 0.46, p < 0.05)

349

and residual concentrations of OMPs (relate analysis, R = 0.54, p < 0.05). Pseudomonas39 and

350

Rhodococcus40

351

Xanthomonadaceae42 synthesise laccases, which can catalytically degrade atrazine43,

352

carbamazepine44, estradiols45 and other OMPs22.

families

(Rhodocyclaceae,

synthesise

cytochrome

Syntrophaceae,

P450,

whereas

Syntrophobacteraceae

Methylophilaceae41

and

353 354

DISCUSSION

355

Our results show that DO perturbations can rapidly alter the functional performance and

356

composition of bacterial communities. These findings are consistent with other studies that 11 ACS Paragon Plus Environment

Environmental Science & Technology

357

exposed microbial communities to shock treatments,46, and reflect the toxic effects of ROS47–

358

49

359

ROS formation and control are illustrated in Figure 5. The arcA and fnr genes play a significant

360

role in regulating the central carbon metabolism51 and are only activated under micro-aerobic

361

or anoxic52–54 conditions. Oxygen modulation between aerobic and anoxic conditions seems to

362

reprogram central carbon metabolism,35, disrupt intracellular electron flow55 and generate O2•-,

363

one the primary ROS molecules within cells.

364

SoxRS genes to produce superoxide dismutase (SOD), which converts O2•- to H2O2. The

365

formation of H2O2 activates OxyR genes, which produce catalases and peroxidases to convert

366

H2O2 to H2O.13 Enzymes such as catalase and peroxidases form activated enzyme-oxygen

367

complexes during H2O2 attenuation that are not completely selective to H2O2 and incidentally

368

degrade contaminants like OMPs reduction of H2O2 to H2O.56 Peroxidases, laccases and

369

cytochromes have been reported to hydroxylate and oxidise aromatic organic substrates via

370

one electron oxidations.22,23,24 In addition, the metal centres of the peroxidases and laccases

371

have reactive affinities for various functional groups in OMPs.57

372

Microbes that are better able to cope with dynamic changes to central carbon metabolism and

373

at scavenging toxic ROS species are preferentially able to survive and grow under oxygen

374

modulation47. Constant aeration, continuous cyclic DO perturbation, and intermittent cyclic

375

DO perturbation will challenge microbial survival in different ways: under constant aeration,

376

high and low DO levels will lead to different rates of carbon metabolism; for continuous cyclic

377

DO perturbations, the frequency will affect the rate of transitioning between aerobic, micro-

378

aerobic and anoxic regimes, thereby inducing shifts in carbon metabolism and altering ROS

379

formation.The introduction of anoxic phase in intermittent cyclic modulations alters carbon

380

metabolism (e.g., potentially inactivating the TCA cycle) and potentially activates anaerobic

381

pathways (e.g. nitrate reduction) along with impacts to ROS formation55. The most common

382

bacterial families surviving the ROS stress caused by constant aeration, continuous cyclic DO

383

perturbation, and intermittent cyclic DO perturbation belong to the proteobacterial phylum

384

(Pseudomonadacaeae, Xanthomonadacaeae, Rhodocyclacea, Comamonadaceace) concordant

385

with previous studies58,59.

386

DO perturbations led to the production of a more diverse protein pool. The significant

387

formation of superoxide dismutase only under constant aeration suggests the activation of Sox

388

genes due to the formation of more O2•-. The production of cytochrome c, cbb3 cytochrome

389

oxidase and NADH-quinone oxidoreductase under cyclic DO perturbations is indicative of a

along with their role as signalling molecules50. Key regulatory mechanisms of intracellular

47

Microbes control O2•- levels by activating

12 ACS Paragon Plus Environment

Page 12 of 26

Page 13 of 26

Environmental Science & Technology

390

change in carbon metabolism and the activation of FNR and ArcAB regulatory systems.51 The

391

formation of putative oxidoreductases occurred under all conditions, but their concentrations

392

were higher for cyclic DO perturbations, with the intermittent cyclic DO perturbation showing

393

the largest enzyme concentration. Furthermore, catalase and peroxidase formation were only

394

seen for intermittent cyclic DO perturbation. The enzyme activity assay results were consistent

395

with proteomics data, with DO perturbations exhibiting more enzyme activity relative to

396

constant aeration conditions. All DO perturbations, however, showed a similar level of enzyme

397

activity, indicating that enzyme assays have reduced sensitivity than proteomics analysis in

398

discriminating between DO treatments

399

The observed OMP data illustrate that cyclic DO perturbations typically resulted in more OMP

400

removal as compared to constant aeration. Additionally, intermittent cyclic DO perturbations

401

showed larger OMP removal than continuous cyclic DO perturbation. DO perturbation

402

characteristics affected OMP removal, but the relationship is complicated and requires further

403

characterization60. The frequency and range of DO perturbation also influenced OMP removal.

404

For example, for continuous cyclic high DO perturbations, the 0.25 cycles/hr treatment (in 0-

405

10 mg/L DO range) showed more OMP removal than 1 cycle/hr (in 5-7 mg/L DO range) and

406

2 cycles/hr (in 6-8 mg/L DO range). Similarly, for the 0.25 cycles/hr treatments, more OMP

407

removal occurred at low DO compared to the high DO condition. Continuous OMP addition in

408

our fed batch mode of reaction operation (i.e., continuous OMP addition) seemed to result in a

409

plateau of OMP removal at ~70%. However, using the fitted first-order rate constants (Table

410

S6) to predict OMP removal in batch mode (i.e., OMPs only added at the start of study) yields

411

an expected 80-100% removal for most OMPs under DO perturbation conditions (Figure S1).

412

Practical Implications

413

This study suggests that exposing wastewater microbial communities to oscillating oxygen

414

concentrations in different aeration regimes can affect bacterial species selection, induce de

415

novo synthesis of oxidative biocatalysts and subsequently increase OMP removal. The use of

416

DO perturbations to improve treatment performance implies the possibility of new and

417

underutilized mechanisms for wastewater treatment plants to use existing infrastructure to

418

improve their treatment efficiencies for problematic organic contaminants.11 With a potential

419

to catalyse OMP biotransformations, the use of induced oxidoreductases for directed

420

transformations in wastewater treatment has just started to be investigated at the laboratory

421

scale.23 In the future, the biodegradation of OMPs may be achieved by improving the

422

performance of existing wastewater treatment systems under conditions that provide directed 13 ACS Paragon Plus Environment

Environmental Science & Technology

423

environmental stresses that result in over-production of beneficial OMP-degrading biocatalyst

424

enzymes. We recommend further investigation of these possibilities.

425 426

ASSOCIATED CONTENT

427

Supporting Information

428

Detailed materials and methods, tables of chemicals used and residual OMPs concentrations

429

and figures of OMPs removal efficiency, microbial speciation and bacterial regulatory system.

430 431

AUTHOR INFORMATION

432

Corresponding Authors

433

*Phone: +64 9 923 4512. Fax: +64 9 373 7462. Email: [email protected] (Naresh

434

Singhal).

435

*Phone: +64 9 923 4512. Fax: +64 9 373 7462. Email: [email protected] (Amrita

436

Bains).

437 438

ORCID

439

Amrita Bains: 0000-0003-2286-6596

440 441

Notes

442

The authors declare no competing financial interest.

443 444 445 446

ACKNOWLEDGEMENTS The study was funded by FRDF grant 3707510/9572 from the Faculty of Engineering to

447

Singhal, Lear and Greenwood. Bains thanks the University of Auckland for a Doctoral

448

Scholarship during her PhD study. Mabel Smith and Matthew Fung, students visiting our lab

449

from the University of Alberta, assisted with the development of analytical methods.

450 451

REFERENCES

452

(1)

453

Kim, S. D.; Cho, J.; Kim, I. S.; Vanderford, B. J.; Snyder, S. A. Occurrence and Removal of Pharmaceuticals and Endocrine Disruptors in South Korean Surface, Drinking, and 14 ACS Paragon Plus Environment

Page 14 of 26

Page 15 of 26

Environmental Science & Technology

454 455

Waste Waters. Water Res. 2007, 41 (5), 1013–1021. (2)

Luo, Y.; Guo, W.; Ngo, H. H.; Nghiem, L. D.; Hai, F. I.; Zhang, J.; Liang, S.; Wang, X.

456

C. A Review on the Occurrence of Micropollutants in the Aquatic Environment and

457

Their Fate and Removal during Wastewater Treatment. Sci. Total Environ. 2014, 473–

458

474, 619–641.

459

(3)

Mezzelani, M.; Gorbi, S.; Regoli, F. Pharmaceuticals in the Aquatic Environments:

460

Evidence of Emerged Threat and Future Challenges for Marine Organisms. Mar.

461

Environ. Res. 2018, 140 , 41–60.

462

(4)

Routti, H.; Berg, M. K.; Lille-Langøy, R.; Øygarden, L.; Harju, M.; Dietz, R.; Sonne,

463

C.; Goksøyr, A. Environmental Contaminants Modulate the Transcriptional Activity of

464

Polar Bear (Ursus Maritimus) and Human Peroxisome Proliferator-Activated Receptor

465

Alpha (PPARA). Sci. Rep. 2019, 9 (1), 1–14.

466

(5)

Richmond, E. K.; Rosi, E. J.; Walters, D. M.; Fick, J.; Hamilton, S. K.; Brodin, T.;

467

Sundelin, A.; Grace, M. R. A Diverse Suite of Pharmaceuticals Contaminates Stream

468

and Riparian Food Webs. Nat. Commun. 2018, 9 (1), 1–9.

469

(6)

Kidd, K. A.; Blanchfield, P. J.; Mills, K. H.; Palace, V. P.; Evans, R. E.; Lazorchak, J.

470

M.; Flick, R. W. Collapse of a Fish Population after Exposure to a Synthetic Estrogen.

471

Proc. Nat. Acad. Sci. 2007, 104 (21), 8897–8901.

472

(7)

473 474

Sumpter, J. P.; Donnachie, R. L.; Johnson, A. C. The Apparently Very Variable Potency of the Anti-Depressant Fluoxetine. Aquat. Toxicol. 2014, 151, 57–60.

(8)

Galus, M.; Jeyaranjaan, J.; Smith, E.; Li, H.; Metcalfe, C.; Wilson, J. Y. Chronic Effects

475

of Exposure to a Pharmaceutical Mixture and Municipal Wastewater in Zebrafish.

476

Aquat. Toxicol. 2013, 132–133, 212–222..

477

(9)

Joss, A.; Zabczynski, S.; Göbel, A.; Hoffmann, B.; Löffler, D.; McArdell, C. S.; Ternes,

478

T. A.; Thomsen, A.; Siegrist, H. Biological Degradation of Pharmaceuticals in

479

Municipal Wastewater Treatment: Proposing a Classification Scheme. Water Res. 2006,

480

40 (8), 1686–1696.

481

(10)

Stackelberg, P. E.; Gibs, J.; Furlong, E. T.; Meyer, M. T.; Zaugg, S. D.; Lippincott, R.

482

L. Efficiency of Conventional Drinking-Water-Treatment Processes in Removal of

483

Pharmaceuticals and Other Organic Compounds. Sci. Total Environ. 2007, 377 (2–3),

484

255–272.

485

(11)

Falås, P.; Wick, A.; Castronovo, S.; Habermacher, J.; Ternes, T. A.; Joss, A. Tracing the

486

Limits of Organic Micropollutant Removal in Biological Wastewater Treatment. Water

487

Res. 2016, 95, 240–249. 15 ACS Paragon Plus Environment

Environmental Science & Technology

488

(12)

489 490

Systems in Bacteria. Nat. Rev. Genet. 2004, 5 (3), 169–178. (13)

491 492

(14) (15) (16) (17)

Gianfreda, L.; Xu, F.; Bollag, J.-M. Laccases: A Useful Group of Oxidoreductive Enzymes. Bioremed J. 1999, 3 (1), 1–26.

(18)

501 502

Torres, E.; Bustos-Jaimes, I.; Le Borgne, S. Potential Use of Oxidative Enzymes for the Detoxification of Organic Pollutants. Appl Catal B: Environmental 2003, 46 (1), 1–15.

499 500

Rao, M. A.; Scelza, R.; Acevedo, F.; Diez, M. C.; Gianfreda, L. Enzymes as Useful Tools for Environmental Purposes. Chemosphere 2014, 107, 145–162.

497 498

Imlay, J. A. The Molecular Mechanisms and Physiological Consequences of Oxidative Stress: Lessons from a Model Bacterium. Nat. Rev. Microbiol. 2013, 11 (7), 443–454.

495 496

Imlay, J. A. The Molecular Mechanisms and Physiological Consequences of Oxidative Stress: Lessons from a Model Bacterium. Nat. Rev. Microbiol. 2013, 11 (7), 443–454..

493 494

McAdams, H. H.; Srinivasan, B.; Arkin, A. P. The Evolution of Genetic Regulatory

Dȩbska, J.; Kot-Wasik, A.; Namieśnik, J. Fate and Analysis of Pharmaceutical Residues in the Aquatic Environment. Crit. Revi. Anal. Chem. 2004, 34 (1), 51–67.

(19)

Popova, I. E.; Deng, S. A High-Throughput Microplate Assay for Simultaneous

503

Colorimetric Quantification of Multiple Enzyme Activities in Soil. Appl. Soil. Ecol.

504

2010, 45 (3), 315–318.

505

(20)

Brissos, V.; Pereira, L.; Munteanu, F. D.; Cavaco-Paulo, A.; Martins, L. O. Expression

506

System of CotA-Laccase for Directed Evolution and High-Throughput Screenings for

507

the Oxidation of High-Redox Potential Dyes. Biotechnol. J .2009, 4 (4), 558–563.

508

(21)

Johnston, W. A.; Weiliang Huang; De Voss, J. J.; Hayes, M. A.; Gillam, E. M. J.

509

Quantitative Whole-Cell Cytochrome P450 Measurement Suitable for High-Throughput

510

Application. J. Biomol. Screen 2008, 13 (2), 135–141.

511

(22)

Khatoon, N.; Jamal, A.; Ali, M. I. Polymeric Pollutant Biodegradation through

512

Microbial Oxidoreductase: A Better Strategy to Safe Environment. Int. J.Biol.

513

Macromol. 2017, 105, 9–16.

514

(23)

Alneyadi, A. H.; Rauf, M. A.; Ashraf, S. S. Oxidoreductases for the Remediation of

515

Organic Pollutants in Water–a Critical Review. Crit. Rev. Biotechnol. 2018, 38 (7), 971–

516

988.

517

(24)

Martínez, A. T.; Ruiz-Dueñas, F. J.; Camarero, S.; Serrano, A.; Linde, D.; Lund, H.;

518

Vind, J.; Tovborg, M.; Herold-Majumdar, O. M.; Hofrichter, M.; et al. Oxidoreductases

519

on Their Way to Industrial Biotransformations. Biotechnol Adv 2017, 35 (6), 815–831.

520

(25)

Allison, S.; Sinsabaugh, B. Enzyme Assays for Fresh Litter and Soil. 2008, 1–4.

521

(26)

German, D. P.; Weintraub, M. N.; Grandy, A. S.; Lauber, C. L.; Rinkes, Z. L.; Allison, 16 ACS Paragon Plus Environment

Page 16 of 26

Page 17 of 26

Environmental Science & Technology

522

S. D. Optimization of Hydrolytic and Oxidative Enzyme Methods for Ecosystem

523

Studies. Soil Biol Biochem 2011, 43 (7), 1387–1397.

524

(27)

Vanderford, B. J.; Pearson, R. A.; Rexing, D. J.; Snyder, S. A. Analysis of Endocrine

525

Disruptors, Pharmaceuticals, and Personal Care Products in Water Using Liquid

526

Chromatography/Tandem Mass Spectrometry. Anal. Chem. 2003, 75 (22), 6265–6274.

527

(28)

Ferrer, I.; Thurman, M. EPA Method 1694: Agilent's 6410A LC/MS/MS Solution for

528

Pharmaceuticals and Personal Care Products in Water, Soil, Sediment, and Biosolids

529

by HPLC/MS/MS: Application Note (Center for Environmental Mass Spectrometry

530

University of Colorado Civil, Environmental, and Architectural Engineering, Boulder,

531

CO). 2008.

532

(29)

Klindworth, A.; Pruesse, E.; Schweer, T.; Peplies, J.; Quast, C.;Horn, M.; Glöckner, F.

533

O. Evaluation of General 16S Ribosomal RNA Gene PCR Primers for Classical and

534

Next-Generation Sequencing-Based Diversity Studies. Nucleic Acids Res. 2013, 41 (1),

535

1–11.

536

(30)

537 538

Edgar, R. C. UPARSE : Highly Accurate OTU Sequences from Microbial Amplicon Reads. Nat Methods. 2013, 10 (10), doi:10.1038/nmeth.2604.

(31)

Ramirez, K. S.; Leff, J. W.; Barbera, A.; Bates, S. T.; Betley, J.; Crowther, T. W.; Kelly,

539

E. F.; Oldfield, E. E.; Shaw, E. A.; Steenbock, C.; Bradford, M. A.; Wall, D. H; Fierer

540

N. Biogeographic Patterns in Below-Ground Diversity in New York City ’ s Central

541

Park Are Similar to Those Observed Globally. Proc Biol Sci. 2014, 281 (1795), pii:

542

20141988. doi: 10.1098/rspb.2014.1988.

543

(32)

Wang, Q.; Garrity, G. M.; Tiedje, J. M.; Cole, J. R. Naı¨ve Bayesian Classifier for Rapid

544

Assignment of RRNA Sequences into New Bacterila Taxonomy. Appl. Environ.

545

Microbial. 2007, 73 (16), 5261–5267.

546

(33)

Caporaso, J. G.; Kuczynski, J.; Stombaugh, J.; Bittinger, K.; Bushman, F. D.; Costello,

547

E. K.; Fierer, N.; Peña, A. G.; Goodrich, J. K.; Gordon, J. I.; Huttley, G. A.; Kelley, S.

548

T.; Knights, D.; Koenig, J. E.; Ley, R. E.; Lozupne, C. A.; McDonald, D.; Muegge, B.

549

D.; Pirrung, M.; Reeder, J.; Sevinsky, J. R.; Turnbaugh, P. J.; Waletrs, W. A.; Widmann,

550

J.; Yatsunenko, T.; Zaneveld, J.; Knight, R. CorrespondEnce QIIME Allows Analysis

551

of High- Throughput Community Sequencing Data Intensity Normalization Improves

552

Color Calling in SOLiD Sequencing. Nat Methods. 2010, 7 (5), 335–336.

553

(34)

Mcdonald, D.; Price, M. N.; Goodrich, J.; Nawrocki, E. P.; Desantis, T. Z.; Probst, A.;

554

Andersen, G. L.; Knight, R.; Hugenholtz, P. An Improved Greengenes Taxonomy with

555

Explicit Ranks for Ecological and Evolutionary Analyses of Bacteria and Archaea. The 17 ACS Paragon Plus Environment

Environmental Science & Technology

556 557

ISME J, 2011, 6 (3), 610–618. (35)

Partridge, J. D.; Sanguinetti, G.; Dibden, D. P.; Roberts, R. E.; Poole, R. K.; Green, J.

558

Transition of Escherichia Coli from Aerobic to Micro-Aerobic Conditions Involves Fast

559

and Slow Reacting Regulatory Components. J Biol Chem 2007, 282 (15), 11230–11237.

560

(36)

Kaur, A.; Van, P. T.; Busch, C. R.; Robinson, C. K.; Pan, M.; Pang, W. L.; Reiss, D. J.;

561

Diruggiero, J.; Baliga, N. S. Coordination of Frontline Defense Mechanisms under

562

Severe Oxidative Stress. Mol Syst Biol 2010, 6 (393), 1–16.

563

(37)

Gonzalez-martinez, A.; Sihvonen, M.; Muñoz-palazon, B.; Mikola, A.; Vahala, R.

564

Microbial Ecology of Full-Scale Wastewater Treatment Systems in the Polar Arctic

565

Circle : Archaea , Bacteria and Fungi. Sci. Rep. 2018, , 8 (1), 2208, doi: 10.1038/s41598-

566

018-20633-5.

567

(38)

Isazadeh, S.; Jauffur, S.; Frigon, D. Bacterial Community Assembly in Activated

568

Sludge: Mapping Beta Diversity across Environmental Variables. Microbiology Open

569

2016, 5 (6), 1050–1060.

570

(39)

571 572

Kumar, S. Engineering Cytochrome P450 Biocatalysts for Biotechnology, Medicine and Bioremediation. Expert Opin Drug Metab Toxicol 2010, 6 (2), 115–131.

(40)

Malandain, C.; Fayolle-Guichard, F.; Vogel, T. M. Cytochromes P450-Mediated

573

Degradation of Fuel Oxygenates by Environmental Isolates. FEMS Microbiol Ecol

574

2010, 72 (2), 289–296.

575

(41)

Seki, M.; Oikawa, J. I.; Taguchi, T.; Ohnuki, T.; Muramatsu, Y.; Sakamoto, K.; Amachi,

576

S. Laccase-Catalyzed Oxidation of Iodide and Formation of Organically Bound Iodine

577

in Soils. Environ Sci Technol 2013, 47 (1), 390–397.

578

(42)

579 580

Majeau, J. A.; Brar, S. K.; Tyagi, R. D. Laccases for Removal of Recalcitrant and Emerging Pollutants. Bioresour Technol 2010, 101 (7), 2331–2350.

(43)

Wackett, L.; Sadowsky, M.; Martinez, B.; Shapir, N. Biodegradation of Atrazine and

581

Related S-Triazine Compounds: From Enzymes to Field Studies. Appl. Microbiol.

582

Biotechnol. 2002, 58 (1), 39–45.

583

(44)

Naghdi, M.; Taheran, M.; Brar, S. K.; Kermanshahi-pour, A.; Verma, M.; Surampalli,

584

R. Y. Biotransformation of Carbamazepine by Laccase-Mediator System: Kinetics, by-

585

Products and Toxicity Assessment. Process Biochem 2018, 67 (January), 147–154.

586

(45)

Yang, J.; Li, W.; Bun Ng, T.; Deng, X.; Lin, J.; Ye, X. Laccases: Production, Expression

587

Regulation, and Applications in Pharmaceutical Biodegradation. Front. Microbiol.

588

2017, 8 (MAY).

589

(46)

Allison, G. The Influence of Species Diversity and Stress Intensity on Community 18 ACS Paragon Plus Environment

Page 18 of 26

Page 19 of 26

Environmental Science & Technology

590 591

Resistance and Resilience. Ecol Monogr 2004, 74 (1), 117–134. (47)

592 593

Green, J.; Paget, M. S. Bacterial Redox Sensors. Nat Rev Microbiol 2004, 2 (12), 954– 966.

(48)

594

Johnson, L. A.; Hug, L. A. Distribution of Reactive Oxygen Species Defense Mechanisms across Domain Bacteria. Free Radic Biol Med 2019, No. October 2018.

595

(49)

Imlay, J. A. Oxidative Stress : Lessons From a Model Bacterium; 2014; Vol. 11.

596

(50)

Schieber, M.; Chandel, N. S. ROS Function in Redox Signaling and Oxidative Stress.

597 598

Curr Biol 2014, 24 (10), R453–R462. (51)

Levanon, S. S.; San, K. Y.; Bennett, G. N. Effect of Oxygen on the Escherichia Coli

599

ArcA and FNR Regulation Systems and Metabolic Responses. Biotechnol Bioeng 2005,

600

89 (5), 556–564.

601

(52)

Bekker, M.; Alexeeva, S.; Laan, W.; Sawers, G.; De Mattos, J. T.; Hellingwerf, K. The

602

ArcBA Two-Component System of Escherichia Coli Is Regulated by the Redox State

603

of Both the Ubiquinone and the Menaquinone Pool. J. Bacteriol. 2010, 192 (3), 746–

604

754.

605

(53)

Loui, C.; Chang, A. C.; Lu, S. Role of the ArcAB Two-Component System in the

606

Resistance of Escherichia Coli to Reactive Oxygen Stress. BMC Microbiol 2009, 9, 1–

607

14.

608

(54)

Ibrahim, S. A.; Crack, J. C.; Rolfe, M. D.; Borrero-De Acuña, J. M.; Thomson, A. J.; Le

609

Brun, N. E.; Schobert, M.; Stapleton, M. R.; Green, J. Three Pseudomonas Putida FNR

610

Family Proteins with Different Sensitivities to O2. J Biol Chem 2015, 290

611

(27), 16812–16823.

612

(55)

Bueno, E.; Mesa, S.; Bedmar, E. J.; Richardson, D. J.; Delgado, M. J. Bacterial

613

Adaptation of Respiration from Oxic to Microoxic and Anoxic Conditions: Redox

614

Control. Antioxid Redox Signal 2012, 16 (8), 819–852.

615

(56)

616 617

Imlay, J. A. Pathways of Oxidative Damage. Annu Rev Microbiol 2004, 57 (1), 395– 418.

(57)

Gonzalez-Gil, L.; Carballa, M.; Lema, J. M. Cometabolic Enzymatic Transformation of

618

Organic Micropollutants under Methanogenic Conditions. Environ Sci Technol 2017,

619

51 (5), 2963–2971.

620

(58)

Svenningsen, N. B.; Pérez-Pantoja, D.; Nikel, P. I.; Nicolaisen, M. H.; De Lorenzo, V.;

621

Nybroe, O. Pseudomonas Putida Mt-2 Tolerates Reactive Oxygen Species Generated

622

during Matric Stress by Inducing a Major Oxidative Defense Response. BMC Microbiol

623

2015, 15 (1), 1–11. 19 ACS Paragon Plus Environment

Environmental Science & Technology

624

(59)

Sigler, K.; Chaloupka, J.; Brozmanová, J.; Stadler, N.; Höfer, M. Oxidative Stress in

625

Microorganisms - I: Microbial vs. Higher Cells - Damage and Defenses in Relation to

626

Cell Aging and Death. Folia Microbiol 1999, 44 (6), 587–624.

627

(60)

Grandclément, C.; Seyssiecq, I.; Piram, A.; Wong-Wah-Chung, P.; Vanot, G.; Tiliacos,

628

N.; Roche, N.; Doumenq, P. From the Conventional Biological Wastewater Treatment

629

to Hybrid Processes, the Evaluation of Organic Micropollutant Removal: A Review.

630

Water Res 2017, 111, 297–317.

631

20 ACS Paragon Plus Environment

Page 20 of 26

Page 21 of 26

632

Environmental Science & Technology

Table 1: Different temporal dissolved oxygen (DO) supply regimes Operational condition

High-DO Aerobic

Low-DO Aerobic

Dissolved Oxygen range (mg/L) 8-9 0-10 4-6 6-8 2-3 0-3 0-4 2-5

Perturbation frequency (cycles/hr)

Dissolved Oxygen shock durations (min.) Oxygenation duration (min.) Constant 30 10 8 Constant 6 6 12

0 0.25 1 2 0 0.16 0.25 0.5

Deoxygenation duration (min.) Constant 210 50 22 Constant 354 234 106

a)

b)

633 634 635 636 637

Figure 1. Dissolved oxygen profiles of microbial mixed cultures under: (a) high-DO and (b) low DO perturbation in comparison to non-perturbed conditions (constant). Stable DO profiles were achieved in the reactors that were provided constant aerobic (high and low) DO conditions, while aerobic-anaerobic DO transitions were induced with perturbed conditions.

638 639

(Note: DO Conditions; CHA- constant high-DO aerobic, CLA- constant low-DO aerobic, PHA -perturbed high-DO aerobic, PLA- perturbed low-DO aerobic). 21 ACS Paragon Plus Environment

Environmental Science & Technology

a)

b) 640

641 642 643 644 645 646 647 648 649 650

Figure 2. Autoscaled enzyme activity at 48 hour and quantification in cultures; (a) Oxidoreductase activities in DO perturbed and non-perturbed samples. Different rows represent DO perturbed frequencies in high and low-DO regimes; different columns represent expressed oxidoreductase activities (Lacc- laccase, HRP- horseradish peroxidase, LiPlignin peroxidase, Cyt P450- cytochrome P450, β-glcNAc- beta-glucosaminidase, β-glubeta-glucosidase), each detected with more than one dye. The color intensity in each panel shows the auto-scaled enzyme activity (µmol mL-1 min.-1) across the gradient from red (highest enzyme activity) to green (lowest enzyme activity), (b) Relative peptide abundances of expressed oxidoreductases under non-perturbed and perturbed DO conditions.

651 652

(Note: DO Conditions; CHA- constant high-DO aerobic, CLA- constant low-DO aerobic, PHA -perturbed high-DO aerobic, PLA- perturbed low-DO aerobic).

653

22 ACS Paragon Plus Environment

Page 22 of 26

Page 23 of 26

Environmental Science & Technology

654 655 656 657 658 659 660 661 662

Figure 3. Observed removal efficiency of OMPs under different DO conditions in mixed culture fed-batch bioreactors; constant non-perturbed (controls at high aerobic (8-9 mg/L) and low aerobic (2-3 mg/L)) and perturbed (different frequencies (2,1,0.5,0.25,0.16)). Different letters above bars in the same graph indicate significant statistical differences between datasets via ANOVA (p = 0.05) and post hoc Tukey tests; bars indicated with the same letter are not significantly different. The error bars represent standard deviations from four measurements (duplicate analytical results for two biological replicates; n=4).

663 664 665

(Note: DO Conditions; AB- autoclaved biomass control, CHA- constant high-DO aerobic, CLA- constant low-DO aerobic, PHA -perturbed high-DO aerobic, PLA- perturbed low-DO aerobic)

23 ACS Paragon Plus Environment

Environmental Science & Technology

666 667 668 669 670 671 672 673 674

Figure 4. Variation in bacterial community composition (Bray-Curtis similarity) as revealed by distance-based redundancy analysis (dbRDA) of the data, illustrating relationships between bacterial families, enzyme activity and OMPs concentrations remaining in DO perturbed and non-perturbed control cultures at different frequencies. The dbRDA was constrained by the best fitting explanatory variables from multivariate multiple regression analysis (DISTLM) and vectors are shown for predictor variables explaining significant proportions of bacterial community variation (p = 0.001) at 97% OTU richness. Dotted black ellipses indicate 95% confidence interval for clusters of each DO condition (perturbed and non-perturbed).

675 676

(Note: DO Conditions; CHA- constant high-DO aerobic, CLA- constant low-DO aerobic, PHA -perturbed high-DO aerobic, PLA- perturbed low-DO aerobic).

24 ACS Paragon Plus Environment

Page 24 of 26

Page 25 of 26

Environmental Science & Technology

677 678 679 680 681 682 683

Figure 5. Cellular mechanism for ROS formation, gene activation, enzyme production and pollutant degradation. Three sections explains the effect of cyclic oxygen perturbations between aerobic and anaerobic conditions has on: (1) carbon metabolism, NADPH and O2●– formation; (2) gene activation and production of enzymes (SOD, catalase, peroxidase); (3) lowering of hydrogen peroxide, formation enzyme-oxidant complex; and (4) pollutant degradation by enzymes.

25 ACS Paragon Plus Environment

Environmental Science & Technology

84x27mm (150 x 150 DPI)

ACS Paragon Plus Environment

Page 26 of 26