1.Introduction

Excessive nitrogen input,often in the forms of nitrate and ammonia,causes problems for water quality,resulting in water eutrophication.There is still some debate over whether nitrogen alone is the main driver of these problems,but there is no question that the increase in nitrogen loading degrades the water quality.In the past few years,more and more cities have used reservoirs as their main drinking water source,due to the pollution of groundwater and the cone of groundwater depression(Su et al.,2013).However,nutrient enrichment can cause eutrophication in reservoirs,leading to harmful algal blooms,development of hypoxic zones,and an alteration in ecosystem structure and function(Carpenter et al.,1998;Zhu et al.,2013;Shen et al.,2014).Reservoir water quality has also been affected by endogenous pollutants released into overlying water under anoxic conditions.A number of studies have previously demonstrated that anoxia can re-introduce nitrogen,phosphorus,iron,and manganese from sediments into overlying water layers(Gantzer et al.,2009;Chai et al.,2011).To effectively control the pollutants released from sediments,aerobic conditions must be maintained in the reservoir's hypolimnetic ecosystem(Beutel and Horne,1999).Over the past few decades,water-lifting aeration(WLA)technology has been developed and used effectively to increase dissolved oxygen concentration and to improve the quality of micro-polluted drinking water reservoirs(Cong et al.,2006,2009;Bryant et al.,2011;Gerling et al.,2014).As we all know,increasing nitrogen,which is the main pollutant in reservoir systems,has resulted in water eutrophication in drinking water reservoirs.Therefore,improving the self-purification of nitrogen under aerobic conditions in reservoir systems is a scientific problem that needs to be solved urgently.It is known that nitrogen removal is mainly focused on denitrification under anaerobic conditions,which makes it impractical in natural ecosystems,especially in drinking water reservoirs.However,identification of the first aerobic denitrifying bacteria by Robertson and Kuenen(1983)revealed simultaneous nitrification and denitrification,which made nitrogen removal in open water bodies easier.Subsequently,several aerobic denitrifying bacteria have been isolated from soil(Kim et al.,2008;Miyahara et al.,2010),wastewater(Neufeld et al.,2007;Jiang et al.,2012;Kim et al.,2012),and activated sludge(Su et al.,2001;Heylen et al.,2008),with most of these strains being used to treat water with high carbon and nitrogen loads.In the past few years,there have been reports on aerobic denitrification in the natural environment.These include the study by Gao et al.(2010),which showed that aerobic denitrification occurred in the surface sediment of the permeable Wadden Sea and that by Coban et al.(2015),in which the rates of aerobic denitrification in constructed wetlands were first quantified.However,little research has been conducted on the use of aerobic denitrifying bacteria to bioremediate the reservoir system.

The Zhoucun Reservoir is a large shallow eutrophic reservoir,located in northern China,with a surface area of 8.54 km2,a maximum depth of about 20 m,and a corresponding storage capacity of 84.29×106 m3,which presents typical seasonal stratification characteristics.The thermal stratification is monomictic,and lasts from April to October.Fish breeding in the reservoir started in the 1980s,with the scale expanding rapidly,and the fish breeding area soon reaching 20%of the reservoir area.Net-case fish breeding has thus caused serious deterioration of the water quality in the reservoir.After the clearing up of the net-case,the water quality of the reservoir has improved,but the sediments of the reservoir have been polluted as well.Thermal stratification plays an important role in the changes of the water environment.In particular,in the stratification period(summer),the hypoxia in the bottom water caused by thermal stratification leads to the release of nutrients and reductants from the sediment.The reservoir is a typical large,shallow,eutrophic,and seasonal stratification reservoir of northern China,but there is not yet enough information available about denitrification in the surface sediments of drinking water reservoirs.In recent years,our research group has carried out a series of studies on oligotrophic aerobic denitrifying bacteria.We have not only isolated oligotrophic aerobic denitrifying bacteria from water and surface sediments,but have also explored the nitrogen removal performance in oligotrophic niches and reservoir water.Huang et al.(2012)carried out a detailed study using WLA combined with oligotrophic biofilm,and exhibited fitting nitrogen removal performances.Huang et al.(2015b)isolated and domesticated three strains of oligotrophic aerobic denitrifiers(N299,G107,and 81Y),and explored their nitrogen removal performance in a reservoir system.Zhou et al.(2016a)studied the nitrogen removal characteristics of indigenous aerobic denitrifying bacteria via in situ oxygen enhancement using WLA technology.Meanwhile,we quantified the abundance of aerobic denitrification bacteria through the functional gene(Zhou et al.,2016b)and aerobic denitrification screening medium(Zhou et al.,2016a).Moreover,when the enhancement was completed,in order to explore the nitrogen removal performance of enhanced indigenous aerobic denitrifying bacteria,we continued the in situ enclosure experiment through addition of the nitrogen source(nitrate)(Zhou et al.,2016b).However,the variations in microbial community(especially those involved in nitrogen function)in enhanced and control systems were not clear throughout the enclosure experiment period.

Therefore,it is necessary to explore the changes in microbial community during the nitrogen removal process based on the MiSeq high-throughput sequencing technology.A previous study found perfect nitrogen removal performance(Zhou et al.,2016b).The objectives of this study were(1)to investigate the diversity and composition of the microbial community by using the MiSeq high-throughput sequencing technology;(2)to determine spatial and temporal changes of the microbial community(especially,denitrifying bacteria)during the enclosure experiment period;and(3)to investigate the relationship between the microbial community structure and environmental driving factors throughout the experimental period through redundancy analysis(RDA).The results of this study can enrich our understanding of the relationship between the bacterial community(especially,denitrifying bacteria)structure and environmental factors,and provide references for research into,protection of,and pollution control of drinking water reservoir environments.

2.Materials and methods

