Insight into c-di-GMP Regulation in Anammox ... - ACS Publications

Jul 24, 2017 - Environmental Science & Technology .... Based on our findings, we believe that c-di-GMP regulation and EPS production of the anammox ...
0 downloads 3 Views 1MB Size
Subscriber access provided by PEPPERDINE UNIV

Article

Insight into c-di-GMP Regulation in Anammox Aggregation in response to Alternating Feed Loadings Yongzhao Guo, Sitong Liu, Xi Tang, Chao Wang, Zhao Niu, and Ying Feng Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.6b06396 • Publication Date (Web): 24 Jul 2017 Downloaded from http://pubs.acs.org on July 25, 2017

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 free 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 accessible to all readers and 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.

Environmental Science & Technology 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 38

Environmental Science & Technology

1

Insight into c-di-GMP Regulation in Anammox Aggregation in

2

response to Alternating Feed Loadings

3

4

Yongzhao Guo 1,2, Sitong Liu 1, 2, *, Xi Tang 1, Chao Wang 3, Zhao Niu 1,2, Ying Feng 1

5

6

1

7

University, Beijing 100871, China

8

2

9

518055, China

Key Laboratory of Water and Sediment Sciences, Ministry of Education of China, Peking

School of Environment and Energy, Peking University Shenzhen Graduate School, Shenzhen

10

3

11

Science and Technology, Dalian University of Technology, Dalian 116024, China

12

*Corresponding author: Sitong Liu

13

Address: College of Environmental Science and Engineering, Peking University, Yiheyuan Road,

14

No.5, Haidian District, Beijing 100871, China.

15

E-mail: [email protected]

16

Tel/Fax: 0086-10-62754290.

Key Laboratory of Industrial Ecology and Environmental Engineering, School of Environmental

1

ACS Paragon Plus Environment

Environmental Science & Technology

Page 2 of 38

17

Abstract

18

Substrate concentrations generally fluctuate in wastewaters. However, how anammox

19

biomass behaves to overcome the stress of alternating feed loadings remains unclear.

20

Here, we combined long-term reactor operation, batch tests, 16S rRNA transcript

21

sequencing, and metabolomics analysis to investigate the aggregation of anammox

22

biomass under the regulation of c-di-GMP, a key second messenger, in response to

23

alternating feed loadings. We demonstrated that the aggregation process was

24

significantly faster under alternating loadings and was significantly correlated with

25

higher levels of c-di-GMP and extracellular polymeric substances (EPS) production.

26

The increase in c-di-GMP was positively correlated with a higher relative transcript

27

expression level in the c-di-GMP pathway-dependent community. The targeted

28

metabolomics results indicated that the increased production of fructose 6-phosphate

29

and UDP-N-acetyl-D-glucosamine, the precursor substances for the synthesis of

30

exopolysaccharides, was induced by higher levels of c-di-GMP. Consequently, the

31

granulation process was accelerated via EPS production. Higher levels of intracellular

32

hydrophobic amino acids were also positively correlated with increased extracellular

33

protein levels, considering the significant increase in peptides under alternating

34

loadings. Based on our findings, we believe that c-di-GMP regulation and EPS

35

production of the anammox biomass are important mechanisms to enhance its

36

tolerance against unfavorable feed stress. These results highlight the role of c-di-GMP

37

in anammox biomass as it works to survive in unfavorable niches.

38 39

Key words: Anammox; c-di-GMP; 16S rRNA transcripts; Metabolomics; Alternating loadings; Granulation 2

ACS Paragon Plus Environment

Page 3 of 38

Environmental Science & Technology

40

1 Introduction

41

The anaerobic ammonium oxidation (anammox) process is a shortcut in the nitrogen

42

cycle that directly converts nitrite and ammonium to nitrogen gas.1 Many applications

43

of the anammox process have emerged worldwide as alternatives for wastewater

44

treatment because of its high nitrogen removal, low energy consumption, no

45

requirement for external organic carbon, and low sludge yield.2 Despite these

46

significant advantages, application of the anammox process has been limited

47

primarily because of a slow growth rate (doubling time of 11-13 d)3 and the high

48

sensitivity to operational conditions such as dissolved oxygen, organic content,

49

temperature, pH, and especially substrates.4, 5 The fluctuation of substrate

50

concentrations occurs commonly and naturally both in actual industrial and domestic

51

wastewater and in natural environments. Therefore, gaining further insight into the

52

phenotype and metabolic response to the nutrient changes is a subject of great interest

53

for applications involving the anammox process.

54

Anammox bacteria are not currently available as pure cultures;6 they exist as

55

matrix-encased species-rich communities with other species as surface-associated

56

biofilms or surface-independent aggregates in natural or bioreactor habitats.7 Thus,

57

the characteristics of an anammox biofilm will be closely related to the performance

58

of the anammox process. Nutrient concentration or availability, as an important

59

environmental factor, can impact biofilm growth, development, and dispersal behavior.

60

8-10

61

dissolved organic carbon, glucose or nitrogen concentrations decrease 11, and biofilms

It has been demonstrated that Pseudomonas fluorescens biofilms detach when

3

ACS Paragon Plus Environment

Environmental Science & Technology

Page 4 of 38

62

of Pseudomonas putida rapidly disperse in response to carbon starvation.12 For E. coli,

63

when the carbon/nitrogen ratio in the nutrient supply is increased, the extracellular

64

polysaccharide/protein ratio also increases which has a strong effect on its biofilm

65

formation.13

66

For the molecular mechanisms triggered by nutrient changes, the ubiquitous

67

intracellular second messenger c-di-GMP (Bis-(3’-5’)-cyclic dimeric guanosine

68

monophosphate) represents an intriguing area of research, particularly with respect to

69

the behavior and ecology of microbial assemblages.14 It has been shown that an

70

increase in c-di-GMP facilitates the biofilm mode, whereas a decrease results in a

71

switch to dispersal and the planktonic mode of existence.10 For instance, carbon

72

starvation of Pseudomonas putida results in a decrease in c-di-GMP levels, which will

73

induce the protease LapG-mediated cleavage of the surface adhesion LapA and lead

74

to dispersal.15 Similarly, biofilms of Pseudomonas aeruginosa undergo dispersal in

75

response to a sudden decrease or increase in carbon-dependent nutrients by

76

modulating c-di-GMP levels;16 17 the sensor regulator for this process in P.

77

aeruginosa is BdlA,18 a chemotaxis regulator that is affected by c-di-GMP levels. For

78

complex microbial communities, a biofilm composed of species-rich communities

79

represents a completely different organism status compared with a pure culture,

80

including different living styles, metabolic behaviors, and resistance capability. Yang

81

et al. 19 used strategies involving different organic loading rates to achieve accelerated

82

aerobic granulation under the regulation of c-di-GMP. When nutrients were increased

83

suddenly, the constituent cells were stimulated to secrete higher levels of c-di-GMP to 4

ACS Paragon Plus Environment

Page 5 of 38

Environmental Science & Technology

84

produce alginate-like exopolysaccharides (ALE), which served as the precursor for

85

aerobic granules. This study produced meaningful results for wastewater treatment,

86

but the mechanism by which c-di-GMP affects community assembly and the

87

metabolic properties of the community members are still worth being explored.

88

16S rRNA sequencing is a widely used method for exploring bacterial community

89

compositions. Unlike the 16S rRNA gene set from genomic DNA that includes

90

dormant or dead bacteria, the 16S rRNA transcript set from mRNA can be used to

