Subscriber access provided by Kaohsiung Medical University
Remediation and Control Technologies
Stable-isotope probing reveals activity and function of autotrophic and heterotrophic denitrifiers in nitrate removal from organic-limited wastewater Wei Xing, Jin-Long Li, Desheng Li, Jincui Hu, Shi-Hai Deng, Yuwei Cui, and Hong Yao Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.8b01993 • Publication Date (Web): 14 Jun 2018 Downloaded from http://pubs.acs.org on June 15, 2018
Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.
is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.
Page 1 of 33
Environmental Science & Technology
1
Stable-isotope probing reveals activity and function of
2
autotrophic and heterotrophic denitrifiers in nitrate
3
removal from organic-limited wastewater
4
Wei Xing †,‡, 1, Jinlong Li †,1, Desheng Li †,‡,*, Jincui Hu †, Shihai Deng †, Yuwei Cui †,
5
Hong Yao †,‡
6 7
†
School of Civil Engineering, Beijing Jiaotong University, Beijing 100044, China
8
‡
Beijing Key Laboratory of Aqueous Typical Pollutants Control and Water Quality
9
Safeguard, Beijing 100044, China
10 11 12
*
13
University, No. 3 Shangyuancun, Haidian District, Beijing 100044, PR China.
14
E-mail:
[email protected]; Tel/Fax: 86-10-5168-4986
15
1
16
Jiaotong University, No. 3 Shangyuancun, Haidian District, Beijing 100044, PR China.
17
E-mail:
[email protected] (Wei Xing);
[email protected] (Jinlong Li)
Corresponding author: Desheng Li, School of Civil Engineering, Beijing Jiaotong
These authors contributed equally to this work. School of Civil Engineering, Beijing
1
ACS Paragon Plus Environment
Environmental Science & Technology
18
ABSTRACT
19
Combined heterotrophic and autotrophic denitrification (HAD) is a sustainable and
20
practical method for removing nitrate from organic-limited wastewater. However, the
21
active microorganisms responsible for denitrification in wastewater treatment have not
22
been clearly identified. In this study, a combined micro-electrolysis, heterotrophic, and
23
autotrophic denitrification (CEHAD) process was established. DNA-based stable isotope
24
probing was employed to identify the active denitrifiers in reactors fed either 13C-labeled
25
inorganic or organic carbon source. The total nitrogen removal efficiencies reached 87.2–
26
92.8% at a low organic carbon concentration (20 mg/L COD). Real-time PCR analysis of
27
the nirS gene as a function of the DNA buoyant density following ultracentrifugation of
28
the total DNA indicated marked 13C-labeling of active denitrifiers. High-throughput
29
sequencing of the fractionated DNA in H13CO3−/12CH312COO−-fed and
30
H12CO3−/13CH313COO−-fed reactors revealed that Thermomonas-like phylotypes were
31
labeled by 13C-bicarbonate, while Thauera-like and Comamonas-like phylotypes were
32
labeled by 13C-acetate. Meanwhile, Arenimonas-like and Rubellimicrobium-like
33
phylotypes were recovered in the “heavy” DNA fraction from both reactors. These results
34
suggest that nitrate removal in CEHAD is catalyzed by various active microorganisms,
35
including autotrophs, heterotrophs, and mixotrophs. Our findings provide a better
36
understanding of the mechanism of nitrogen removal from organic-limited water and
37
wastewater and can be applied to further optimize such processes.
2
ACS Paragon Plus Environment
Page 2 of 33
Page 3 of 33
Environmental Science & Technology
38
TOC/Abstract art
39 40
3
ACS Paragon Plus Environment
Environmental Science & Technology
41
42
INTRODUCTION Nitrate removal is one of the central objectives in water supply and wastewater
43
treatment.1, 2 Biological denitrification, including heterotrophic and autotrophic
44
denitrification, is a common method for nitrate removal.3 In particular, heterotrophic
45
denitrification has been widely reported for its highly efficient denitrification in case of
46
sufficient organic carbon substrate. Conversely, in autotrophic denitrification, microbes
47
reduce nitrate by utilizing inorganic elements or compounds as electron donors and
48
consuming inorganic carbon (bicarbonate) as a carbon source. Notably, autotrophic
49
denitrification is considered an economically affordable and environmentally friendly
50
process for organic-limited water and wastewater, because the dosing of additional
51
organics (e.g., methanol, acetate) can increase the cost and risk of organic residue
52
contamination.4 In recent years, sulfur-, hydrogen-, and biocathode-based autotrophic
53
denitrification processes have been rapidly developed and have attracted widespread
54
attention.5-7 In a previous study, we established a novel combined micro-electrolysis and
55
autotrophic denitrification (CEAD) process based on iron-carbon micro-electrolysis
56
carriers (MECs), which produced H2 and Fe2+ for autotrophic denitrifying bacteria via
57
numerous microscopic galvanic cells.8, 9 With this system, total nitrogen (TN) removal
58
efficiencies of 65–76% were achieved under organic-free conditions, without the need for
59
external H2 or electricity.
60
Combined heterotrophic and autotrophic denitrification (HAD) is considered a
61
practical method for nitrogen removal from nitrogen-contaminated groundwater and
62
drinking water, because natural water is usually organic-limited, rather than organic-free
63
or organic-abundant.10, 11 HAD is also a promising method for wastewater treatment. For 4
ACS Paragon Plus Environment
Page 4 of 33
Page 5 of 33
Environmental Science & Technology
64
instance, secondary effluents from municipal wastewater treatment plants (WWTPs)
65
contain insufficient organic carbon for tertiary denitrification; the application of HAD in
66
advanced treatment can reduce TN emission by utilizing both the small amount of
67
organics in the secondary effluent and inorganic electron donors, without the need for
68
organic dosing. In previous studies, higher denitrification rates were achieved when
69
limited organic carbon sources (liquid carbon, woodchips, and sludge) were subjected to
70
sulfur- and ferrous-based autotrophic denitrification.12, 13 Therefore, to resolve
71
complications with TN removal from organic-limited water and wastewater, it is essential
72
to study the effects of organic carbon on autotrophic denitrification processes and to
73
improve TN removal performance via the HAD process.
74
In the HAD process, heterotrophic and autotrophic denitrifiers are supposed to
75
coexist under organic-limited conditions, but it has been difficult to reveal their actual
76
functions in this process. Since denitrification performance is closely associated with the
77
activity of denitrifiers, identification of the functional denitrifiers is essential to reveal the
78
microbiological mechanisms and for denitrification process optimization. Multiple
79
studies have demonstrated that strains in genera such as Pseudomonas, Thauera, and
80
Thiobacillus reduce nitrate.14, 15 Advanced molecular techniques (e.g., high-throughput
81
sequencing) that can provide comprehensive information without cultivation have been
82
successfully applied for microbial community analysis of denitrification processes.16, 17
83
However, the characterization of the autotrophic and heterotrophic lifestyles of
84
denitrifiers on the basis of 16S rRNA gene analysis alone has been challenging. For
85
instance, it is difficult to determine whether uncultured bacteria of the genus Thauera are
86
autotrophic or heterotrophic denitrifiers, because Thauera can predominate in all of
5
ACS Paragon Plus Environment
Environmental Science & Technology
87
autotrophic, heterotrophic, and HAD systems.13, 15, 18 DNA-based stable-isotope probing (DNA-SIP) is a promising technique for the
88 89
identification of metabolically active bacteria that relies on the cell division of active
90
microorganisms assimilating isotope-labeled substrates.19 Combined with fingerprinting
91
or high-throughput sequencing, DNA-SIP has been successfully applied to investigate the
92
functions and phylogenetic relationships of bacteria in natural environments,20, 21 and has
93
been recently used in studies on drinking water purification and wastewater treatment.22,
94
23
95
dominant autotrophic denitrifier in the CEAD process under organic-free condition.9
96
By using DNA-SIP, we previously identified a Thiobacillus-like bacterium as a
In the present study, to further enhance nitrate removal from wastewater with
97
subsistent but insufficient organic carbon, a combined micro-electrolysis, heterotrophic,
98
and autotrophic denitrification (CEHAD) process was established by using synthetic
99
wastewater with a low total organic carbon to nitrogen ratio (TOC/N; 7.5 mgTOC/L, 40
100
mgNO3−-N/L). The DNA-SIP-based approach was applied to identify active autotrophic,
101
heterotrophic, and mixotrophic denitrifiers on the basis of their different carbon
102
assimilation lifestyles. We expected autotrophic, heterotrophic, and mixotrophic
103
denitrifying phylotypes to be labeled by H13CO3−, 13CH313COO−, and both, respectively.
104
Combined results of taxonomic identification and trophic characterization of denitrifiers
105
provided insights into the microbiological mechanisms involved in nitrogen removal
106
from organic-limited water and wastewater.
107
108
6
ACS Paragon Plus Environment
Page 6 of 33
Page 7 of 33
Environmental Science & Technology
109
MATERIALS AND METHODS
110
CEHAD principles and experimental operation
111
MECs were prepared using powdered iron (17.5% volume), scrap iron (25.0%
112
volume), powdered activated carbon (35.0% volume, passed through a 200-mesh screen),
113
three types of catalysts (each of 2.5% volume), adhesive X (10.0% volume), and foaming
114
agent Y (5.0% volume), and were roasted under an oxygen-free atmosphere at 900–
115
1000°C for 3 h.24 Numerous galvanic cell reactions occur once zero-valent iron and
116
activated carbon in carriers coexist in wastewater. Fe2+ and H2 generated from the anode
117
and cathode, respectively, were utilized as electron donors for the autotrophic denitrifiers,
118
which have been described in detail in our previous studies.8, 25 Both autotrophic and
119
heterotrophic denitrifiers, which use Fe2+/H2 and acetate as electron donors, respectively,
120
contribute to nitrogen removal in the CEHAD process under low TOC/N conditions
121
(Figure 1). Details of these principles are described in the Supporting Information.
122
Three paralleled batch reactors were operated using three 250-mL serum bottles with
123
tight caps, each containing 250 g (wet weight) of MECs and 170 mL of synthetic
124
wastewater. All three reactors were placed in a shaking incubator (HZQ-F100; Donglian
125
Instrument Co. Ltd., Harbin, China) at constant parameters of 30 °C and 90 rpm.
126
Inoculated activated sludge (IAS) was collected from the GBD WWTP (Beijing, China).
127
After centrifugation (3000 rpm for 3 min), 3 g of IAS was weighed and suspended in
128
synthetic substrate, and then incubated with MECs in each reactor for the first three days,
129
without substrate replacement. Starting on the fourth day, all of the liquid and excess
130
suspended sludge (170 mL) in the reactors was replaced with the same volume of
131
synthetic wastewater every 24 h (HRT = 24 h) for supplementation. After each liquid 7
ACS Paragon Plus Environment
Environmental Science & Technology
132
replacement, N2-stripping was performed, after which the reactors were sealed to
133
establish an anoxic condition and to avoid atmospheric CO2 existing in the headspace.
134
Synthetic wastewater was distilled water containing 40 mg/L NO3−-N, 40 mg/L inorganic
135
carbon (IC; equal to 280 mg/L NaHCO3), 7.5 mg/L total organic carbon (TOC;
136
CH3COONa, equal to 20 mg/L COD), and 1 mL/L trace element solution, the latter of
137
which was prepared as described by Till et al.26 HCl solution (10%) was added to the
138
reactors thrice a day to maintain the initial pH of 7.0–7.3.
139
DNA-SIP incubations
140
To distinguish active autotrophic and heterotrophic denitrifiers, the three reactors
141
were fed synthetic wastewater labeled with 13C-bicarbonate or 13C-acetate in different
142
serum bottles (Figure 1). The three reactors were fed the following substrates: (1)
143
unlabeled H12CO3− and 12CH312COO−, with NO3− (control reactor, RC); (2) 13C-labeled
144
H13CO3− and 12CH312COO−, with NO3− (13C-labeled inorganic carbon reactor, RIC); and
145
(3) 13C-labeled 13CH313COO− and H12CO3−, with NO3− (13C-labeled organic carbon
146
reactor, ROC). 13C-labeled H13CO3− and 13CH313COO− (99.9%) were purchased from
147
Cambridge Isotope Laboratories, Inc. (Tewksbury, MA, USA). The reactors were
148
operated for 46 days. To reduce the effect of cross-feeding from decayed microbes, 13C-
149
labeled substrates were employed for 14 days only. During the first 32 days, 12C-
150
unlabeled carbon substrates were used in all three reactors; as of the 33rd day, when the
151
nitrogen removal efficiencies were stabilized, 13C-labeled carbon substrates were
152
employed in RIC and ROC. At the end of the 46-day operation, MECs (100 g, wet
153
weight) with attached biofilms were collected from each reactor for molecular analysis.
154
Because most nitrogen in the reactors would be converted to N2 rather than being 8
ACS Paragon Plus Environment
Page 8 of 33
Page 9 of 33
Environmental Science & Technology
155
assimilated by the microbes, 15N labeling is not practical for evaluating denitrification.
156
However, bacteria that can utilize 13C-substrate but are not involved in denitrification
157
should be excluded. Hence, two supplementary DNA-SIP treatments (T1 and T2) were
158
established together with two supplementary controls (C1 and C2) using the following
159
substrates: (1) only H13CO3−, without NO3− (treatment T1); (2) only 13CH313COO−,
160
without NO3− (treatment T2); (3) unlabeled H12CO3− and NO3− (control C1); (4)
161
unlabeled 13CH313COO− and NO3− (control C2). The concentration of each substrate (if
162
applicable) was identical to that used for the synthetic wastewater treatment in the above-
163
described three reactors. More details are described in the Supporting Information.
164
Analytical methods
165
The concentrations of nitrogen, including NO3−-N, NO2−-N, and NH4+-N, were
166
measured using an ultraviolet spectrophotometer (2102C; UNICO Co. Ltd., Shanghai,
167
China), according to the APHA standard methods.27 The sum of values for NO3−-N,
168
NO2−-N, and NH4+-N was used to estimate TN. TOC and IC were determined in the
169
stable phase using a TOC analyzer (TOC 5000A; Shimatzu, Kyoto, Japan). pH was
170
measured using a pH electrode (pHG-7685A; Tian-yi Co. Ltd., Shanghai, China).
171
DNA extraction, SIP gradient fractionation, and microbial community analysis
172
At the end of operation, seven biofilm samples were collected from the surfaces of
173
the MECs of each of the three reactors (RC, RIC, and ROC) and the four supplementary
174
treatments (T1 and T2) and controls (C1 and C2). Genomic DNA was extracted from the
175
samples using a FastDNA Spin Kit for Soil (MP Biomedicals, Irvine, CA, USA). SIP
176
fractionation was then performed by buoyant density gradient centrifugation according to
177
previously described methods. 9, 28 Buoyant densities of each DNA fraction from all 9
ACS Paragon Plus Environment
Environmental Science & Technology
178
samples are shown in Table S1. Genomic DNA from RC, RIC, and ROC was used for
179
real-time quantitative PCR (qPCR) analysis of the nirS and nirK genes (both encoding
180
nitrite reductase). For all seven samples, the nirS gene was further quantified in the
181
different DNA fractions across the buoyant density gradient. Each analysis was
182
performed in triplicate. qPCR revealed no detectable labeling in the “heavy” DNA
183
fractions of T1 and T2; therefore, only fractionated DNA from the three reactors (RC,
184
RIC, and ROC) was used for high-throughput sequencing and clone analysis. High-
185
throughput sequencing was performed on the Illumina MiSeq platform (Majorbio Bio-
186
pharm Technology Co. Ltd., Shanghai, China) using a set of barcode-containing primers
187
338F/806R for 16S rRNA gene amplification. In total, 40 samples were sequenced, which
188
included genomic DNA samples from the IAS, RC, RIC, and ROC, as well as 12
189
fractionated DNA samples from each reactor. After pre-processing using Trimmomatic
190
software and FLASH, 20,761–41,629 effective sequences were obtained for the different
191
samples. With random and even subsampling, 20,761 sequences in each sample were
192
analyzed, and the average sequence length was 444 bp. Sequences were divided and
193
clustered into operational taxonomic units (OTUs) with a 97% similarity threshold by
194
Uparse. A heatmap was generated using R (version 3.2.3). Bacterial nirS and nirK clone
195
libraries were constructed using the “heavy” DNA fractions of RIC and ROC (the 11th
196
DNA fractions with buoyant densities of 1.745 and 1.746 g/mL, respectively) and the
197
“light” DNA fractions of RC, RIC, and ROC (the 7th DNA fractions with buoyant
198
densities of 1.729, 1.730, and 1.730 g/mL, respectively). For each sample, 30 clone
199
sequences were obtained and analyzed. Details of the PCR conditions and primers (Table
200
S2) are described in the Supporting Information.
10
ACS Paragon Plus Environment
Page 10 of 33
Page 11 of 33
201
202
Environmental Science & Technology
Data deposition Sequences obtained by high-throughput sequencing of genomic DNA (four samples)
203
have been deposited in the NCBI short-read archive under accession numbers
204
SRR4457989–4457992; sequences obtained from fractionated DNA (36 samples) were
205
archived under accession numbers SRR4457945–4458007. All non-identical sequences
206
from clone libraries also have been deposited in the NCBI archive under accession
207
numbers MH432369–432507 (139 sequences for nirS) and MH432265–432359 (95
208
sequences for nirK).
209
RESULTS AND DISCUSSION
210
Nitrogen removal performance in the CEHAD process
211
As shown in Figure 2, nitrogen removal showed a similar trend in each reactor.
212
Nitrate was rapidly reduced initially, likely because the decay of the seed sludge provided
213
electron donors for denitrification. The average concentration of nitrate in the effluent
214
was 9.8 mg/L between days 8 and 17, and decreased to less than 2.7 mg/L after day 21,
215
indicating the establishment of biofilm culture. The average nitrate removal efficiencies
216
during the stable phase in the three reactors ranged from 98.7% to 100%. Notably, NO3−-
217
N can be converted to NH4+-N by chemical reactions (Eq. 6 in Supporting Information)
218
or via the dissimilatory nitrate reduction to ammonium (DNRA) pathway in the presence
219
of iron.29 The concentration of NH4+-N decreased gradually from 20.5–23.8 mg/L (in the
220
first 17 days) to 0.6–1.5 mg/L (in the last three days), demonstrating that biological
221
denitrification dominated over ammonium generation. At the end of operation, TN
222
removal efficiencies of 87.2–92.8% were noted, indicating that highly efficient
11
ACS Paragon Plus Environment
Environmental Science & Technology
223
224
denitrification was achieved in all three reactors. Although small amounts of NO2−-N (1.9–4.4 mg/L) and NH4+-N (0.6–1.5 mg/L)
225
existed in the effluent, CEHAD achieved higher TN removal efficiencies than those
226
reported for Fe(0)-supported autotrophic denitrification, sulfur autotrophic denitrification,
227
and electrochemical denitrification.12, 30, 31 Notably, the CEHAD process yielded lower
228
levels of NH4+-N than the CEAD process under organic-free conditions, despite the use
229
of the same carriers. As shown in Table S3, at the end of operation, the concentrations of
230
NH4+-N in the effluents of the CEAD reactors were 5.6 ± 2.1 mg/L and 10.5 ± 2.2 mg/L
231
(in duplicate) under similar conditions, without the addition of organic carbon.9
232
Comparison of the CEHAD and CEAD performances indicated that a small amount of
233
organic carbon, which was insufficient for complete heterotrophic denitrification,8 could
234
enhance biological denitrification and improve TN removal efficiency through combined
235
heterotrophic and autotrophic denitrification processes. These results were in accordance
236
with those of other studies.12, 13 Thus, the CEHAD process is an efficient and promising
237
method for organic-limited wastewater treatment, such as advanced nitrogen removal
238
from secondary effluent in WWTPs.
239
Because of precipitation and continuous consumption, Fe2+ and H2 supplied to
240
autotrophic denitrifiers could not be accurately quantified. Instead, IC was determined to
241
assess autotrophic denitrification activity from the process perspective. The IC and TOC
242
concentrations in the influent and effluent during the last three days (stable phase) of
243
activity indicated that 6.6 ± 3.6 mg IC/L and 6.3 ± 0.6 mg TOC/L were consumed in the
244
reactors, respectively. These data suggested the existence of autotrophic and
245
heterotrophic denitrification in the CEHAD reactors. However, part of the IC in the
12
ACS Paragon Plus Environment
Page 12 of 33
Page 13 of 33
Environmental Science & Technology
246
reactors was consumed by the autotrophic ammonia-oxidizing bacteria (AOB, with
247
Nitrosomonas-like bacteria being dominant) and nitrite-oxidizing bacteria (NOB, with
248
Nitrotoga-like bacteria being dominant), as evidenced by the high-throughput sequencing
249
results which are discussed in the next subsection. The quantitative contributions of IC
250
and TOC to autotrophic and heterotrophic denitrifiers were difficult to determine and
251
require further study.
252
Characterization of the bacterial communities present in the CEHAD reactor
253
Bacterial community composition was analyzed by high-throughput sequencing of
254
the 16S rRNA gene. Diversity indexes are listed in Table S4. The high coverage index
255
values indicated that the sequence data were sufficient to effectively characterize
256
bacterial communities. Figure S1 shows that the bacterial communities obtained from all
257
three reactors were different from those observed in the IAS. Although they were not
258
identical in specific OTUs (e.g., OTU181) due to random errors, the bacterial
259
communities in the three reactors were generally similar. In the CEHAD reactors, the
260
phylum Proteobacteria accounted for >80% of total bacteria. In total, 355 OTUs were
261
classified at 97% similarity. As shown in Figure 3, OTU143, OTU1092, OTU181,
262
OTU352, OTU170, and OTU174, which accounted for 21.7% ± 3.4%, 11.1% ± 2.7%, 8.1%
263
± 6.7%, 7.7% ± 4.7%, 6.8% ± 2.7%, and 5.7% ± 2.8% of total bacteria, respectively,
264
were the most predominant, and were affiliated with the genera Thauera, Arenimonas,
265
Comamonas, Rubellimicrobium, Nitrotoga, and Thermomonas, respectively.
266
Many Proteobacteria, particularly those of the α, β, and γ subclasses, are involved in
267
global nitrogen cycling and pollutant removal.32, 33 Many strains in the most predominant
268
genus, Thauera, have been reported to mediate nitrate or nitrite reduction under either 13
ACS Paragon Plus Environment
Environmental Science & Technology
269
autotrophic or heterotrophic conditions using hydrogen or organic carbon as the electron
270
donor.13, 15 However, some strains of Thauera can only reduce nitrate to nitrite, leading to
271
nitrite accumulation, which can be resolved by maintaining a pH close to neutral.34 In
272
addition, Arenimonas, Comamonas, Thermomonas, and Rubellimicrobium are typical
273
denitrifiers that are capable of reducing nitrate to nitrogen gas.17, 35 In addition to
274
denitrifiers, one of the dominant genera, Nitrotoga, reportedly belongs to NOB.36 The
275
genus Nitrosomonas, typically comprising AOB, was also detected in the reactors. This
276
might be because NO3−-N can be converted to NH4+-N by chemical reduction or via the
277
DNRA pathway, thereby providing substrates for AOB and NOB. Furthermore, although
278
N2-stripping was performed to establish an anoxic condition, a strict oxygen-free
279
condition is hard to maintain, and both Nitrotoga and Nitrosomonas can reportedly grow
280
under low-oxygen conditions.37, 38
281
Notably, Thiobacillus, which was identified as the most dominant autotrophic
282
denitrifier in the CEAD process in our previous study,9 was not detected as a dominant
283
genus in the current study. This discrepancy can be explained by the presence of organic
284
carbon in the CEHAD process, which might reduce the competitiveness of obligate
285
autotrophic denitrifiers. Moreover, the biofilms in the CEHAD reactors were cultivated
286
directly from activated sludge of the WWTP, rather than from previously operating
287
reactors; thus, other random factors might be responsible for the absence of Thiobacillus
288
from the reactors used in this study.
289
Identification of autotrophic and heterotrophic denitrifiers by DNA-SIP
290
qPCR of the nirS and nirK genes from genomic DNA and fractionated DNA
291
Both nirS and nirK genes were quantified from genomic DNA of RC, RIC, and ROC 14
ACS Paragon Plus Environment
Page 14 of 33
Page 15 of 33
Environmental Science & Technology
292
by qPCR (Figure S2). The results showed that nirK copy numbers in the reactors were
293
much lower than those of nirS. After incubation for 14 days in the presence of 13C-
294
bicarbonate or 13C-acetate as carbon substrates, respectively, autotrophic and
295
heterotrophic denitrifiers were identified by DNA-SIP.
296
Figure 4 shows the variation in nirS gene copy numbers in different fractionated
297
DNA of the seven biofilm samples. In the unlabeled treatment (RC), all nirS gene copies
298
were concentrated in the 5th–7th fractions (“light” fractions with buoyant densities of
299
1.720–1.728 g/mL). As 13C-DNA is heavier than 12C-DNA, the large peak of nirS copies
300
in the ROC shifted to the 11th fraction (“heavy” DNA fraction with a buoyant density of
301
1.746 g/mL), indicating that 13C-acetate was strongly assimilated by heterotrophic
302
denitrifiers. In the RIC, peak nirS copy number was observed in the 7th (“light”) DNA
303
fraction; however, labeled 13C-DNA was also observed in the 9th–11th fractions (“heavy”
304
fractions with buoyant densities of 1.737–1.745 g/mL), demonstrating that autotrophic
305
denitrifiers that assimilated 13C-bicarbonate coexisted with heterotrophic denitrifiers. The
306
nirS copy number in the “heavy” fractions of the ROC was higher than that in the “heavy”
307
fractions of the RIC, indicating that the denitrifiers in this reactor were inclined to utilize
308
acetate as a carbon source.
309
As shown in Figure 4a, in the absence of nitrate, treatment T1, employing only
310
H13CO3−, resulted in no detectable labeling in the 9th–11th (“heavy”) DNA fractions. In
311
contrast, RIC fed H13CO3−, 12CH312COO−, and NO3− showed labeling in the “heavy”
312
fractions. In addition, remarkable labeling in the “heavy” fractions was detected in the
313
previously reported CEAD reactor fed H13CO3− and NO3−.9 Similar results were observed
314
when comparing ROC and T2 (Figure 4b). These results suggest that the 13C-labeled
15
ACS Paragon Plus Environment
Environmental Science & Technology
315
DNA shown in the “heavy” fractions of RIC and ROC was mainly derived from bacteria
316
actively involved in N-cycling. Before inoculation of the 13C-substrate, the biofilms in the
317
CEHAD and CEAD reactors were enriched from specific organic-limited wastewater
318
containing nitrate; thus, denitrifiers are expected to have been dominant. Therefore, in
319
this study, bacteria that can utilize 13C-substrate but are not involved in N-cycling were
320
not abundant and could not have greatly affected the identification of denitrifiers.
321
However, for environmental samples in which denitrifiers are not highly enriched, 13C-
322
labeled DNA in the “heavy” fraction in treatments T1 and T2 is supposed to originate
323
from bacteria utilizing 13C-substrate without NO3−. In such case, microbial community
324
analysis in the “heavy” DNA fractions of T1 and T2 should be performed to exclude
325
potentially interfering communities.
326
High-throughput sequencing of the fractionated DNA
327
To analyze the carbon utilization patterns of the predominant OTUs, each DNA
328
gradient fraction in RC, RIC, and ROC was subjected to high-throughput 16S rRNA gene
329
sequencing. Based on differences in the relative abundances of OTUs, three typical
330
patterns were observed (Figure 5). The first pattern consisted a predominant peak in the
331
“heavy” fraction of the ROC (demonstrating that 13C-acetate was assimilated) and in the
332
“light” fraction of the RC and RIC (demonstrating that 13C-bicarbonate was not utilized).
333
The bacteria affiliated with these OTUs, such as OTU143 (Figure 5a, affiliated with
334
Thauera) and OTU181 (Figure 5b, affiliated with Comamonas), were identified as
335
heterotrophic denitrifiers. Thauera and Comamonas have been shown to exhibit nitrate
336
removal capacity using organic carbon compounds as electron donors.34, 39 However,
337
Thauera reportedly also is dominant in hydrogen-based autotrophic denitrification.15
16
ACS Paragon Plus Environment
Page 16 of 33
Page 17 of 33
Environmental Science & Technology
338
Therefore, it is difficult to determine whether Thauera is a heterotrophic or autotrophic
339
denitrifier based solely on phylogenetic analysis of the 16S rRNA gene. In the current
340
study, DNA-SIP analysis provided clear molecular evidence that the OTU affiliated with
341
Thauera in the CEHAD reactors represented heterotrophic denitrifiers assimilating
342
acetate.
343
The second pattern consisted of a predominant peak within the “heavy” fraction of
344
the RIC (demonstrating that 13C-bicarbonate was assimilated) and in the “light” fraction
345
of the RC and ROC (demonstrating that 13C-acetate was not utilized). OTU174 (Figure
346
5c), which was affiliated with Thermomonas, displayed this second pattern and
347
represented an autotrophic denitrifier that uses inorganic carbon as the carbon source. In
348
our previous study, Thermomonas was identified as a dominant genus (>20%) in the
349
CEAD process under organic-free conditions.8
350
The third pattern revealed labeled 13C-DNA in the “heavy” fractions of both RIC and
351
ROC, demonstrating that both 13C-bicarbonate and 13C-acetate were assimilated.
352
OTU1092 (affiliated with Arenimonas) and OTU 352 (affiliated with Rubellimicrobium),
353
presented this mixotrophic pattern, evidenced by the shifts to heavier DNA fractions in
354
both RIC and ROC, as shown in Figure 5d and 5e. In previous studies, Arenimonas has
355
been observed in both heterotrophic and bioelectrochemical autotrophic denitrification
356
processes,35, 40 and Rubellimicrobium mesophilum was reported to have a denitrification
357
pathway based on genome sequence data.41 OTUs revealing a mixotrophic pattern likely
358
show partial labeling, since individual species affiliated with a certain OTU can
359
incorporate inorganic carbon while others representing the same OTU can incorporate
360
organic carbon. In this study, OTU1092 showed a more legible peak in the “heavy”
17
ACS Paragon Plus Environment
Environmental Science & Technology
361
fractions of the RIC than that of the ROC, indicating that it was preferentially labeled by
362
13
363
C-bicarbonate. Although multiple studies have reported various genera as denitrifiers, their
364
heterotrophic or autotrophic lifestyle is primarily deduced on the basis of whether organic
365
carbons are employed in the process.5, 13, 15 It is difficult to distinguish heterotrophic and
366
autotrophic denitrifiers in HAD processes because organic and inorganic electron donors
367
coexist in such systems. By combining DNA-SIP with qPCR and high-throughput
368
sequencing, we identified heterotrophic, autotrophic, and mixotrophic denitrifiers, using
369
labeled organic and inorganic carbon compounds as substrates.
370
Clone libraries and phylogenetic analysis of the nirS and nirK genes
371
As denitrifiers are difficult to be identified solely by 16S rRNA gene analysis, clone
372
libraries of the nirS and nirK genes were constructed using the “heavy” and “light” DNA
373
fractions. In nirS clone libraries, 21 OTUs were obtained from 150 clone sequences by
374
clustering at 95% similarity (Table S5). A phylogenetic tree generated on the basis of the
375
nirS clones is shown in Figure S3. NirS-Clone-OTUs 1 to 4 represented four dominant
376
OTUs, which accounted for 42.0%, 9.3%, 10.7%, and 8.7% of the 150 clone sequences.
377
NirS-Clone-OTU3 and nirS-Clone-OTU4 were both highly similar to Thauera sp. 27 in
378
GenBank. They existed in the “heavy” fraction of the ROC and the “light” fraction of the
379
RIC, but were negligible in the “light” fraction of the ROC and the “heavy” fraction of
380
the RIC, which indicated that they could reduce nitrate using acetate rather than
381
bicarbonate. Thus, these two clone OTUs affiliated with Thauera related to heterotrophic
382
denitrifiers. These findings are in agreement with the above-mentioned high-throughput
383
sequencing results, which demonstrated that Thauera-like phylotypes were predominant 18
ACS Paragon Plus Environment
Page 18 of 33
Page 19 of 33
Environmental Science & Technology
384
in the reactors and could be identified as heterotrophic denitrifiers. In contrast, nirS-
385
Clone-OTU1 was present in similar proportions (13/30, 16/30, 10/30, 12/30, 12/30,
386
respectively) in each DNA fraction, indicating that it represented a mixotrophic
387
denitrifier. Furthermore, it was found to be similar to Arenimonas donghaensis.
388
Consistent with this finding, the OTU affiliated with Arenimonas in the high-throughput
389
sequencing analysis also showed a mixotrophic lifestyle. Finally, nirS-Clone-OTU2
390
existed only in the “heavy” fraction of the RIC and the “light” fraction of the ROC and
391
thus, represented an autotrophic denitrifier. The nirS sequence in OTU2 was found to be
392
similar to that of an uncultured bacterium in GenBank, and in all nirS genes of known
393
species, the closest one was from Thiohalobacter thiocyanaticus.
394
In nirK clone libraries, nine OTUs were obtained from 150 clone sequences by
395
clustering at 95% similarity (Table S6). A phylogenetic tree constructed on basis of the
396
nirK clones is shown in Figure S4. More than 80% of clone sequences were clustered
397
with nirK-Clone-OTU1, which showed high sequence similarity with Rubellimicrobium
398
mesophilum, indicating that most denitrifiers in CEHAD reactors detected by the nirK
399
gene primers were affiliated with Rubellimicrobium.
400
Intriguingly, the genera Comamonas and Thermomonas, which were detected by
401
high-throughput sequencing, were not found in the clone libraries. In a primer coverage
402
assessment using NCBI Primer-BLAST, we found that Comamonas serinivorans strain
403
DSM 26136 (CP021455) cannot be detected with the nirS1F/nirS6R pair used in this
404
study, though it was reported to harbor a nirS gene. Meanwhile, both Thermomonas fusca
405
(WP_028838532) and Thermomonas hydrothermalis (WP_072755949) harbor nirK
406
genes, but cannot be detected with the primer pair F1a/R3Cu. Moreover, some strains in
19
ACS Paragon Plus Environment
Environmental Science & Technology
407
the genus Thauera, such as Thauera sp. 28 (WP_037984032), cannot be detected with
408
nirS1F/nirS6R, while others, such as Thauera humireducens strain SgZ-1 (CP014646)
409
can, although both contain a nirS gene. Because bacterial nirS and nirK genes might be
410
more labile than the 16S rRNA gene, it might be difficult to find an appropriate primer
411
pair targeting the nirS or nirK genes in all denitrifiers. It should be emphasized that the
412
primers for nirS and nirK used in this study need to be optimized for more
413
comprehensive studies.
414
Implications and perspectives In this study, autotrophic, heterotrophic, and mixotrophic denitrifiers were identified
415 416
by DNA-SIP combined with 16S rRNA, nirS, and nirK gene analysis. However, our
417
study had some limitations. Technical errors, e.g., in centrifugation efficiency, and
418
fractionation operation are inherent to the SIP technology. Part of the 12C-unlabeled
419
genomic DNA, which is supposed to be distributed to the “light” fractions, may have
420
been spun down during the isopycnic ultracentrifugation of the total DNA extract and
421
therefore be detected in the “heavy” DNA fractions, especially for the bacteria with
422
relatively high GC content or highly abundant bacteria. Therefore, in general, the
423
difference in the microbial abundances in the same DNA fraction of different treatments
424
(RC, RIC, and ROC) was reliable for identifying whether the bacteria assimilated
425
carbonate or acetate, rather than the absolute abundances in a certain treatment.
426
Furthermore, cross-feeding may have occurred in the reactors. After assimilating 13C-
427
labeled bicarbonate, autotrophic denitrifiers would decay and release 13C-labeled cell
428
substances, which could be secondarily assimilated by heterotrophic bacteria. Meanwhile,
429
13
C-labeled acetate could be decomposed by heterotrophic denitrifiers and generate
20
ACS Paragon Plus Environment
Page 20 of 33
Page 21 of 33
Environmental Science & Technology
430
inorganic carbon, which could subsequently be used as inorganic carbon in autotrophic
431
denitrification. To avoid erroneous assumptions regarding carbon utilization patterns, the
432
following three strategies were applied. (1) 13C-labeled substrates were employed only
433
for 14 days, after the reactors had achieved stable performance. Thus, the incubation time
434
for 13C-labeled substrates was shortened, and the possibility of cross-feeding by decayed
435
microbes could be reduced. (2) The liquid and excess suspended sludge was replaced
436
with fresh synthetic wastewater daily, so that the potential for recycling of cell substances
437
was minimized. (3) Because TOC and IC were maintained at 7.5 mg/L and 40 mg/L,
438
respectively, in the synthetic wastewater, inorganic carbon generated by acetate
439
decomposition (a small portion of the TOC, considering that most organic carbon could
440
be assimilated) was much less than the inorganic carbon content in the influent.
441
Therefore, the effects of inorganic carbon generated by acetate decomposition on
442
autotrophic denitrifiers would not be significant. However, time-course DNA-SIP
443
analysis should be carried out in future, to characterize the optimal duration for
444
incubation with the 13C-substrate to separate 13C-labeled DNA and distinguish
445
heterotrophic, autotrophic, and mixotrophic denitrifiers more precisely.
446
This study showed that nitrate removal in CEHAD was catalyzed by a variety of
447
active microorganisms, including autotrophs (Thermomonas-like phylotypes),
448
heterotrophs (Thauera-like and Comamonas-like phylotypes), and mixotrophs
449
(Arenimonas-like and Rubellimicrobium-like phylotypes). The results lend strong support
450
for the performance of denitrification activity in the reactors, which indicated that a small
451
amount of organic carbon enhanced biological denitrification and benefited higher TN
452
removal efficiency. However, real wastewater has a more complex and variable
21
ACS Paragon Plus Environment
Environmental Science & Technology
453
components, and thus, different microbial communities might be existed in systems fed
454
with real wastewater. Therefore, further studies using real wastewater are desired to
455
analyze the microbial communities in situ, by using the real wastewater as the influent
456
and adding the target 13C-labeled substrates in it. Further, studies with different TOC/N
457
and TOC/IC ratios should be conducted to illuminate the effect of component variation
458
on microbial species composition. The microbiological mechanisms of denitrification in
459
real wastewater treatment could thus be profiled comprehensively.
460
ACKNOWLEDGEMENTS
461
This study was financially supported by the National Natural Science Foundation of
462
China (grant nos. 51408028 and 51278034). We would like to express our deep gratitude
463
to Prof. Zhongjun Jia at the Institute of Soil Science, Chinese Academy of Sciences, for
464
his technical support with the DNA-SIP experiments.
465
ASSOCIATED CONTENT
466
Supporting Information
467
Detailed methods for the establishment of the CEHAD process, supplementary
468
experiments, and molecular biological methods. Additional Tables S1 to S6 showing
469
buoyant density of each DNA fraction, PCR primers, CEAD performance, diversity
470
indices, and the results of nirS and nirK clone libraries. Additional Figures S1 to S4
471
related to abundances of bacterial populations, qPCR results of genomic DNA, and
472
phylogenetic trees on the basis of nirS and nirK clones.
473
This information is available free of charge via the Internet at http://pubs.acs.org.
22
ACS Paragon Plus Environment
Page 22 of 33
Page 23 of 33
Environmental Science & Technology
474
AUTHOR INFORMATION
475
Corresponding Author:
476
Desheng Li
477
First Authors:
478
Wei Xing
Tel: +86-10-5168-5917; E-mail:
[email protected] 479
Jinlong Li
Tel: +86-10-5168-5917; E-mail:
[email protected] 480
Notes
481
The authors declare no competing financial interest.
482
REFERENCES
483
1. Buttiglieri, G.; Malpei, F.; Daverio, E.; Melchiori, M.; Nieman, H.; Ligthart, J. Denitrification
484
of drinking water sources by advanced biological treatment using a membrane bioreactor.
485
Desalination. 2005, 178 (1-3), 211-218; DOI:10.1016/j.desal.2004.11.038.
486
Tel: +86-10-5168-4986. E-mail:
[email protected] 2. Chung, J.; Amin, K.; Kim, S.; Yoon, S.; Kwon, K.; Bae, W. Autotrophic denitrification of
487
nitrate and nitrite using thiosulfate as an electron donor. Water Res. 2014, 58, 169-178;
488
DOI:10.1016/j.watres.2014.03.071.
489
3. Wang, A.; Liu, C.; Han, H.; Ren, N.; Lee, D. J. Modeling denitrifying sulfide removal process
490
using artificial neural networks. J. Hazard. Mater. 2009, 168 (2-3), 1274-1279;
491
DOI:10.1016/j.jhazmat.2009.03.006.
492
4. Sahinkaya, E.; Dursun, N. Sulfur-oxidizing autotrophic and mixotrophic denitrification
493
processes for drinking water treatment: elimination of excess sulfate production and alkalinity
494
requirement. Chemosphere. 2012, 89 (2), 144-149; DOI:10.1016/j.chemosphere.2012.05.029. 23
ACS Paragon Plus Environment
Environmental Science & Technology
495
5. Liu, C.; Zhao, D.; Yan, L.; Wang, A.; Gu, Y.; Lee, D. J. Elemental sulfur formation and
496
nitrogen removal from wastewaters by autotrophic denitrifiers and anammox bacteria.
497
Bioresour. Technol. 2015, 191, 332-336; DOI:10.1016/j.biortech.2015.05.027.
498
6. Chen, D.; Yang, K.; Wang, H. Effects of important factors on hydrogen-based autotrophic
499
denitrification in a bioreactor. Desalin. Water Treat. 2016, 57 (8), 3482-3488;
500
DOI:10.1080/19443994.2014.986533.
501
7. Xiao, Y.; Zheng, Y.; Wu, S.; Yang, Z. H.; Zhao, F. Bacterial community structure of autotrophic
502
denitrification biocathode by 454 pyrosequencing of the 16S rRNA gene. Microb. Ecol. 2015,
503
69 (3), 492-499; DOI:10.1007/s00248-014-0492-4.
504
8. Xing, W.; Li, D.; Li, J.; Hu, Q.; Deng, S. Nitrate removal and microbial analysis by combined
505
micro-electrolysis and autotrophic denitrification. Bioresour. Technol. 2016, 211, 240-247;
506
DOI:10.1016/j.biortech.2016.03.044.
507
9. Xing, W.; Li, J.; Cong, Y.; Gao, W.; Jia, Z.; Li, D. Identification of the autotrophic denitrifying
508
community in nitrate removal reactors by DNA-stable isotope probing. Bioresour. Technol.
509
2017, 229, 134-142; DOI:10.1016/j.biortech.2017.01.010.
510
10. Sahinkaya, E.; Dursun, N.; Kilic, A.; Demirel, S.; Uyanik, S.; Cinar, O. Simultaneous
511
heterotrophic and sulfur-oxidizing autotrophic denitrification process for drinking water
512
treatment: control of sulfate production. Water Res. 2011, 45 (20), 6661-6667;
513
DOI:10.1016/j.watres.2011.09.056.
514
11. Qambrani, N. A.; Jung, S. H.; Ok, Y. S.; Kim, Y. S.; Oh, S. E. Nitrate-contaminated
515
groundwater remediation by combined autotrophic and heterotrophic denitrification for sulfate
516
and pH control: batch tests. Environ. Sci. Pollut. R. 2013, 20 (12), 9084-9091;
517
DOI:10.1007/s11356-013-1623-z.
518
12. Hao, R.; Meng, C.; Li, J. Impact of operating condition on the denitrifying bacterial
24
ACS Paragon Plus Environment
Page 24 of 33
Page 25 of 33
Environmental Science & Technology
519
community structure in a 3DBER-SAD reactor. J. Ind. Microbiol. Biotechnol. 2017, 44 (1), 9-
520
21; DOI:10.1007/s10295-016-1853-4.
521
13. Wang, R.; Zheng, P.; Zhang, M.; Zhao, H. P.; Ji, J. Y.; Zhou, X. X.; Li, W. Bioaugmentation of
522
nitrate-dependent anaerobic ferrous oxidation by heterotrophic denitrifying sludge addition: A
523
promising way for promotion of chemoautotrophic denitrification. Bioresour. Technol. 2015,
524
197, 410-415; DOI:10.1016/j.biortech.2015.08.135.
525
14. Hosono, T.; Alvarez, K.; Lin, I. T.; Shimada, J. Nitrogen, carbon, and sulfur isotopic change
526
during heterotrophic (Pseudomonas aureofaciens) and autotrophic (Thiobacillus denitrificans)
527
denitrification reactions. J. Contam. Hydrol. 2015, 183, 72-81;
528
DOI:10.1016/j.jconhyd.2015.10.009.
529
15. Mao, Y.; Xia, Y.; Zhang, T. Characterization of Thauera-dominated hydrogen-oxidizing
530
autotrophic denitrifying microbial communities by using high-throughput sequencing.
531
Bioresour. Technol. 2013, 128, 703-710; DOI:10.1016/j.biortech.2012.10.106.
532
16. Zhou, S.; Huang, T.; Zhang, C.; Fang, K.; Xia, C.; Bai, S.; Zeng, M.; Qiu, X. Illumina MiSeq
533
sequencing reveals the community composition of NirS-Type and NirK-Type denitrifiers in
534
Zhoucun reservoir - a large shallow eutrophic reservoir in northern China. RSC Adv. 2016, 6
535
(94), 91517-91528; DOI:10.1039/c6ra18017e.
536
17. Zhong, F.; Wu, J.; Dai, Y.; Yang, L.; Zhang, Z.; Cheng, S.; Zhang, Q. Bacterial community
537
analysis by PCR-DGGE and 454-pyrosequencing of horizontal subsurface flow constructed
538
wetlands with front aeration. Appl. Microbiol. Biotechnol. 2015, 99 (3), 1499-1512;
539
DOI:10.1007/s00253-014-6063-2.
540
18. Xu, G.; Peng, J.; Feng, C.; Fang, F.; Chen, S.; Xu, Y.; Wang, X. Evaluation of simultaneous
541
autotrophic and heterotrophic denitrification processes and bacterial community structure
542
analysis. Appl. Microbiol. Biotechnol. 2015, 99 (15), 6527-6536; DOI:10.1007/s00253-015-
25
ACS Paragon Plus Environment
Environmental Science & Technology
543 544 545 546
6532-2. 19. Radajewski, S.; Ineson, P.; Parekh, N. R.; Murrell, J. C. Stable-isotope probing as a tool in microbial ecology. Nature. 2000, 403 (6770), 646-649; DOI:10.1038/35001054. 20. Friedrich, M. W. Stable-isotope probing of DNA: insights into the function of uncultivated
547
microorganisms from isotopically labeled metagenomes. Curr. Opin. Biotechnol. 2006, 17 (1),
548
59-66; DOI:10.1016/j.copbio.2005.12.003.
549 550 551
21. Chen, Y.; Murrell, J. C. When metagenomics meets stable-isotope probing: progress and perspectives. Trends Microbiol. 2010, 18 (4), 157-163; DOI:10.1016/j.tim.2010.02.002. 22. Niu, J.; Kasuga, I.; Kurisu, F.; Furumai, H.; Shigeeda, T. Evaluation of autotrophic growth of
552
ammonia-oxidizers associated with granular activated carbon used for drinking water
553
purification by DNA-stable isotope probing. Water Res. 2013, 47 (19), 7053-7065;
554
DOI:10.1016/j.watres.2013.07.056.
555
23. McIlroy, S. J.; Starnawska, A.; Starnawski, P.; Saunders, A. M.; Nierychlo, M.; Nielsen, P. H.;
556
Nielsen, J. L. Identification of active denitrifiers in full-scale nutrient removal wastewater
557
treatment systems. Environ. Microbiol. 2016, 18 (1), 50-64; DOI:10.1111/1462-2920.12614.
558
24. Deng, S.; Li, D.; Yang, X.; Zhu, S.; Li, J. Process of nitrogen transformation and microbial
559
community structure in the Fe(0)-carbon-based bio-carrier filled in biological aerated filter.
560
Environ. Sci. Pollut. R. 2016, 23 (7), 6621-6630; DOI:10.1007/s11356-015-5892-6.
561
25. Deng, S.; Li, D.; Yang, X.; Xing, W.; Li, J.; Zhang, Q. Biological denitrification process based
562
on the Fe(0)-carbon micro-electrolysis for simultaneous ammonia and nitrate removal from
563
low organic carbon water under a microaerobic condition. Bioresour. Technol. 2016, 219, 677-
564
686; DOI:10.1016/j.biortech.2016.08.014.
565 566
26. Till, B. A.; Weathers, L. J.; Alvarez, P. J. Fe(0)-Supported Autotrophic Denitrification. Environ. Sci. Technol. 1998, 32 (5), 634-639; DOI:10.1021/es9707769. 26
ACS Paragon Plus Environment
Page 26 of 33
Page 27 of 33
567 568 569
Environmental Science & Technology
27. APHA. Standard methods for the examination of water and wastewater. American Public Health Association (APHA): Washington, DC, USA, 2005. 28. Xia, W.; Zhang, C.; Zeng, X.; Feng, Y.; Weng, J.; Lin, X.; Zhu, J.; Xiong, Z.; Xu, J.; Cai, Z.;
570
Jia, Z. Autotrophic growth of nitrifying community in an agricultural soil. ISME J. 2011, 5 (7),
571
1226-1236; DOI:10.1038/ismej.2011.5.
572
29. Robertson, E. K.; Roberts, K. L.; Burdorf, L. D.; Cook, P.; Thamdrup, B. Dissimilatory nitrate
573
reduction to ammonium coupled to Fe(II) oxidation in sediments of a periodically hypoxic
574
estuary. Limnol. Oceanogr. 2016, 61 (1), 365-381; DOI:10.1002/lno.10220.
575 576
30. Dash, B. P.; Chaudhari, S. Electrochemical denitrificaton of simulated ground water. Water Res. 2005, 39 (17), 4065-4072; DOI:10.1016/j.watres.2005.07.032.
577
31. An, Y.; Li, T.; Jin, Z.; Dong, M.; Li, Q.; Wang, S. Decreasing ammonium generation using
578
hydrogenotrophic bacteria in the process of nitrate reduction by nanoscale zero-valent iron.
579
Sci. Total Environ. 2009, 407 (21), 5465-5470; DOI:10.1016/j.scitotenv.2009.06.046.
580
32. Saarenheimo, J.; Tiirola, M. A.; Rissanen, A. J. Functional gene pyrosequencing reveals core
581
proteobacterial denitrifiers in boreal lakes. Front Microbiol. 2015, 6, 1-12; DOI:Artn
582
67410.3389/Fmicb.2015.00674.
583
33. Shen, Z.; Zhou, Y.; Wang, J. Comparison of denitrification performance and microbial
584
diversity using starch/polylactic acid blends and ethanol as electron donor for nitrate removal.
585
Bioresour. Technol. 2013, 131, 33-39; DOI:10.1016/j.biortech.2012.12.169.
586
34. Li, W.; Shan, X. Y.; Wang, Z. Y.; Lin, X. Y.; Li, C. X.; Cai, C. Y.; Abbas, G.; Zhang, M.; Shen,
587
L. D.; Hu, Z. Q.; Zhao, H. P.; Zheng, P. Effect of self-alkalization on nitrite accumulation in a
588
high-rate denitrification system: Performance, microflora and enzymatic activities. Water Res.
589
2016, 88, 758-765; DOI:10.1016/j.watres.2015.11.003.
590
35. Liu, H.; Yan, Q.; Shen, W. Biohydrogen facilitated denitrification at biocathode in 27
ACS Paragon Plus Environment
Environmental Science & Technology
591
bioelectrochemical system (BES). Bioresour. Technol. 2014, 171, 187-192;
592
DOI:10.1016/j.biortech.2014.08.056.
593
36. Lucker, S.; Schwarz, J.; Gruber-Dorninger, C.; Spieck, E.; Wagner, M.; Daims, H. Nitrotoga-
594
like bacteria are previously unrecognized key nitrite oxidizers in full-scale wastewater
595
treatment plants. ISME J. 2015, 9 (3), 708-720; DOI:10.1038/ismej.2014.158.
596
37. Park, H. D.; Noguera, D. R. Characterization of two ammonia-oxidizing bacteria isolated
597
from reactors operated with low dissolved oxygen concentrations. J. Appl. Microbiol. 2010,
598
102 (5), 1401-1417; DOI: 10.1111/j.1365-2672.2006.03176.x.
599
38. Park, H. D.; Noguera, D. R. Nitrospira community composition in nitrifying reactors operated
600
with two different dissolved oxygen levels. J. Microbiol. Biotechnol. 2008, 18 (8), 1470-1474;
601
DOI:10.1111/j.1365-2672.2006.03176.x.
602
39. Lemmer, H.; Zaglauer, A.; Neef, A.; Meier, H.; Amann, R., Denitrification in a methanol-fed
603
fixed-bed reactor. Part 2: Composition and ecology of the bacterial community in the biofilms.
604
Water Res. 1997, 31 (8), 1903-1908; DOI:10.1016/S0043-1354(97)00027-4.
605
40. Liu, C.; Zhao, C.; Wang, A.; Guo, Y.; Lee, D. J. Denitrifying sulfide removal process on high-
606
salinity wastewaters. Appl. Microbiol. Biotechnol. 2015, 99 (15), 6463-6469;
607
DOI:10.1007/s00253-015-6505-5.
608
41. Riedel, T.; Spring, S.; Fiebig, A.; Petersen, J.; Göker, M.; Klenk, H. P. Genome sequence of
609
the pink to light reddish-pigmented Rubellimicrobium mesophilum type strain (DSM
610
19309T), a representative of the Roseobacter group isolated from soil, and emended
611
description of the species. Standard. Genomic Sci. 2014, 9 (3), 902-913;
612
DOI:10.4056/sigs.5621012.
613
28
ACS Paragon Plus Environment
Page 28 of 33
Page 29 of 33
Environmental Science & Technology
614 615
Figure 1. Schematic of the experimental procedures. Three paralleled reactors were
616
operated for 46 days, including 32-day unlabeled incubation and 14-day DNA-based
617
stable-isotope probing (DNA-SIP) incubation periods. Unlabeled 12C-substances, 13C-
618
labeled inorganic carbon, and 13C-labeled organic carbon were employed in the RC
619
(control), RIC (reactor with inorganic carbon), and ROC (reactor with organic carbon),
620
respectively.
621 622
29
ACS Paragon Plus Environment
Environmental Science & Technology
623 624
Figure 2. Performance of the combined micro-electrolysis, heterotrophic, and
625
autotrophic denitrification (CEHAD) process in the presence of low total carbon
626
to nitrogen (TOC/N) ratios. (a–c) Results obtained from the (a) control reactor with
627
12
628
with 13C-labeled organic carbon (ROC). The dashed line indicates the time point at
629
which the labeled substrates were added; it divides the periods of unlabeled incubation
630
(first 32 days) and DNA-based stable-isotope probing (DNA-SIP) incubation
631
(subsequent 14 days). NRE represents total nitrogen (TN) removal efficiency.
C-substrate (RC), (b) reactor with 13C-labeled inorganic carbon (RIC), and (c) reactor
30
ACS Paragon Plus Environment
Page 30 of 33
Page 31 of 33
Environmental Science & Technology
632 633
Figure 3. Hierarchical heat map analysis of the genomic DNA obtained from the
634
IAS, RC, RIC, and ROC, according to high-throughput sequencing of 16S rRNA at
635
the OTU level. OTUs were clustered with 97% similarity. The relative abundance of
636
OTUs is indicated by color intensity, which represents the percentage of the targeted
637
OTU reads in the total bacterial 16S rRNA gene reads in each genomic DNA. All OTUs
638
accounting for >5% in the RC, RIC, and ROC are listed. OTU, operational taxonomic
639
unit; IAS, inoculated active sludge; RC, control reactor with 12C-substrate; RIC, reactor
640
with 13C-labeled inorganic carbon; ROC, reactor with 13C-labeled organic carbon.
641 642
31
ACS Paragon Plus Environment
Environmental Science & Technology
643 644
Figure 4. Quantitative distribution of the nirS gene copies ratio in different
645
fractionated DNA in a) RC, RIC, C1, and T1; b) RC, ROC, C2, and T2. The nirS
646
copies ratios were defined by copy numbers in the present fraction to that in total
647
fractions. RC, control reactor with 12C-substrate (H12CO3− and 12CH312COO−, with NO3−);
648
RIC, reactor with 13C-labeled inorganic carbon (H13CO3− and 12CH312COO−, with NO3−);
649
ROC, reactor with 13C-labeled organic carbon (H12CO3− and 13CH313COO−, with NO3−);
650
supplementary T1 (only contains H13CO3−, without NO3−); control C1 (H12CO3− and
651
NO3−); supplementary T2 (only contains 13CH313COO−, without NO3−); control C2
652
(12CH312COO− and NO3−).
653
32
ACS Paragon Plus Environment
Page 32 of 33
Page 33 of 33
Environmental Science & Technology
654 655
Figure 5. Relative abundance of 16S rRNA gene sequences across the entire buoyant
656
gradient of the DNA fractions for three reactors, based on high-throughput
657
sequencing. (a–e) Results of (a) OTU143, (b) OTU181, (c) OTU174, (d) OTU1092, and
658
(e) OTU352, affiliated with the genera Thauera, Comamonas, Thermomonas,
659
Arenimonas, and Rubellimicrobium, respectively. The relative abundance is expressed as
660
the percentage of the targeted OTU reads in the total bacterial 16S rRNA gene reads in
661
each DNA gradient fraction. OTU, operational taxonomic unit; RC, control reactor with
662
12
663
labeled organic carbon.
C-substrate; RIC, reactor with 13C-labeled inorganic carbon; ROC, reactor with 13C-
33
ACS Paragon Plus Environment