2.1.Experimental system

In order to explore the nitrogen removal characteristics of indigenous aerobic denitrifying bacteria via in situ oxygen enhancement and the changes of the microbial community during the nitrogen removal process,we implemented the enclosure experiment system in the Zhoucun Reservoir(34°56′43′N,117°40′54′E)of Zaozhuang City,in Shandong Province.The details of the enhanced and control systems are as follows:

Enhanced system:The enhanced system was used to simulate the operational environment of WLA,and was placed in the Zhoucun Reservoir with a steel pipe(with an inner diameter of 1.0 m).The enclosure system was filled with approximately 11 m3 of raw water in a polyethylene body(with a height of 14.0 m)to simulate natural reservoir conditions.Compressed air was released into the bottom of the enclosure system in the form of small bubbles,which maintained low levels of dissolved oxygen(DO)concentration(4-5 mg/L)through direct mixing and oxygenation(Zhou et al.,2016b).The low DO conditions of the enhanced system lasted one month;after that the indigenous aerobic denitrifying bacteria had been enhanced completely.Next,nitrate was added to the enclosure system to maintain the nitrogen concentration at 2-2.5 mg/L.After two days,the enclosure system reached a stable state.The nitrogen removal performance and the changes in the microbial community of water and surface sediment systems in the enclosure experiment were investigated.

干预实施1学年后，儿童吃水果的频率明显增加。选择“每天都吃水果”的比例从干预前的18.6%提高到50.0%，而选择“基本不吃”的比例从干预前的5.7%下降至0。差异具有显著性(P=0.003)。

2.2.Sample handling and collection

In this study,water samples(2.0 L for water depths of 0.5 m,7.5 m,and 13 m in the enhanced and control systems)and surface sediment samples(50 mL at different periods in the enhanced and control systems)were collected throughout the enclosure experiment.Surface sediment samples were collected at a layer of 0-2 cm using a sterilized Petersen stainless steel grab sampler(Huang et al.,2015a).All samples were kept in an icebox and shipped to the water research laboratory within 6 h.

2.3.DNA extraction

3.4.2.Surface sediment system

2.4.Quantification of nirS and nirK genes abundance

Rarefaction curves,Shannon-Wiener curves,and Rankabundance curves for water and surface sediment systems at 97% similarity revealed that MiSeq high-throughput sequencing could supply enough information to explore microbial community variations.All OTUs with 97%similarity were assigned to eight phyla each for water and surface sediment systems.

2.5.Microbial community analysis

To explore the microbial community composition of water and surface sediment systems in the enhanced and control systems,MiSeq high-throughput sequencing was performed at Shanghai Majorbio Bio-pharm Technology Co.,Ltd.DNA extracted from surface sediment samples(as described above)was amplified by PCR using primers 27F(5′-AGAGTTTGATCCTGGCTCAG-3′) and 338R (5′-TGCTGCCTC CCGTAGGAGT-3′)targeting the V2 regions of bacterial 16S rRNA genes(Ravel et al.,2011).The purified amplicons were sequenced on an Illumina MiSeq platform.The amplicons were sequences shorter than 200 base pairs and low-quality sequences(with a quality score lower than 25)were removed(Quince et al.,2011).The taxonomic classification of effective sequences was determined using the Ribosomal Database Project(RDP).

2.6.Instrumental analysis

The environment parameters of the sampling site were measured in situ at 0.5-m increments using a multi-parameter water quality analyzer(Hydrolab DS5,HACH Company,USA).Specifically,parameters including temperature,dissolved oxygen(DO),pH,oxidation-reduction potential(ORP),electrical conductivity(EC),and Chlorophyll-a(Chl-a)were measured.The water parameters of enclosure experiments were measured using a spectrophotometer(DR6000,HACH Company,USA)(CNEPA,2002).Total nitrogen(TN)and total phosphorus(TP)of surface sediment were determined by persulfate(Zhou et al.,2016a).The moisture content(MC)was determined by weighing,following the method described by Huang et al.(2015a).The aerobic denitrifier concentration(a),in colony-forming units per mL(cfu/mL),was calculated with a plate count(Zhou et al.,2016b).

2.7.Statistical analysis

The rarefaction curves(RC),abundance-based coverage estimators(ACE),and community diversity indices(the Chao richness estimator,the Shannon index,and the Simpson index)were calculated by MOTHUR(Loman et al.,2012).The difference in composition of species was compared with this method(Fouts et al.,2012).Principal component analysis(PCA)and RDA were performed using the Canoco(5.0)project.In our study,the length of the maximum gradient was less than 3.0,and thus RDA was chosen to construct the reˇlationship between species and environmental factors(Smilauer and Lepˇs,2014).One-way analysis of variance(ANOVA)was used to determine the significance of difference(P＜0.05)using the software statistical product and service solutions(SPSS)20.

后来再跟她遇到，是一个颁奖典礼，林建岳办的。那时候他叫我去颁奖给邓丽君，我心里想的是，不要，因为我知道她应该还在恨我。林建岳知道我和邓丽君之间的事，我就问他，她知道我会颁奖给她吗？林建岳说，知道啊。我说，真的吗？如果是这样，那我就去给她颁奖，就当是有个机会可以表达抱歉，弥补一点我的不懂事。林建岳还说，你可以想办法在颁奖环节上给她一个惊喜。

你看毕加索，不同时期的作品变化有多大？黄宾虹前后期绘画，风格迥异，人称白宾虹与黑宾虹，八十岁以后居然还搞变法。

3.Results and discussion

3.1.Variations in environmental factors and functional genes

As shown in Fig.1 and Table 1,Zhou et al.(2016a)studied the nitrogen removal characteristics of indigenous aerobic denitrifying bacteria via in situ oxygen enhancement using WLA technology.The values of temperature and DO concentrations revealed similar variations for enhanced and control systems throughout the enclosure experiment period.

