Relationship of Microbiota and Cyanobacterial ... - ACS Publications

Mar 27, 2017 - Centre for Environmental Sensing and Modelling, Singapore-MIT ... National University of Singapore, 1 Engineering Drive 2, E1A 07-03, ...
1 downloads 0 Views 1MB Size
Subscriber access provided by University of Newcastle, Australia

Article

Relationship of microbiota and cyanobacterial secondary metabolites in Planktothricoides-dominated bloom Shu Harn Te, Boon Fei Tan, Janelle Renee Thompson, and Karina Yew-Hoong Gin Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.6b05767 • Publication Date (Web): 27 Mar 2017 Downloaded from http://pubs.acs.org on March 28, 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 35

Environmental Science & Technology

1

Relationship of microbiota and cyanobacterial

2

secondary metabolites in Planktothricoides-

3

dominated bloom

4

Shu Harn, Te1, Boon Fei, Tan2, Janelle, R. Thompson2, 3and Karina Yew-Hoong, Gin1,4*

5 6

1

7

Drive 1, #02-01 T-Lab Building, Singapore 117411

8

2

9

Technology Centre, 1 CREATE Way, #09-03 CREATE Tower, Singapore 138602

NUS Environmental Research Institute, National University of Singapore, 5A Engineering

Centre for Environmental Sensing and Modelling, Singapore-MIT Alliance for Research and

10

3

11

Massachusetts Avenue, Cambridge, MA 02139-4307

12

4

13

Engineering Drive 2, E1A 07-03, Singapore 117576

Department of Civil and Environmental Engineering, Massachusetts Institute of Technology, 77

Department of Civil and Environmental Engineering, National University of Singapore, 1

14 15

* corresponding author: Karina Yew-Hoong, Gin; [email protected]; +656516 8104 (Tel)

16 1 ACS Paragon Plus Environment

Environmental Science & Technology

17

Page 2 of 35

ABSTRACT GRAPHICS

18

19 20 21 22

KEYWORDS

23

Cyanobacterial bloom, 16S rRNA OTU (operational taxonomy unit), microbial community,

24

Planktothricoides, off-flavours.

25

2 ACS Paragon Plus Environment

Page 3 of 35

Environmental Science & Technology

26

ABSTRACT

27

Identification of phytoplankton species and microbial biodiversity is necessary to assess water

28

ecosystem health and the quality of water resources. We investigated the short-term (2 days)

29

vertical and diel variations in bacterial community structure and microbially-derived secondary

30

metabolites during a cyanobacterial bloom that emerged in a highly urbanized tropical reservoir.

31

The waterbody was largely dominated by the cyanobacteria Planktothricoides spp., together with

32

the Synechococcus, Pseudanabaena, Prochlorothrix and Limnothrix. Spatial differences (i.e.

33

water depth) rather than temporal differences (i.e day vs night) better explained the short-term

34

variability in water quality parameters and bacterial community composition. Difference in

35

bacterial structure suggested a resource driven distribution pattern for the community. We found

36

that the freshwater bacterial community associated with cyanobacterial blooms is largely

37

conserved at the phylum level, with Proteobacteria (β-proteobateria), Bacteroidetes and

38

Actinobacteria as the main taxa despite the cyanobacterial species present and geographical

39

(Asia, Europe, Australia and North America) or climatic distinctions. Through multivariate

40

statistical analyses of the bacterial community, environmental parameters and secondary

41

metabolite concentrations, we observed positive relationships between the occurrences of

42

cyanobacterial groups and off-flavour compounds (2-methyisoborneol and β-ionone), suggesting

43

a cyanobacterial origin. This study demonstrates the potential of 16S rRNA gene amplicon

44

sequencing as a supporting tool in algal bloom monitoring or water resource management.

45 46

3 ACS Paragon Plus Environment

Environmental Science & Technology

Page 4 of 35

47

1.

48

Proliferation of cyanobacteria in freshwater systems has increased in scale and frequency in past

49

decades as a consequence of eutrophication and climate change, becoming one of the major

50

challenges in protecting aquatic ecosystem health and maintaining safety of water resources 1.

51

The production of secondary metabolites and exudates that are toxic to organisms and humans or

52

having unpleasant smells is one of the key concerns related to cyanobacterial blooms 2. These

53

compounds can be synthesized by different taxa or genera 3, 4 which complicates the water

54

management strategy in controlling the sources of these pollutants.

55

Extensive effort has been carried out to explain bloom formation with regards to nutrients and

56

physiochemical factor 5. Recently, more studies have focused on the bacterial community that

57

thrives alongside the cyanobacterial population 6. Indeed, heterotrophic bacteria can be

58

influenced by biotic and abiotic factors induced by blooms such as nutrients, temperature, pH,

59

alkalinity and organic carbon 7, 8. Similarly, heterotrophic bacteria can also alter the water

60

environments 9 leading to positive or negative impacts on blooms. For instance, Flavobacterium,

61

Streptomyces and Oxalobacteraceae can enhance cyanobacterial growth, Pedobacter and

62

Arthrobacter may cause negative impacts 6, 10, while Limnobacter, Rhizobium and Lacibacter

63

develop symbiotic relationships with cyanobacteria 11. Previous studies demonstrated that

64

bloom-associated microbial communities are highly varied at phylum level 8, 12. However, the

65

functional potential revealed by genes could be highly conserved despite differences in

66

community assemblage 12, such as nitrogen fixing / assimilation 8, organic matter degradation 6,

67

11

68

Riding on the development of high-throughput sequencing (HTS), information regarding

69

bacterial groups and microbial diversity in bloom-affected water bodies is fast expanding but the

INTRODUCTION

and iron uptake 13.

4 ACS Paragon Plus Environment

Page 5 of 35

Environmental Science & Technology

70

existing literatures report only a few genera such as Microcystis, Anabaena and

71

Cylindrospermopsis 8, 11, 14, 15. In addition, most of the relevant studies focused on temperate or

72

subtropical regions with very limited investigations done in the tropical climate zone 16. In the

73

current study, we carried out an intensive short-term sampling in one of the reservoirs in an

74

urbanized area of Singapore, during which a rarely reported filamentous cyanobacterium,

75

Planktothricoides, dominated the waterbody. This study aims to: (i) determine the bacterial

76

community associated with Planktothricoides bloom (ii) investigate diel and vertical variations

77

of physiochemical variables, cyanobacterial metabolites and microbial community, and (iii)

78

identify potential metabolite producers using 16S rRNA amplicon sequencing combined with

79

correlation analysis.

80 81 82

2.

MATERIALS AND METHODS

83

2.1.

Site description and sampling

84

This study was conducted in a freshwater reservoir located in a fully urbanized catchment

85

covered mainly by commercial, residential and industrial areas. Having a total catchment area of

86

100 km2, water body surface area of 2.4 km2 and average depth of nearly 4 m, the water body is

87

the largest multi-functional reservoir in Singapore providing recreation, flood control and water

88

storage to the country. Sampling was carried out over two days from 7 to 9 January 2014 at the

89

confluence of three main tributaries feeding water to the reservoir. To study possible diel vertical

90

variations caused by thermal or nutrient stratifications 17, 18, samples were collected at different

91