91

analyze the metabolically active community composition.20 In addition, metabolomics

92

has emerged as a technique for defining abundant small molecules from a complex

93

network of chemical and biochemical pathways; the resulting spectrum of masses

94

offers a chemical fingerprint of the microbiota functional status.

95

To improve our understanding of the aggregation mechanism of anammox sludge

96

under nutrient stress, two SBRs were operated for 19 weeks. Particle size, c-di-GMP

97

levels, and EPS contents were subsequently determined, along with reactor

98

performance. Furthermore, short-term batch assays were combined with 16S rRNA

99

transcript sequencing and metabolomics analysis to identify the community assembly

100

mechanism regulated by c-di-GMP.

5

ACS Paragon Plus Environment

Environmental Science & Technology

Page 6 of 38

101

2 Materials and Methods

102

2.1 Bioreactor Operation.

103

To study how c-di-GMP levels and EPS change along with anammox sludge

104

aggregation under alternating feed loadings, two identical sequencing batch reactors

105

(SBRs) with a 1.5-L working volume were operated for 135 days (19 weeks). R1, as

106

control reactor, was operated at conventional stepwise nitrogen loading rates (NLRs)

107

when the concentrations of effluent NO2--N and NH4+-N were both lower than 20

108

mg-N/L (Figure 1A). In accordance with this condition, R2 was operated at a NLR of

109

0.13 kg-N/m3/d, alternating weekly to be higher or lower than that of R1 (Figure 1A).

110

(Figure 1A). Both reactors were fed with the same medium solution mainly consisting

111

of ammonium and nitrite in the form of (NH4)2SO4 and NaNO2 as the nitrogen source.

112

The details of the medium compositions are presented in the Supporting Information21.

113

The reactors were maintained at 37 ± 1 °C, pH of 7.5 to 8.0, and under anaerobic

114

conditions. The concentrations of ammonium, nitrite, and nitrate in the SBRs were

115

measured every two days using American Public Health Association (APHA) standard

116

engineering methods.22 The particle size was determined weekly using a laser particle

117

analysis system (Malvern Mastersizer 2000, UK). Sludge samples were taken from

118

the reactor once a week and immediately stored at -80 °C for the determination of

119

extracellular polymeric substances (EPS) and intracellular c-di-GMP content.

120

2.2 Batch Assays.

121

Batch assays were used to explore the mechanism of aggregation regulated by

122

c-di-GMP. The anammox sludge was first pulverized and sieved through a 100-mesh 6

ACS Paragon Plus Environment

Page 7 of 38

Environmental Science & Technology

123

screen to form a uniform particle size and then harvested by centrifugation at 6,000

124

rpm for 10 min. Batch assays were performed in 100-mL serum bottles with 50 mL of

125

anammox sludge suspensions at a final concentration of 1.45 gVSS/L (VSS: volatile

126

suspended solids). The normal load test (B1) contained 25 mg NO2--N/L and 25 mg

127

NH4+-N/L and the other components shown in Text S1. B2-low and B2-high

128

contained NO2--N/L and NH4+-N/L at concentrations 8 mg/L lower and higher,

129

respectively, than that of B1. Every 30 min, 1/3 of the supernatant was replaced with

130

the same volume of fresh medium to obtain alternating loadings (Figure S2A). The

131

vials were made anoxic by flushing with a gas mixture of N2/CO2 (95%/5%, v/v) for

132

10 min, sealed tightly with rubber caps, and shaken at 150 rpm and 37 ± 1 °C. The

133

assays were performed for 2 h. Supernatant samples were taken every 10 min for

134

nitrogen (NO2-, NH4+, and NO3-) detection. Biomass samples were taken at the

135

beginning, 1 h and 2 h and immediately stored at -80 °C for c-di-GMP extraction, EPS

136

determination, 16S rRNA transcripts sequencing and metabolomics analysis. Each

137

loading pattern was conducted in triplicate.

138

2.3 EPS Extraction and Determination.

139

EPS were extracted using the cation exchange resin (CER) method proposed by Hou

140

et al.23 CER was added at a dosage of 70 g/g VSS. The suspensions were stirred for 3

141

h at 200 rpm and 4 °C. More details are provided in Supporting Information Text S2.

142

After the final centrifugation, the bacteria without EPS were collected for intracellular

143

c-di-GMP extraction. The cellular proteins were extracted using a Bacterial Protein

144

Extraction Kit (Sigma, USA). The extracellular polysaccharides and proteins were 7

ACS Paragon Plus Environment

Environmental Science & Technology

145

determined using the Anthrone method with glucose as a standard 24and the

146

Bicinchoninic Acid (BCA) assay 25, respectively. Each was performed in triplicate.

147

More detailed information is provided in Texts S3 and S4.

148

2.4 Intracellular c-di-GMP Detection.

149

Intracellular c-di-GMP was extracted using ultrasonication and acetonitrile/methanol

150

(50/50, v/v) extraction following the method of Christian Spangler et al. 26 with some

151

modifications. Quantitative analysis of intracellular c-di-GMP concentrations was

152

performed using xanthosine 3’,5’-cyclic monophosphate (cXMP) (Sigma-Aldrich,

153

Steinheim, Germany) as the internal standard. LC-MS/MS analysis was performed

154

using a Dionex Ultimate 3000 UPLC system coupled to a TSQ Quantiva Ultra

155

triple-quadrupole mass spectrometer (Thermo Fisher, CA, USA) and equipped with a

156

heated electrospray ionization (HESI) probe. Detailed procedures are described in

157

Supporting Information Text S5.

158

2.5 16S rRNA-transcript Sequencing Analysis.

159

Sequencing of 16S rRNA transcripts was conducted to analyze the metabolically

160

active community composition under alternating feed loadings. Total RNA was

161

extracted from anammox sludge using a RiboPureTM RNA Purification Kit (Ambion,

162

Life Technologies, Lithuania) according to the manufacturer’s guidelines. Then, the

163

extracted RNA was subjected to DNA removal and complimentary DNA (cDNA)

164

synthesis (Tiangen, Beijing, China). The quality and quantity of RNA were checked

165

using a Nanodrop 2000 spectrophotometer (Thermo Fisher Scientific, USA); DNA

166

contamination was detected by performing PCR reactions using the isolated RNA

Page 8 of 38

8

ACS Paragon Plus Environment

Page 9 of 38

Environmental Science & Technology

167

samples. No PCR products were detected after 27 cycles of PCR (Figure S3). (Figure

168

S3). The cDNA with bacterial 16S rRNA transcripts was amplified with the primers

169

338F and 806R targeting the V3–V4 region.27 High-throughput sequencing was

170

conducted at Majorbio Co., Ltd. (Shanghai, China) using the Illumina MiSeq platform

171

according to the manufacturer’s instructions. The raw 16S rRNA sequences

172

determined by MiSeq sequencing were quality-filtered using QIIME v.1.7.0 with the

173

following criteria: (1) the sequences with Ns, lengths shorter than 200 bp, or average

174

quality scores less than 25 were discarded; (2) two nucleotide mismatches in primer

175

matching and reads containing ambiguous characters were removed; and (3) only

176

sequences that overlapped >10 bp were assembled according to their overlap

177

sequence. The sequences were clustered into operational taxonomy units (OTUs) with

178

97% similarity cutoff using the UPARSE pipeline. 29 The detailed sequencing

179

methods are presented in Text S6. All 16S rRNA sequences from Miseq sequencing

180