目前，我国建筑行业中可以单独处理建筑智能化系统中出现的问题的人员较少，很多人专业素养不够，少部分人虽然掌握了智能化技术，但在实践中对于具体的智能产品又了解甚少，无法熟练操控，导致即使熟悉设计图纸也无法达到预期要求。除此之外，还存在很多设计人员虽然可以上岗，却未经过专业的培训，无法完全理解智能系统工程中的子项目，设计出的图纸不能十分到位地协调系统之间的运作，智能系统的作用无法得到充分利用，进而无法为建筑施工带来便捷[1]。

The parameters for the entire water layer(0.5 m,5 m,7.5 m,10 m,and 13 m)were analyzed.The temperature showed a decrease from 11.85 ± 0.02°C and 11.92 ± 0.02°C(Nov.24)to 4.10 ± 0.04°C and 4.09 ± 0.11°C(Dec.30),respectively.The DO concentration in the enhanced system ranged from 4.00±0.44 mg/L(Nov.24)to 10.51±0.26 mg/L(Dec.30),owing to the mixed state of the reservoir,whereas that in the control system increased from 8.06±0.14 mg/L(Nov.24)to 10.46±0.65 mg/L(Dec.30)(Fig.1(a)and(b)).The TN concentration of the enhanced system decreased from 2.27±0.05 mg/L(Nov.24)to 1.40±0.05 mg/L(Dec.30)with a 38.24%removal rate,whereas in the control system it increased from 1.29±0.05 mg/L(Nov.24)to 1.42±0.03 mg/L(Dec.30)(Fig.2(a)and(b)).The nitrate concentration decreased from 1.70±0.01 mg/L(Nov.24)to 0.80±0.06 mg/L(Dec.30)in the enhanced system,whereas in the control system it increased from 0.74±0.01 mg/L(Nov.21)to 0.90±0.00 mg/L(Dec.30)(Fig.3(a)and(b)).There was no nitrite accumulation in the enclosure system(Zhou et al.,2016b).Both aerobic denitrifying bacteria(Fig.4(a)and(b))and denitrification functional genes(nirS and nirK)of the enhanced system(Table 1)increased and were higher than those in the controlsystem during the enclosure experiment.Meanwhile,compared with the control system(from 1468 to 2721 μg/L),the sediment of the enhanced system also showed significant nitrogen removal(from 2751 to 2095μg/L)(Zhou etal.,2016b).The total organic carbon(TOC)of the enhanced system also fit with regard to removal performances.

3.2.MiSeq high-throughput sequencing overview of bacteria

The aerobic denitrification microbial compositions for the water and surface sediment systems showed significant differences in the enhanced and control systems during the enclosure experiment(as shown in Table 3).There were several bacteria in the enclosure experiment that exhibited evident spatial-temporal variations.For example,Yan et al.(2015)showed in detail that a large amount of Proteobacteria(such as Hydrogenophaga and Acidovorax)were involved in nitrogen cycling;Heylen et al.(2008)isolated Acidovorax caeni(denitrifying bacteria)from activated sludge in Belgium;Heylen et al.(2006)isolated Dechloromonas(denitrifying bacteria)from activated sludge in a municipal wastewater treatment plant;Rhodocyclaceae(denitrifying bacteria,Thauera and Methylophaga)was isolated from wastewater(Neufeld et al.,2007;Jiang et al.,2012;Kim et al.,2012);Zoogloea sp.299(aerobic denitrifying bacteria)was isolated from reservoirs with perfect nitrogen performances(Huang et al.,2015a);and Patureau et al.(2000)isolated Thiobacillus sp.from environmental samples.

Quantitative PCR was used to estimate the numbers of nirS and nirK copies in water and sediment systems collected during the enclosure experiment.Primers used for nirK(Hallin and Lindgren,1999) quantification were F1aCu,5′-ATYGGCGGVCAYGGCGA-3′,and R3Cu,5′-GCCTCGATCAGRTTRTGGTT-3'.Primers used for nirS(Throb¨ack et al.,2004)quantification were cd3aF,5′-GTSAACGTSAAGGARACSGG-3′, and R3cd, 5′-GASTTCGGRTGSGT CTTGA-3'.Real-time PCR(qPCR)was performed on an Applied Biosystems(ABI)7500 real-time system (Life Technologies,USA)using a qPCR kit,following the method described by Huang et al.(2015a).

3.3.Spatial and temporal changes of microbial community in water and surface sediment systems

As shown in Fig.5 and Fig.6,(X_Y_Z refers to the system_period_location;X=T and C are the enhanced and control systems,respectively;Y=1,2,and 3 are Nov.30,Dec.4,and Dec.30 during the enclosure experiment,respectively;and Z=0.5,7.5,13,and Nare the water depths of0.5 m,7.5 m,13 m,and surface sediment,respectively).In the water system(Fig.5(a)),the main phyla included Proteobacteria(dominant 1,43.99%±11.17%,and P=0.00 for the enhanced system;28.96%±1.67%,and P=0.00 for the control system),Actinobacteria(dominant 2,22.16%±3.83%,and P=0.00 for the enhanced system;24.72%±5.32%,and P=0.00 for the control system),Cyanobacteria(dominant3,19.78%±11.73%for the enhanced system,and 27.34%±3.95%for the control system),Bacteroidetes(dominant 4,4.88%±1.35%for the enhanced system,and 6.71%±2.71%for the control system),and Planctomycetes(dominant 5,3.22%±1.64%for the enhanced system,and 5.77%±1.70%for the control system).

