ORIGINAL RESEARCH article
Shifts in the pelagic ammonia-oxidizing microbial communities along the eutrophic estuary of Yong River in Ningbo City, China
- 1Faculty of Architectural Civil Engineering and Environment, Ningbo University, Ningbo, China
- 2Department of Microbial Wetland Ecology, Netherlands Institute of Ecology (NIOO-KNAW), Wageningen, Netherlands
- 3Institute of Environmental Biology, Utrecht University, Utrecht, Netherlands
Aerobic ammonia oxidation plays a key role in the nitrogen cycle, and the diversity of the responsible microorganisms is regulated by environmental factors. Abundance and composition of ammonia-oxidizing archaea (AOA) and ammonia-oxidizing bacteria (AOB) were investigated in the surface waters along an environmental gradient of the Yong River in Ningbo, East China. Water samples were collected from three pelagic zones: (1) freshwaters in the urban canals of Ningbo, (2) brackish waters in the downstream Yong River, and (3) coastal marine water of Hangzhou Bay. Shifts in activity and diversity of the ammonia-oxidizing microorganisms occurred simultaneously with changes in environmental factors, among which salinity and the availabilities of ammonium and oxygen. The AOA abundance was always higher than that of AOB and was related to the ammonia oxidation activity. The ratios of AOA/AOB in the brackish and marine waters were significantly higher than those found in freshwaters. Both AOA and AOB showed similar community compositions in brackish and marine waters, but only 31 and 35% similarity, respectively, between these waters and the urban inland freshwaters. Most of AOA-amoA sequences from freshwater were affiliated with sequences obtained from terrestrial environments and those collected from brackish and coastal areas were ubiquitous in marine, coastal, and terrestrial ecosystems. All AOB from freshwaters belonged to Nitrosomonas, and the AOB from brackish and marine waters mainly belonged to Nitrosospira.
Excessive inorganic nitrogen is often the cause of eutrophication, which is a key important factor contributing to habitat variation and to catastrophic algal bloom expansion in aquatic ecosystems (Glibert et al., 2005). At the same time, complex relations between ammonia oxidizers and canonical denitrifiers are supposed to contribute to about 70% of the global N2O emissions, which has detrimental effects on global warming and on the Earth’s ozone layer (IPCC, 2007; Hatzenpichler, 2012). Being the first step of a coupled nitrification/denitrification process, which enables the removal of a large part of anthropogenically produced reactive nitrogen, aerobic ammonia oxidation is of pivotal significance in aquatic ecosystems (Seitzinger, 1988; Francis et al., 2007). Ammonia oxidation is controlled by the enzyme ammonia monooxygenase, which is partly encoded by the amoA gene (Rotthauwe et al., 1997; Kowalchuk and Stephen, 2001). Hence, amoA gene is a useful molecular marker for ammonia-oxidizing microorganisms in the environment (Francis et al., 2005, 2007).
Besides ammonia-oxidizing bacteria (AOB), ammonia-oxidizing archaea (AOA) are also potential contributors to nitrification (Francis et al., 2007; Prosser and Nicol, 2012). Previous studies have shown that the AOA and AOB communities were sensitive to many environmental factors (Erguder et al., 2009; Fernàndez-Guerra and Casamayor, 2012; Hatzenpichler, 2012; Cao et al., 2013). Estuaries are considered as major players in the transformation and removal of reactive nitrogen before it reaches the sea (Seitzinger et al., 2006). The steep physicochemical gradient in an estuary provides an ideal natural laboratory for studying ammonia-oxidizing archaea and bacteria (Bernhard and Bollmann, 2010; Santoro, 2010). Many worldwide studies in estuaries have reported that spatial variation in environmental factors drive shifts in abundance and community composition of ammonia-oxidizing microorganisms in sediments, i.e., in the Sacramento San Joaquin Delta (Damashek et al., 2015), the Scheldt estuary (Sahan and Muyzer, 2008), the Chesapeake Bay (Francis et al., 2003), the Plum Island Sound (Bernhard et al., 2005, 2007, 2010), the Ythan Estuary (Freitag et al., 2006), the Tokyo Bay (Urakawa et al., 2006), the Bahía del Tóbari estuary (Beman and Francis, 2006), the San Francisco Bay (Mosier and Francis, 2008), the Huntington Beach (Santoro et al., 2008), the Barn Island salt marsh (Moin et al., 2009), the Jiaozhou Bay (Dang et al., 2010), the Pearl River estuary (Cao et al., 2011; Jin et al., 2011; Xie et al., 2014), the Elkhorn Slough estuary (Wankel et al., 2011), and the Qiantang River (Liu et al., 2013). The studies of de Bie et al. (2001) and of Bollmann and Laanbroek (2002) are the only studies performed in the water column of an estuary, i.e., the Scheldt estuary. The series of studies have indicated that environmental factors such as salinity, trophic status, oxygen and temperature are likely to be potential selective factors shaping ammonia-oxidizing microbial communities.
The Yong River that passes through Ningbo City in eastern China is a typical example of a eutrophic coastal river system. Due to a limited runoff in relation to the tidal flows, pollutants emitted by anthropogenic activities have caused serious eutrophication problem in the aquatic ecosystem of Ningbo City, because most canals of the water system are constrained by dams built to defend against high tides, before they converge into the mainstream Yong River (Zhao and Yang, 2009). Therefore, to relieve problems with reactive nitrogen in the urban ecosystem, more should be known about the composition of the ammonia-oxidizing microbial communities and the dynamics of the ammonia oxidation process in this water system.
We hypothesized that the activity, abundance and species composition of AOA and AOB in the surface waters of the Yong River will be affected by the physicochemical water properties at the different zones of the estuary. The characteristics of the ammonia-oxidizing microorganisms were investigated in fresh, brackish and marine waters. The composition and abundance of amoA genes for both AOB and AOA in surface waters were compared between the fresh, brackish and marine areas. The relationship between ammonia oxidizers diversity and ammonia oxidation rate was also determined.
Materials and Methods
Sampling Description and Water Properties
As a coastal river system, the Yong River consists mainly of the Yao River, the Fenghua River, and the Yong River. The mainstream Yong River is about 26 km in length, from Sanjiangkou to Zhenhai, which is the entrance to the East China Sea (Figure 1). The runoff of the Yong River is only 5–23% of the tidal flows (Zhao and Yang, 2009). Sampling zones were divided into three geographic regions along the Yong River spanning freshwater, brackish water, and marine coastal water sites: (1) two freshwater sites in urban canals at Jiangdong (W1; N29°51′50.7″, E121°35′43.6″) and Haishu (W2; N29°53′55.4″, E121°31′45.9″), respectively; (2) two brackish sites near Ningbo University (W3; N29°54′32.1″, E121°37′35.4″) and in the Beilun estuary (W4; N29°58′14.2″, E121°46′36.8″), respectively; and (3) two coastal water sites near the island of Zhoushan (W5; N29°59′04.2″, E121°50′52.3″) and (W6; N30°02′08.1″, E121°50′ 29.4″) both in the Hangzhou Bay remote from residential and large-scale industrial areas and hence waste and freshwater were absent. Triplicate water samples were collected within a 5 m radius at each sampling site from a depth of 30 cm under the water surface. Water sampling was carried out at low tide on November 23, 2012.
FIGURE 1. Locations of the sampling sites in the Ningbo City urban canals (W1 and W2), in the Yong River estuary (W3 and W4) and in the Hangzhou Bay (W5 and W6).
The water parameters of salinity, pH, and dissolved oxygen (DO) were determined in situ by a portable meter (Hach Company, Loveland, CO, USA), and the in situ temperature was about 12°C. The triplicate water samples were separated in two portions after removal of larger particles and transferred into sterilized bottles; one portion was placed on ice in a cool box for extracting DNA and the other portion was kept at 4°C for determining water properties and potential ammonia oxidation activity (PAA). Samples were transported to the laboratory within 2 h. The portion of 500 mL waters kept at 0°C were filtered through 0.22 μm pore-sized filters (47 mm diameter, Durapore GVWP, Millipore) to collect all microorganisms. The filters were stored and frozen at -80°C until further analysis. The ammonium (NH4+) and nitrite (NO2-) concentrations were determined by a flow injection analyzer (Lachat-QC8500, Lachat Instruments, Loveland, CO, USA), total nitrogen (TN) and nitrate (NO3-) were also determined by the flow injection analyzer after alkaline persulfate digestion and cadmium-copper reduction, respectively. Total organic matter (TOC) was determined by a total organic carbon analyzer (TOC-VCPH, Shimadzu).
Potential Ammonia Oxidation Rates
Potential ammonia oxidation activity of water samples was measured according to a previously described method with minor modification (Kasuga et al., 2010). Samples were taken at 0, 6, 12, 24, and 48 h for determining the ammonium concentrations. Generally, ammonium consumption accelerated after 6 h, but slowed down again after 12 h of incubation. The decrease in ammonium concentrations between 0 and 6 h was used to determine the potential activity as this reflected mostly the potential ammonium consumption rate in the water column of the estuary. The rate determinations were carried out by incubating 400 mL in the dark at 25°C at 100 rpm. Inorganic medium per liter of water was composed of (NH4)2SO4 0.18 mM, MgSO4⋅7H2O 0.81 mM, CaCl2 0.18 mM, K2HPO4 0.76 mM, NaH2PO4 1.30 mM, Fe-EDTA 0.02 mM, CaCO3 2.00 mM, Na2MoO4⋅2H2O 4 nM, MnCl2⋅4H2O 10 nM, CoCl2⋅6H2O 8 nM, CuSO4⋅5H2O 80 nM and ZnSO4⋅7H2O 300 nM. The ammonium concentrations were determined according to the method described above.
Extraction of the nucleic acids from the filters containing the microorganisms was performed using the Fast DNA® spin kit for soil (MP Biomedicals, LLC, Illkirch, France). The filters were cut into small pieces and transferred into the Lysing Matrix E tubes containing a matrix designed for the lysis of most of the cell types, using ethanol flamed forceps and surgical scissors. The mixture was homogenized in a FastPrep®-24 (Bio 101, CA) for 15 s at 4 m⋅s-1 (Sekiguchi et al., 2002). The extraction procedure was performed as described in the manufacturer’s instructions. DNA extracts were suspended in 100 μL of DES solution (Qbiogene, Vista, CA, USA) and quantified by a Nanodrop spectrophotometer. The DNA extracts were stored frozen at -80°C until further analysis.
CrenamoA23f/CrenamoA616r (Tourna et al., 2008) and amoA1F-GC/amoA2R (Rotthauwe et al., 1997) were primers for the archaeal and bacterial amoA genes, respectively (Table 1). Amplifications were performed for AOA and AOB in 50 μL reaction mixtures including 1× PCR buffer, 3.0 mM MgCl2, 400 μM of each dNTP, 2.5 U Taq DNA polymerase (Takara, Dalian, China) and 0.2 mg mL-1 bovine serum albumin (BSA) plus 0.2 mM of each primer, 2 and 4 μL genomic DNA (80–120 ng μL-1) were added as template for AOA and AOB, respectively.
DGGE analyses were carried out using a Dcode Universal Mutation Detection System (Bio-Rad Laboratories, Hercules, CA, USA) under the conditions presented in Table 1. 8 μL of the PCR products of either AOA-amoA or AOB-amoA genes were loaded into polyacrylamide gel wells. The gels were stained with 1:10,000 SYBRGreen I (Sigma, Taufkirchen, Germany) for 30 min and scanned by a GelDoc XR (BIO-RAD, Hercules, CA, USA).
Cloning and Sequencing
Selected DGGE bands were excised for cloning and sequencing according to the method described in a previous study (Zhang et al., 2011). Excised bands were crushed and immersed overnight at 4°C in 30 μL of sterilized water to release DNA, and subsequently amplified with the primers CrenamoA 23f/CrenamoA616r (Tourna et al., 2008) and amoA1F/amoA2R (Rotthauwe et al., 1997), respectively. PCR products were purified and linked into the PMDTM 19-T Vector (TaKaRa Code: D102A, Dalian). The resulting ligation mixture was transformed into Escherichia coli DH5α competent cells following the instructions of the manufacturer. Positive clones were amplified using the above-mentioned primers with a GC clamp and cross-checked by DGGE for migration behavior and sent to Invitrogen (Shanghai, China) for sequencing.
Real-time PCR (qPCR)
qPCR assays of amoA gene of all samples were performed in triplicates using an Eppendorf Mastercycler® ep realplex2 machine (Eppendorf, Hamburg, Germany). Each reaction was performed in a 25 μL volume containing 1 μL of DNA extract (80–120 ng μL-1) as template, 0.2 mg mL-1 BSA, 0.4 mM of each primer and 12.5 μL Premix Ex taqTM (Takara, Dalian). A total of 40 cycles was run with the procedure and primers listed in Table 1. Product specificity was confirmed by melting curve analysis and visualization on 1.2% agarose gels. Specific products were observed at the expected sizes of ca. 491 and 629 bp for the bacterial and archaeal amoA genes, respectively.
Standard curves for qPCR were developed as previously described (Zhang et al., 2011). Archaeal and bacterial amoA genes fragments were cloned as described in the section of cloning and sequencing. The qPCR standard DNA was amplified from representative plasmids containing AOA-amoA or AOB-amoA gene linearized by EcoR I, and then purified by a Purification Kit (Tiangen Biotech Co., China). The concentration of the purified DNA was measured with a spectrophotometer (Nanodrop). Ten-fold serial dilutions of a known copy number of the amoA gene were generated to produce the standard curve over seven orders of magnitude (from 2.23 × 102 to 2.23 × 108 copies for AOA-amoA, and from 6.51 × 102 to 6.51 × 108 copies for AOB-amoA) per assay. High amplification efficiencies of 95.0–93.6% were obtained for archaeal and bacterial amoA quantification with R2 values 0.990 and 1.000, slopes were -3.35 and -3.38, respectively. A negative control without template DNA (NTC) was present in every qPCR assay. No evident inhibition was found with any of the samples. Inhibition of DNA extraction was tested using a previously described method (Dumonceaux et al., 2006). No significant inhibition was observed at any dilution (1–120 ng μL-1 DNA extract).
Based on DGGE banding patterns, similarity of water samples was elucidated with an unweighted pair group method with mathematical averages (UPGMA) using the software package Quantity One (Bio-Rad Laboratories, Hercules, CA, USA). The diversity indices of Shannon (H), Evenness (J) and Simpson (D) were calculated for DGGE band pattern of amoA genes with the PRIMER5.0 software. The one-way ANOVA and Pearson’s coefficient analysis were carried out with SPSS software (version 17.0; SPSS Statistics Inc., Chicago, IL, USA). To explore which environmental factors contributed to the changes among sites, principal-component analysis (PCA) was performed with PAST ver. 2.17c (Hammer et al., 2001), based on the water properties, including water salinity, pH, DO, TN, ammonium (NH4+-N), nitrite (NO2-–N), nitrate (NO3-–N) and total organic carbon (TOC) for all six stations.
Sequences were compared with GenBank database using BLAST (http://www.ncbi.nlm.nih.gov/BLAST/), and the closest matches were included in the alignment and phylogenetic analysis. The neighbor-joining trees were constructed using MEGA 4 and bootstrapped 1,000 times to calculate linear distances, and alignments were edited using the Bioedit sequence alignment Editor 184.108.40.206 (Tamura et al., 2007).
Sequence Accession Numbers
All sequences from DGGE bands have been deposited in the EMBL nucleotide sequence database with accession numbers from LN845817 to LN845837 for the AOA-amoA genes, and from LN845795 to LN845813 for the AOB-amoA genes.
Water Properties and Potential Nitrification Rates
Water properties varied among the different sampling sites (Table 2). Water salinity ranged from 0.20 to 24.60%. The highest salinities were found at Zhoushan 1 and 2 (W5 and W6), followed by Beilun estuary (W4: 7.60%) and Yong River next to Ningbo University (W3: 4.07%), and the lowest value of 0.20% was observed at Jiangdong (W1) and Haishu (W2). The salinities in the coastal waters at W5 and W6 were significantly more saline than those at the other sites, and the salinities in the Yong River downstream of W3 and W4 were also significantly higher than those of the urban canals at W1 and W2. Water pH at the urban canal sites W1 and W2 was lower than at W3, W5, and W6, but similar to the pH of W4. DO concentrations at W1 and W2 were significantly lower than at other sites, whereas the highest DO was found at W4 in the Beilun estuary. The values of TN, NH4+-N, NO2--N, NO3--N and TOC measured at W1, W2, and W3 were significantly lower than those determined at W4, W5, and W6. According to the “Environmental quality standard for surface water” (GB3838-2002, China; Zhao and Yang, 2009), water samples collected at the urban canal sites W1 and W2 were hyper eutrophic.
A Principal Component Analyses was performed on the water properties at the different sampling locations. The first two principal components explained a large amount (85.7%) of the variation in the environmental factors and strong geographical clusters were formed (Figure 2). The first principal component (PC 1), which explained 74.3% of the variation, separated DO, salinity, and pH from the other parameters (TN, NH4+-N, NO2--N, NO3--N, and TOC), while the second principal component (PC 2) explaining 11.4% of the variation separated mainly DO from pH, salinity and TOC. The nitrogen-related data (e.g., nitrate, nitrite, ammonium, and TN) were not separated by these two principal components. The freshwater sampling locations W1 and W2 clustered together and were mostly characterized by higher concentrations of TN, NH4+-N, NO2--N, NO3--N, and TOC. The marine water locations W5 and W6 grouped also together at higher values of salinity, and pH. The brackish water sampling locations W3 and W4 did not cluster together. W3 was positioned between the freshwater and marine locations, whereas W4 separated from the all other locations by high oxygen concentrations.
FIGURE 2. Ordinate plot from principal-component analysis (PCA) base on water properties from six different sites. The values in parentheses are percentages of the total variances of PCA derived from water parameters.
Potential ammonia oxidation activities were also different among the fresh, brackish and marine water sites (Table 2). PAA was higher at the brackish and marine sites (W3, W4, W5, and W6) than at the freshwater sites (W1 and W2). PAA was positively correlated with DO (R = 0.539, n = 18, P < 0.05) and salinity (R = 0.768, n = 18, P < 0.01). However, PAA was negatively related to TN (R = -0.684, n = 18, P < 0.01), NO3--N (R = -0.659, n = 18, P < 0.01), NO2--N (R = -0.572, n = 18, P < 0.05), NH4+-N (R = -0.655, n = 18, P < 0.01) and TOC (R = -0.636, n = 18, P < 0.01). There was no significant correlation between PAA and pH (R = 0.439, n = 18, P = 0.068).
Spatial Variation in amoA Abundance
The abundance of AOA-amoA and AOB-amoA genes was different in the three water areas (Table 3). For AOA, amoA gene copy numbers per liter of water varied between 2.31 × 106 and 7.33 × 107. The lowest number of AOA-amoA was found in water of the Yong River next to Ningbo University (W3). AOA-amoA gene numbers at W4, W5, and W6 were significantly higher than those observed at W1, W2, and W3. AOA-amoA gene copy numbers positively related to salinity (R = 0.818, n = 18, P < 0.01), DO (R = 0.716, n = 18, P < 0.01), and PAA (R = 0.697, n = 18, P < 0.01), respectively. AOA-amoA gene copy numbers were negatively correlated with TN (R = -0.794, n = 18, P < 0.01), NO3--N (R = -0.941, n = 18, P < 0.01), NO2--N (R = -0.698, n = 18, P < 0.01), NH4+-N (R = -0.810, n = 18, P < 0.01) and TOC (R = -0.732, n = 18, P < 0.01), respectively.
qPCR assays showed that bacterial amoA gene copy numbers per liter of water ranged from 2.97 × 105 to 6.09 × 105. The highest number of AOB-amoA was found in water of the Yong River next to Ningbo University (W3). AOB-amoA gene numbers in the brackish and marine areas (W3, W4, W5, and W6) were significantly higher than those observed in the freshwater sites W1 and W2. Positive correlations were found between AOB-amoA gene copy numbers and salinity (R = 0.895, n = 18, P < 0.01), pH (R = 0.515, n = 18, P < 0.01) and DO (R = 0.752, n = 18, P < 0.01), respectively, but AOB-amoA gene copy numbers negatively correlated with TN (R = -0.904, n = 18, P < 0.01), NO2--N (R = -0.819, n = 18, P < 0.01), and NH4+-N (R = -0.883, n = 18, P < 0.01). No significant correlation was observed between AOB-amoA gene copy numbers and PAA (R = 0.326, n = 18, P = 0.187) and TOC (R = -0.360, n = 18, P = 0.142), respectively.
Ammonia-oxidizing archaea-amoA gene numbers were always more abundant than those of AOB in all water samples (Table 3). Especially, the ratios of AOA/AOB at the more marine sites (W4, W6, and W5) were much higher than those at the other sites.
Composition of AOA-amoA Genotypes
The AOA community composition was analyzed by PCR-DGGE (Figures 3 and 4). The triplicate banding profiles of each sampling sites clustered together. The patterns of AOA at the freshwater sites W1 and W2 were distinct from those observed at the brackish sites W3 and W4, and at the coastal sites W5 and W6. The similarity between the freshwater sites and the other sites was only 31%, whereas one of the marine sites (W5) showed 64% similarity with a cluster of sites that was comprised of the other marine site (W6) and the brackish sampling sites (W3 and W4), which showed 76% mutual similarity.
FIGURE 3. Denaturing gradient gel electrophoresis profiles of archaeal amoA genes in surface water collected from the sampling sites W1, W2, W3, W4, W5, and W6 along the Yong River. Bands used for sequencing and phylogenetic analysis are highlighted.
FIGURE 4. Cluster analysis unweighted pair group method with mathematical averages (UPGMA) of ammonia-oxidizing archaea (AOA) community composition based on the DGGE band pattern in surface water from the different sampling sites W1, W2, W3, W4, W5, and W6. The numbers in the figure indicate the clustering similarities.
The Shannon, Evenness and Simpson indices were used to establish the ammonia oxidizers diversity (Table 4). Of the AOA communities, the Shannon and Simpson indices were significantly higher at the freshwater site W1 than at all the other sites. In contrast to this difference in similarity indices, the Evenness index was significantly higher at the marine site W6 than at the other sites. The Shannon index was significantly and positively correlated with TN (R = 0.646, n = 18, P < 0.01), NO2--N (R = 0.530, n = 18, P < 0.05) and NH4+-N (R = 0.622, n = 18, P < 0.01), but negatively with pH (R = -0.491, n = 18, P < 0.05) and DO (R = -0.552, n = 18, P < 0.05). A significantly positive correlation was found between the Evenness index and salinity (R = 0.498, n = 18, P < 0.05), and a negative correlation between this index and TN (R = -0.475, n = 18, P < 0.05) and NH4+-N (R = -0.475, n = 18, P < 0.05).
TABLE 4. Diversity properties of AOA and AOB communities in the water layer at the different sampling locations based on DGGE band pattern data.
Twenty-one bands of AOA-amoA gene were excised from the DGGE gels (Figure 3). The recovered sequences (629 bp) were used to build a phylogenetic tree (Figure 5). The sequences were affiliated with sequences found in marine and coastal habitats, and in terrestrial environments. Most of sequences collected from freshwater belong to the cluster from terrestrial environments, only the sequence of band 5 fell into the marine and coastal environments. The sequences from brackish and coastal areas were mainly grouped in the cluster from marine and coastal environments, which also included the cultures of Candidatus Nitrosoarchaeum koreensis MY1 and Candidatus Nitrosoarchaeum limnia, and two Crenarchaeote enrichment culture clones collected from the San Francisco Bay estuary. But some of the sequences from brackish and coastal areas were also distributed in the cluster from terrestrial environments, overlapping with those from freshwater samples.
FIGURE 5. Phylogenetic relationships among archaeal amoA gene sequences retrieved from surface waters. Designation of the clones in bold includes the following information: accession number in the EMBL with DGGE band label. Colored points next to each retrieved gene sequences indicate the sampling station. Bootstrap values ≥ 50% (1000 neighbor-joining bootstraps) are shown at the branch points. The scale bar represents 5% estimated sequence divergence. The tree was rooted with Nitrosospira briensis C-128 (NBU76553).
Composition of AOB-amoA Genotypes
Similar to the AOA patterns, the triplicate banding profiles of each sampling sites clustered together and the AOB community compositions in at the freshwater sites W1 and W2 were again distinct from those of the brackish and marine sites W3, W4, W5, and W6 (Figures 6 and 7). The similarity in DGGE patterns between the freshwater sites and the other sites was only 35%. With a similarity of 47%, both freshwater sites were more dissimilar from each other than the estuarine and marine sites that shared a similarity of higher than 74%.
FIGURE 6. Denaturing gradient gel electrophoresis profiles of bacterial amoA in surface water collected from the sampling sites W1, W2, W3, W4, W5, and W6 along the Yong River. Bands used for sequencing and phylogenetic analysis are highlighted.
FIGURE 7. Cluster analysis (UPGMA) of AOB community composition based on the DGGE band pattern in surface waters from the sampling sites W1, W2, W3, W4, W5, and W6. The numbers in the figure indicate the clustering similarity.
The diversity and evenness indices of bacterial ammonia oxidizers diversity in freshwater (W1 and W2) were obviously different from brackish (W3 and W4) and coastal areas (W5 and W6; Table 4). In contrast to AOA, the AOB Shannon index was lower at W1 and W2 than at the other sites, and the Evenness index was higher at W1 and W2 compared with the other sites. The Simpson index was significantly lower at W1 than at the other sites. The Shannon was significantly and positively correlated with salinity (R = 0.606, n = 18, P < 0.01), pH (R = 0.504, n = 18, P < 0.05) and DO (R = 0.741, n = 18, P < 0.01), and negatively with TN (R = -0.846, n = 18, P < 0.01), NO2--N (R = -0.699, n = 18, P < 0.01) and NH4+-N (R = -0.815, n = 18, P < 0.01), respectively. A significantly positive correlation was found between the Evenness index and salinity (R = 0.498, n = 18, P < 0.05), and a negative correlation between this index and TN (R = -0.475, n = 18, P < 0.05) and NH4+(R = -0.475, n = 18, P < 0.05).
A total of nineteen bands of the AOB-amoA gene were excised from the DGGE gels (Figure 6). The recovered sequences (491 bp) were used to build a phylogenetic tree (Figure 8). A neighbor-joining analysis indicated that the AOB-amoA sequences were grouped into a Nitrosomonas and a Nitrosospira lineage. The sequences from the freshwater sites fell into the Nitrosomonas lineage. Furthermore, the Nitrosomonas lineage could be divided into two clusters (I and II). Some of the sequences fell into cluster I, including also Nitrosomonas sp. Nm41, Nitrosomonas communis, and Nitrosomonas sp. Nm33. Sequences distantly related to Nitrosomonas sp. ENI-11 and Nitrosomonas sp. GH22 grouped into cluster II. With the exception of the sequences of the marine bands 17, 18, and 19, which clustered into Nitrosomonas cluster II, the majority of the AOB-amoA sequences recovered from W3, W4, W5, and W6 was affiliated with the Nitrosospira lineage.
FIGURE 8. Phylogenetic relationship among bacterial amoA gene sequences retrieved from the surface waters along Yong River. Designation of the clones in bold includes the following information: accession number in the EMBL with DGGE band label. Colored points next to each retrieved gene sequences indicate the sampling station. Bootstrap values ≥ 50% (1000 neighbor-joining bootstrap) are shown at the branch points. The scale bar represents 5% estimated sequence divergence. The tree was rooted with Nitrosococcus sp. C-113 (AF153344).
The abundance and composition of AOA-amoA and AOB-amoA genes in surface waters were determined along the steep physicochemical freshwater to marine water gradient in Ningbo City, China. Clear distinctions among dominant ammonia oxidizers as observed from the DGGE patterns are revealed between freshwater, brackish, and coastal surface waters. Despite the relatively low resolution of DGGE banding patterns, the technique is still ideal as an initial screening method for comparing multiple samples with its low costs and short processing time (Green et al., 2010), provided that the quality of DGGE profile is high enough. Furthermore, when using marker genes DGGE fingerprinting may be an appropriate method to indicate qualitatively and semi-quantitatively the predominant functional microbial groups in ecosystems. It has successfully been applied to study the responses of ammonia monooxygenase (amoA) genes to environmental factors in many studies. In our study, clear shifts in predominant AOA and AOB were observed using the PCR-DGGE method and the detected sequences have been analyzed and combined with qPCR data and potential nitrification rates. All these data together showed the responses of the ammonia-oxidizing community on the environmental gradient in eutrophic waters.
Relative Abundance of AOA and AOB in Relation to Different Environmental Factors
Compared with samples from the brackish (W3 and W4) and coastal (W5 and W6) surface waters, high concentrations of NH4+-N, NO2--N, NO3--N, and TOC, and low amounts of DO were observed in samples from the hyper eutrophic waters of W1 and W2 that are situated in long-term stagnant canals (Zhao and Yang, 2009; Table 2). The lower amoA gene copy numbers of both AOA and AOB in these hyper eutrophic surface waters may be due to one or a combination of these environmental factors.
Ammonia-oxidizing archaea and AOB in surface waters were strongly impacted by salinity. In contrast to the AOB-amoA gene abundances, AOA-amoA gene numbers in Yong River varied markedly between sampling areas. From freshwater to marine sites AOA increased as much as 13–22 times, whereas AOB only increased 1–2 times. Opposite to our observation, Santoro et al. (2008) observed a nearly constant abundance of AOA across an estuarine gradient at Huntington Beach, while the abundance of AOB was dramatically lower in the freshwater stations compared with saline stations. Bernhard et al. (2007) detected lowest numbers of AOB at the low salinity site of the Plum Island Sound estuary. In a subsequent study of AOA in the same estuary, Bernhard et al. (2010) observed the highest abundance of AOA at the high salinity site in April, but at intermediate salinities in August/September. Such differences between AOA and AOB abundances are reflected in the AOA/AOB ratios. In our study, the ratios of AOA/AOB in W1 (14.74), W2 (9.98), and W3 (3.80) were significantly lower than in W4 (191.88), W5 (143.71), and W6 (153.43) (Table 3), which is in contrast with the AOA/AOB ratios that were reported for Plum Island Sound (Bernhard et al., 2010) San Francisco Bay (Mosier and Francis, 2008) and Huntington Beach (Santoro et al., 2008). At these other locations, the relative numbers of AOA decreased with increasing salinity. The difference between the Yong River and the other estuaries may be related to the physical environment, which was the water column for Yong River and sediments for the other studies.
The ratio of AOA/AOB declined with increasing ammonium concentrations. In general AOB dominate at higher ammonium concentrations as shown in a simulated creek ecosystem (Herrmann et al., 2011), lakes (Vissers et al., 2013; Bollmann et al., 2014), wetlands (Sims et al., 2012) and soils (Verhamme et al., 2011), and in the ammonia-rich area of Sacramento River (Damashek et al., 2015). Therefore, our findings support the idea that AOA are more adapted than AOB to oligotrophic environments (Martens-Habbena et al., 2009).
Potential Niche Differentiation between AOA and AOB
To estimate the ammonia oxidation capacity of the ammonia-oxidizing microorganisms, the potential ammonia-oxidizing activity was determined as described before (Kasuga et al., 2010). Since the method of Kasuga et al. (2010) uses ammonium consumption rates and not nitrite production rates for measuring potential ammonia oxidation, some of the ammonium may have been used by heterotrophic bacteria. However, in activated sludge systems only 10–20% of the ammonium was allocated to ammonium removal by heterotrophic organisms (cited from Kasuga et al., 2010). Potential ammonia oxidation rates do not represent in situ rates, but give rather sizes of the active ammonia-oxidizing microbial community and hence can be used for comparison between sites (see also Damashek et al., 2015). The data indicated that PAAs in the freshwater samples (W1 and W2) were lower than those measured in the brackish and coastal water samples, which was different from previous studies that showed increased nitrification activities in freshwater samples compared to more marine environments (Seitzinger, 1988; Bollmann and Laanbroek, 2002; Bernhard et al., 2005). This can be due to the hyper eutrophic conditions and lower DO concentrations (Rysgaard et al., 1994; Liikanen and Martikainen, 2003) in the freshwater part of the Yong River estuary that may inhibit the ammonia oxidation.
The PAA was positively correlated with AOA abundance, but not with AOB abundance, which was consistent with previous studies (Caffrey et al., 2007; Lam et al., 2007). Such a correlation between AOA abundance and PAA, could suggest that AOA were the most active contributors of ammonia oxidation in the surface waters in Yong River. The potential of AOA involved in nitrogen removal was already reported by previous studies (Herrmann et al., 2008; You et al., 2009). However, some studies supported that AOB were the major contributors to ammonia oxidation (Jia and Conrad, 2009; Zhang et al., 2010), especially in high-ammonium habitats where AOB were quantitatively dominant (Wells et al., 2009; Limpiyakorn et al., 2011), and where the potential nitrification activity was positively correlated with AOB, but not with AOA (Herrmann et al., 2011).
Composition of AOA and AOB Communities in Relation to Environmental Factors
A distinct difference in the composition of the AOA and AOB communities was observed between freshwater on one site and brackish or coastal areas on the other site (Figures 5 and 8). It suggested that ammonia oxidizers community composition is shaped by environmental factors. A relationship between amoA composition and ecological niches was described in previous studies (Biller et al., 2012; Fernàndez-Guerra and Casamayor, 2012; Hatzenpichler, 2012; Cao et al., 2013).
Within the estuary, a relatively high similarity of AOA community composition was observed between the brackish (W3 and W4) and coastal (W5 and W6) areas (Figures 3 and 4). The similarity in community composition between the brackish and marine areas can be due to the intrusion of coastal water, creating similar estuarine environments. A comparable AOA community composition was found at the mouth of the Changjiang Estuary, which flows into the East China Sea approximately 150 km north of Ningbo City (Dang et al., 2008).
However, a clear shift of AOA composition was observed from freshwater to brackish and marine waters. Salinity was considered the most important factor affecting AOA community structure based on amoA sequence clustering in aquatic habitats (Cao et al., 2013; Xie et al., 2014). Other environmental factors such as oxygen (Molina et al., 2010; Biller et al., 2012, ammonium concentration (Verhamme et al., 2011) and pH (Nicol et al., 2008) could also contribute to the change in AOA community composition.
The AOA-amoA DNA sequences derived from Bands 10 to 21 of DGGE gel, which were obtained from brackish and coastal water samples were highly similar to those from an aquarium system (Sauder et al., 2011), the Changjiang estuary (Dang et al., 2008), the East China Sea (Hu et al., 2011), a wastewater treatment (Bai et al., 2012), hot springs (Zhang et al., 2008; Jiang et al., 2010), the Elkhorn Slough estuary (Wankel et al., 2011), the Pearl River estuary (Jin et al., 2011). The AOA-amoA sequences selected from Bands 1 to 9 obtained from freshwater had a high similarity to those collected from a wastewater treatment (Bai et al., 2012), hot spring (Zhang et al., 2008; Jiang et al., 2010), the Elkhorn Slough estuary (Wankel et al., 2011), fresh-water aquaria (Sauder et al., 2011), the Changjiang estuary (Dang et al., 2008), the paddy soil (Fujii et al., 2010), temperate forest soil (Onadera et al., 2010) and a littoral wetland zone at Baiyangdian Lake (Wang et al., 2012).
A relatively high similarity in AOB community was also found between the brackish and coastal areas (Figures 6 and 7). However, the Shannon diversity indices of the AOB communities in the brackish and coastal areas were significantly higher than those at calculated from the freshwater stations, which was in contrast with a previous study that showed the loss of diversity of AOB with increasing salinity in Plum Island Sound estuary (Bernhard et al., 2005). The collected sequences of Nitrosomonas clustered to sequences obtained from the Yangtze Estuary (Zheng et al., 2014), the Chongming eastern tidal flat (Zheng et al., 2013), a wastewater treatment plant (Gao et al., 2014), a paddy field soil (GenBank unpublished), and Dianchi Lake (GenBank unpublished). One of the Nitrosospira-related sequences from the marine area showed 100% identity with sequences from the Chongming eastern tidal flat (Zheng et al., 2013) and the Jiaozhou Bay (Dang et al., 2010), and 99% identity with sequences from the Chongming eastern tidal (Zheng et al., 2013), the Yangtze Estuary (Zheng et al., 2014), red soil (GenBank unpublished) and the Changjiang Estuary (GenBank unpublished), respectively.
Nitrosomonas-related sequences were predominant in the freshwater part of the Yong River, whereas Nitrosospira-associated sequences came to the fore in the more saline parts of the estuary (Figure 8). A shift from a Nitrosomonas-dominated community to a Nitrosospira-dominated community has been observed before in the Scheldt River when the estuary changes from fresh to brackish to marine conditions (de Bie et al., 2001). From a study of ammonia oxidizers at Huntington Beach, it was suggested that changes in the amount of DO caused the shift in the ammonia oxidizer community (Santoro et al., 2008). However, the ammonia-oxidizing bacterial community in the river Scheldt was impacted by salt rather than by oxygen (Bollmann and Laanbroek, 2002). Organic and nutrient loads of wastewater shaped significantly the AOB community structure by in the eutrophic Tokyo Bay (Urakawa et al., 2006). In our study on the hyper eutrophic Yong River, the composition of the AOB communities at freshwater sites was dissimilar from those at the brackish and marine sampling locations. The freshwater sites were characterized by relatively high concentrations of nitrogen and low salinities and concentrations of oxygen. Hence, it seems likely that these factors govern the structure of the AOB communities leading from predominance of Nitrosomonas-related sequences in the nitrogen-rich, but oxygen-limited freshwater parts to predominance Nitrosospira-dominated communities in the nitrogen-poor, but more oxygen-rich, saline waters. Therefore, our study suggests that salinity is not the main contributing factor in determining AOB diversity in the low saline regions. Perhaps in more oligotrophic waters, salinity is the most important determinant for diversity of the AOB community composition, but in the hypereutrophic waters, other factors, such as DO, or levels of nitrogen and organic matter become more important, especially in the urban canals of the Yong River system where the surface waters are retained for a long period.
A clear shift in ammonia-oxidizing activity was observed along the estuary of the River Yong with lower activities in the freshwater, urban zones. Simultaneously with the shift in ammonia-oxidizing activity, shifts in the abundance and species composition of both the archaeal and bacterial ammonia oxidizers were observed together with an increase in the AOA/AOB ratio with increasing salinity. The shifts in activity, abundance and species composition were accompanied by shifts in ammonium and oxygen availabilities, which suggest that these factors together with salinity might impact the characteristics of the ammonia-oxidizing microbial communities. The difference in behavior between AOA and AOB along the urban, estuarine gradient might underlines the physiological differences between these groups. Hence, a better understanding of the nitrogen-converting processes in such highly eutrophic ecosystems is depending on a further exploration of the physiological abilities of AOA and AOB.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
This study was financially supported by National Natural Science Foundation of China (grant no. 41571252), Science Technology Department of Zhejiang and Ningbo (grant nos. 2012C23065, 2011C50050) and Natural Science Foundation of Zhejiang and Ningbo (grant nos. Y5110331, 2011A610085).
Bai, Y., Sun, Q., Wen, D., and Tang, X. (2012). Abundance of ammonia-oxidizing bacteria and archaea in industrial and domestic wastewater treatment systems. FEMS Microbiol. Ecol. 80, 323–330. doi: 10.1111/j.1574-6941.2012.01296.x
Beman, J. M., and Francis, C. A. (2006). Diversity of ammonia-oxidizing archaea and bacteria in the sediments of a hypernutrified subtropical estuary: Bahia del Tobari, Mexico. Appl. Environ. Microbiol. 72, 7767–7777. doi: 10.1128/AEM.00946-06
Bernhard, A. E., Donn, T., Giblin, A. E., and Stahl, D. A. (2005). Loss of diversity of ammonia-oxidizing bacteria correlates with increasing salinity in an estuary system. Environ. Microbiol. 7, 1289–1297. doi: 10.1111/j.1462-2920.2005.00808.x
Bernhard, A. E., Landry, Z. C., Blevins, A., José, R., Giblin, A. E., and Stahl, D. A. (2010). Abundance of ammonia-oxidizing archaea and bacteria along an estuarine salinity gradient in relation to potential nitrification rates. Appl. Environ. Microbiol. 76, 1285–1289. doi: 10.1128/AEM.02018-09
Bernhard, A. E., Tucker, J., Giblin, A. E., and Stahl, D. A. (2007). Functionally distinct communities of ammonia-oxidizing bacteria along an estuarine salinity gradient. Environ. Microbiol. 6, 1439–1447. doi: 10.1111/j.1462-2920.2007.01260.x
Biller, S. J., Mosier, A. C., Wells, G. F., and Francis, C. A. (2012). Global biodiversity of aquatic ammonia-oxidizing archaea is partitioned by habitat. Front. Microbiol. 3:252. doi: 10.3389/fmicb.2012.00252
Bollmann, A., Bullerjahn, G. S., and McKay, R. M. (2014). Abundance and diversity of ammonia-oxidizing archaea and bacteria in sediments of trophic end members of the Laurentian Great Lakes. Erie and Superior. PLoS ONE 9:e97068. doi: 10.1371/journal.pone.0097068
Bollmann, A., and Laanbroek, H. J. (2002). Influence of oxygen partial pressure and salinity on the community composition of ammonia-oxidizing bacteria in the Schelde estuary. Aquat. Microbiol. Ecol. 28, 239–247. doi: 10.3354/ame028239
Caffrey, J. M., Bano, N., Kalanetra, K., and Hollibaugh, J. T. (2007). Ammonia oxidation and ammonia-oxidizing bacteria and archaea from estuaries with differing histories of hypoxia. ISME J. 1, 660–662. doi: 10.1038/ismej.2007.79
Cao, H., Hong, Y., Li, M., and Gu, J. D. (2011). Diversity and abundance of ammonia-oxidizing prokaryotes in sediments from the coastal Pearl River estuary to the South China Sea. Anton. Leeuw. 100, 545–556. doi: 10.1007/s10482-011-9610-1
Damashek, J., Smith, J. M., Mosier, A. C., and Francis, C. A. (2015). Benthic ammonia oxidizers differ in community structure and biogeochemical potential across a riverine delta. Front. Microbiol. 5:743. doi: 10.3389/fmicb.2014.00743
Dang, H., Li, J., Chen, R., Wang, L., Guo, L., Zhang, Z., et al. (2010). Diversity, abundance, and spatial distribution of sediment ammonia-oxidizing beta proteobacteria in response to environmental gradients and coastal eutrophication in Jiaozhou Bay. China Appl. Environ. Microbiol. 76, 4691–4702. doi: 10.1128/AEM.02563-09
Dang, H., Zhang, X., Sun, J., Li, T., Zhang, Z., and Yang, G. (2008). Diversity and spatial distribution of sediment ammonia-oxidizing crenarchaeota in response to estuarine and environmental gradients in the Changjiang Estuary and East China Sea. Microbiology 154, 2084–2095. doi: 10.1099/mic.0.2007/013581-0
de Bie, M. J., Speksnijder, A. G., Kowalchuk, G. A., Schuurman, T., Zwart, G., Stephen, J. R., et al. (2001). Shifts in the dominant populations of ammonia-oxidizing β-subclass Proteobacteria along the eutrophic Schelde estuary. Aquat. Microbiol. Ecol. 23, 225–236. doi: 10.3354/ame023225
Dumonceaux, T. J., Hill, J. E., Briggs, S. A., Amoako, K. K., Hemmingsen, S. M., and Kessel, A. G. V. (2006). Enumeration of specific bacterial populations in complex intestinal communities using quantitative PCR based on the chaperonin-60 target. J. Microbiol. Methods 64, 46–62. doi: 10.1016/j.mimet.2005.04.006
Erguder, T. H., Boon, N., Wittebolle, L., Marzorati, M., and Verstraete, W. (2009). Environmental factors shaping the ecological niches of ammonia-oxidizing archaea. FEMS Microbiol. Rev. 33, 855–869. doi: 10.1111/j.1574-6976.2009.00179.x
Francis, C. A., Beman, J. M., and Kuypers, M. M. (2007). New processes and players in the nitrogen cycle: the microbial ecology of anaerobic and archaeal ammonia oxidation. ISME J. 1, 19–27. doi: 10.1038/ismej.2007.8
Francis, C. A., O’Mullan, G. D., and Ward, B. B. (2003). Diversity of ammonia monooxygenase (amoA) genes across environmental gradients in Chesapeake Bay sediments. Geobiology 1, 129–140. doi: 10.1046/j.1472-4669.2003.00010.x
Francis, C. A., Roberts, K. J., Beman, J. M., Santoro, A. E., and Oakley, B. B. (2005). Ubiquity and diversity of ammonia-oxidizing archaea in water columns and sediments of the ocean. Proc. Nat. Acad. Sci. U.S.A. 102, 14683–14688. doi: 10.1073/pnas.0506625102
Freitag, T. E., Chang, L., and Prosser, J. I. (2006). Changes in the community structure and activity of beta proteobacterial ammonia-oxidizing sediment bacteria along a freshwater-marine gradient. Environ. Microbiol. 8, 684–696. doi: 10.1111/j.1462-2920.2005.00947.x
Fujii, C., Nakagawa, T., Onodera, Y., Matsutani, N., Sasada, K., Takahashi, R., et al. (2010). Succession and community composition of ammonia-oxidizing archaea and bacteria in bulk soil of a Japanese paddy field. Soil Sci. Plant Nutr. 56, 212–219. doi: 10.1111/j.1747-0765.2010.00449.x
Gao, J., Luo, X., Wu, G., Li, T., and Peng, Y. (2014). Abundance and diversity based on amoA genes of ammonia-oxidizing archaea and bacteria in ten wastewater treatment systems. Appl. Microbiol. Biotechnol. 98, 3339–3354. doi: 10.1007/s00253-013-5428-2
Green, S. J., Leigh, M. B., and Neufeld, J. D. (2010). “Denaturing gradient gel electrophoresis (DGGE) for microbial community analysis,” in Handbook of Hydrocarbon and Lipid Microbiology, eds K. N. Timmis, T. J. McGenity, J. R. van der Meer, and V. de Lorenzo (Berlin Heidelberg: Springer), 4137–4158.
Hammer, Ø., Harper, D. A. T., and Ryan, P. D. (2001). PAST: Paleontological statistics software package for education and data analysis. Palaeontol. Electron. 4:9. Available at: http://palaeo-electronica.org/2001_1/past/issue1_01.htm
Herrmann, M., Saunders, A. M., and Schramm, A. (2008). Archaea dominate the ammonia-oxidizing community in the rhizosphere of the freshwater macrophyte Littorella uniflora. Appl. Environ. Microbiol. 74, 3279–3283. doi: 10.1128/AEM.02802-07
Herrmann, M., Scheibe, A., Avrahami, S., and Küsel, K. (2011). Ammonium availability affects the ratio of ammonia-oxidizing bacteria to ammonia-oxidizing archaea in simulated creek ecosystems. Appl. Environ. Microbiol. 77, 1896–1899. doi: 10.1128/AEM.02879-10
Hu, A., Jiao, N., Zhang, R., and Yang, Z. (2011). Niche partitioning of marine group I Crenarchaeota in the euphotic and upper mesopelagic zones of the East China Sea. Appl. Environ. Microbiol. 77, 7469–7478. doi: 10.1128/AEM.00294-11
Jiang, H., Huang, Q., Dong, H., Wang, P., Wang, F., Li, W., et al. (2010). RNA-based investigation of ammonia-oxidizing archaea in hot springs of Yunnan Province. China Appl. Environ. Microbiol. 76, 4538–4541. doi: 10.1128/AEM.00143-10
Jin, T., Zhang, T., Ye, L., Lee, O. O., Wong, Y. H., and Qian, P. Y. (2011). Diversity and quantity of ammonia-oxidizing archaea and bacteria in sediment of the Pearl River Estuary. China Appl. Microbiol. Biotechnol. 90, 1137–1145. doi: 10.1007/s00253-011-3107-8
Kasuga, I., Nakagaki, H., Kurisu, F., and Furumai, H. (2010). Predominance of ammonia-oxidizing archaea on granular activated carbon used in a full-scale advanced drinking water treatment plant. Water Res. 44, 5039–5049. doi: 10.1016/j.watres.2010.07.015
Lam, P., Jensen, M. M., Lavik, G., McGinnis, D. F., Müller, B., Schubert, C. J., et al. (2007). Linking crenarchaeal and bacterial nitrification to anammox in the Black Sea. Proc. Natl. Acad. Sci. U.S.A. 104, 7104–7109. doi: 10.1073/pnas.0611081104
Liikanen, A., and Martikainen, P. J. (2003). Effect of ammonium and oxygen on methane and nitrous oxide fluxes across sediment–water interface in a eutrophic lake. Chemosphere 52, 1287–1293. doi: 10.1016/S0045-6535(03)00224-8
Limpiyakorn, T., Sonthiphand, P., Rongsayamanont, C., and Polprasert, C. (2011). Abundance of amoA genes of ammonia-oxidizing archaea and bacteria in activated sludge of full-scale wastewater treatment plants. Bioresource Technol. 102, 3694–3701. doi: 10.1016/j.biortech.2010.11.085
Liu, S., Shen, L., Lou, L., Tian, G., Zheng, P., and Hu, B. (2013). Spatial distribution and factors shaping the niche segregation of ammonia-oxidizing microorganisms in the Qiantang River. China Appl. Environ. Microbiol. 79, 4065–4071. doi: 10.1128/AEM.00543-13
Martens-Habbena, W., Berube, P. M., Urakawa, H., José, R., and Stahl, D. A. (2009). Ammonia oxidation kinetics determine niche separation of nitrifying Archaea and Bacteria. Nature 461, 976–979. doi: 10.1038/nature08465
Moin, N. S., Nelson, K. A., Bush, A., and Bernhard, A. E. (2009). Distribution and diversity of archaeal and bacterial ammonia oxidizers in salt marsh sediments. Appl. Environ. Microbiol. 75, 7461–7468. doi: 10.1128/AEM.01001-09
Molina, V., Belmar, L., and Ulloa, O. (2010). High diversity of ammonia-oxidizing archaea in permanent and seasonal oxygen-deficient waters of the eastern South Pacific. Environ. Microbiol. 12, 2450–2465. doi: 10.1111/j.1462-2920.2010.02218.x
Mosier, A. C., and Francis, C. A. (2008). Relative abundance and diversity of ammonia-oxidizing archaea and bacteria in the San Francisco Bay estuary. Environ. Microbiol. 10, 3002–3016. doi: 10.1111/j.1462-2920.2008.01764.x
Nicol, G. W., Leininger, S., Schleper, C., and Prosser, J. I. (2008). The influence of soil pH on the diversity, abundance and transcriptional activity of ammonia oxidizing archaea and bacteria. Environ. Microbiol. 10, 2966–2978. doi: 10.1111/j.1462-2920.2008.01701.x
Onadera, Y., Nakagawa, T., Takahashi, R., and Tokuyama, T. (2010). Seasonal change in vertical distribution of ammonia-oxidizing archaea and bacteria and their nitrification in temperate forest soil. Microbes Environ. 25, 28–35. doi: 10.1264/jsme2.ME09179
Prosser, J. I., and Nicol, G. W. (2012). Archaeal and bacterial ammonia-oxidisers in soil: the quest for niche specialisation and differentiation. Trends. Microbiol. 20, 523–531. doi: 10.1016/j.tim.2012.08.001
Rotthauwe, J. H., Witzel, K. P., and Liesack, W. (1997). The ammonia monooxygenase structural gene amoA as a functional marker: molecular fine-scale analysis of natural ammonia-oxidizing populations. Appl. Environ. Microbiol. 63, 4704–4712. doi: 10.4014/jmb.1204.04047
Rysgaard, S., Risgaard-Petersen, N., Sloth, N. P., Kim, J., and Nielsen, L. P. (1994). Oxygen regulation of nitrification and denitrification in sediments. Limnol. Oceanogr. 39, 1643–1652. doi: 10.4319/lo.19220.127.116.113
Sahan, E., and Muyzer, G. (2008). Diversity and spatio-temporal distribution of ammonia-oxidizing archaea and bacteria in sediments of the Westerschelde estuary. FEMS Microbiol. Ecol. 64, 175–186. doi: 10.1111/j.1574-6941.2008.00462.x
Santoro, A. E., Francis, C. A., De Sieyes, N. R., and Boehm, A. B. (2008). Shifts in the relative abundance of ammonia- oxidizing bacteria and archaea across physicochemical gradients in a subterranean estuary. Environ. Microbiol. 10, 1068–1079. doi: 10.1111/j.1462-2920.2007.01547.x
Sauder, L. A., Engel, K., Stearns, J. C., Masella, A. P., Pawliszyn, R., and Neufeld, J. D. (2011). Aquarium nitrification revisited: thaumarchaeota are the dominant ammonia oxidizers in freshwater aquarium biofilters. PLoS ONE 6:e23281. doi: 10.1371/journal.pone.0023281
Seitzinger, S., Harrison, J. A., Böhlke, J. K., Bouwman, A. F., Lowrance, R., Peterson, B., et al. (2006). Denitrification across landscapes and waterscapes: a synthesis. Ecol. Appl. 16, 2064–2090. doi: 10.1890/1051-0761(2006)016[2064:DALAWA]2.0.CO;2
Sekiguchi, H., Watanabe, M., Nakahara, T., Xu, B., and Uchiyama, H. (2002). Succession of bacterial community structure along the Changjiang River determined by denaturing gradient gel electrophoresis and clone library analysis. Appl. Environ. Microbiol. 68, 5142–5150. doi: 10.1128/AEM.68.10.5142-5150.2002
Sims, A., Horton, J., Gajaraj, S., McIntosh, S., Miles, R. J., Mueller, R., et al. (2012). Temporal and spatial distributions of ammonia-oxidizing archaea and bacteria and their ratio as an indicator of oligotrophic conditions in natural wetlands. Water Res. 46, 4121–4129. doi: 10.1016/j.watres.2012.05.007
Tourna, M., Freitag, T. E., Nicol, G. W., and Prosser, J. I. (2008). Growth, activity and temperature responses of ammonia-oxidizing archaea and bacteria in soil microcosms. Environ. Microbiol. 10, 1357–1364. doi: 10.1111/j.1462-2920.2007.01563.x
Urakawa, H., Kurata, S., Fujiwara, T., Kuroiwa, D., Maki, H., Kawabata, S., et al. (2006). Characterization and quantification of ammonia-oxidizing bacteria in eutrophic coastal marine sediments using polyphasic molecular approaches and immunofluorescence staining. Environ. Microbiol. 8, 787–803. doi: 10.1111/j.1462-2920.2005.00962.x
Verhamme, D. T., Prosser, J. I., and Nicol, G. W. (2011). Ammonia concentration determines differential growth of ammonia-oxidizing archaea and bacteria in soil microcosms. ISME J. 15, 1067–1071. doi: 10.1038/ismej.2010.191
Vissers, E. W., Blaga, C. I., Bodelier, P. L., Muyzer, G., Schleper, C., Sinninghe Damsté, J. S., et al. (2013). Seasonal and vertical distribution of putative ammonia-oxidizing thaumarchaeotal communities in an oligotrophic lake. FEMS Microbiol. Ecol. 83, 515–526. doi: 10.1111/1574-6941.12013
Wang, Y., Zhu, G., Ye, L., Feng, X., den Camp, H. J. O., and Yin, C. (2012). Spatial distribution of archaeal and bacterial ammonia oxidizers in the littoral buffer zone of a nitrogen-rich lake. J. Environ. Sci. 24, 790–799. doi: 10.1016/S1001-0742(11)60861-9
Wankel, S. D., Mosier, A. C., Hansel, C. M., Paytan, A., and Francis, C. A. (2011). Spatial variability in nitrification rates and ammonia-oxidizing microbial communities in the agriculturally impacted Elkhorn Slough estuary. California. Appl. Environ. Microbiol. 77, 269–280. doi: 10.1128/AEM.01318-10
Wells, G. F., Park, H. D., Yeung, C. H., Eggleston, B., Francis, C. A., and Criddle, C. S. (2009). Ammonia-oxidizing communities in a highly aerated full-scale activated sludge bioreactor: beta proteobacterial dynamics and low relative abundance of Crenarchaea. Environ. Microbiol. 11, 2310–2328. doi: 10.1111/j.1462-2920.2009.01958.x
Xie, W., Zhang, C., Zhou, X., and Wang, P. (2014). Salinity-dominated change in community structure and ecological function of Archaea from the lower Pearl River to coastal South China Sea. Appl. Microbiol. Biotechnol. 98, 7971–7982. doi: 10.1007/s00253-014-5838-9
Zhang, C. L., Ye, Q., Huang, Z., Li, W., Chen, J., Song, Z., et al. (2008). Global occurrence of archaeal amoA genes in terrestrial hot springs. Appl. Environ. Microbiol. 74, 6417–6426. doi: 10.1128/AEM.00843-08
Zhang, L. M., Offre, P. R., He, J. Z., Verhamme, D. T., Nicol, G. W., and Prosser, J. I. (2010). Autotrophic ammonia oxidation by soil thaumarchaea. Proc. Natl. Acad. Sci. U.S.A. 107, 17240–17245. doi: 10.1073/pnas.1004947107
Zhang, Q. F., Peng, J. J., Chen, Q., Li, X. F., Xu, C. Y., Yin, H. B., et al. (2011). Impacts of Spartina alterniflora invasion on abundance and composition of ammonia oxidizers in estuarine sediment. J. Soil. Sediment. 11, 1020–1031. doi: 10.1007/s11368-011-0369-9
Zhao, Y. W., and Yang, Z. F. (2009). Integrative fuzzy hierarchical model for river health assessment: a case study of Yong River in Ningbo City. China Commun. Nonlinear. Sci. 14, 1729–1736. doi: 10.1016/j.cnsns.2007.09.019
Zheng, Y., Hou, L., Liu, M., Lu, M., Zhao, H., Yin, G., et al. (2013). Diversity, abundance, and activity of ammonia-oxidizing bacteria and archaea in Chongming eastern intertidal sediments. Appl. Microbiol. Biotechnol. 97, 8351–8363. doi: 10.1007/s00253-012-4512-3
Zheng, Y., Hou, L., Newell, S., Liu, M., Zhou, J., Zhao, H., et al. (2014). Community dynamics and activity of ammonia-oxidizing prokaryotes in intertidal sediments of the Yangtze Estuary. Appl. Environ. Microbiol. 80, 408–419. doi: 10.1128/AEM.03035-13
Keywords: ammonia-oxidizing archaea (AOA), ammonia-oxidizing bacteria (AOB), eutrophic status, salinity, surface water
Citation: Zhang Q, Tang F, Zhou Y, Xu J, Chen H, Wang M and Laanbroek HJ (2015) Shifts in the pelagic ammonia-oxidizing microbial communities along the eutrophic estuary of Yong River in Ningbo City, China. Front. Microbiol. 6:1180. doi: 10.3389/fmicb.2015.01180
Received: 19 July 2015; Accepted: 12 October 2015;
Published: 27 October 2015.
Edited by:Hongyue Dang, Xiamen University, China
Reviewed by:Anne Bernhard, Connecticut College, USA
Annika C. Mosier, University of Colorado Denver, USA
Copyright © 2015 Zhang, Tang, Zhou, Xu, Chen, Wang and Laanbroek. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Qiufang Zhang, email@example.com