times of the day (early morning, late morning, afternoon, evening and night, i.e. 6am, 10am,

5 ACS Paragon Plus Environment

Environmental Science & Technology

Page 6 of 35

92

2pm, 6pm, 10pm, 2am) from 2 depths (surface, 0.1 m; bottom, 3.2 m) using a subsurface grab

93

sampler. A 4-hr sampling interval was used to improve detection of fast-degrading off-flavour

94

compounds which may have half-lives on the order of hours to days 19. All samples were kept in

95

sterile polypropylene bottles and transported to the lab immediately for processing or storage.

96

2.2.

97

Temperature, pH, conductivity, total dissolved solid (TDS), salinity and dissolved oxygen (DO)

98

were measured with a YSI 6600 V2 multi-parameter water quality sonde (YSI Inc.). Other water

99

quality including chlorophyll-a (chl-a), total phosphorus (TP), orthophosphate (PO4), total

100

nitrogen (TN), ammonium (NH4), nitrate (NO3), nitrite (NO2), total organic carbon (TOC),

101

dissolved organic carbon (DOC), calcium (Ca) and magnesium (Mg) were determined with

102

APHA standard methods (10200H, 4500-P (B,E), 5310B, 4110B, 3120B) and JIS method (K

103

0102). Light intensity was measured with a LI-250A light meter and LI-192 sensor (LICOR).

104

2.3.

105

Water samples for cyanotoxin analysis were stored at -20 °C upon arrival to the lab.

106

Concentrations of two bloom-associated cyanotoxins, i.e. microcystin and cylindrospermopsin,

107

were measured using the Microcystin/Nodularin (ADDA) and Cylindrospermopsin ELISA kits

108

(Abraxis LLC), respectively. Water samples were sonicated without heat for 15 min to release

109

cell-bound toxin and final development of colour was measured with a microtiter plate reader.

110

Four off-flavour compounds, i.e. 2-methylisoborneol (2-MIB), geosmin, β-cyclocitral and β-

111

ionone, were determined with HS-SPME/GC-MS/MS method using an Agilent GC 7890A

112

coupled with a 7000B Triple Quad series mass spectral detector as described previously 20. All

113

secondary metabolite measurements were run in duplicate.

Physical and chemical water quality variables

Cyanobacterial secondary metabolites

6 ACS Paragon Plus Environment

Page 7 of 35

Environmental Science & Technology

114

2.4.

115

Microbial biomass was harvested within 3 hours after each sampling by filtering reservoir water

116

through a 0.45 µm cellulose nitrate filter (Sartorius) and frozen at -20 °C. DNA was extracted

117

from the filter using a PowerWater ® DNA Isolation Kit (MOBIO) following manufacturer’s

118

instructions.

119

The abundance of bacteria and cyanobacteria in water were inferred from the 16S rRNA gene

120

copy numbers (GNC) in the samples. Total bacterial 16S rRNA gene was quantified using a

121

droplet digital PCR system (ddPCR) (QX200, Bio-Rad), while cyanobacterial 16S rRNA gene

122

was quantified with a StepOnePlus real-time PCR (qPCR) system. Two most commonly reported

123

cyanobacterial bloom genera – Microcystis (MIC) and Cylindrospermopsis (CYL) were also

124

quantified using qPCR assays developed previously 21. Detailed methods are provided in

125

supporting information SI (A). To investigated the major off-flavor producers, we amplified and

126

sequenced partial 2-MIB biosynthesis genes using primer sets reported previously 22, and

127

searched for sequence homologues against the database of NCBI by using BLAST 23.

128

2.5.

129

Isolated genomic DNA (gDNA) was amplified using the universal pyrotaq primer set (926wF

130

and 1392R) targeting the small subunit 16 rRNA variable regions 6 – 8 of bacterial communities

131

24

132

rRNA sequences have been deposited at the National Centre for Biotechnology Information

133

(NCBI) Sequence Read Archive (SRA) under accession numbers SRR3996072-3996095.

134

The standardized biological data (relative abundance of bacterial OTUs and cyanobacterial

135

genera) were square-root-transformed to reduce the dominance of OTUs in high abundances in

Quantitative analysis of bacterial and cyanobacteria abundance

16S rRNA amplicon sequencing and data analyses

. Details on Illumina Miseq sequencing and sequence pre-processing are in SI (B). All 16S

7 ACS Paragon Plus Environment

Environmental Science & Technology

Page 8 of 35

136

analysis, while all environmental data sets were normalized using means and standard deviations

137

of the variables. The similarity (resemblance) of each pair of samples in terms of biological and

138

physico-chemical characteristics was calculated using the Bray-Curtis similarity and Euclidean

139

distance, respectively. Relationships between OTUs and environmental variables were analyzed

140

using RELATE-BEST and a simple model indicating variance explained by selected

141

environmental variables was calculated using the distance-based linear model (DistLM)

142

implemented in Primer-7 25. To detect spatial and temporal differences in the microbial

143

community between samples, permutational multivariate analysis of variance (PERMANOVA)

144

and Analysis of Similarity test (ANOSIM) were performed. Samples were separated into 4

145

groups and compared according to the sampling time and depth – day versus night; surface

146

versus bottom. This was followed by similarity percentage calculations (SIMPER) to indicate

147

OTUs that contributed the most to the variation.

148

Diel and vertical differences in water quality variables and biological variables (qPCR and

149

ddPCR) were determined using paired-T test and general linear model (GML) in the statistical

150

software package PASW18 (SPSS Inc.). Non-normally distributed data sets were transformed

151

with a 2-steps transformation approach 26 prior to the tests. Day and night samplings were

152

determined from the ambient light intensity of which q quantum measurement > 8 umol /m2.s

153

was considered as day time.

154

To investigate changes of bacterial community composition, non-metric multi-dimensional

155

scaling (nMDS) was conducted for sample comparisons based on OTUs, followed by Spearman

156

correlation to determine environmental variables that correlated most highly with the first and

157

second ordination axes. Results were visualized using nMDS and vector plots.

8 ACS Paragon Plus Environment

Page 9 of 35

Environmental Science & Technology

158 159

3.

RESULTS

160

3.1.

Water quality characteristics, total bacterial and cyanobacterial abundance

161

The reservoir was ranked between mesotrophic and eutrophic (18-month chl-a: 7.6 – 99.7 µg/L)

162

states and susceptible to occasional algal bloom outbreaks driven by monsoon changes and

163

nutrient loadings. During the 2-day sampling duration, the weather was cloudy with temperature

164

recorded (26.6) lower than average (28.4 oC), while ambient light intensity was below 1000

165

µmol photon/m2.s around noon time. In general, windy conditions were observed in the day

166

while calm conditions were observed at night.

167

Physicochemical water characteristics measured in this study are summarized in Figure 1, S1 and

168

Table S1. During the study, the waterbody experienced a bloom attributed to a prolonged period

169

of drought lasting for three months from December 2013 to February 2014. Compared to non-

170

bloom periods (Chl-a < 20 µg/L; phytoplankton count < 60,000 cell/mL), an increase in

171

phytoplankton biomass was observed in these three consecutive months (Chl-a from 14.1 to 73

172

µg/L; phytoplankton count from 43,000 to 540,000 cell/mL). For nutrient analyses, the

173

geometric mean concentration of total nitrogen (TN) was consistent at 0.78±1.053 mg/L, while

174