Meanwhile,for the phylum Proteobacteria,β-proteobacteria(dominant 1,25.63%±11.13%for the enhanced system,and 14.76%±0.48%for the control system)was the largest class,with α-proteobacteria at 15.60% ± 1.92%for the enhanced system and 11.27%±2.39%for the control system(Fig.5(b));for the phylum Actinobacteria level(Fig.5(c)),the main classes were Acidimicrobiia(from 5.36%to 1.87%for the enhanced system,and from 8.09%to 1.69%for the control system)and Actinobacteria(18.16%±2.81%for the enhanced system,and 19.80%±2.59%for the control system).Based on the analysis of microbial communities,in the enhanced system,unclassified Comamonadaceae decreased from 17.79%±4.32%(at the beginning)to 2.07%±0.67%(at the end);hgcI clade remained at 15.44%±3.23%(for the whole experiment),CL500-29 marine group decreased from 7.45%±0.53%(at the beginning)to 1.44%±0.40%(at the end),and Lactococcus ranged from 2.03%±0.35%(at the beginning)to 0.71%±0.21%(at the end).

Control system:All the reservoir water except for that in the enclosure system was considered the control system(Zhou et al.,2016b).

2013年，“联盟”成立之前，在领导班子支持下，医院信息化系统被进一步转型，有了进一步的拆分和整合。“并非零散打碎后再重组，核心思想是一体化拆分，将各系统平台拆出。”孙岩国表示，这些动作基于医院班子医疗和质控管理的新思路，“拆完LIS拆PACS，再到电子病历，各自提升后，整体融合入新的、大的系统中运行。”他补充，这也正是相关流程改造后，医院从数年前日门诊量两三千人次，至而今八千人次仍然运行流畅的重要原因。

In the surface sediment system(Fig.6(a)),the main phyla included Proteobacteria(dominant 1,44.22%±2.17%for the enhanced system,which was consistent with the change of the water system,and 31.02%±8.38%for the control system),Bacteroidetes(dominant 2,11.37%±0.34%for the enhanced system,and 17.30%±2.18%for the control system),Chloroflexi(dominant 3,6.61%±1.52%for the enhanced system,and 5.25%±0.03%for the control system),Firmicutes(dominant 4,4.71%±0.20%for the enhanced system,and 4.61%±0.29%for the control system),and Actinobacteria(dominant 5,2.62%±0.27%for the enhanced system,and 3.88%±1.78%for the control system).

Meanwhile,for phylum Proteobacteria(Fig.6(b)),β-proteobacteria was the main class(21.40%±2.53%for the enhanced system,and 5.81%-11.70%for the control system),which was consistent with the water system.The δ-proteobacteria reached 13.51%and 13.41%for the enhanced and control systems,respectively.In the enhanced system,vadin HA17 norank(dominant 1)remained at 6.15%±0.61%,GIF9 norank decreased from 5.29%to 2.29%,and Ferribacterium increased from 4.74%to 9.35%.In the control system,the main bacteria both showed a decrease(vadin HA17 norank,from 14.81%to 9.63%,and Candidate division OP8 norank,from 8.07%to 2.66%).The main classes were Acidimicrobiia(Fig.6(c))(from 0.24%to 41%for the enhanced system,and from 0.21%to 0.69%for the control system),and Actinobacteria(from 0.94%to 0.81%for the enhanced system,and from 0.84%to 1.73%for the control system).

自2014年下半学期开始,对水处理综合实验教学大纲进行修订,纳入探究式实验,努力适应解决复杂工程问题的工程教育认证标准。传统实验的数量从5次减少到3次,并要求必须有一个探究式试验。对于探究式实验,要求每一组确定一个问题,并对此问题设计解决方案和开展实验,记录结果、提交个人的实验报告。探究式实验成绩占水处理综合实验课程最终成绩的40%,是终极的实验室训练。

①我国地下水、地表水、湖泊水体都受到了不同程度的有机氯农药的污染，但都远低于地表水和地下水环境质量标准的限值，其中地下水的污染程度相对较低。水体中DDTs和HCHs是主要的OCPs污染物。

3.4.Changes of denitrification microbial diversity

The microbial diversity and community structure of water and surface sediment systems were investigated using MiSeq high-throughput sequencing.After removing low-quality sequences,a total of 184024 sequences with an average length of 315 base pairs and 65945 sequences with an average length of 327 base pairs for water and surface sediment systems,respectively,were obtained during the enclosure experiment.After this,quality trimming for 20 samples(15 samples for the water system and five samples for the surface sediment system)was carried out in triplicate.Each sample generated 334 to 723 and 1262 to 2504 operational taxonomic units(OTUs)for water and surface sediment systems,respectively.The coverage ranged from 0.9782 to 0.9875 and from 0.9186 to 0.9518 for water and sediment systems,respectively,at 97%similarity(Table 2).The largest number of OTUs was found in the enhanced system for water and surface sediment systems throughout the experimental period.The abundance-based coverage estimator(ACE)and Chao diversity estimators for water and surface sediment systems presented significant seasonal variation throughout the experiment period,and the diversity of water and surface sediment systems in the enhanced system were evidently higher than those in the control system.However,the Shannon richness and Simpson richness of water and surface sediment systems exhibited a similar trend throughout the experiment period.

3.4.1.Water system

The denitrifying bacteria of the enhanced system(Acidovorax,Bacillus,Dechloromonas,Hydrogenophaga,Pseudomonas,Rhodocyclaceae,Sphingomonas,Thiobacillus,and Zoogloea)mostly decreased during the enclosure experiment.Sphingomonas was seen as an exception,displaying an increase from 0.11%to 0.43%(Tables 3 and 4)during the enclosure experiment.The denitrifying bacteria showed a decrease(Bacillus:from 0.49%to 0.17%;Pseudomonas:from 0.47%to 0.21%)in the control system(Tables 3 and 4).

