Configuration of biological wastewater treatment line and influent composition as the main factors driving bacterial community structure of activated sludge

The structure of microbial consortia in wastewater treatment facilities is a resultant of environmental conditions created by the operational parameters of the purification process. In the research, activated sludge from nine Polish wastewater treatment plants (WWTPs) was investigated at a molecular level to determine the impact of the complexity of biological treatment line and the influent composition on the species structure and the diversity of bacterial consortia. The community fingerprints and technological data were subjected to the canonical correspondence and correlation analyses. The number of separated biological processes realized in the treatment line and the presence of industrial wastewater in the influent were the key factors determining the species structure of total and ammonia-oxidizing bacteria in biomass. The N2O-reducers community composition depended significantly on the design of the facility; the highest species richness of denitrifiers was noted in the WWTPs with separated denitrification tanks. The contribution of industrial streams to the inflow affected the diversity of total and denitrifying bacterial consortia and diminished the diversity of ammonia oxidizers. The obtained data are valuable for engineers since they revealed the main factors, including the design of wastewater treatment plant, influencing the microbial groups critical for the stability of purification processes. Electronic supplementary material The online version of this article (doi:10.1007/s11274-013-1273-9) contains supplementary material, which is available to authorized users.


Introduction
In wastewater treatment plants (WWTPs), the removal of nutrients should be maintained at a level ensuring the concentrations of pollutants in the effluent in accordance with the legal requirements. For this purpose, the alternating oxygen conditions must occur in the technological systems to enable the efficient nitrogen and phosphorus elimination from wastewater. In most WWTPs, in order to ensure the high efficiency of processes associated with the nutrient removal, anaerobic, anoxic and aerobic tanks are deliberately introduced into the technological line. In the WWTPs with a low flow rate, the biological treatment is performed only in aeration tanks, however, alternating oxic conditions may occur spontaneously as a result of limited oxygen diffusion into large flocs or thick biofilm.
The species structure of microbial communities in WWTPs decides the efficiency of the wastewater treatment. The knowledge about the ecology of the microbial communities responsible for the pollutant removal and the influence of the operational parameters and environmental variables on their structure and dynamics is significant for the optimization of the treatment systems (Gray et al. 2002;Molina-Muñoz et al. 2009). It was proved that the stability of wastewater treatment depends on the diversity of microbial groups responsible for the particular processes. The high species richness of microbial community ensures conservation of a given functionality and a quick recovery of the consortia after the stress conditions because of the functional redundancy and the alternative ways to use the flow of energy (Fernandez et al. 1999).
Electronic supplementary material The online version of this article (doi:10.1007/s11274-013-1273-9) contains supplementary material, which is available to authorized users.
The threat posed by nitrogen and phosphorus to aquatic ecosystems requires the reduction of their load discharged into the environment. In most WWTPs, phosphorus is removed by chemical precipitation, so it is relatively easy to deal with its excessive load. Nitrogen is only removed as a result of the biological conversions, therefore it is necessary to ensure the high diversity and number of nitrifiers and denitrifiers in biomass. Nitrification is a two-phase process in which ammonium is oxidized to nitrite by ammonia-oxidizing bacteria (AOB) and then nitrite is oxidized to nitrate by nitrite-oxidizing bacteria (NOB). The path of nitrogen removal from wastewater is mostly via heterotrophic denitrification involving the reduction of nitrate and nitrite to gaseous products (NO, N 2 O, N 2 ) under the conditions of the limited oxygen availability. Denitrifying bacteria belong to a broad variety of groups and encompass a wide range of physiological traits (Zumft 1997).
Since most of the microorganisms from the environmental samples cannot be cultivated, the monitoring of their diversity and number must involve molecular tools. In particular, polymerase chain reaction-denaturing gradient gel electrophoresis (PCR-DGGE) is reckoned as a rapid and reliable method for the relative comparison of the different bacterial communities (Muyzer et al. 1993). This fingerprinting method allows the direct measurement of the genetic biodiversity of a sample and coupled with sequencing and phylogenetic analyses can provide an overview of the microbial composition and diversity in the system (Sanz and Köchling 2007).
Many studies have focused on the analyzing the relations between the bacterial communities, especially nitrifiers, in biomass, and a type of wastewater (Limpiyakorn et al. 2011;Whang et al. 2009) in different bioreactors (Wan et al. 2011;Ye et al. 2011). The literature scarce the data on the influence of the full-scale wastewater treatment plant organization on the bacterial consortia in biomass. Therefore, the goal of the study was a robust statistical analysis of the dependence between the presence of the particular tanks and processes in the biological treatment line of WWTPs, the influent characteristics and the community structure of total, nitrifying and denitrifying bacteria in the activated sludge. The relations between the microbial assemblages and the technological data were investigated using the canonical correspondence analysis and correlation matrix.