nitrate (NO3) was higher at the surface than bottom waters (paired T-test and GLM, P < 0.05)

175

and the opposite was true for ammonia (NH4). Total phosphorus (TP) showed greater fluctuation

176

between 0.005 and 0.07 mg/L, but orthophosphate was consistently below the detection limit of

177

5 µg/L. An average TN:TP ratio of 32.5 ± 2.23 (geometric) indicated that the reservoir was

178

phosphorus-limited, and trophic status based on Chl-a and TP concentrations still classified the

179

system as eutrophic 27. The water temperature was higher at the surface than bottom throughout

180

the day-night cycle (0.6 °C, paired T-test and GLM, P < 0.01). Light intensity at the surface 9 ACS Paragon Plus Environment

Environmental Science & Technology

Page 10 of 35

181

during mid-day was 35.8 – 184.1 µmol photon/m2.s and decreased with depth with virtually no

182

light at the bottom (3.4 m depth) at mid-day (i.e. 0 – 0.4 µmol photon/m2.s). Physical parameters

183

showing significant vertical variation were turbidity and pH (paired T-test and GLM, P < 0.05).

184

Notable levels of 2-MIB and β-ionone were found in surface and bottom samples (2-MIB: 12-39

185

ng/L, β-ionone: 10-48 ng/L), exceeding their olfactory threshold concentrations of 10 and 7 ng/L,

186

respectively. 2-MIB showed significant depth variation where its concentration at the surface

187

was 32% higher than the bottom.

188

The total bacterial 16S rRNA GCN (including cyanobacteria) showed a stable community

189

abundance ranging from 1.3 – 2.6 × 106 GCN/mL. Cyanobacterial population abundance

190

fluctuated at a slightly greater magnitude that varied from 3.7 × 105 – 1.7 × 106 GCN/mL (Figure

191

2). Significant diurnal differences in phytoplankton and bacterial abundance were observed, with

192

daytime samples showing higher Chl-a, Cyan 16S rRNA and Bact 16S rRNA gene copies

193

(paired T-test and GLM, P < 0.05), consistent with the result from an online multi-parameter

194

sensor probe located near the sampling site (Figure S2). Cyanobacteria, including the targeted

195

populations of Cylindrospermopsis and Microcystis, formed denser assemblages at the top layer

196

(paired T-test and GLM, P < 0.05) but did not exhibit a diurnal trend in abundance. With

197

concentrations ranging from 36 – 1356 and 377 – 3630 GCN/mL respectively,

198

Cylindrospermopsis and Microcystis spp. accounted for not more than 0.5% of the total

199

cyanobacterial population (16S rRNA equivalent), indicating that they were not the main

200

phytoplankton taxa.

201

The raw sequencing output showed good evenness across samples based on the number of reads

202

where all samples had a minimum of 40k raw paired-end reads (Table S2). The average copy of

203

artefactual OTU was 1.39 ± 0.89 (Table S2b) calculated from three blanks sequenced together 10 ACS Paragon Plus Environment

Page 11 of 35

Environmental Science & Technology

204

with samples. Proportions of sequence removed by quality filtering for length & ambiguity,

205

alignment and predicted chimeras were 30%, 5% and 30%, respectively, generating a total of

206

680,581 quality-controlled sequences for 24 samples (Table S2a). We obtained a total of 27,975

207

unique OTUs at 95% cut-off distance, however, approximately 60% of the OTUs were rare taxa

208

contributing less than 5 reads in the whole dataset. Owing to a high percentage of chimeric

209

sequences and large amount of rare OTUs, singletons and doubletons were removed to avoid

210

overestimation of species, following the average copies of artefactual OTU obtained from the

211

blanks. The sample standardization (subsampled to 17,679 reads) resulted in a loss of no more

212

than 1% of reads, which is in agreement with the rarefaction curves (Figure S3) that plateaued

213

after 15,000 reads. A total of 3,577 final OTUs were obtained, of which 3,447 were bacterial

214

OTUs and 15.6% of them were unclassified. OTUs of archaeal and eukaryotic origin were not

215

included in the analysis. 361 out of 3,447 (10.5%) bacterial OTUs were assigned as

216

Cyanobacteria based on Silva reference.

217

Of the 37 phyla detected in the system, phyla with highest relative sequence abundance were

218

Cyanobacteria (39%), Proteobacteria (29%), Bacteroidetes (15%) and Actinobacteria (10%)

219

(Figure 3) comparable to other studies of freshwater systems 28, 29. Other bacterial phyla detected

220

including Planctomycetes, Verrucomicrobia, Chloroflexi, Acidobacteria and Armatimonadetes,

221

in sum contributing 6.6% of the total sequences. Classification at family level revealed nine

222

major non-cyanobacterial taxa (≥ 1% of the total non-cyanobacteria community) belonging to

223

Alphaproteobacteria (Sphingomonadaceae, Rhodobacteraceae), Betaproteobacteria

224

(Comamonadaceae, Burkholderiaceae, Oxalobacteraceae), Bacteroidetes (Chitinophagaceae,

225

Cytophagaceae) and Actinobacteria (Sporichthyaceae, Acidimicrobiaceae). Cyanobacteria was

226

the most abundant phylum in the upper water layers (39%) and co-dominated the lower layers 11 ACS Paragon Plus Environment

Environmental Science & Technology

Page 12 of 35

227

with Proteobacteria, with similar proportions for both. Planktothricoides was the key genus

228

within cyanobacterial populations varying from 39 to 92% in relative abundance. This was

229

followed by Synechococcus, Pseudanabaena, Prochlorothrix and Limnothrix, making up 96.6 to

230

99.2% of the entire phylum.

231

3.2.

232

Richness and diversity (Chao1 and Inverse-Simpson’s diversity indices) calculated from

233

normalized OTUs showed that samples from all four groups (day vs night, surface vs bottom)

234

shared similar richness while higher alpha diversity was found for bottom samples (T-test, P
0.1 mg/L; TN:TP = 1:10) were strongly associated with higher cell

362

biomass and relative abundance of Planktothricoides spp. 4, 8. It is not surprising that our

Planktothricoides bloom

17 ACS Paragon Plus Environment

Environmental Science & Technology

Page 18 of 35

363

sampling site which is located in a tropical urban catchment where water temperatures are warm

364

throughout the year, is of favourable range for Planktothricoides spp. 52. However, a significant

365

connection between phosphorus and total cyanobacteria / Planktothricoides was absent in our

366

study, likely due to low levels of phosphorus measured during the short sampling period.

367

Interestingly, the bloom detected in Brazil coincided with a La Nina period with low

368

precipitation, while a prolonged drought with slightly lower temperature was recorded when

369

samples were taken for this study (rainfall < 50% of average; monthly temperature 26.6°C

370

compared to annual average of 28.4°C). Our findings suggest that Planktothricoides may share

371

features akin to Oscillatoriales which can explain their ecological succession and disappearance

372

in waters.

373 374

4.4.

375

We present here one of the few studies describing freshwater microbial community structure in

376

an equatorial region. To our knowledge, this is also the first study on a bacterial community

377

associated with a Planktothricoides bloom. Given the highly variable freshwater

378

bacterioplankton composition responding to a diverse range of physiochemical and biological

379

variables, a question arises as to whether cyanobacterial blooms would shape the non-