In order to obtain the total DNA,2 L of water sample,filtered with a 0.22-μm cellulose-acetate filter,was collected,and approximately 50 mL of surface sediment(with a depth of 0-2 cm)was collected as well.Microbial DNA was extracted from sediment samples using the E.Z.N.A.®Soil DNA Kit(Omega Bio-tek,Norcross,GA,USA)according to the manufacturer's protocols.DNA was purified using the AxyPrep DNA Gel Extraction Kit(Axygen Biosciences,Union City,CA,USA)according to the manufacturer's instructions.The extracted DNA was stored at-80°C for polymerase chain reaction(PCR)amplification analysis.

The changes in denitrifying bacteria have an important influence on the nitrogen cycle. However, aerobic denitrification functional bacteria are rarely cultured and studied in the natural ecosystem,especially in water and sediment systems of a drinking water reservoir.Therefore,studies(such as those based on high-throughput GeoChip functional gene microarray analysis,denitrification activity,denitrifying enzymes,and nitrogen metabolic genes)involved in the mechanism of nitrogen removal are necessary to further investigate the community structure and function of denitrifying bacteria.

In the enhanced system(Tables 3 and 4),Acidovorax(from 0.05%to 0.12%),Dechloromonas(from 1.50%to 2.76%),Hydrogenophaga(from 0.01%to 0.28%),Rhodocyclaceae(from 0.66%to 0.87%),and Zoogloea(from 0 to 0.03%)presented an increase throughout the enclosure experiment period,whereas the abundance of Thiobacillus showed minimal change(4.20%±0.39%).In the control system(Tables 3 and 4),variations in denitrifying bacteria presented the same trend.Dechloromonas increased from 0.17%to 0.91%,Bacillus increased from 0.13%to 0.26%,and Thiobacillus increased from 0.32%to 1.51%.

3.5.Microbial community structure comparison

In order to explore the bacterial community relationship between water and surface sediment samples during the enclosure experiment period in the Zhoucun Reservoir,community structure comparisons were carried out with PCA for different samples during different periods based on Bray-Curtis dissimilarity(Ligi et al.,2014).The PCA results for water and surface sediment systems are shown in Fig.7(a)and(b),respectively.Results revealed that the first two principle components(PCA axis 1 and PCA axis 2)explained 63.13%and 64.16%of the variability for water and surface sediment systems,respectively.The accumulated contribution ratios of PCA axis 1 and PCA axis 2 for the water system achieved 41.62%and 21.51%,respectively(Fig.7(a)),whereas for the surface sediment system,the accumulated contribution ratios of PCA axis 1 and PCA axis 2 achieved 43.74%and 20.42%,respectively(Fig.7(b)).Microflora of water and sediment systems was well separated in different samples throughout the experiment period.

Particularly in the water system,samples from the same system showed tighter clustering,whereas others were relatively distributed.The samples of the enhanced system were mostly located on the negative side of PCA axis 2(from quadrant 3 to quadrant 4),whereas the samples of the control system were located on the positive side of PCA axis 2(from quadrant 2 to quadrant 1)(Fig.7(a)).Moreover,the variations of all the water samples throughout the enclosure experiment could be further explored using the PCA axis 1 and PCA axis 2 values of the water samples in enhanced and control systems.In the sediment system,the samples of the enhanced system were all located on the negative side of PCA axis 1(from quadrant 2 to quadrant 3),whereas the samples of the control system were located on the positive side of PCA axis 1(from quadrant 1 to quadrant 4)(Fig.7(b)).The low PCA axis 2 and high PCA axis 1 values suggest that samples of different systems throughout the experiment period in the Zhoucun Reservoir may be mainly dominated by PCA axis 1.Meanwhile,the samples of sediment systems exhibited evident distribution in the enclosure experiment period based on the values of PCA axis 1 and PCA axis 2.

The results of Venn analysis(Fouts et al.,2012)in the enhanced and control systems suggested different and rapidly changing bacterial communities throughout the experiment period.The results showed that the numbers of shared genera for the enhanced and control systems showed a decrease in the water system(from 361 OTUs to 277 OTUs in the enhanced system,and from 323 OTUs to 202 OTUs in the control system).Meanwhile,in the sediment system,the number of shared genera for the enhanced system ranged from 1505 OTUs(at the beginning)to 1451 OTUs(in the middle)to 1087 OTUs(at the end),whereas for the control system they increased from 1452 OTUs(at the beginning)to 1962 OTUs(at the end).Hierarchical clustering based on OTU information was also generated.The samples of water and surface sediment for the enhanced and control systems from the same period were tightly grouped,and could be well separated in different periods.

3.6.Environmental effect on microbial community

RDA was used to study the relationship between the microbial community and environmental variables based on genus level(Fig.8),which can provide the missing link between diversity and activity.In detail,the separation locations of water samples and surface sediment samples of the enhanced and control systems had a better resolution at the genus.Therefore,the different bacterial communities for the enhanced and control systems were well discriminated at the genus level(Fig.8(a)and(b)).

For the water system(Fig.8(a)),the first two RDA dimensions according to the 16 parameters explained 62.69%of the microbial community variation.The accumulated contribution ratios of RDA axis 1 and RDA axis 2 could explain 41.85%and 20.84%,respectively.The lower RDA axis 2 and higher RDA axis 1 scores suggest that samples of different systems throughout the experiment period may be mainly dominated by RDA axis 1.As shown in Table 5,physical and chemical parameters including temperature,pH,DO,nitrate,TN,total dissolved nitrogen(TDN),dissolved organic carbon(DOC),Mn,TP,and nirS significantly influenced water bacterial community composition(Fig.8(a)).For the surface sediment system(Fig.8(b)),the first two RDA dimensions explained 72.55%of the bacterial community variation(RDA axis 1:52.63%;RDA axis 2:19.92%;Fig.8(b)),and RDA analysis revealed that temperature,sediment total nitrogen(STN),sediment total phosphorus(STP),and nirK were the critical environmental factors influencing the spatial and temporal variations of the bacterial community.