WWTPs
The WWTPs were chosen to obtain a broad spectrum of the possible technological solutions. We analyzed nine facilities operated in the activated sludge method that differed in the influent composition and the characteristics of the biological treatment line. The People Equivalent (PE) of analyzed WWTPs varied between 2,400 and 177,000. Five of the WWTPs, Tyrowo (TY), Jędrychówko (JE), Elbląg (EL), Nowe Gizewo (NG) and Olsztyn (OL) received industry and domestic wastewater. The remaining WWTPs, namely Gietrzwałd (GT), Jonkowo (JO), Rakowiec (RA) and Łyse (LY) treated domestic wastewater. The detailed characteristics of the facilities is given in Table 1. The technological results were obtained from WWTP operators, the influents and effluents from RA and GT were analyzed in our laboratory, according to APHA (1992).

PCR-DGGE
The biomass from each WWTP was sampled twice from aerobic tanks in a period from January to March 2012 and was frozen in -20°C prior to molecular analysis. DNA was extracted from approximately 400 mg of centrifuged sample using a FastDNA Ò SPIN Ò Kit (Q-BIOgene, Canada). The working solutions with the DNA concentration of 50 ng/lL were prepared and the concentration of the DNA was measured spectrophotometrically using BioPhotometer (Eppendorf, Germany). The PCRs were performed in an Eppendorf Ò Mastercycler Gradient (Eppendorf). The determination of the total bacteria diversity was based on 16S rDNA analysis with a primer set 341F/515R. The first-phase nitrifiers' diversity was based on the analyses of amoA gene (301F/302R and amoA-1F/amoA-2R primer sets) that codes for the ammonium monooxygenase involved in the ammonium oxidation to nitrite. The denitrifying bacteria diversity was assessed based on the presence of nosZ gene (NosZ1/NosZ2 primer set), which codes for the nitrous oxide reductase responsible for the last step of denitrification (the gene is present in bacteria conducting full denitrification). The primer sequences are given in Table 2.
The PCR mixture contained 1.7 ng/lL of extracted DNA, 0.5 lM of each primer, 100 lM of deoxynucleoside triphosphate mixture (Promega, Madison, USA), 1.5 U of GoTaq Ò DNA Polymerase (Promega), 6 lL of 10 9 reaction buffer supplied with polymerase, 1.5 mM MgCl 2 and sterile water to a final volume of 30 lL. The amplification of the amoA gene was performed as a nested-PCR (Cydzik-Kwiatkowska and Wojnowska-Baryła 2011). The thermal profile for the 16S rDNA amplification was: 94°C for 5 min, 35 cycles of: denaturation at 94°C for 45 s, annealing at 62°C for 45 s, extension at 72°C for 1 min, and a final elongation at 72°C for 5 min. The thermal profile for the nosZ gene amplification was: 94°C for 5 min, 6 cycles of touchdown PCR (denaturation 94°C for 30 s, annealing for 1 min with an 1°C for two cycles decrement at temperature 61°C, extension at 72°C for 1 min), followed by 25 cycles of 94°C for 30 s, 58°C for 1 min, 72°C for 1 min and a final elongation for 10 min at 72°C. The presence of the PCR products was confirmed by agarose electrophoresis. The amplified products were resolved on DGGE gels using a dCode System (Bio-Rad, USA), the electrophoresis conditions are presented in Table 2.
Denaturing gradient gel electrophoresis gels were stained with SYBR Gold (Invitrogen, USA) and digitalized using Kodak 1D 3.6 Image Analysis Software (Eastman Kodak Company, USA). The amoA amplicons that were clear and had a high intensity were excised from the DGGE gel, reamplified and sequenced in the Institute of Biochemistry and Biophysic, Polish Academy of Science (http://www.oligo.ibb.waw.pl). The nucleotide sequences were compared with the sequences in the GenBank using the BLASTn program (Altschul et al. 1997) and deposited in the GenBank under the accession No. JX546276-79. The sequences determined in this study were aligned and the genetic relationships were determined (the Maximum Likelihood method) using the MEGA5 software (Tamura et al. 2011).

Calculation methods
In all tests, the significant effects were those with p value \0.05. For the calculations, the numerical values were assigned to express the number of separate biological processes realized in the biological treatment line of WWTPs and the presence of the tanks favoring denitrification (Table 3).
The canonical correspondence analysis (CCA) was performed on the relative DGGE band intensities with the Monte Carlo permutation testing (499 permutations). To the CCA analysis, next to metadata obtained by DGGE, the technological data such as the COD/N and BOD/COD ratios, TKN and COD in the influent, the presence (IN ?) or absence (IN-) of industrial wastewater in the influent, the presence of the denitrification tanks (DT) and the number of processes designed to occur in the biological treatment line of WWTP (PR) were taken. The flow rate was removed from the analysis since it was strongly correlated (R = -0.98) with the COD concentration in the influent (ter Braak 1986). The values of DT and PR for each WWTP are presented in Table 4.
The analyses were carried out using the CANOCO for Windows ver. 4.51 and CANODRAW (ter Braak, Smilauer, Biometris, Wageningen, The Netherlands). The correlations between the bacterial diversity and the technological parameters were analyzed using the correlation matrix in the Statistica 10.0 (StatSoft, USA). The statistical analysis assumed that the values from 0.9 to 1.0 point out to almost full correlation, very high correlation Table 1 The treatment lines and characteristics of influents and effluents of investigated WWTPs  reaches the values from 0.7 to 0.9, strong-from 0.5 to 0.7, the medium correlation from 0.3 to 0.5, while weak from 0.1 to 0.3 (Stanisz 2000).

Results
In the present research, to analyze the bacterial assemblages, the DNA isolated from the biomass was amplified using the specific primer sets and the obtained products were separated in DGGE. In general, a higher number of bands indicates greater diversity of the analyzed microbial consortia. The DGGE separation of the amoA PCR products is presented in Fig. 1 showed that most of the sequenced bands (A, B, C, D) were related to the Nitrosospira sp. The band E was closely related to Nitrosomonas eutropha. This band was in all the plants except EL, we also detected only a weak band E in GT.
The canonical correspondence analysis is designed for relating the species composition of communities to their environment and can provide an insight into the impact of wastewater treatment plant design and the operational parameters on the bacterial assemblages. The data analysis can be used in an explanatory way and it leads to an ordination diagram of samples, species and environmental variables, which optimally displays how the community composition varies with the environment. When used in a confirmatory way, it leads to statistical tests of the effects of particular environmental variables on the community    Fig. 3a, the points represented the total bacterial consortia in different WWTPs while the environmental variables were represented by arrows. The bacterial community structure differed between the WWTPs and no separate clusters were observed in the diagram. The first axis of a biplot explained the 24 % of the species-environment relation, while the summary variation explained by the two-dimensional diagram was 42 %. The first and the second eigenvalues equaled 0.31 and 0.23, respectively. The species-environment correlations of the first two axes were very high (above 99 %) showing that the measured environmental variables were sufficient to explain the major variations among the analyzed WWTPs. The correlation coefficients showed that the first axis was a presence or an absence of industrial wastewater in the WWTP influent (R = ± 0.71). The discrete variable IN had the highest power to explain the patterns in the species data and the significance of the explanatory effect was statistically important (Monte Carlo permutation test, F = 1.48, p = 0.05). The correlations of the second axis showed a contrast between the WWTPs with the different number of processes realized in the treatment line (R = 0.50) and Fig. 1 The DGGE analysis of PCR amplifications of the partial amoA gene. The abbreviations above each lane represent the WWTP the activated sludge samples were taken from. The bands that were sequenced are marked with capital letters Fig. 2 Phylogenetic tree showing the relationships of the partial amoA gene sequences to reference the sequences from GeneBank database (accession numbers given in parentheses). The sequence of band A is given in Supplementary material. The evolutionary history was inferred by using the Maximum Likelihood method based on the Tamura-Nei model (Tamura and Nei 1993). Initial tree(s) for the heuristic search were obtained automatically as follows. When the number of common sites was \100 or less than one fourth of the total number of sites, the maximum parsimony method was used; otherwise the BIONJ method with MCL distance matrix was applied. The tree is drawn to scale, with branch lengths measured in the number of substitutions per site. Evolutionary analyses were conducted in the MEGA5 (Tamura et al. 2011) World J Microbiol Biotechnol (2013) 29:1145-11531149 different TKN and COD concentrations in the influent (R = 0.48). The length of an arrow representing an environmental variable is a measure of how much the total bacteria community structure differs along that environmental variable. Since the PR arrow was the longest one, the number of processes realized in the WWTPs was the most important continuous environmental variable influencing the total bacteria community structure in activated sludge from all analyzed. From Fig. 3b it can be concluded that the points representing the N 2 O-reducing bacteria communities in JE, JO, OL and EL were similar and grouped together in the diagram. The first and the second axes explained 36 and 20 % of the species-environment relation with the respective eigenvalues of 0.41 and 0.24. The species-environment correlations for the first axis was 0.98 while for the second one it equaled 0.99. The correlation coefficients showed that the first axis is DT in the WWTP treatment line (R = -0.73). This environmental continuous variable had the highest power to explain the DGGE patterns of denitrifiers, and the significance of the explanatory effect was statistically important (Monte Carlo permutation test, F = 1.97, p = 0.02). The correlations of the second axis showed that the BOD/COD ratio of the influent was the second environmental variable influencing mostly the N 2 O-reducing bacteria communities in activated sludge from the analyzed WWTPs (R = -0.58).
The CCA analysis of the DGGE patterns characterizing the AOB communities showed that the most similar assemblages were in GT, EL, TY, JE and OL (Fig. 3c). The first axis of a biplot explained the 30 % of the species-environment relation, while the summary variation explained by the two-dimensional diagram was high and equaled 53 %. The first and the second eigenvalues equaled 0.57 and 0.45, respectively. The species-environment correlations of the first two axes were above 99 %. The correlation coefficients showed that the first axis is concentration of TKN in the influent (R = 0.67) and that this factor was the most important continuous environmental variable. The second axis was the discrete variable namely the presence/absence of industrial wastewater in the influent (IN; R = ± 0.65). The significance of the explanatory effect of IN was statistically important (Monte Carlo permutation test, F = 1.82, p \ 0.00). Fig. 4 presents the microbial diversity expressed as the average number of amplicons in the DGGE patterns. The number of 16S rDNA bands varied depending on the WWTP. The highest number of amplicons (47) was in TY, the WWTP with a technological line ensuring the removal of carbon, nitrogen and phosphorus with the respective efficiency of 96, 86, 95 %. The lowest number of 16S rDNA amplicons was obtained in GT. Three facilities (JE, RA, GT) operated without the separated denitrification tanks had the lowest diversity of denitrifiers (13, 16 and 15 nosZ bands, respectively) in activated sludge. The highest diversity of denitrifiers (from 26 to 32 nosZ bands) was obtained in activated sludge from EL, TY and NG. The highest diversity of AOB (ca. 27 bands) was noted for the biomass from JO and LY. The lowest ammonia-oxidizers diversity characterized activated sludge from GT.
The correlation matrix was constructed to find the dependences between the diversity of particular groups of bacteria (total, N 2 O-reducing and ammonia-oxidizing) and the investigated WWTP design and the influent composition. It was observed that the total bacteria diversity was strongly correlated with the COD/N ratio (R = 0.84) and the number of processes realized in the facility (R = 0.79). We also obtained a positive strong correlation (R = 0.50) between the 16S rDNA band number and the presence of denitrification tank. As to the facility design, the presence of anoxic tanks had the greatest impact on the diversity of N 2 O-utilizing denitrifiers (R = 0.57). Analyzing the wastewater characteristics, both the presence of the industrial stream and the high COD/N ratio in the WWTP influent influenced the species richness of denitrifiers (R = 0.51 and R = 0.48, respectively). The AOB diversity was weakly positively correlated with the TKN concentration in the influent (R = 0.25) and negatively influenced by the presence of industrial wastewater in the influent (R = -0.25).

Discussion
In the research, activated sludge from nine WWTPs was investigated to determine the influence of the complexity of the biological treatment line and the influent characteristics on the structure and diversity of microbial consortia. For the statistical analysis, the CCA and correlation matrix were applied.
The species composition of the total bacterial communities in activated sludge was influenced by both the design of the treatment line and the wastewater characteristics. The CCA analysis showed that the number of processes realized in the WWTPs was the most important of all continuous environmental variables deciding about the structure of total bacteria assemblages. The complexity of the treatment line positively affected also the overall diversity of bacteria in the biomass. The highest number of different bacterial species was noted for TY with a complex, many-stage technological line. Reversely, the lowest number of 16S rDNA bands characterized activated sludge community in GT, the smallest of all analyzed WWTPs (PE = 2,405) with the simplest biological treatment line consisting only of the aeration tank. It can be concluded that the presence of many tanks with the different oxic conditions in the treatment line favors the growth of multispecies microbial consortia thus promoting the stability of the purification processes. As to the wastewater composition, the statistically significant influence of the industrial wastewater presence on both the structure and the diversity of the total bacteria communities was proven. Since the origin of industrial wastewater in the analyzed treatment facilities was broad, it was difficult to claim which components of the sewage promoted the biodiversity. Our results, however, pointed out that the presence of the industrial stream favored species richness of total bacteria in activated sludge due to the accessibility of a broader range of substrates as compared to that present in typical domestic wastewater.
The WWTPs are one of the major sources of nitrogen oxides, especially N 2 O. Many parameters such as a substrate concentration, a C/N ratio, a type of carbon source, a nitrite accumulation and an NO concentration influence the N 2 O production (Adouani et al. 2010). The knowledge about the links between the wastewater treatment line design and the structure of N 2 O-utilizing bacteria communities in the biomass would allow engineers to apply the purification strategy favoring the effective reduction of the greenhouse gas emission. The research proved that the number of separated denitrification tanks in the WWTPs was the major factor deciding about the species composition and the diversity of the N 2 Oreducing microorganisms in activated sludge. In three facilities (JE, RA, GT) operated without the separated denitrification tanks, the lowest species richness of denitrifiers in the biomass was noted. In JE, the tank with the simultaneous nitrification/ denitrification was exploited meaning that nitrifiers and denitrifiers co-existed in activated sludge, however, the alternating oxic conditions did not favor the full denitrification to N 2 since the activity of bacteria that possess the nosZ genes is inhibited in the presence of even low oxygen concentrations (Hochstein et al. 1984). For comparison, in OL, nitrogen removal was realized by SND in a carrousel-type bioreactor, however there was an additional tank for pre-denitrification favoring the growth of more diverse assemblage of the N 2 O-utilizing denitrifiers. The highest diversity of denitrifiers was obtained in activated sludge from EL, TY and NG. These WWTPs were characterized by both the presence of separate anoxic tanks and the highest ratio of COD/N (12-14) that was also proven to positively influence the species richness of investigated bacteria (R = 0.48).
The microorganisms capable of denitrification belong to a broad variety of groups and encompass a wide range of the physiological traits. Most denitrifiers are the aerobic heterotrophic organisms that transfer redox equivalents  Fig. 4 The average number of 16S rDNA, nosZ, amoA bands in the DGGE patterns characterizing activated sludge communities from analyzed WWTPs; the averages of two different measurements, standard deviations are given World J Microbiol Biotechnol (2013) 29:1145-11531151 from the oxidation of a carbon source to N oxides under the anaerobic conditions. The ability to reduce N 2 O by the nitrifying and phosphorus accumulating bacteria is also documented (Zumft 1997). This can explain the tendency observed in the current research that more processes realized separately in the treatment line favored the diversity of N 2 O-reducers. Siripong and Rittman (2007) investigated the AOB communities in activated sludge in seven typical singlestage municipal plants. Among the analyzed factors (flow rate, influent and effluent BOD and TKN, effluent ammonia, nitrite and nitrate, pH, and sewage temperature), only the seasonal temperature variations seemed to change the nitrifying community, especially the balance between Nitrosospira sp. and Nitrosomonas sp., although both genera coexisted in winter and summer samples. Lydmark et al. (2007) showed that ammonium concentration was an important structuring factor for an AOB community. In the present research, the CCA proved that the main variable influencing the AOB consortia in activated sludge was the presence of the industrial stream in the WWTPs influent. In fact, the IN variable was the only one that had a statistically significant impact on the species composition of the AOB. These results are well depicted in Fig. 3c. The DGGE patterns characterizing activated sludge from GT, EL, TY, JE and OL facilities with the lower values of COD and TKN in the influent (except from TY) and the presence of industrial influent (except from GT) were grouped together in a separate cluster.
The coexistence of various nitrifiers in WWTP is an evidence of a functional redundancy, a feature that may help in the maintaining the stability of the system for nitrification. Wang et al. (2010) investigated the communities of microorganisms in activated sludge of eight wastewater treatment systems and suggested the negative impact of the presence of industrial wastewater in the influent on the nitrifiers' diversity. Our observations pointed out to a weak negative correlation between the industrial wastewater presence and the species richness of AOB (R = -0.25), nevertheless, the highest diversity of these bacteria was noted in biomass from JO and LY that received only domestic wastewater. The lowest AOB diversity characterized activated sludge from GT, the WWTP with the simplest technological line consisting of only an aeration tank. It can be concluded that in this system, under a stable oxygen concentration, the microorganisms underwent a strong selective pressure favoring the growth of only a few best adapted species. Whang et al. (2009) evaluated the nitrifying community and the nitrification performance of the full-scale municipal (20 mg N/L) and swine (220 mg N/L) WWTPs. Authors observed dissimilar nitrifying populations prevailing in these two plants and related this fact to different input nitrogen concentrations.
In our research, in the analyzed range of the influent TKN (62 ± 14-129 ± 18 mg/L), only a weak correlation between this parameter and the AOB diversity was proven.
In general, the members of the Nitrosospira spp. or/and the Nitrosomonas oligotropha clusters are the dominant AOB in the ammonia-low environments, whereas the members of the N. europaea-Nitrosococcus mobilis cluster comprise the majority of AOB in the ammonia-rich environments (Limpiyakorn et al. 2005). The sequencing of the chosen amoA bands and the phylogenetic analysis showed that the AOB in the analyzed activated sludge samples were related to both Nitrosomonas sp. and Nitrosospira sp. (Fig. 2). An interesting observation can be made on band E, obtained from N. eutropha. This species was present in all the WWTPs except from EL, in GT only a weak band was observed. N. eutropha is commonly found in strongly eutrophic environments such as the municipal and industrial sewage disposal systems. This microorganism has a high tolerance for both elevated ammonia concentrations and the fluctuating conditions (especially oxic/anoxic cycles) (Koops et al. 1991;Stein et al. 2007). The low number of N. eutropha in GT (DGGE is a semi-quantitative technique so such a statement is justified) can be explained by the simplicity of the biological treatment line in this WWTP-it consisted of only an aeration tank. The stable oxic conditions and a lack of environmental fluctuations did not favor the growth of the analyzed species. The EL WWTP, on the other hand, received wastewater from a wood industry that did not use water as an input for the manufacturing processes but generated several wastewater streams after washing/cleaning procedures. In the production of wood-based floors and wood-laminates, small volumes of highly polluted wastewaters with the high contents of formaldehyde and COD are generated (e.g. cleaning of machines that are used to apply urea-formaldehyde resins onto wood-fiber boards) (Kaczala et al. 2010). Wastewater with formaldehyde could have negatively influenced the AOB community in EL, especially N. eutropha, since its genome lacks genes for urease metabolism (Stein et al. 2007).
Open Access This article is distributed under the terms of the Creative Commons Attribution License which permits any use, distribution, and reproduction in any medium, provided the original author(s) and the source are credited.