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