In order to explore the effects of environmental factors on nitrogen-function(N-function)microbial community,RDA was carried out during the enclosure experiment period.The sample-environmental variations of water and surface sediment systems(Fig.9(a)and(b))were explored,respectively.The first two RDA dimensions could explain 84.31%and 93.73%for the water and surface sediment systems of the total fitting cumulative variable,respectively,indicating that these results are credible.Meanwhile,the samples of the same system showed tighter clustering,whereas the others were relatively distributed.The lower RDA axis 2 and higher RDA axis 1 scores suggest that samples of different systems throughout the experiment period may be mainly dominated by RDA axis 1.As shown in Fig.9(a)and(b)(sample-environmental variation),the enhanced system and control system were located on the positive and negative sides of RDA axis 1(water)and the negative and positive sides of RDA axis 1(sediment),respectively.It is clear that the environment variations(TN,DOC,nirK,and nirS)were the most important factors affecting the N-function bacterial community and composition for water and sediment systems,respectively.

4.Conclusions

The results showed perfect nitrogen removal and pollutant inhibition in the enhanced system.The microbial diversity and community structure of the water and surface sediment systems were investigated using MiSeq high-throughput sequencing,respectively.The microbial community of the water and surface sediment systems presented spatial and temporal changes in the enhanced and control systems,respectively.RDA was carried out throughout the enclosure experiment period.Based on the results,temperature,DO,TN,nitrate,nirK,nirS,and DOC were shown to play key roles in affecting bacterial community function and composition.The changes in the microbial community during the nitrogen removal process were investigated successfully with the MiSeq high-throughput sequencing technology.This can provide a theoretical basis and technical support for bioremediation of micro-polluted reservoir water.

References

Beutel,M.W.,Horne,A.J.,1999.A review of the effects of hypolimnetic oxygenation on lake and reservoir water quality.Lake Reserv.Manage.15(4),285-297.https://doi.org/10.1080/07438149909354124.

Bryant,L.D.,Hsu-Kim,H.,Gantzer,P.A.,Little,J.C.,2011.Solving the problem at the source:Controlling Mn release at the sediment-water interface via hypolimnetic oxygenation.Water Res.45(19),6381-6392.https://doi.org/10.1016/j.watres.2011.09.030.

Carpenter,S.R.,Caraco,N.F.,Correll,D.L.,Howarth,R.W.,Sharpley,A.N.,Smith,V.H.,1998.Nonpoint pollution of surface waters with phosphorus and nitrogen.Ecol.Appl.8(3),559-568.https://doi.org/10.1890/1051-0761(1998)008[0559:NPOSWW]2.0.CO;2.

Chai,B.B.,Huang,T.L.,Zhu,W.H.,Yang,F.Y.,2011.A new method of inhibiting pollutant release from source water reservoir sediment by adding chemical stabilization agents combined with water-lifting aerator.J.Environ.Sci.23(12),1977-1982.https://doi.org/10.1016/S1001-0742(10)60661-4.

Chen,F.,Xia,Q.,Ju,L.K.,2003.Aerobic denitrification of Pseudomonas aeruginosa monitored by online NAD(P)H fluorescence.Appl.Environ.Micro.69(11),6715-6722.https://doi.org/10.1128/AEM.69.11.6715-6722.2003.

Chinese National Environmental Protection Agency(CNEPA),2002.Water and Wastewater Monitoring Methods.China Environmental Science Press,Beijing(in Chinese).

Coban,O.,Kuschk,P.,Kappelmeyer,U.,Spott,O.,Martienssen,M.,Jetten,M.S.M.,Knoeller,K.,2015.Nitrogen transforming community in a horizontal subsurface- flow constructed wetland.Water Res.74(1),203-212.https://doi.org/10.1016/j.watres.2015.02.018.

Cong,H.B.,Huang,T.L.,Zhao,J.W.,Zhou,Z.M.,He,W.J.,Han,H.D.,2006.Application of the technology of lifting water and aeration for improving water quality.Environ.Pollut.Contr.28(3),215-218.https://doi.org/10.15985/j.cnki.1001-3865.2006.03.016.

Cong,H.B.,Huang,T.L.,Chai,B.B.,Zhao,J.W.,2009.A new mixingoxygenating technology for water quality improvement of urban water source and its implication in a reservoir.Renew.Energy 34(9),2054-2060.https://doi.org/10.1016/j.renene.2009.02.007.

Fouts,D.E.,Szpakowski,S.,Purushe,J.,Torralba,M.,Waterman,R.C.,MacNeil,M.D.,Alexander,L.J.,Nelson,K.E.,2012.Next generation sequencing to define prokaryotic and fungal diversity in the bovine rumen.PLOS One 7(11),e48289.https://doi.org/10.1371/journal.pone.0048289.

Gantzer,P.A.,Bryant,L.D.,Little,J.C.,2009.Controlling soluble iron and manganese in a water-supply reservoir using hypolimnetic oxygenation.Water Res.43(5),1285-1294.https://doi.org/10.1016/j.watres.2008.12.019.

Gao,H.,Schreiber,F.,Collins,G.,Jensen,M.M.,Kostka,J.E.,Lavik,G.,de Beer,D.,Zhou,H.,Kuypers,M.M.,2010.Aerobic denitrification in permeable Wadden Sea sediments.ISME J.4(3),417-426.https://doi.org/10.1038/ismej.2010.166.