have been deposited into the NCBI Sequence Read Archive database under the

181

accession number SRP095288.

182

2.6 LC-MS-based Metabolomic Profiling and Quantification Analysis.

183

Metabolomics analysis was applied to provide deeper insights into the metabolic

184

responses of anammox sludge to environmental alterations. Intracellular metabolites

185

were extracted using ultrasonication and methanol (80%) extraction. Untargeted

186

metabolomic analysis was performed using the Q Exactive Orbitrap (Thermo, CA).

187

Targeted metabolomics profiles were analyzed by TSQ Quantiva (Thermo, CA, USA).

188

Details are presented in Text S7.

9

ACS Paragon Plus Environment

Environmental Science & Technology

Page 10 of 38

189

2.7 Data Analysis.

190

Means and standard deviations were calculated using Microsoft Excel. ANOVA tests

191

30

192

determine the significance and differences between experimental and control groups,

193

and p < 0.05 was considered statistically significant. Pearson correlation coefficients

194

were determined using SPSS software. Prior to analysis, data sets of abundance were

195

normalized by log10 transformation for 16S rRNA transcript sequencing and log2

196

transformation for metabolomics profiling, followed by mean centering.32 Differences

197

in active microbial community composition were visualized through unconstrained

198

non-metric multidimensional scaling (NMDS) ordination, based on Bray-Curtis

199

distance in R with the Vegan package. All the heatmaps in this study were generated

200

in R (Package “pheatmap”) using the Pearson correlation as the distance measure.

201

Principal component analysis (PCA) was performed using R (Package

202

“pcaMethods”).

and Bonferroni corrections 31 were conducted using PASW (SPSS version 18.0) to

203

10

ACS Paragon Plus Environment

Page 11 of 38

Environmental Science & Technology

204

3 Results

205

3.1 Anammox Sludge Granulation Process in SBRs.

206

Two SBRs (R1 and R2) were operated for 135 d (19 weeks) with two different

207

loading strategies (Figure 1A). Three distinct phases were defined according to the

208

developmental trends for particle size (Figure 1B). At the point of inoculation, the

209

biomass had a mean particle diameter of 100 µm (50th percentile distribution). This

210

was followed by different granulation processes due to different nitrogen loading rate

211

(NLR) strategies. For R2 as shown in Figure 1B, weeks 0-4 (phase I) were

212

characterized by a dramatic increase in the mean particle size from 100 to 1100 µm.

213

Granule size remained steady in the range of 1000 to 1100 µm for the subsequent 15

214

weeks (phase II, weeks 4-13; phase III, weeks 13-19).

215

More obvious granulation occurred under alternating loadings. Compared to R2,

216

R1 had an even slower granulation process from 100 to 800 µm (phases I-II, weeks

217

0-13), followed by a maintenance phase (granule size: 772 ± 31 µm, phase III, weeks

218

13-19) (Figure 1B).

219

3.2 c-di-GMP and EPS Production in response to Alternating Feed Loadings.

220

The amount of c-di-GMP in the sludge of R2 increased markedly during phase I,

221

which was normalized to 11.07 ± 0.42 to 36.70 ± 1.39 pmol/mg bacterial protein. For

222

R1, in contrast, there was a smooth uptrend within phases I-II, from 11.07 ± 0.42 to

223

21.96 ± 1.11 pmol/mg bacterial protein (Figure 2). The Pearson analysis results

224

showed that c-di-GMP levels in both reactors were strongly and positively correlated 11

ACS Paragon Plus Environment

Environmental Science & Technology

Page 12 of 38

225

with anammox sludge particle size (as a granulation measure) during the granulation

226

processes (R1: r = 0.8907, p < 0.01 for phases I-II; R2: r = 0.9797, p < 0.01 for phase

227

I) (Table S1). The maintenance phases (no obvious changes in particle size) in R1

228

(phases II-III) and R2 (phase III) were characterized by granule maintenance with a

229

steady trend of particle size development, during which c-di-GMP levels showed a

230

weak correlation and no statistical difference (Table S1).

231

EPS contents (PN, PS, and ratio of PN/PS) increased during the granulation

232

process in R1 and R2, and showed a strong positive Pearson correlation with particle

233

size (Figure 1B and Figure 3, Table S1). Interestingly, EPS contents showed no

234

significant correlation with particle size during the maintenance phases in R1 and R2

235

(r < 0.3469, p > 0.0615; minimum absolute value ‘r’ and maximum ‘p’ values for

236

most PNs, PSs, and PN/PS ratios).

237

Additionally, statistical analysis indicated that the EPS contents correlated with

238

the c-di-GMP levels of anammox biomass in the SBRs (Figure 2 and Figure 3, Table

239

S2). For R2, especially for phase I, PN, PS, and the PN/PS ratio showed strong

240

positive correlations with the c-di-GMP levels (PN: r = 0.9894, p < 0.01; PS: r =

241

0.8975, p < 0.01; PN/PS ratio: r = 0.8289, p < 0.01). However, the c-di-GMP levels

242

and EPS contents were not significantly correlated (maximum ‘r’ absolute value=

243

0.5487) or statistically correlated (minimum ‘p’ = 0.2547) during phases II-III. For R1,

244

phases I-II, as expected, were emphasized. We found that PN and the PN/PS ratio

245

both showed strong positive correlations with the c-di-GMP levels (PN: r = 0.9098, p

246

< 0.01; PN/PS ratio: r = 0.8530, p < 0.01). Even the Pearson correlation coefficient 12

ACS Paragon Plus Environment

Page 13 of 38

Environmental Science & Technology

247

between PS and c-di-GMP was statistically significant (p < 0.01), but not significantly

248

correlated (r = 0.3407).

249

Furthermore, we made attempt to explore the cause and effect between c-di-GMP

250

and EPS. The short-term batch test showed that there was a significant increase in

251

c-di-GMP was produced at the 1 h time point under alternating loadings; EPS

252

production did not increase until the 2 h time point (Figure S4). Although the increase

253

in the amount of EPS in the batch culture was not as high as that observed in the

254

natural granulation course, the EPS increase was statistically significant (B2-high or

255

B2-low vs B1; PS: p < 0.0006; PN: p < 0.0021) and likely to be of biological

256

importance.

257

3.3 Metabolically Active Community Composition Analysis.

258

Sequencing of 16S rRNA transcripts offers support in uncovering the composition of

259

the metabolically active community. A total of 330,314 raw sequences were obtained

260

initially. After filtering for low quality sequences, an average of 35,212 sequences

261

were yielded per sample (n = 9, SD = 2,461). Individual samples contained OTU

262

number from 248 to 307. The NMDS plot, based on sequences at an OTU level with >

263

97% similarity, showed a clear separation of the community composition between

264

normal loading (B1) and alternating loadings (B2-low or B2-high) (Figure 4),

265

revealing obvious changes in the community structure under different loadings. A

266

close examination of the community changes based on 16S rRNA transcript

267

sequencing revealed a total of 46 community members that significantly changed in

13

ACS Paragon Plus Environment

Environmental Science & Technology

Page 14 of 38

268

abundance (p < 0.05) under alternating loadings, as shown in Figure 5. The

269

metabolically active members were divided into two clusters based on their

270

abundance under different feed loading patterns (Figure 5). Under alternating loadings,

271

the members in cluster 1 decreased significantly in abundance compared with those

272

under normal loading (B1). Importantly, it was obvious that 39 out of the 46

273

community members (cluster 2) (over 80%) exhibited an increase in abundance under

