Stable-Isotope Probing Reveals the Activity and Function of

Jun 14, 2018 - I&EC Analytical Edition .... (1,2) Biological denitrification, including heterotrophic and autotrophic ... efficiencies of 65–76% wer...
1 downloads 0 Views 1MB Size
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