Gerling,A.B.,Browne,R.G.,Gantzer,P.A.,Mobley,M.H.,Little,J.C.,Carey,C.C.,2014.First report of the successful operation of a side stream supersaturation hypolimnetic oxygenation system in a eutrophic,shallow reservoir.Water Res.67(15),129-143.https://doi.org/10.1016/j.watres.2014.09.002.

Hallin,S.,Lindgren,P.-E.,1999.PCR detection of genes encoding nitrite reductase in denitrifying bacteria.Appl.Environ.Microbiol.65(4),1652-1657.https://doi.org/10.2478/v10006-011-0029-7.

Heylen,K.,Vanparys,B.,Wittebolle,L.,Verstraete,W.,Boon,N.,de Vos,P.,2006.Cultivation of denitrifying bacteria:Optimization of isolation conditions and diversity study.Appl.Environ.Micro.72(4),2637-2643.https://doi.org/10.1128/AEM.72.4.2637-2643.2006.

Heylen,K.,Lebbe,L.,de Vos,P.,2008.Acidovorax caeni sp.nov.,a denitrifying species with genetically diverse isolates from activated sludge.Int.J.Syst.Evol.Microbiol.58(1),73-77.https://doi.org/10.1099/ijs.0.65387-0.

Huang,T.L.,Wei,W.,Su,J.F.,Zhang,H.H.,Li,N.,2012.Denitrification performance and microbial community structure of a combined WLAOBCO system.PLoS One 7(11),e48339.https://doi.org/10.1371/journal.pone.0048339.

Huang,T.L.,2015.Water Pollution and Water Quality Control of Selected Chinese Reservoir Basins.Springer,Belin.https://doi.org/10.1007/978-3-319-20391-1.

Huang,T.L.,Zhou,S.L.,Zhang,H.H.,Bai,S.Y.,He,X.X.,Yang,X.,2015a.Nitrogen removal characteristics of a newly isolated indigenous aerobic denitrifier from oligotrophic drinking water reservoir,Zoogloea sp.N299.Int. J. Mol. Sci. 16(5), 10038-10060. https://doi.org/10.3390/ijms160510038.

Huang,T.L.,Zhou,S.L.,Zhang,H.H.,Zhou,N.,Guo,L.,Di,S.,Zhou,Z.,2015b.Nitrogen removal from micro-polluted reservoir water by indigenous aerobic denitrifiers.Int.J.Mol.Sci.16(4),8008-8026.https://doi.org/10.3390/ijms16048008.

Jiang,K.,Sanseverino,J.,Chauhan,A.,Lucas,S.,Copeland,A.,Lapidus,A.,Del Rio,T.G.,Dalin,E.,Tice,H.,Bruce,D.,2012.Complete genome sequence of Thauera aminoaromatica strain MZ1T.Stand.Genomic Sci.6(3),325-335.https://doi.org/10.4056/sigs.2696029.

Kim,H.G.,Han,G.H.,Kim,D.,Choi,J.S.,Kim,S.W.,2012.Comparative analysis of two types of methanol dehydrogenase from Methylophaga aminisul fidivorans MPT grown on methanol.J.Basic Microbiol.52(2),141-149.https://doi.org/10.1002/jobm.201000479.

Kim,J.K.,Park,K.J.,Cho,K.S.,Nam,S.-W.,Park,T.-J.,Bajpai,R.,2005.Aerobic nitrification-denitrification by heterotrophic Bacillus strains.Bioresour.Technol.96(17),1897-1906.https://doi.org/10.1016/j.biortech.2005.01.040.

Kim,M.,Jeong,S.-Y.,Yoon,S.J.,Cho,S.J.,Kim,Y.H.,Kim,M.J.,Ryu,E.Y.,Lee,S.-J.,2008.Aerobic denitrification of Pseudomonas putida AD-21 at different C/N ratios.J.Biosci.Bioeng.106(5),498-502.https://doi.org/10.1263/jbb.106.498.

Kornaros,M.,Zafiri,C.,Lyberatos,G.,1996.Kinetics of denitrification by Pseudomonas denitrificans under growth conditions limited by carbon and/or nitrate or nitrite.Water Environ.Res.68(5),934-945.https://doi.org/10.2175/106143096X127947.

Ligi,T.,Oopkaup,K.,Truu,M.,Preem,J.-K.,N~olvak,H.,Mitsch,W.J.,Mander,U¨.,Truu,J.,2014.Characterization of bacterial communities in soil and sediment of a created riverine wetland complex using highthroughput 16S rRNA amplicon sequencing.Ecol.Eng.72,56-66.https://doi.org/10.1016/j.ecoleng.2013.09.007.

Loman,N.J.,Misra,R.V.,Dallman,T.J.,Constantinidou,C.,Gharbia,S.E.,Wain,J.,Pallen,M.J.,2012.Performance comparison of benchtop high throughput sequencing platforms.Nat.Biotechnol.30(5),434-439.https://doi.org/10.1038/nbt2198.

Lukow,T.,Diekmann,H.,1997.Aerobic denitrification by a newly isolated heterotrophic bacterium strain TL1.Biotechnol.Lett.19(11),1157-1159.https://doi.org/10.1023/A:1018465232392.

Miyahara,M.,Kim,S.-W.,Fushinobu,S.,Takaki,K.,Yamada,T.,Watanabe,A.,Miyauchi,K.,Endo,G.,Wakagi,T.,Shoun,H.,2010.Potential of aerobic denitrification by Pseudomonas stutzeri TR2 to reduce nitrous oxide emissions from wastewater treatment plants.Appl.Environ.Micro.76(14),4619-4625.https://doi.org/10.1128/AEM.01983-09.

Neufeld,J.D.,Sch¨afer,H.,Cox,M.J.,Boden,R.,McDonald,I.R.,Murrell,J.C.,2007.Stable-isotope probing implicates Methylophaga spp and novel Gammaproteobacteria in marine methanol and methylamine metabolism.ISME J.1(6),480-491.https://doi.org/10.1038/ismej.2007.65.