274

alternating loadings. From these clusters, we found that metabolically-active

275

community structural composition shifts occurred at the phylum level mainly for

276

Proteobacteria (up regulation, approximately 57% in cluster 2) and Planctomycetes

277

(down regulation) (Table S3).

278

Up-regulation of the proteobacterial branch bacteria, Nitrosomonas europaea, an

279

aerobic ammonium oxidizer (AOB) belonging to the phylum Proteobacteria (Tag 4 in

280

Figure 4), was found to be significantly higher under alternating loadings [2.7-fold

281

change, p = 0.0052 (B2-low) and 2.9-fold change, p = 0.0034 (B2-high)]. Moreover,

282

within cluster 2 as shown in Figure 5, the nitrite-oxidizing bacteria (NOB)

283

Nitrospira_sp. species (Tag 38) [2.5-fold change, p = 0.012 (B2-high)] and

284

Denitrifiers Denitratisoma genus (Tag 7 and Tag 32) presented significantly to be

285

higher under alternating loadings [Tag 7: 2.8-fold change, p = 0.0091 (B2-low) and

286

2.9-fold change, p = 0.027 (B2-high); Tag 32: 3.6-fold change, p = 0.0011 (B2-low)

287

and 3.2-fold change, p = 0.012 (B2-high)].

288

3.4 Metabolomics Response to Alternating Feed Loadings.

14

ACS Paragon Plus Environment

Page 15 of 38

Environmental Science & Technology

289

An untargeted approach was used for the global profiling of metabolites, and 331

290

metabolites belonging to the following Kyoto Encyclopedia of Genes and Genomes

291

(KEGG) metabolic pathways were identified: amino acids, carbohydrates, lipids,

292

peptides, cofactors and vitamins, nucleotides, and xenobiotics (Supporting Data S1).

293

The PCA scatter plot presented a sufficient separation between B1 and B2-low or

294

B2-high (Figure 6A), which indicated that the metabolome associated with alternating

295

loadings (B2-low and B2-high) was drastically different from that with normal

296

loading (B1). In addition, the distance between B2-low and B2-high was close in the

297

PCA scatter plot (Figure 6A), indicating that the two samples harvested from the

298

alternating loading assays (B2-low and B2-high) had similar metabolic profiles.

299

The abundance of amino acids appeared as a polarization state (Figures S5A).

300

The alterations of amino acids in the batch assays coincided with a significant

301

increase in almost all the peptides detected for B2-low, and B2-high (Figures S5B)

302

compared to the levels in B1. The B2-low and B2-high states were characterized by a

303

slight decrease in abundance for the majority of the metabolites in the carbohydrate

304

metabolic pathway. Interestingly, the exceptions were fructose 6-phosphate and

305

UDP-N-acetyl-D-glucosamine, which showed obvious increases in abundance (Figure

306

S5C).

307

To more precisely define the changes in the metabolome identified using the

308

untargeted metabolomics approach, we used a targeted metabolomics approach to

309

measure the levels of amino acids and carbohydrates. The direct measurement of

310

amino acids confirmed these results (Figure 6B). The total amounts of amino acids in 15

ACS Paragon Plus Environment

Environmental Science & Technology

311

B1, B2-low, and B2-high were calculated, and the results showed that significantly

312

higher levels of amino acids were produced in B2-low and B2-high than in B1

313

(B2-low > B1, p = 0.00036; B2-high > B1, p = 0.00032) (Table S4).

314

Page 16 of 38

For the carbohydrate metabolic pathway, fructose 6-phosphate and

315

poly-β-1,6-N-acetylglucosamine (PNAG) were emphasized. Fructose 6-phosphate

316

was found to be up regulated by 21.55% in B2-low and by 50.91% in B3 compared to

317

the level in B1. Similar to fructose 6-phosphate, the amount of

318

UDP-N-acetyl-D-glucosamine also increased by 29.36 % and 19.38 % in B2-low and

319

B2-high, respectively (Figure 6C).

16

ACS Paragon Plus Environment

Page 17 of 38

Environmental Science & Technology

320

4 Discussion

321

4.1 Response of c-di-GMP Levels to Alternating Feed Loadings.

322

In this study, two long-term SBRs with different feed loading strategies were operated,

323

and the intracellular c-di-GMP levels were determined for anammox sludge. Two

324

dramatically different aggregation processes were found under different operational

325

strategies. Moreover, this treatment has been successfully replicated within the first

326

month of SBR operation, indicative of its repeatability (Figure S8). Additionally, the

327

discrepancy within the initial stage of reactor operation appeared once the nitrogen

328

loading changed, indicating that the feed variability was in fact the main factor

329

responsible for inducing the observed differences between the two treatments (Figures

330

2, 3 and S8). It should be emphasized that under alternating feed loadings, anammox

331

biomass formed larger particles in a shorter time interval to reach a stable level

332

(almost 4 weeks vs 13 weeks, R2 vs R1) (Figure 1). The c-di-GMP levels of R2 kept

333

pace with its granulation process, especially during the quick maturation period

334

(phase I) (Figures 1 and 2). Additionally, c-di-GMP levels were weakly correlated in

335

phase III of R1 and in phases II-III of R2. The reason for this result could be that

336

anammox biomass reached a relatively stable granule status to adapt to the

337

environment after long-term development, during which c-di-GMP probably no

338

longer played a major role. Based on this, we hypothesized that the anammox sludge

339

community made use of the second messenger c-di-GMP to regulate its aggregation

340

under alternating loadings.

17

ACS Paragon Plus Environment

Environmental Science & Technology

341

Page 18 of 38

Here the community levels of c-di-GMP signaling and responses to alternating

342

loadings were analyzed using anammox batch assays. An 16S rRNA data set was

343

generated from reverse-transcribed RNA and used for sequencing to obtain insights

344

into the changes in the metabolically active community composition in response to

345

alternating feed loading.20 Within cluster 2 of the community members with a

346

statistically significant change, we found that the Proteobacteria phylum was

347

obviously up regulated (Figure 5). Such shifts could indicate that metabolically active

348

microbial species were most likely to adapt or respond to the nutrient stress. Bacteria,

349

especially those in the proteobacterial branch, contain numerous enzymes involved in

350

c-di-GMP turnover that are used to monitor various environmental and intracellular

351

inputs and adjust c-di-GMP levels in a precise manner 33, 34. In addition, the receptor

352

PilZ domain, which has been well studied, is extensively involved in the

353

proteobacterial branch bacteria 33-35. Hence, it was likely that c-di-GMP levels were

354

increased under alternating loadings by the regulation of certain community

355

compositions, especially of the Proteobacteria. To better understand the up-regulation

356

of proteobacterial branch bacteria, an AOB species, N. europaea, was evaluated.

357

Hydroxylamine, an intermediate in the oxidation of ammonia to nitrite for AOB,36

358

was found to be consumed under alternating loadings (Figure S6). Additionally, N.

359

europaea is capable of growing mixotropically on ammonia and hydroxylamine.

360

Under anoxic conditions, hydroxylamine is oxidized with nitrite as an electron

361

acceptor, and nitrous oxide is produced.37 Therefore, during alternating loadings in the

362

absence of oxygen, AOB are likely to be activated when using hydroxylamine and 18

ACS Paragon Plus Environment

Page 19 of 38

Environmental Science & Technology

363

ammonia as substrates. The hydroxylamine might be released from anammox species

364

through the ladderane lipids that are involved.38 Hydroxylamine-dependent metabolic