380

cyanobacterial bacterial community. A substantial portion of bacterial taxa recovered from this

381

study are known as typical freshwater or river taxa covering Proteobacteria, Cyanobacteria,

382

Actinobacteria, Bacteriodetes and Verrucomicrobia 28, 29. Generally, bacterial taxa discovered

383

were comparable to those reported in other bloom studies with different dominating

384

cyanobacterial species (Table 3). In those studies, bacterial clones, DGGE bands or OTUs

385

isolated were mainly affiliated with Proteobacteria, Bacteriodetes, Actinobacteria,

Microbial community

18 ACS Paragon Plus Environment

Page 19 of 35

Environmental Science & Technology

386

Verrucomicrobia and Planctomycetes of varied orders or dominance. Comparing our findings

387

with other studies, we conclude that the freshwater bacterial community associated with

388

cyanobacterial blooms are largely conserved at phylum level with Proteobacteria

389

(Betaproteobacteria), Bacteroidetes and Actinobacteria as the main taxa despite differences in the

390

phytoplankton species, geography and climate. We speculate that this could be due to the core

391

functions of the cyanobacterial bloom associated community being conserved across a wide

392

variety of taxa, as demonstrated in several studies 12, 16. Further investigation on the

393

metagenomes of cyanobacterial blooms of other dominant species could shed light on the

394

functional association of the bacterial consortium in cyanobacteria-dominant ecology.

395

Vertical stratification of the bacterial community was also demonstrated, although not to a large

396

extent. Moving down the water column, relative abundance of Sphingomonadaceae

397

(Proteobacteria), Chitinophagaceae, Cytophagaceae (Bacteriodetes) and Acidimicrobiaceae

398

(Actinobacteria) were enriched among non-Cyanobacterial taxa in the community composition

399

(paired T-test, P< 0.05), which is in agreement with their natural habitats of soil, sediment or

400

deep water layers. Further examination revealed that the aquatic microbial community

401

comprised of only several major families (>1%) (Table 3) which made up 65% of the non-

402

cyanobacterial microbial portion. The alphaproteobacteria, Sphingomonadaceae, are common

403

bacteria having widespread distribution in the soil and aquatic systems and even drinking water

404

57

405

cyanobacterial infected waters 11 and are capable of breaking down complex organic compounds

406

and recalcitrant by products released by cyanobacteria, such as microcystin 7. Comamonadaceae

407

was the most abundant Betaproteobacteria (10.1%) in the system. Little is known about this

408

group of bacteria but they exhibit diverse physiological mechanisms including denitrification 59

, and are typically connected to anthropogenic influence 58. They are frequently detected in

19 ACS Paragon Plus Environment

Environmental Science & Technology

Page 20 of 35

409

and degradation of aromatic compounds 60. Although evidence indicates that they co-occur with

410

several filamentous cyanobacteria (e.g. Anabaena circinalis & Phormidium) 15, and are likely

411

involved in the decomposition of Microcystis biomass and production of volatile organic

412

compounds 61, their implication in the current tropical lake is unclear given the scope of study.

413

Chitanophagaceae (Bacteroidetes) was the second most abundant family comprising of major

414

genera classified as sediment or deep water taxa such as Sediminibacterium 62, Lacibacter 63 and

415

Ferruginibacter 64, comparable to the relative abundance found for Comamonadaceae (major

416

Proteobacteria). Members of the phylum, Bacteroidetes, are specialized in decomposition of

417

organic matter in food chains 29. Along with Cytophagaceae and Flavobacteriaceae of the same

418

family group, Chitanophagaceae are chemoorganotrophs that can lyse complex organic

419

compounds such as cellulose 65 and high molecular weight organic matter in water 66,

420

emphasizing their role in regulating the carbon-cycle in aquatic environments. This feature

421

implies a carbon-rich layer near the bottom, which is favourable to these bacterial groups,

422

emphasizing a resource-driven distribution pattern for the bacteria. Actinobacteria is a typical

423

freshwater phylum 29 consisting of globally ubiquitous surface dwellers 67 although some of its

424

members live in bottom waters 68. Association of Actinobacteria with phytoplankton seems to be

425

complicated where their abundance could either be enhanced or decreased by a bloom 10, 69. The

426

Actinobacterial clades, hgcl and CL-500_29, belonging to the most prevalent lineages (acI and

427

acIV) among freshwater Actinobacteria 29, were also the two most abundant groups detected in

428

the current system.

429

It is worth noting that several major bacterial taxa, i.e. hgcl-clade (Actinobacteria),

430

Polynucleobacter (Betaproteobacteria), and Sediminibacterium (Bacteroidetes), detected in our

431

study matched the members of core microbiomes identified in three constructed water bodies of 20 ACS Paragon Plus Environment

Page 21 of 35

Environmental Science & Technology

432

different trophic levels spanning from oligotrophic to eutrophic 7. This supports the cosmopolitan

433

nature of these taxa 29 and suggest that they are less susceptible to eutrophication pressure.

434 435

Funding Sources

436

This research grant is supported by the Singapore National Research Foundation under its

437

Environmental & Water Technologies Strategic Research Programme and administered by PUB,

438

the Singapore’s National Water Agency (Grant number: 1102-IRIS-14-02). BFT and JRT were

439

supported by the National Research Foundation Singapore through the Singapore MIT Alliance

440

for Research and Technology’s (SMART) Center for Environmental Sensing and Modeling

441

(CENSAM) research program.

442

ACKNOWLEDGMENT

443

We thank Yueat Tin Wong and Grace Chan for logistics support and Jin Zhi Lim for off-flavour

444

compound analysis.

445

21 ACS Paragon Plus Environment

Environmental Science & Technology

Page 22 of 35

446

SUPPORTING INFORMATION

447

SI (A). qPCR and ddPCR methods

448

SI (B). 16S rRNA sequence pre-processing

449

Table S1. Water physicochemical characteristics.

450

Table S2a. Summary on numbers sequence and OTU.

451

Table S2b. Artefacts from sequencing blanks.

452

Table S3. PERMANOVA analysis.

453

Table S4. Marginal test for DISTLM analyses.

454

Table S5. Primer and probe sequences.

455

Figure S1. Bloxplot of water physicochemical characteristics measured at day and night times.

456

Figure S2. Cyanobacterial concentration estimated with sensor.

457

Figure S3. Rarefaction curves of samples.

458

Figure S4. Regression between the relative abundances of Microcystis and Cylindrospermopsis

459

estimated by qPCR and 16S rRNA sequence.

460 461

22 ACS Paragon Plus Environment

Page 23 of 35

Environmental Science & Technology

462

REFERENCES

463 464 465 466 467 468 469 470 471 472 473 474 475 476 477 478 479 480 481 482 483 484 485 486 487 488 489 490 491 492 493 494 495 496 497 498 499 500 501 502 503 504 505 506 507 508