Patureau,D.,Zumstein,E.,Delgenes,J.,Moletta,R.,2000.Aerobic denitrifiers isolated from diverse natural and managed ecosystems.Microb.Ecol.39(2),145-152.https://doi.org/10.1007/s002480000009.

Quince,C.,Lanzen,A.,Davenport,R.J.,Turnbaugh,P.J.,2011.Removing noise from pyrosequenced amplicons.BMC Bioinf.12(1),1-18.https://doi.org/10.1186/1471-2105-12-38.

Ravel,J.,Gajer,P.,Abdo,Z.,Schneider,G.M.,Koenig,S.S.,McCulle,S.L.,Karlebach,S.,Gorle,R.,Russell,J.,Tacket,C.O.,2011.Vaginal microbiome of reproductive-age women.P.Natl.A.Sci.108(s1),4680-4687.https://doi.org/10.1073/pnas.1002611107.

Robertson,L.A.,Kuenen,J.G.,1983.Thiosphaera pantotropha gen.nov.sp.nov.,a facultatively anaerobic,facultatively autotrophic sulphur bacterium.Microbiology 129(9),2847-2855.https://doi.org/10.1099/00221287-129-9-2847.

Shen,Y.Z.,Qiu,J.L.,Hong,Q.,Chen,L.,2014.Simulation of spatial and temporal distributions of non-point source pollution load in the Three Gorges Reservoir Region.Sci.Total Environ.493,138-146.https://doi.org/10.1016/j.scitotenv.2014.05.109.

ˇSmilauer,P.,Lepˇs,J.,2014.Multivariate Analysis of Ecological Data Using CANOCO 5.Cambridge University Press,Cambridge.https://doi.org/10.2136/vzj2004.1057.

Su,J.-Y.,Liu,B.-Y.,Liu,C.-Y.,2001.Comparison of aerobic denitrification under high oxygen atmosphere by Thiosphaera pantotropha ATCC 35512 and Pseudomonas stutzeri SU2 newly isolated from the activated sludge of a piggery wastewater treatment system.J.Appl.Microbiol.90(3),457-462.https://doi.org/10.1046/j.1365-2672.2001.01265.x.

Su,X.X.,Wang,H.,Zhang,Y.L.,2013.Health risk assessment of nitrate contamination in groundwater:A case study of an agricultural area in northeast China.Water Resour.Manag.27(8),3025-3034.https://doi.org/10.1007/s11269-013-0330-3.

Takaya,N.,Catalan-Sakairi,M.A.B.,Sakaguchi,Y.,Kato,I.,Zhou,Z.,Shoun,H.,2003.Aerobic denitrifying bacteria that produce low levels of nitrous oxide.Appl.Environ.Micro 69(6),3152-3157.https://doi.org/10.1128/AEM.69.6.3152-3157.2003.

Throb¨ack,I.N.,Enwall,K.,Jarvis, Å.,Hallin,S.,2004.Reassessing PCR primers targeting nirS,nirK and nosZ genes for community surveys of denitrifying bacteria with DGGE.FEMS Microbiol.Ecol.49(3),401-417.https://doi.org/10.1016/j.femsec.2004.04.011.

Willems,A.,Falsen,E.,Pot,B.,Jantzen,E.,Hoste,B.,Vandamme,P.,Gillis,M.,Kersters,K.,de Ley,J.,1990.Acidovorax,a new genus for Pseudomonas facilis,Pseudomonas delafieldii,E.Falsen(EF)group 13,EF group 16,and several clinical isolates,with the species Acidovorax facilis comb.nov.,Acidovorax delafieldii comb.nov.,and Acidovorax temperans sp.nov.Int.J.Syst.Evol.Microbiol.40(4),384-398.https://doi.org/10.1099/00207713-40-4-384.

Yan,Q.,Bi,Y.,Deng,Y.,He,Z.,Wu,L.,van Nostrand,J.D.,Shi,Z.,Li,J.,Wang,X.,Hu,Z.,2015.Impacts of the Three Gorges Dam on microbial structure and potential function.Sci.Rep.UK 5,8605.https://doi.org/10.1038/srep08605.

Zhao,W.,Hao,R.,Wang,R.,Du,P.,2015.Denitrification of composite carbon filler and character of microbial community.China Environ.Sci.35(10),3003-3009(in Chinese).

Zhou,S.L.,Huang,T.L.,Ngo,H.H.,Zhang,H.H.,Liu,F.,Zeng,M.Z.,Shi,J.C.,Qiu,X.P.,2016a.Nitrogen removal characteristics of indigenous aerobic denitrifiers and changes in the microbial community of a reservoir enclosure system via in situ oxygen enhancement using water lifting and aeration technology.Bioresour.Technol.214,63-73.https://doi.org/10.1016/j.biortech.2016.04.071.

Zhou,S.L.,Huang,T.L.,Zhang,H.H.,Zeng,M.Z.,Liu,F.,Bai,S.Y.,Shi,J.C.,Qiu,X.P.,Yang,X.,2016b.Nitrogen removal characteristics of enhanced in situ indigenous aerobic denitrification bacteria for micro-polluted reservoir source water.Bioresour.Technol.201,195-207.https://doi.org/10.1016/j.biortech.2015.11.041.

Zhu,M.,Zhu,G.,Zhao,L.,Yao,X.,Zhang,Y.,Gao,G.,Qin,B.,2013.Influence of algal bloom degradation on nutrient release at the sediment water interface in Lake Taihu,China.Environ.Sci.Pollut.R.20(3),1803-1811.https://doi.org/10.1007/s11356-012-1084-9.