365

pathways were demonstrated within Candidatus Brocadia sinica,39 and hydrazine, an

366

intermediate in Candidatus Kuenenia species and Candidatus Jettenia species, was

367

found to be able to diffuse through ladderane lipids though its unusual density.40 Thus,

368

some c-di-GMP-dependent community members could be activated and induce higher

369

levels of c-di-GMP when responding to alternating loadings.

370

Moreover, we also found that some NOB and denitrifiers were presented at

371

significantly higher levels under alternating loadings (Figure 5). The intra-taxon

372

relation between AOB and NOB was syntrophic in which nitrite was released by the

373

former and used by the latter. In addition, the NO3- that accumulated from nitrifiers

374

(AOB and NOB) can be further used by denitrifiers. The coexisting nitrifying and

375

putative heterotrophic bacteria in the anammox biofilm might consume a trace

376

amount of O2 or organic by-products of anammox bacteria, which would

377

consequently establish suitable microenvironments for anammox bacteria to resist the

378

unfavorable conditions.41

379

4.2 Manipulation of EPS Production in response to Alternating Feed Loadings.

380

EPS as a key component of biofilms and granules were characterized and quantified

381

based on extracellular protein (PN) and polysaccharide (PS) contents.42 In this study,

382

EPS production was determined to have a strong positive correlation with granulation,

383

especially for the extracellular protein content during the granular formation phases

19

ACS Paragon Plus Environment

Environmental Science & Technology

Page 20 of 38

384

(phases I-II for R1 and phase I for R2) (Figures 1 and 3, Table S1). It has been

385

proposed that extracellular proteins are likely to be strongly associated with the

386

aggregation of suspended flocs into granules by affecting the surface properties of

387

sludge. 23, 43 Thus, there was a strong correlation during the rapid granulation phases,

388

and weak correlations between EPS contents and granulation were found in phase III

389

of R1 and phases II-III of R2 when the granulation process stopped and reached its

390

maintenance stage.

391

Importantly, it can be concluded that EPS production increased significantly

392

faster under alternating loadings than during normal anammox SBR operation (Figure

393

3). When comparing different feed loadings, it could be concluded that PN production

394

was markedly more rapid and higher in amount at the granulation stage during phase I

395

of R2 (10-fold in amount; week 4 vs week 0) than during phases I-II of R1 (6-fold in

396

amount; week 13 vs week 0). These observations are consistent with a previous report

397

showing that excess EPS form the interiors of aerobic granules that are secreted by

398

bacteria under environmental stress.44 Hence, we hypothesized that the increased

399

production of EPS was likely to induce faster granular formation under unfavorable

400

conditions (that is, alternating loadings in this study).

401

To further examine the levels of the metabolome using the targeted metabolomics

402

approach, amino acids were grouped by their hydrophobicity. The ratio of

403

hydrophobic amino acids to hydrophilic amino acids was determined for each loading

404

pattern (Table S4). The average ratios of B2-low and B2-high were significantly

405

higher than that of B1, and the p values were 0.00018 and 0.00022 relative to B1, 20

ACS Paragon Plus Environment

Page 21 of 38

Environmental Science & Technology

406

respectively. More importantly, as indicated by the metabolomics results for amino

407

acids and peptides (Figures S5), an increasing trend was not found for all amino acids,

408

but there was a polarization state as mentioned above, which indicated that bacterial

409

protein was likely not degraded on a large scale but utilized for synthesis considering

410

the up-regulation of peptides under alternating loadings. Therefore, the levels of

411

cellular proteins, especially hydrophobic amino acid-rich proteins, would increase.

412

Anammox bacteria have a slow growth rate, with a doubling time of 11-13 d.3

413

Thus, after a very short-term batch cultivation of anammox sludge (2 h), the biomass

414

would barely increase. Additionally, the growth might even be slower under

415

alternating loadings because of the unfavorable stress. Thus, intracellular proteins

416

may remain stable due to the negligible cellular growth in the batch assays. Based on

417

that conclusion, a higher level of cellular protein synthesis likely implied that more

418

extracellular proteins were produced. To verify this, the naked bacteria after EPS

419

extraction were subjected to cellular protein detection. Significantly more

420

extracellular proteins for the same amount of total bacterial proteins were produced

421

(Figure S7), indicating that more hydrophobic amino acids were probably used for the

422

synthesis of extracellular proteins. This result was also in accordance with the results

423

of Hou et al. 23 in which higher levels of hydrophobic amino acids in extracellular

424

proteins significantly contributed to the high aggregation ability of anammox sludge.

425

Previous studies have shown that the transcription of genes encoding matrix

426

proteins in Vibrio cholerae was increased when intracellular levels of c-di-GMP were

427

elevated 45-47. In addition, through the manipulation of c-di-GMP metabolism, cAMP 21

ACS Paragon Plus Environment

Environmental Science & Technology

Page 22 of 38

428

signaling was found to affect biofilm formation in V. cholera by directly affecting

429

polysaccharide biosynthesis genes and matrix proteins 48. However, it should be noted

430

that there is no direct evidence that cAMP signaling also participated in the process of

431

anammox sludge. Further studies are needed to clarify whether extracellular proteins

432

were also regulated by c-di-GMP in the anammox community.

433

4.3 Functional Pathways of c-di-GMP for Anammox Biomass Aggregation.

434

Previous studies have shown that c-di-GMP signaling is important for biofilm

435

development, and in some species, this is mediated partly through the regulation of

436

EPS production, 49-51 although, to the best of our knowledge, most previous studies

437

used pure culture laboratory systems. Based on these results, we used the short-term

438

batch test as described above to clarify the potential cause and effect from these data.

439

The c-di-GMP levels under nutrient stress increased significantly first at the 1 h time

440

point of the batch test, followed by an increase in EPS (Figure S4). Hence, the EPS

441

induction observed in the batch experiments was likely to be partly a consequence of

442

the regulation of c-di-GMP.

443

The decrease for most carbohydrates under alternating loadings was a logical

444

response to the nutrient stress and unfavorable growth environment. In addition,

445

fructose 6-phosphate and UDP-N-acetyl-D-glucosamine were emphasized here

446

because they have been found to participate in the production of alginate and PNAG,

447

which are two important exopolysaccharide (PS) components. 52-54 The synchronous

448

increased levels of fructose 6-phosphate and c-di-GMP were confirmed at the end of

22

ACS Paragon Plus Environment

Page 23 of 38

Environmental Science & Technology

449

the batch assays (Figure 6C and Figure S4). Fructose 6-phosphate is a substrate used

450

for the formation of GDP-mannuronic acid, which is the precursor of alginate, an

451

important component of EPS. 53 Notably, the process of the conversion of

452

GDP-mannuronic acid into alginate is regulated by c-di-GMP. 52, 53 At the same time,

453

alginate concentrations also increased significantly under alternating feed loadings

454

(Figure 6D). Similar trends were also determined for UDP-N-acetyl-D-glucosamine.

455

Under the regulation of c-di-GMP, it can act as a monomer in the assembly of PNAG,

456

which is another important exopolysaccharide and an essential intercellular

457

polysaccharide adhesin for adherence and biofilm formation.54 In parallel, PNAG

458

production was about 1.5-fold for B2-low and 2.0-fold for B2-high higher than that of

459

B1 (Figure 6E). Overall, the results above may support the conclusion that c-di-GMP

460

plays an important role in EPS production for anammox biomass.

461

4.4 Significance of this Study and Prospects for the Future.