1. Paerl, H. W.; Paul, V. J., Climate change: Links to global expansion of harmful cyanobacteria. Water Res. 2012, 46, (5), 1349-1363. 2. Leão, P. N.; Vasconcelos, M. T. S. D.; Vasconcelos, V. M., Allelopathy in freshwater cyanobacteria. Crit. Rev. Microbiol. 2009, 35, (4), 271-282. 3. Sivonen, K.; Jones, G., Cyanobacterial Toxins. In Toxic Cyanobacteria in Water a Guide to their Public Health Consequences, Monitoring and Management, Chorus, I.; Bartram, J., Eds. World Health Organization: 1999; pp 41-111. 4. Jüttner, F.; Watson, S. B., Biochemical and Ecological Control of Geosmin and 2Methylisoborneol in Source Waters. Appl. Environ. Microbiol. 2007, 73, (14), 4395-4406. 5. Paerl, H. W.; Fulton, R. S., Ecology of Harmful Cyanobacteria. In Ecology of Harmful Algae, Granéli, E.; Turner, J. T., Eds. Springer Berlin Heidelberg: Berlin, Heidelberg, 2006; pp 95-109. 6. Berg, K. A.; Lyra, C.; Sivonen, K.; Paulin, L.; Suomalainen, S.; Tuomi, P.; Rapala, J., High diversity of cultivable heterotrophic bacteria in association with cyanobacterial water blooms. ISME J 2008, 3, (3), 314-325. 7. Llirós, M.; Inceo?lu, Ö.; García-Armisen, T.; Anzil, A.; Leporcq, B.; Pigneur, L.-M.; Viroux, L.; Darchambeau, F.; Descy, J.-P.; Servais, P., Bacterial Community Composition in Three Freshwater Reservoirs of Different Alkalinity and Trophic Status. PLoS ONE 2014, 9, (12), e116145. 8. Louati, I.; Pascault, N.; Debroas, D.; Bernard, C.; Humbert, J.-F.; Leloup, J., Structural Diversity of Bacterial Communities Associated with Bloom-Forming Freshwater Cyanobacteria Differs According to the Cyanobacterial Genus. PLoS ONE 2015, 10, (11), e0140614. 9. Havens, K. E., Cyanobacteria blooms: effects on aquatic ecosystems. In Cyanobacterial Harmful Algal Blooms: State of the Science and Research, Hudnell, H. K., Ed. Springer: New York, 2007; Vol. 619, pp 675–732. 10. Liu, F. H.; Lin, G. H.; Gao, G.; Qin, B. Q.; Zhang, J. S.; Zhao, G. P.; Zhou, Z. H.; Shen, J. H., Bacterial and archaeal assemblages in sediments of a large shallow freshwater lake, Lake Taihu, as revealed by denaturing gradient gel electrophoresis. J. Appl. Microbiol. 2009, 106, (3), 1022-1032. 11. Xie, M.; Ren, M.; Yang, C.; Yi, H.; Li, Z.; Li, T.; Zhao, J., Metagenomic analysis reveals symbiotic relationship among bacteria in Microcystis-dominated community. Front. Microbiol. 2016, 7, 1-10. 12. Steffen, M. M.; Li, Z.; Effler, T. C.; Hauser, L. J.; Boyer, G. L.; Wilhelm, S. W., Comparative Metagenomics of Toxic Freshwater Cyanobacteria Bloom Communities on Two Continents. PLoS ONE 2012, 7, (8), e44002. 13. Ploug, H.; Adam, B.; Musat, N.; Kalvelage, T.; Lavik, G.; Wolf-Gladrow, D.; Kuypers, M. M. M., Carbon, nitrogen and O(2) fluxes associated with the cyanobacterium Nodularia spumigena in the Baltic Sea. The ISME Journal 2011, 5, (9), 1549-1558. 14. Eiler, A.; Bertilsson, S., Composition of freshwater bacterial communities associated with cyanobacterial blooms in four Swedish lakes. Environ. Microbiol. 2004, 6, (12), 1228-1243. 15. Woodhouse, J. N.; Kinsela, A. S.; Collins, R. N.; Bowling, L. C.; Honeyman, G. L.; Holliday, J. K.; Neilan, B. A., Microbial communities reflect temporal changes in cyanobacterial composition in a shallow ephemeral freshwater lake. ISME J 2016, 10, (6), 1337-1351. 16. Penn, K.; Wang, J.; Fernando, S. C.; Thompson, J. R., Secondary metabolite gene expression and interplay of bacterial functions in a tropical freshwater cyanobacterial bloom. ISME J 2014, 8, (9), 18661878. 17. Reynolds, C. S., Phytoplankton assemblages and their periodicity in stratifying lake systems. Ecography 1980, 3, (3), 141-159. 18. Mellard, J. P.; Yoshiyama, K.; Litchman, E.; Klausmeier, C. A., The vertical distribution of phytoplankton in stratified water columns. J. Theor. Biol. 2011, 269, (1), 16-30.

23 ACS Paragon Plus Environment

Environmental Science & Technology

509 510 511 512 513 514 515 516 517 518 519 520 521 522 523 524 525 526 527 528 529 530 531 532 533 534 535 536 537 538 539 540 541 542 543 544 545 546 547 548 549 550 551 552 553 554 555

Page 24 of 35