462

The long-term operation of bioreactors and short-term batch assays combined with

463

16S rRNA transcript sequencing and metabolomics analysis were used for the first

464

time to investigate the aggregation behaviors of anammox sludge under alternating

465

feed loadings. This study provides strong evidence that a more active aggregation

466

process is strongly and positively correlated with c-di-GMP levels and EPS content

467

under alternating loadings. When subjected to nutrient stress, anammox sludge tended

468

to produce higher levels of c-di-GMP through shifts in the metabolically active

469

community, according to the 16S rRNA transcript sequencing results. Metabolomics

470

analysis revealed that hydrophobic amino acids were significantly up regulated and 23

ACS Paragon Plus Environment

Environmental Science & Technology

471

likely to be used for the synthesis of more extracellular proteins in accordance with

472

the characteristics of proteins within anammox sludge EPS. In addition, the EPS

473

production process was probably partly under the regulation of c-di-GMP.

474

Page 24 of 38

c-di-GMP, a ubiquitous second messenger in bacteria, can play an important role

475

in the regulation of the behaviors of anammox sludge in response to nutrient stress, as

476

discussed here. Because anammox is a very promising wastewater treatment process,

477

further studies are needed to improve its performance via c-di-GMP regulation

478

especially under unfavorable conditions.

24

ACS Paragon Plus Environment

Page 25 of 38

Environmental Science & Technology

479

Acknowledgments

480

The authors are grateful to the National Natural Science Foundations of China (No.

481

51308007 and No. 51478006) for financial support. The financial support from

482

Shenzhen Science and Technology Innovation Committee (No.

483

JSGG20160429162015597) should also be highly appreciated.

484

Supporting Information

485

The supporting information contains (1) tables of Pearson correlation analysis results

486

and relative abundance of 16S rRNA samples at phylum level; (2) details of EPS

487

content detection procedures, c-di-GMP detection methods, 16S rRNA transcript

488

sequencing, metabolites detection, hydroxylamine determination, alginate and PNAG

489

detection, and replicated reactor operation; (3) figures illustrating the details of SBRs

490

operation and batch assays, heatmaps and bar graph of metabolites with different

491

pathways, and the results of the replicated reactor operation within the first month; (4)

492

one file (excel format) about the details of untargeted metabolomics results.

25

ACS Paragon Plus Environment

Environmental Science & Technology

Page 26 of 38

493

References

494 495

1.

496 497

2.

498 499 500 501

3.

502 503

4.

504 505 506

5.

507 508 509

6.

510 511

7.

512 513

8.

514 515

9.

516 517 518

10. Mcdougald, D.; Rice, S. A.; Barraud, N.; Steinberg, P. D.; Kjelleberg, S. Should we stay or should

519 520 521

11. Delaquis, P. J.; Caldwell, D. E.; Lawrence, J. R.; Mccurdy, A. R. Detachment ofPseudomonas

522 523

12. Rochex, A.; Lebeault, J. M. Effects of nutrients on biofilm formation and detachment of a

524 525

13. Huang, C. T.; Peretti, S. W.; Bryers, J. D. Effects of medium carbon-to-nitrogen ratio on biofilm

526 527

14. Wood, T. K.; Hong, S. H.; Ma, Q. Engineering biofilm formation and dispersal. Trends Biotechnol.

528 529 530

15. Gjermansen, M.; Nilsson, M.; Yang, L.; Tolker-Nielsen, T. Characterization of starvation-induced

531 532

16. Sauer, K.; Cullen, M.; Rickard, A.; Zeef, L.; Davies, D.; Gilbert, P. Characterization of

Kartal, B.; Kuenen, J. G.; van Loosdrecht, M. C. M. Sewage Treatment with Anammox. Science

2010, 328 (5979), 702-703. Lancet, T. Nitrogen removal with the anaerobic ammonium oxidation process. Biotechnol. Lett.

2013, 35 (8), 1145-1154. van der Star, W. R. L.; Abma, W. R.; Blommers, D.; Mulder, J. W.; Tokutomi, T.; Strous, M.;

Picioreanu, C.; van Loosdrecht, M. C. M. Startup of reactors for anoxic ammonium oxidation: Experiences from the first full-scale anammox reactor in Rotterdam. Water Res. 2007, 41 (18), 4149-63. Jin, R. C.; Yang, G. F.; Yu, J. J.; Zheng, P. The inhibition of the Anammox process: A review.

Chem. Eng. J. 2012, 197, 67-79. Carvajal-Arroyo, J. M.; Sun, W.; Sierra-Alvarez, R.; Field, J. A. Inhibition of anaerobic

ammonium oxidizing (anammox) enrichment cultures by substrates, metabolites and common wastewater constituents. Chemosphere 2013, 91 (1), 22-7. Kindaichi, T.; Awata, T.; Mugimoto, Y.; Rathnayake, R. M.; Kasahara, S.; Satoh, H. Effects of

organic matter in livestock manure digester liquid on microbial community structure and in situ activity of anammox granules. Chemosphere 2016, 159, 300-307. Chen, T. T.; Zheng, P.; Shen, L. D. Growth and metabolism characteristics of anaerobic

ammonium-oxidizing bacteria aggregates. Appl. Microbiol. Biotechnol. 2013, 97 (12), 5575-83. Telgmann, U.; Horn, H.; Morgenroth, E. Influence of growth history on sloughing and erosion

from biofilms. Water Res. 2004, 38 (17), 3671-3684. Wijeyekoon, S.; Mino, T.; Satoh, H.; Matsuo, T. Effects of substrate loading rate on biofilm

structure. Water Res. 2004, 38 (10), 2479-2488.

we go: mechanisms and ecological consequences for biofilm dispersal. Nat. Rev. Microbiol. 2011, 10 (1), 39-50.

fluorescens from biofilms on glass surfaces in response to nutrient stress. Microb. Ecol. 1989, 18 (3), 199-210.

Pseudomonas putida strain isolated from a paper machine. Water Res. 2007, 41 (13), 2885-92.

formation and plasmid stability. Biotechnol. Bioeng. 1994, 44 (3), 329-36.

2011, 29 (2), 87-94.

dispersion in Pseudomonas putida biofilms: genetic elements and molecular mechanisms. Mol. Microbiol. 2010, 75 (4), 815-826.

nutrient-induced dispersion in Pseudomonas aeruginosa PAO1 biofilm. J. Bacteriol. 2004, 186 (21), 26

ACS Paragon Plus Environment

Page 27 of 38

Environmental Science & Technology

533

7312-7326.

534 535

17. Hunt, S. M.; Werner, E. M.; Huang, B.; Hamilton, M. A.; Stewart, P. S. Hypothesis for the role of

536 537

18. Morgan, R.; Kohn, S.; Hwang, S.-H.; Hassett, D. J.; Sauer, K. BdlA, a chemotaxis regulator

538 539

19. Yang, Y. C.; Xiang, L.; Wan, C.; Sun, S.; Lee, D. J. Accelerated aerobic granulation using

540 541 542 543

20. Twomey, K. B.; Alston, M.; An, S. Q.; O'Connell, O. J.; Mccarthy, Y.; Swarbreck, D.; Febrer, M.;

544 545 546

21. van de Graaf, A. A.; de Bruijn, P.; Robertson, L. A.; Jetten, M. S. M.; Kuenen, J. G. Autotrophic

547 548

22. APHA. Standard Methods for the Examination of Water and Wastewater, 20th ed. American

549 550

23. Hou, X.; Liu, S.; Zhang, Z. Role of extracellular polymeric substance in determining the high

551 552

24. Loewus, F. A. Improvement in Anthrone Method for Determination of Carbohydrates. Anal. Chem.

553 554 555

25. Osnes, T.; Sandstad, O.; Skar, V.; Osnes, M.; Kierulf, P. Total protein in common duct bile

556 557 558

26. Spangler, C.; Böhm, A.; Jenal, U. A liquid chromatography-coupled tandem mass spectrometry

559 560 561

27. Derakhshani, H.; Tun, H. M.; Khafipour, E. An extended single-index multiplexed 16S rRNA

562 563 564

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

565 566

29. Edgar, R. C. UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat.

567 568

30. Quinn, G. P.; Keough, M. J., Experimental design and data analysis for biologists. Cambridge

569 570

31. Goeman, J. J.; Solari, A. Multiple hypothesis testing in genomics. Stat. Med. 2014, 33 (11),

571 572

32. van den Berg, R. A.; Hoefsloot, H. C.; Westerhuis, J. A.; Smilde, A. K.; van der Werf, M. J.

nutrient starvation in biofilm detachment. Appl. Environ. Microbiol. 2004, 70 (12), 7418-7425.

essential for biofilm dispersion in Pseudomonas aeruginosa. J. Bacteriol. 2006, 188 (21), 7335-7343.

alternating feed loadings: Alginate-like exopolysaccharides. Bioresour. Technol. 2014, 171, 360-366.

Dow, J. M.; Plant, B. J.; Ryan, R. P. Microbiota and metabolite profiling reveal specific alterations in bacterial community structure and environment in the cystic fibrosis airway during exacerbation. Plos One 2013, 8 (12), e82432.

growth of anaerobic ammonium-oxidizing micro-organisms in a fluidized bed reactor. Microbiology 1996, 142 (8), 2187-2196.

Public Health Association: Washington, DC, USA 1998.

aggregation ability ofanammox sludge. Water Res. 2015, 75, 51-62.

2002, 24 (1), 219-219.

measured by acetonitrile precipitation and a micro bicinchoninic acid (BCA) method. Scand. J. Clin. Lab. Invest. 2009, 53 (7), 757-63.

method for quantitation of cyclic di-guanosine monophosphate. J. Microbiol. Methods 2010, 81 (3), 226-31.

sequencing for microbial community analysis on MiSeq illumina platforms. J. Basic Microbiol. 2016, 56, 321-326.

Fierer, N.; Peña, A. G.; Goodrich, J. K.; Gordon, J. I. QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 2010, 7 (5), 335–336.

Methods 2013, 10 (10), 996-998.

University Press: 2002.

1946-1978.

Centering, scaling, and transformations: improving the biological information content of metabolomics 27

ACS Paragon Plus Environment

Environmental Science & Technology

Page 28 of 38

573

data. BMC Genomics 2006, 7, 142.

574 575

33. Hengge, R. Principles of c-di-GMP signalling in bacteria. Nat. Rev. Microbiol. 2009, 7 (7),

576 577

34. Ryjenkov, D. A.; Simm, R.; Römling, U.; Gomelsky, M. The PilZ Domain Is a Receptor for the

578 579

35. Amikam, D.; Galperin, M. Y. PilZ domain is part of the bacterial c-di-GMP binding protein.

580 581

36. Ahn, Y. H. Sustainable nitrogen elimination biotechnologies: A review. Process Biochem. 2006,

582 583

37. Böttcher, B.; Koops, H. P. Growth of lithotrophic ammonia-oxidizing bacteria on hydroxylamine.

584 585 586

38. van de Graaf, A. A.; de Bruijn, P.; Robertson, L. A.; Jetten, M. S. M.; Kuenen, J. G. Metabolic

587 588 589

39. Oshiki, M.; Ali, M.; Shinyako-Hata, K.; Satoh, H.; Okabe, S. Hydroxylamine-dependent

590 591

40. Qiao, S.; Yin, X.; Tian, T.; Jin, R.; Zhou, J. Hydrazine production by anammox biomass with NO

592 593 594

41. Kindaichi, T.; Tsushima, I.; Ogasawara, Y.; Shimokawa, M.; Ozaki, N.; Satoh, H.; Okabe, S. In

595 596 597

42. Tan, C. H.; Koh, K. S.; Xie, C.; Tay, M.; Zhou, Y.; Williams, R.; Ng, W. J.; Rice, S. A.; Kjelleberg,

598 599

43. Zhang, L.; Feng, X.; Zhu, N.; Chen, J. Role of extracellular protein in the formation and stability

600 601 602

44. Mcswain, B. S.; Irvine, R. L.; Hausner, M.; Wilderer, P. A. Composition and Distribution of

603 604

45. Tischler, A.D.; Camilli, A. Cyclic diguanylate (c-di-GMP) regulates Vibrio cholerae biofilm

605 606

46. Beyhan, S.; Tischler, A.D.; Camilli, A.; Yildiz, F.H. Transcriptome and phenotypic responses of

607 608 609

47. Beyhan, S.; Yildiz, F.H. Smooth to rugose phase variation in Vibrio cholerae can be mediated by a

610 611

48. Fong, J.C.; Yildiz, F.H.; Interplay between cyclic AMP-cyclic AMP receptor protein and cyclic

612

49. Jonas, K.; Melefors, O.; Römling, U. Regulation of c-di-GMP metabolism in biofilms. Future

263-73.

Second Messenger c-di-GMP. J. Biol. Chem. 2006, 281 (41), 30310-4.

Bioinformatics 2006, 22 (1), 3-6.

41 (8), 1709-1721.

FEMS Microbiol. Lett. 1994, 122 (3), 263–266. pathway of anaerobic ammonium oxidation on the basis of 15N studies in a fluidized bed reactor. Microbiology 1997, 143(7), 2415-2421.

anaerobic ammonium oxidation (anammox) by “Candidatus Brocadia sinica”. Environ. Microbiol. 2016, 18 (9), 3133-3143.

reversible inhibition effects. Green Chem. 2016, 18, 4908-4915.

situ activity and spatial organization of anaerobic ammonium-oxidizing (anammox) bacteria in biofilms. Appl. Environ. Microbiol. 2007, 73 (15), 4931-9.

S. The role of quorum sensing signalling in EPS production and the assembly of a sludge community into aerobic granules. ISME J. 2014, 8 (6), 1186-97.

of aerobic granules. Enzyme Microb. Technol. 2007, 41 (5), 551-557.

Extracellular Polymeric Substances in Aerobic Flocs and Granular Sludge. Appl. Environ. Microbiol. 2005, 71 (2), 1051-7.

formation. Mol. Microbiol. 2004, 53, 857–869.

Vibrio cholerae to increased cyclic di-GMP level. J. Bacteriol. 2006, 188, 3600–3613.

single nucleotide change that targets c-di-GMP signalling pathway. Mol. Microbiol. 2007, 63, 995–1007.

di-GMP signaling in Vibrio cholerae biofilm formation. J. Bacteriol. 2008, 190 (20), 6646-6659.

28

ACS Paragon Plus Environment

Page 29 of 38

Environmental Science & Technology

613

Microbiol. 2009, 4 (3), 341-58.

614 615

50. Cotter, P. A.; Stibitz, S. c-di-GMP-mediated regulation of virulence and biofilm formation. Curr.

616 617 618