19. Li, Z.; Hobson, P.; An, W.; Burch, M. D.; House, J.; Yang, M., Earthy odor compounds production and loss in three cyanobacterial cultures. Water Res. 2012, 46, (16), 5165-5173. 20. Lim, J. Z. Production of off-flavours and cylindrospermopsin by cyanobacteria in tropical aquatic ecosystem. Unpublished doctoral dissertation, National University of Singapore, Singapore, 2015. 21. Te, S. H.; Chen, E. Y.; Gin, K. Y.-H., Comparison of Quantitative PCR and Droplet Digital PCR Multiplex Assays for Two Genera of Bloom-Forming Cyanobacteria, Cylindrospermopsis and Microcystis. Appl. Environ. Microbiol. 2015, 81, (15), 5203-5211. 22. Wang, Z.; Xu, Y.; Shao, J.; Wang, J.; Li, R., Genes Associated with 2-Methylisoborneol Biosynthesis in Cyanobacteria: Isolation, Characterization, and Expression in Response to Light. PLoS ONE 2011, 6, (4), e18665. 23. Altschul, S. F.; Gish, W.; Miller, W.; Myers, E. W.; Lipman, D. J., Basic local alignment search tool. J. Mol. Biol. 1990, 215, (3), 403-410. 24. Rinke, C.; Lee, J.; Nath, N.; Goudeau, D.; Thompson, B.; Poulton, N.; Dmitrieff, E.; Malmstrom, R.; Stepanauskas, R.; Woyke, T., Obtaining genomes from uncultivated environmental microorganisms using FACS–based single-cell genomics. Nat. Protocols 2014, 9, (5), 1038-1048. 25. Clarke, K.; Gorley, R., PRIMER v7: User Manual/Tutorial. PRIMER-E: Plymouth, 2015. 26. Templeton, G. F., A Two-Step Approach for Transforming Continuous Variables to Normal: Implications and Recommendations for IS Research. Commun. Assoc. Inf. Syst. 2011, 28, 41-58. 27. Carlson, R. E., A trophic state index for lakes. Limnol. Oceanogr. 1977, 22, (2), 361-369. 28. Gabriel, Z.; Byron, C. C.; Miranda, P. K.-v. A.; Ferry, H.; Suk-Kyun, H., Typical freshwater bacteria: an analysis of available 16S rRNA gene sequences from plankton of lakes and rivers. Aquat. Microb. Ecol. 2002, 28, (2), 141-155. 29. Newton, R. J.; Jones, S. E.; Eiler, A.; McMahon, K. D.; Bertilsson, S., A Guide to the Natural History of Freshwater Lake Bacteria. Microbiol. Mol. Biol. Rev. 2011, 75, (1), 14-49. 30. Mur, L. R.; Skulberg, O. M.; Utkilen, H., Cyanobacteria in the Environment. In Toxic Cyanobacteria in Water a Guide to their Public Health Consequences, Monitoring and Management, Chorus, I.; Bartram, J., Eds. World Health Organization: 1999. 31. Chu, Z.; Jin, X.; Yang, B.; Zeng, Q., Buoyancy regulation of Microcystis flos-aquae during phosphorus-limited and nitrogen-limited growth. J. Plankton Res. 2007, 29, (9), 739-745. 32. Sanford, R. A.; Cole, J. R.; Tiedje, J. M., Characterization and Description of Anaeromyxobacter dehalogenans gen. nov., sp. nov., an Aryl-Halorespiring Facultative Anaerobic Myxobacterium. Appl. Environ. Microbiol. 2002, 68, (2), 893-900. 33. Gan, H. M.; Chew, T. H.; Hudson, A. O.; Savka, M. A., Genome Sequence of Novosphingobium sp. Strain Rr 2-17, a Nopaline Crown Gall-Associated Bacterium Isolated from Vitis vinifera L. Grapevine. J. Bacteriol. 2012, 194, (18), 5137-5138. 34. Criminger, J. D.; Hazen, T. H.; Sobecky, P. A.; Lovell, C. R., Nitrogen Fixation by Vibrio parahaemolyticus and Its Implications for a New Ecological Niche. Appl. Environ. Microbiol. 2007, 73, (18), 5959-5961. 35. Hoefman, S.; van der Ha, D.; Boon, N.; Vandamme, P.; De Vos, P.; Heylen, K., Niche differentiation in nitrogen metabolism among methanotrophs within an operational taxonomic unit. BMC Microbiol. 2014, 14, (1), 1-11. 36. Kutvonen, H.; Rajala, P.; Carpén, L.; Bomberg, M., Nitrate and ammonia as nitrogen sources for deep subsurface microorganisms. Front. Microbiol. 2015, 6, 1079. 37. Williams, T. J.; Ertan, H.; Ting, L.; Cavicchioli, R., Carbon and nitrogen substrate utilization in the marine bacterium Sphingopyxis alaskensis strain RB2256. ISME J 2009, 3, (9), 1036-1052. 38. Young, W. F., Taste and odour threshold concentrations of potential potable water contaminants. Water Res. 1996, 30, (2), 331-340.

24 ACS Paragon Plus Environment

Page 25 of 35

556 557 558 559 560 561 562 563 564 565 566 567 568 569 570 571 572 573 574 575 576 577 578 579 580 581 582 583 584 585 586 587 588 589 590 591 592 593 594 595 596 597 598 599 600 601

Environmental Science & Technology

39. Young, C. C.; Suffet, I. H.; Crozes, G.; Bruchet, A., Identification of a Woody-hay Odor-causing Compound in a Drinking Water Supply. Water Sci. Technol. 1999, 40, (6), 273-278. 40. Cotsaris, E.; Bruchet, A.; Mallevialle, J.; Bursill, D. B., The identification of odorous metabolites produced from algal monocultures. Water Sci. Technol. 1995, 31, (11), 251-258. 41. Ikawa, M.; Sasner, J. J.; Haney, J. F., Activity of cyanobacterial and algal odor compounds found in lake waters on green alga Chlorella pyrenoidosa growth. Hydrobiologia 443, (1), 19-22. 42. Ozaki, K.; Ohta, A.; Iwata, C.; Horikawa, A.; Tsuji, K.; Ito, E.; Ikai, Y.; Harada, K.-i., Lysis of cyanobacteria with volatile organic compounds. Chemosphere 2008, 71, (8), 1531-1538. 43. Giglio, S.; Chou, W. K. W.; Ikeda, H.; Cane, D. E.; Monis, P. T., Biosynthesis of 2-Methylisoborneol in Cyanobacteria. Environ. Sci. Technol. 2011, 45, (3), 992-998. 44. Doi, H.; Takahara, T.; Minamoto, T.; Matsuhashi, S.; Uchii, K.; Yamanaka, H., Droplet Digital Polymerase Chain Reaction (PCR) Outperforms Real-Time PCR in the Detection of Environmental DNA from an Invasive Fish Species. Environ. Sci. Technol. 2015, 49, (9), 5601-5608. 45. Bartram, A. K.; Lynch, M. D. J.; Stearns, J. C.; Moreno-Hagelsieb, G.; Neufeld, J. D., Generation of Multimillion-Sequence 16S rRNA Gene Libraries from Complex Microbial Communities by Assembling Paired-End Illumina Reads. Appl. Environ. Microbiol. 2011, 77, (11), 3846-3852. 46. Green, S. J.; Venkatramanan, R.; Naqib, A., Deconstructing the Polymerase Chain Reaction: Understanding and Correcting Bias Associated with Primer Degeneracies and Primer-Template Mismatches. PLoS ONE 2015, 10, (5), e0128122. 47. Sanschagrin, S.; Yergeau, E., Next-generation Sequencing of 16S Ribosomal RNA Gene Amplicons. Journal of Visualized Experiments : JoVE 2014, (90), 51709. 48. Chaudhary, N.; Sharma, A. K.; Agarwal, P.; Gupta, A.; Sharma, V. K., 16S Classifier: A Tool for Fast and Accurate Taxonomic Classification of 16S rRNA Hypervariable Regions in Metagenomic Datasets. PLOS ONE 2015, 10, (2), e0116106. 49. Albertsen, M.; Karst, S. M.; Ziegler, A. S.; Kirkegaard, R. H.; Nielsen, P. H., Back to Basics – The Influence of DNA Extraction and Primer Choice on Phylogenetic Analysis of Activated Sludge Communities. PLOS ONE 2015, 10, (7), e0132783. 50. Yu, Z.; Morrison, M., Comparisons of Different Hypervariable Regions of rrs Genes for Use in Fingerprinting of Microbial Communities by PCR-Denaturing Gradient Gel Electrophoresis. Appl. Environ. Microbiol. 2004, 70, (8), 4800-4806. 51. Rosselli, R.; Romoli, O.; Vitulo, N.; Vezzi, A.; Campanaro, S.; de Pascale, F.; Schiavon, R.; Tiarca, M.; Poletto, F.; Concheri, G.; Valle, G.; Squartini, A., Direct 16S rRNA-seq from bacterial communities: a PCR-independent approach to simultaneously assess microbial diversity and functional activity potential of each taxon. Sci. Rep. 2016, 6, 32165. 52. Suda, S.; Watanabe, M.; Otsuka, S.; Mahakahant, A.; Yongmanitchai, W.; Nopartnaraporn, N.; Liu, Y.; Day, J., Taxonomic revision of water-bloom-forming species of oscillatorioid cyanobacteria. Int. J. Syst. Evol. Microbiol. 2002, 52, 1577-1595. 53. Rohrlack, T.; Christiansen, G.; Kurmayer, R., Putative Antiparasite Defensive System Involving Ribosomal and Nonribosomal Oligopeptides in Cyanobacteria of the Genus Planktothrix. Appl. Environ. Microbiol. 2013, 79, (8), 2642-2647. 54. Guidi-Rontani, C.; Jean, M. R. N.; Gonzalez-Rizzo, S.; Bolte-Kluge, S.; Gros, O., Description of new filamentous toxic Cyanobacteria (Oscillatoriales) colonizing the sulfidic periphyton mat in marine mangroves. FEMS Microbiol. Lett. 2014, 359, (2), 173-181. 55. Andrade, R. d. R.; Giroldo, D., Limnological characterisation and phytoplankton seasonal variation in a subtropical shallow lake (Guaiba Lake, Brazil): a long-term study. Acta Limnol. Brasil. 2014, 26, 442-456.

25 ACS Paragon Plus Environment

Environmental Science & Technology

602 603 604 605 606 607 608 609 610 611 612 613 614 615 616 617 618 619 620 621 622 623 624 625 626 627 628 629 630 631 632 633 634 635 636 637 638 639 640 641 642 643 644 645

Page 26 of 35

56. Wu, X.; Xi, W.; Ye, W.; Yang, H., Bacterial community composition of a shallow hypertrophic freshwater lake in China, revealed by 16S rRNA gene sequences. FEMS Microbiol. Ecol. 2007, 61, (1), 8596. 57. Ren, L.; He, D.; Zeng, J.; Wu, Q. L., Bacterioplankton communities turn unstable and become small under increased temperature and nutrient-enriched conditions. FEMS Microbiol. Ecol. 2013, 84, (3), 614-624. 58. Vaz-Moreira, I.; Nunes, O. C.; Manaia, C. M., Diversity and Antibiotic Resistance Patterns of Sphingomonadaceae Isolates from Drinking Water. Applied and Environmental Microbiology 2011, 77, (16), 5697-5706. 59. Khan, S. T.; Horiba, Y.; Yamamoto, M.; Hiraishi, A., Members of the Family Comamonadaceae as Primary Poly(3-Hydroxybutyrate-co-3-Hydroxyvalerate)-Degrading Denitrifiers in Activated Sludge as Revealed by a Polyphasic Approach. Applied and Environmental Microbiology 2002, 68, (7), 3206-3214. 60. Zhu, D.; Xie, C.; Huang, Y.; Sun, J.; Zhang, W., Description of Comamonas serinivorans sp. nov., isolated from wheat straw compost. International Journal of Systematic and Evolutionary Microbiology 2014, 64, (12), 4141-4146. 61. Li, H.; Xing, P.; Wu, Q. L., Characterization of the bacterial community composition in a hypoxic zone induced by Microcystis blooms in Lake Taihu, China. FEMS Microbiology Ecology 2012, 79, (3), 773784. 62. Qu, J.-H.; Yuan, H.-L., Sediminibacterium salmoneum gen. nov., sp. nov., a member of the phylum Bacteroidetes isolated from sediment of a eutrophic reservoir. Int. J. Syst. Evol. Microbiol. 2008, 58, (9), 2191-2194. 63. Jin, L.; Shin, S.-Y.; Lee, H.-G.; Ahn, C.-Y.; Oh, H.-M., Lacibacter daechungensis sp. nov., isolated from deep freshwater of a reservoir. Int. J. Syst. Evol. Microbiol. 2013, 63, (12), 4519-4523. 64. Jin, L.; Lee, H.-G.; La, H.-J.; Ko, S.-R.; Ahn, C.-Y.; Oh, H.-M., Ferruginibacter profundus sp. nov., a novel member of the family Chitinophagaceae, isolated from freshwater sediment of a reservoir. Antonie van Leeuwenhoek 2014, 106, (2), 319-323. 65. Bailey, V. L.; Fansler, S. J.; Stegen, J. C.; McCue, L. A., Linking microbial community structure to [beta]-glucosidic function in soil aggregates. ISME J 2013, 7, (10), 2044-2053. 66. Kirchman, D. L., The ecology of Cytophaga–Flavobacteria in aquatic environments. FEMS Microbiology Ecology 2002, 39, (2), 91-100. 67. Glöckner, F. O.; Zaichikov, E.; Belkova, N.; Denissova, L.; Pernthaler, J.; Pernthaler, A.; Amann, R., Comparative 16S rRNA Analysis of Lake Bacterioplankton Reveals Globally Distributed Phylogenetic Clusters Including an Abundant Group of Actinobacteria. Appl. Environ. Microbiol. 2000, 66, (11), 50535065. 68. Terkina, I. A.; Drukker, V. V.; Parfenova, V. V.; Kostornova, T. Y., The Biodiversity of Actinomycetes in Lake Baikal. Microbiology 71, (3), 346-349. 69. Allgaier, M.; Grossart, H.-P., Seasonal dynamics and phylogenetic diversity of free-living and particle-associated bacterial communities in four lakes in northeastern Germany. Aquat. Microb. Ecol. 2006, 45, (2), 115-128. 70. Kolmonen, E.; Sivonen, K.; Rapala, J.; Haukka, K., Diversity of cyanobacteria and heterotrophic bacteria in cyanobacterial blooms in Lake Joutikas, Finland. Aquat. Microb. Ecol. 2004, 36, 201-211. 71. Tromas, N.; Fortin, N.; Bedrani, L.; Terrat, Y.; Cardoso, P.; Bird, D.; Greer, C. W.; Shapiro, B. J., Characterizing and predicting cyanobacterial blooms in an 8-year amplicon sequencing time-course. bioRxiv 2016.

646

26 ACS Paragon Plus Environment

Page 27 of 35

647 648

OTU

1 3 4 2 21 51 57 5 9 23 26 649 650

Environmental Science & Technology

Table 1. Bacterial OTUs responsible for dissimilarity between surface and bottom communities calculated using SIMPER. Total average dissimilarity = 34.84 (Bray-Curtis). Only OTUs ≥ 0.3 % contribution were shown.

Classification

With Cyanobacterial OTUs Surface Bottom Diss/SD Av Abund Av Abund

Contrib%

Without Cyanobacterial OTUs Surface Bottom Diss/SD Av Abund Av Abund

Contrib%

Planktothricoides, Cyanobacteria 72.39 60.17 1.4 1.67 Synechococcus, Cyanobacteria 32.61 24 1.5 0.99 hgcl_Clade, Actinobacteria 27.6 27.6 1.4 0.61 27.60 27.6 1.41 0.70 Comamonas serinivorans, Proteobacteria 36.83 38.13 1.49 0.51 36.83 38.13 1.50 0.58 Opitutae, Verrucomicrobia 8.96 12.79 1.45 0.46 8.96 12.79 1.46 0.53 Methylomonas, Proteobacteria 4.14 9.21 2.54 0.45 4.14 9.21 2.53 0.51 Acinetobacter, Proteobacteria 6.6 2.89 0.87 0.42 6.60 2.89 0.87 0.48 Sediminibacterium, Bacteroidetes 20.93 23.91 1.42 0.41 20.93 23.91 1.43 0.47 CL500-29, Actinobacteria 13.85 18.03 1.28 0.39 13.85 18.03 1.28 0.45 Pseudanabaena, Cyanobacteria 12.15 10.36 1.48 0.3 CL500-29, Actinobacteria 9.06 11.64 0.91 0.3 9.06 11.64 0.91 0.34 Abbreviations: Av. – Average; Abund – abundance; Diss – dissimilarity; SD – standard deviation; Contrib% – percentage contribution to the total average dissimilarity.