51. Zhu, B.; Liu, C.; Liu, S.; Cong, H.; Chen, Y.; Gu, L.; Ma, L. Z. Membrane association of SadC

619 620

52. Oglesby, L. L.; Jain, S.; Ohman, D. E. Membrane topology and roles of Pseudomonas aeruginosa

621 622

53. Liang, Z. X. The expanding roles of c-di-GMP in the biosynthesis of exopolysaccharides and

623 624

54. Whiteley, C. G.; Lee, D. J. Bacterial diguanylate cyclases: Structure, function and mechanism in

Opin. Microbiol. 2007, 10 (1), 17-23.

enhances its diguanylate cyclase activity to control exopolysaccharides synthesis and biofilm formation in Pseudomonas aeruginosa. Environ. Microbiol. 2016, 18 (10), 3440-3452.

Alg8 and Alg44 in alginate polymerization. Microbiology 2008, 154 (6), 1605-15.

secondary metabolites. Nat. Prod. Rep. 2015, 32 (5), 663-683.

exopolysaccharide biofilm development. Biotechnol. Adv. 2015, 33 (1), 124–141.

625

29

ACS Paragon Plus Environment

Environmental Science & Technology

Page 30 of 38

626

Figure captions

627

Figure 1. (A) Alternating feed loadings in R1 and R2 during long-term operation; (B)

628

development process of particle size distribution (measured on a volume basis) under

629

different feed loadings. 50th percentile indicates that 50% of total particles were

630

below the corresponding size distribution. The dotted line divides the different

631

developmental phases of granulation, namely phases I-III.

632

Figure 2. HPLC-MS/MS profiling of c-di-GMP extracted from anammox sludge

633

during aggregation under different feed loadings. The data were normalized to the

634

respective sample biomass total cellular proteins. Error bars are indicated as s.e.m. (n

635

= 3, technical replicates). The dotted line separates the different developmental phases

636

of granulation.

637

Figure 3. The concentrations of (A) extracellular proteins (PN) and (B) extracellular

638

polysaccharides components of EPS within different development phases of

639

granulation under different feed loadings. Each EPS component was normalized to the

640

respective sample biomass; (C) the ratio of extracellular proteins to polysaccharides

641

(PN/PS). Error bars are defined as s.e.m. (n = 3, technical replicates). The dotted line

642

separates the different developmental phases of granulation.

643

Figure 4. The NMDS plot based on community members at OUT level with > 97%

644

similarity of normal loading (B1) or alternating loadings (B2-low and B2-high).

645

Figure 5. Unsupervised clustering of community members based on 16S rRNA

646

transcript sequences showing statistically significant changes (p < 0.05) following 30

ACS Paragon Plus Environment

Page 31 of 38

Environmental Science & Technology

647

normal loading (B1) or alternating loadings (B2-low and B2-high). Unsupervised

648

clustering was conducted using Pearson correlation as the distance metric for the

649

purpose of discerning differences between sample classes. The heatmap scale

650

represents the abundance of microbial members normalized by log10 transformation

651

and then mean centering. Tags refer to bacteria at the highest taxonomic resolution

652

identified in the RDP database. Each column represents one biological replicate,

653

which was represented as one colored boxed at the bottom of each column. These

654

community members were colored on the right according to the phylum that they

655

belong to.

656

Figure 6. (A) Principal component analysis (PCA) scatter plot of untargeted

657

metabolomics under different feed loadings in the batch assays. The colored regions

658

in PCA plots indicated 95% confidence region. The abundance of (B) amino acids

659

with statistical significance, (C) fructose 6-phosphate and

660

UDP-N-acetyl-D-glucosamine from targeted metabolomic analysis, (D) alginate

661

concentrations and (E) PNAG detected using WGA-biotin (visualized by

662

chemiluminescence detection) in the batch assays. Error bars are defined as s.e.m. (n

663

= 3, biological replicates). Two-way ANOVA was performed and Bonferroni

664

post-tests were conducted to compare each treatment to normal load (as control)

665

where significant differences are indicated as follows: *p < 0.05 and **p < 0.01.

31

ACS Paragon Plus Environment

Environmental Science & Technology

Page 32 of 38

666 667

Figure 1. (A) Alternating feed loadings in R1 and R2 during long-term operation; (B)

668

development process of particle size distribution (measured on a volume basis) under

669

different feed loadings. 50th percentile indicates that 50% of total particles were

670

below the corresponding size distribution. The dotted line divides the different

671

developmental phases of granulation, namely phases I-III.

32

ACS Paragon Plus Environment

Page 33 of 38

Environmental Science & Technology

672 673

Figure 2. HPLC-MS/MS profiling of c-di-GMP extracted from anammox sludge

674

during aggregation under different feed loadings. The data were normalized to the

675

respective sample biomass total cellular proteins. Error bars are indicated as s.e.m. (n

676

= 3, technical replicates). The dotted line separates the different developmental phases

677

of granulation.

33

ACS Paragon Plus Environment

Environmental Science & Technology

Page 34 of 38

678 679

Figure 3. The concentrations of (A) extracellular proteins (PN) and (B) extracellular

680

polysaccharides components of EPS within different development phases of

681

granulation under different feed loadings. Each EPS component was normalized to the

682

respective sample biomass; (C) the ratio of extracellular proteins to polysaccharides

683

(PN/PS). Error bars are defined as s.e.m. (n = 3, technical replicates). The dotted line

684

separates the different developmental phases of granulation. 34

ACS Paragon Plus Environment

Page 35 of 38

Environmental Science & Technology

685 686

Figure 4. The NMDS plot based on community members at OUT level with > 97%

687

similarity of normal loading (B1) or alternating loadings (B2-low and B2-high).

688

35

ACS Paragon Plus Environment

Environmental Science & Technology

Page 36 of 38

689 690

Figure 5. Unsupervised clustering of community members based on 16S rRNA

691

transcript sequences showing statistically significant changes (p < 0.05) following

692

normal loading (B1) or alternating loadings (B2-low and B2-high). Unsupervised

693

clustering was conducted using Pearson correlation as the distance metric for the

694

purpose of discerning differences between sample classes. The heatmap scale

695

represents the abundance of microbial members normalized by log10 transformation

696

and then mean centering. Tags refer to bacteria at the highest taxonomic resolution

697

identified in the RDP database. Each column represents one biological replicate,

698

which was represented as one colored boxed at the bottom of each column. These

699

community members were colored on the right according to the phylum that they

700

belong to. 36

ACS Paragon Plus Environment

Page 37 of 38

Environmental Science & Technology

701 702

Figure 6. (A) Principal component analysis (PCA) scatter plot of untargeted

703

metabolomics under different feed loadings in the batch assays. The colored regions

704

in PCA plots indicated 95% confidence region. The abundance of (B) amino acids

705

with statistical significance, (C) fructose 6-phosphate and

706

UDP-N-acetyl-D-glucosamine from targeted metabolomic analysis, (D) alginate

707

concentrations and (E) PNAG detected using WGA-biotin (visualized by

708

chemiluminescence detection) in the batch assays. Error bars are defined as s.e.m. (n

709

= 3, biological replicates). Two-way ANOVA was performed and Bonferroni

710

post-tests were conducted to compare each treatment to normal load (as control)

711

where significant differences are indicated as follows: *p < 0.05 and **p < 0.01.

37

ACS Paragon Plus Environment

Environmental Science & Technology

712

Page 38 of 38

TOC/Abstract graphic

713 714

For Table of Contents Only

715

38

ACS Paragon Plus Environment