27 ACS Paragon Plus Environment

Environmental Science & Technology

Page 28 of 35

Table 2. Sequential test results of environmental variables on the composition of bacterial community and cyanobacterial population predicted by distance-based linear model (DistLM) using forward procedure, adjusted R2 criterion and 4999 permutations. Significant values are in bold.

Variable Depth SpCond Ca TP NO3 MIB Turbidity Chl-a Salinity TN NH4 bionone TOC Mcyst Temp Cyr

Bacteria Adj R2 0.0849 0.1516 0.1791 0.1890 0.1978 0.2033 0.2086 0.2137 0.2183 0.2232

SS(trace) 1642.2 1326.8 801.4 577.8 554.5 513.2 505.2 496.9 487.0 484.6

Pseudo-F 3.133 2.730 1.705 1.244 1.207 1.125 1.115 1.103 1.088 1.089

P 0.0002 0.0002 0.0016 0.0854 0.1352 0.2650 0.2954 0.3264 0.3570 0.3674

Prop. 0.1247 0.1007 0.0608 0.0439 0.0421 0.0390 0.0384 0.0377 0.0370 0.0368

Cyanobacteria Adj R2 SS(trace)

Pseudo-F

0.1695

836.5

5.693

0.0050

0.2056

0.5208 0.3960

86.1 337.1

1.015 3.155

0.3820 0.0240

0.0212 0.0829

0.3340 0.4343 0.4639 0.4900 0.5131 0.5203

758.2 235.7 194.3 173.2 155.7 105.3

6.435 2.355 2.048 1.920 1.807 1.240

0.0010 0.0880 0.0990 0.1250 0.1320 0.2930

0.1863 0.0579 0.0477 0.0426 0.0383 0.0259

P

Prop.

Best Solution Adj R2 R2

Best Solution RSS Adj R2 R2 RSS 0.5207 0.22319 0.56093 5784.2 5 0.70828 1187.0 Table 3. Comparison of major bacterial taxa associated with blooms dominated by different cyanobacterial species. 28

ACS Paragon Plus Environment

Page 29 of 35

Environmental Science & Technology

Blooming cyanobacteria taxa

Geographical region, climate

Taxon description

Planktothricoides

Southeast Asia, Proteobacteria tropical/equatorial

Alphaproteobacteria Betaproteobacteria

Bacteroidetes Bacteroidetes Bacteroidetes Actinobacteria Actinobacteria Microcystis

(culture)

Anabaena / Aphanizomenon

Sphingobacteriia Cytophagia Flavobacteriia Actinobacteria Acidimicrobiia

Sphingomonadaceae** Rhodobacteraceae Comamonadaceae Burkholderiaceae* Oxalobacteraceae Chitinophagaceae** Cytophagaceae** Flavobacteriaceae Sporichthyaceae Acidimicrobiaceae**

Proportion Reference in community This study 2.55% 1.52% 10.13% 2.68% 1.51% 9.93% 1.11% 0.97% 5.71% 3.37%

East Asia, Proteobacteria (α,β,γ), Bacteriodetes, Actinobacteria, Firmicutes+ subtropical East Asia, Proteobacteria (β,α,γ,δ), Actinobacteria, Bacteriodetes, Planctomycetes, subtropical Verrucomicrobia+ Southeast Asia, Proteobacteria, Bacteroidetes, Actinobacteria tropical/equatorial Western Europe, Proteobacteria (α,β,γ), Bacteroidetes temperate

11

Northern Europe, cold temperate Northern Europe, cold temperate

70

Verrucomicrobia, Actinobacteria, Chloroflexi, Proteobacteria, Bacteriodetes Bacteroidetes, Proteobacteria (β,α,γ), Actinobacteria, Planctomycetes, Verrucomicrobia, Firmicutes

Aphanizomenon and Western Europe, Proteobacteria , Actinobacteria , Verrucomicrobia, Bacteroidetes,

56

16

8

14

8

29

ACS Paragon Plus Environment

Environmental Science & Technology

Cylindrospermopsis temperate

Page 30 of 35

Planctomycetes, Chloroflexi and Candidate Division OP10

Cylindrospermopsis Western Europe, Proteobacteria (α,β,γ), Bacteroidetes (culture) temperate

8

Microcystis / Oscillatoria / Anabeana

Actinobacteria, Bacteroidetes, Proteobacteria (β,α,γ)

15

Actinobacteria, Proteobacteria (β,α,γ), Verrucomicrobia, Firmicutes

71

Proteobacteria, Bacteroidetes

12

Australia, temperate North America, temperate

Cyanothece / North America, Trichodesmium-like temperate +

Abundance in descending order

30

ACS Paragon Plus Environment

Page 31 of 35

Environmental Science & Technology

Figure 1. Boxplot of water physicochemical characteristics measured at upper and lower water layers (surface and bottom). Unit for parameters: µg/L (Cyn, Mcyst, Chl-a); mg/L (TN, TP, DOC, TOC, NO2, NO3, NH4, Salinity, Ca, Mg, DO); m.S/cm (Cond); g/L (TDS); ng/L (2-MIB, β-cyclocitral, Geosmin, β-ionone); µmol/m2.s (light); °C (Temp), pH unit (pH) and NTU (Turbidity). Asterisks (*) indicate parameters with significant differences determined using paired T-test and GLM, P < 0.05.

31 ACS Paragon Plus Environment

Environmental Science & Technology

Page 32 of 35

Figure 2. Abundance of total bacteria, total cyanobacteria and two common cyanobacteria genera analysed based on molecular quantification methods for different sampling times. First number denotes the day (1 – Day 1, 2 – Day 2), second number denotes the time, D and N represent day and night sampling, respectively.

32 ACS Paragon Plus Environment

Page 33 of 35

Environmental Science & Technology

Figure 3. A) OTU Richness (Chao1) and diversity (Inverse Simpson Index) and B) relative abundance of major bacterial phyla in all samples.

33 ACS Paragon Plus Environment

Environmental Science & Technology

A

Page 34 of 35

B

Figure 4. (A) Non-metric multi-dimensional scaling (nMDS) plot of the bacterial community dissimilarity (Bray-Curtis) between samples. Cyan circles represent surface (0.1 m) samples, while blue squares represent bottom (3.2m) samples, other markers indicate OTUs assigned to different taxa groups:  Planktothricoides (Cyanobacteria); ∅ Other Cyanobacteria;  Proteobacteria;  Bacteroidetes;  Actinobacteria;  Others and | Unclassified. (B) Vector plot of environmental variables correlate (Spearman correlation) significantly to the changes in bacterial community.

34 ACS Paragon Plus Environment

Page 35 of 35

Environmental Science & Technology

16S rRNA Sequencing

80%

60%

40%

20%

y = 0.6497x + 0.0617; R² = 0.8129 0% 0%

20%

40%

60%

80%

qPCR

Figure 5. Regression showing comparison between the relative abundance of cyanobacterial populations in the bacterial community (16S rRNA cyanobacteria / 16S rRNA total bacteria) estimated by qPCR and amplicon sequence classification.

35 ACS Paragon Plus Environment