key: cord-012085-ubdzhkfq
authors: Jin, Tao; Zhang, Tong; Ye, Lin; Lee, On On; Wong, Yue Him; Qian, Pei Yuan
title: Diversity and quantity of ammonia-oxidizing Archaea and Bacteria in sediment of the Pearl River Estuary, China
date: 2011-02-01
journal: Appl Microbiol Biotechnol
DOI: 10.1007/s00253-011-3107-8
sha: 
doc_id: 12085
cord_uid: ubdzhkfq

The diversity and abundance of ammonia-oxidizing archaea (AOA) and ammonia-oxidizing bacteria (AOB) in the sediment of the Pearl River Estuary were investigated by cloning and quantitative real-time polymerase chain reaction (qPCR). From one sediment sample S16, 36 AOA OTUs (3% cutoff) were obtained from three clone libraries constructed using three primer sets for amoA gene. Among the 36 OTUs, six were shared by all three clone libraries, two appeared in two clone libraries, and the other 28 were only recovered in one of the libraries. For AOB, only seven OTUs (based on 16S rRNA gene) and eight OTUs (based on amoA gene) were obtained, showing lower diversity than AOA. The qPCR results revealed that AOA amoA gene copy numbers ranged from 9.6 × 10(6) to 5.1 × 10(7) copies per gram of sediment and AOB amoA gene ranged from 9.5 × 10(4) to 6.2 × 10(5) copies per gram of sediment, indicating that the dominant ammonia-oxidizing microorganisms in the sediment of the Pearl River Estuary were AOA. The terminal restriction fragment length polymorphism results showed that the relative abundance of AOB species in the sediment samples of different salinity were significantly different, indicating that salinity might be a key factor shaping the AOB community composition. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1007/s00253-011-3107-8) contains supplementary material, which is available to authorized users.

Microorganisms in marine sediments contribute significantly to global cycles of organic and inorganic matters (Whitman et al. 1998) , including nitrogen cycle. Biological nitrogen fixation is the largest source of nitrogen to the ocean, whereas anoxic/anaerobic microbial processes are responsible for nitrogen losses (Gruber and Sarmiento 1997) . These two processes are connected by nitrification, the microbial oxidation of ammonia (NH 3 ) to nitrate (NO 3 − ) via nitrite (NO 2 − ). Coupled nitrification/ denitrification can remove a substantial percentage (10~80%) of anthropogenic nitrogen pollution in estuaries (Seitzinger 1988) . Previous studies have focused on the ammonia-oxidizing bacteria (AOB), mainly targeting both the 16S rRNA gene and amoA gene which encodes α subunit of the ammonia monooxygenase (AMO) (Purkhold et al. 2000) . All the known AOB are phylogenetically restricted to the βand γ-proteobacteria. Recently, it was discovered that many archaea in the phylum Crenarchaeota are also capable of performing ammonia oxidation like AOB (Könneke et al. 2005; Venter et al. 2004 ). Both ammonia oxidizing archaea (AOA) and AOB contains the key enzyme AMO which is composed of three subunits encoded by genes of amoA, amoB, and amoC (Klotz et al. 1997) and responsible for the conversion of ammonia to hydroxylamine. By performing the key first step in nitrification, both AOA and AOB play an important role in the global nitrogen cycle. However, estuaries and their intertidal sediments experience large fluctuations in various hydrological and chemical conditions. Therefore, the communities of AOA and AOB might show dynamic changes in response to environmental factors. A few studies demonstrated that diverse AOA/AOB phylotypes and distinct AOA/AOB communities existed in different marine environments both on a large geographical scale and in local estuarine gradients (Beman and Francis 2006; Francis et al. 2005; Sahan and Muyzer 2008) . Studies already showed that archaeal amoA gene copy numbers were about 10 to 1,000 times of those of bacterial amoA in the North Sea and North Atlantic (Wuchter et al. 2006) , indicating dominance of AOA in ocean environments. So far, most of the studies on AOA and AOB in marine water and sediment focused on the coast of the East Pacific Ocean and Atlantic. The diversity and spatial distribution of AOA and AOB are still unknown in other ocean areas of the world, including the Pearl River Estuary.

The Pearl River, the third largest river in China, stretching for 2,214 km and draining an area of 452,000 km 2 , flows into the Pearl River Estuary through eight entrances. The Pearl River Estuary is located in the subtropical region in Guangdong Province, Southern China. It represents one of the most important and complex ecosystems linking the highly developing land area and the South China Sea. In recent years, the Pearl River has a high load of anthropogenic nutrients from increased agricultural activities, fish dike farming, and wastewater runoff due to the increase in population and economic development in Southern China and the Pearl River Delta region (Huang et al. 2003) . Previous study also showed high nitrification rates in summer (ammonia oxidation rate~1.5-33.1 μmol NL −1 day −1 , nitrite oxidation rate~0.6-32.0 μmol NL −1 day −1 ) (Dai et al. 2008 ). However, the microorganisms responsible for the nitrification in this area are little known.

In this study, five sediment samples were collected in the Pearl River Estuary from different sites along the salinity gradient (Fig. 1) . The diversity of AOA and AOB were revealed based on clone libraries. Quantification of AOA and AOB were performed using quantitative PCR (qPCR). For AOA, three primer sets of amoA gene were selected for construction of clone libraries and qPCR to evaluate the specificity and universality of these primers. Terminal restriction fragment length polymorphism (T-RFLP) was used to analyze the spatial distribution of different AOB populations.

Surface sediment samples were collected from five locations along a transect from the Pearl River Estuary to the South China Sea in the summer (June 2009; Fig. 1 ) using a 250-cm 3 Van Veen Grab (KC Denmark, Denmark). Upon collection of the sediment samples, about 2 kg of sediment were transferred to sterile plastic bags and frozen at −80°C. Frozen sediment samples were then transferred back to the laboratory with dry ice. Except for the depth and temperature which were measured on site, other parameters were measured using the pore water of the sediments. Concentrations of ammonium, nitrite, and nitrate were measured according to the Standard Methods (APHA 2005) by Nesslerization method, colorimetric method, and ultraviolet spectrophotometric screening method, respectively. For AOB 16S rRNA and amoA genes, PCR was performed in a total volume of 25 μL containing 20 ng of DNA template, 2.5 μL buffer, 250 mM of each deoxyribonucleotide triphosphate, 1 U of ExTaq (TaKaRa) DNA polymerase, and 0.2 μM of each primer. The major cycling program for each primer set was listed in Table 1 . The presence and sizes of the PCR amplification products were determined by agarose (1%) gel electrophoresis. PCR products were purified by using Quick-spin™ PCR Products Purification Kit (iNtRON Biotechnology, Seongnam-Si, Korea). Then, the PCR products were cloned using InsT/Aclone™ PCR Cloning Kit (Fermentas, Vilnius, Lithuania). White colonies were selected for insert screening by using PCR with primers M13F and M13R. The sequencing of selected clones was performed by using M13F primer on ABI 3730xl capillary sequencers (PE Applied Biosystems, Foster City, USA).

The estimation of diversity index in each library was determined using the ESTIMATES V8.0 (http://viceroy.eeb. uconn.edu/estimates). The Chao 1 estimator of species richness was calculated after 1,000 randomizations of sampling, without replacement. The percentage of coverage was calculated using Good's method (Good 1953) . Operational taxonomic unit (OTU) in this study was defined as a sequence group in which sequences differed by ≤3%.

The nucleotide sequences were compared with those from the GenBank using BLASTn in the National Center for the Biotechnology Information (NCBI) server (http://www. ncbi.nlm.nih.gov). The sequences in this study and the reference sequences retrieved from the GenBank were aligned by ARB (http://www.arb-home.de) (Ludwig et al. 2004 ) to construct phylogenetic trees using the neighborjoining method (based on Jukes-Cantor-corrected distance). Bootstrap value was calculated based on 1,000 replications.

Quantitative PCR was performed using an iCycler IQ System (Bio-Rad, Hercules, USA), three replicates for each sample. For AOA amoA gene, qPCR was conducted in a total volume of 30 μL containing 15 μL of FailSafe™ PCR Premix F, 5 μL of DNA template, 0.2 μM of each primer, 1.5 U of AmpliTaq® DNA polymerase, 15 mM MgCl 2 , and 0.5× SYBR® Green I (Invitrogen, Eugene, USA). For AOB 16S rRNA and amoA genes, PCR was performed in a total volume of 30 μL containing 15 μL of iQ™ SYBR® Green Voytek and Ward 1995 CTO654r

CTAGCYTTGTAGTT TCAAACGC Kowalchuck et al. 1997 The concentration of DNA was determined by NanoDrop® Spectrophotometer ND-1000 (Thermo Fisher Scientific, USA). PCR was performed using FailSafe™ PCR Premix F (EPICENTRE Biotechnologies, Madison, USA) following the cycling program for each primer set listed here a Primers used for T-RFLP analysis were labeled by Hex Super Mix (Bio-Rad, Hercules, USA), 5 μL of DNA template, and 0.2 μM of each primer. The qPCR thermocycling steps were set as follows: 95°C for 4 min and 45 cycles at 95°C for 45 s, 55°C for 45 s, and 72°C for 45 s. Cycling was completed by a final elongation step at 72°C for 10 min. The negative control containing no DNA was subjected to the same procedure to exclude or detect any possible contamination. After qPCR assay, the specificity of amplification was verified by melting curve analysis and checking with agarose gel electrophoresis.

For T-RFLP, PCR was performed by using fluorescently labeled primer amoA-2R-Hex for AOB. PCR conditions were as the same as those described above for clone library construction. AluI was selected to digest PCR products of AOB amoA gene. The digestion mixture, containing 17 μL of purified PCR product (about 200 ng DNA), 2 μL of buffer, and 1 μL of restriction enzyme (10 U/μL), was incubated at 37°C for 5 h. The fluorescently labeled T-RFs were run though an ABI 3730xl capillary sequencers in the GeneScan mode. T-RFLP data was analyzed using Gene-Marker V1.6 (SoftGenetics LLC, Pennsylvania, USA). Because of the detection range of internal marker GS500, T-RFs smaller than 50 bp and larger than 500 bp were excluded from further analysis. The peaks were first selected by default parameters setting of the software GeneMarker with the threshold cutoff of peak intensity of 100. The relative abundance of each T-RF was then determined by calculating the ratio between the area of each peak and the total area of all peaks in one sample. The peaks with relative abundance less than 1% were neglected in this study.

The sequences reported in this study were deposited in GenBank under accession numbers GU988755-GU988900.

Physical and chemical properties of the sediment samples Diversity of AOA amoA gene sequences

Archaeal amoA clone libraries were generated for the sediment sample S16 by using three primer sets (Table 1) . A total of 145 clones were sequenced and grouped into 36 OTUs. In details, three clone libraries (A1, A2, and A3) recovered 19, 16, and 15 OTUs using primer sets I, II, and III, respectively (Table S2 ). Six OTUs (Fig. 2a, OTU1-6) were shared by the three clone libraries and accounted for 57%, 64%, and 70% of total clones in the three libraries, respectively. Two OTUs (Fig. 2a, OTU7 and OTU8) appeared in two clone libraries. The other 28 OTUs were only recovered in one of the three libraries.

For AOB 16S rRNA gene clone library, as shown in Fig. S1 , 26 clones were sequenced and grouped into seven OTUs, i.e., one OTU in the genus Nitrosospira and six OTUs in the genus Nitrosomonas. Similarly, 26 AOB amoA gene sequences were generated from another clone library and grouped into eight OTUs (Fig. 2b) : three OTUs in the genus Nitrosospira and five OTUs in the genus Nitrosomonas. No γ-AOB sequences were found in the amoA gene clone library indicating that the major AOB species in the samples were β-AOB.

To quantify the AOA abundance, three primer sets of AOA amoA gene were used for qPCR. The results showed that qPCR could be successfully conducted by using the primer sets I and III. The qPCR results showed that the AOA abundance was significantly different by using two primer sets for the samples S9, S16, S19, and S21. Generally, as shown in Fig. 3 , the AOA abundance calculated from qPCR using primer set I was about 10 times higher than that using primer set III. Only for the sample S4, two primer sets resulted in the same level of AOA abundance. AOA amoA gene copy numbers based on qPCR by using primer set I ranged from 6.8×10 3 to 2.2×10 4 copies per nanogram DNA and from 9.6×10 6 to 5.1×10 7 copies per gram of sediment (wet weight). Over all, there was no significantly difference among five sediment samples considering the variation of qPCR itself. AOB abundance was quantified by qPCR using both AOB 16S rRNA and amoA genes. As shown in Fig. 3 , AOB 16S rRNA gene ranged from 15 to 126 copies per nanogram DNA and 2.7×10 4 to 1.4×10 5 copies per gram of sediment. Similarly, AOB amoA gene resulted in 67 to 333 copies per nanogram DNA and 9.5×10 4 to 6.2×10 5 copies per gram of sediment. If assuming that there is one copy of 16S rRNA gene per cell for both Nitrospira and Nitrosomonas (Dionisi et al. 2002) , the mixed AOB community in this study should contain one copy 16S rRNA gene per cell, according to the clone library results. Assuming that there are 2.5 copies of amoA gene per AOB cell (Norton et al. 2002) , thus the copy number of AOB amoA gene in one sample should be 2.5 times that of AOB 16S rRNA gene. For five sites, the average copy numbers per gram of sediment were 2.9×10 5 and 1.2×10 5 for AOB and 16S rRNA gene, respectively. The results of cells of Fig. 2 Phylogenetic trees of AOA (a) and AOB (b). The neighbor-joining trees were constructed using ARB based on Jukes-Cantor-corrected DNA distances, showing the phylogenetic affiliation of amoA gene sequences from the Pearl River Estuary sediment S16 (printed in bold) and from other environmental samples or pure strain retrieved from the GeneBank. In (a), except for the singletons, one clone was selected from each OTU for phylogenetic analysis and the total number of clones grouped to that OTU was listed after the clone name. The bootstraps (based on 1,000 replications) larger than 60% were indicated by the black dots AOB based on amoA gene copy number were generally consistent with those based on 16S rRNA gene (Fig. 3) .

T-RFLP analysis clearly showed the dynamic changes of AOB community among five sampling sites. There were total of four T-RFs appearing in profiles. Based on the relative abundance of each one, T-RFs of 95 and 264 bp represented the dominant AOB and other two T-RFs of 143 and 486 bp only occupied very low abundance, less than 6% individually. All the four T-RFs could be correlated to the clones of AOB amoA gene clone library, as shown in Table S3 . The T-RF of 95 bp corresponded to AOB in the genus Nitrosomonas, while T-RF of 264 bp was mainly derived from AOB in the genus Nitrosospira, although several Nitrosomonas clones also resulted in this T-RF. The clone library was constructed using the sediment sample S16, and here, it was assumed that the same T-RFs in other four samples represented similar AOB species as in the sample S16. Obviously, as shown in Fig. 4 , the species compositions of AOB among five samples were dramatically different. Gradient changes showed along the sample collection sites. From the sites S4 to S21, the relative abundance of Nitrosomonas, indicated by T-RF of 95 bp, gradually decreased from 58% to the level under detection by T-RFLP. In contrast, the relative abundance of AOB in the genus Nitrosospira, indicated by T-RF of 264 bp, significantly increased to double. It was also possible that the same T-RF in different samples came from different AOB species, or the T-RF of 264 bp in other samples was composed by AOB in two genera in different ratios, thus, in fact the diversity changes of AOB community among the sampling sites should be even more complicated than that observed here.

Phylogenetic analysis indicated that most of AOA amoA gene sequences in this study were closely related to uncultured sediment AOA found at other places, except Fig. 2 (continued) three small clusters (clusters I, II, and III). Clusters I and III were grouped with sequences from mixed environmental samples, including soil (AB529723, GQ481040), hot spring (GQ226087, GQ226097), sediment (FJ951725, EU885576), active sludge of bioreactor (EU860263), and the identified pure AOA species Nitrosopumilus maritimus (EU239959). In cluster II, the two sequences were affiliated to those obtained from freshwater samples. Phylogenetic tree also showed that the most closely related sequences of five shared OTUs, except OTU6, were all from the San Francisco Bay estuary sediments (Mosier and Francis 2008) . In total, these five OTUs reached 57.2% abundance (83 out of 145 clones) in the sediment of this study. These results indicated that the AOA represented by OTU1 to OTU5 in this study, and their relatives from other estuaries were possibly major AOA population in the estuary sediments, although so far, these AOA are still quite difficult to be cultured and identified. Among them, OTU1 represented the most dominant AOA group in this study, occupying 25.5% relative abundance (37 of 145 clones).

For AOB, generally, the results from two clone libraries were comparable and consistent. It was also consistent with previous reports that AOB in two genera Nitrosomonas and Nitrosospira were the dominant AOB in estuarine sediments (Mosier and Francis 2008) . Compared to AOA amoA gene, AOB showed lower diversity for the same sample in this study.

The clone library and agarose gel electrophoresis results indicated that both primer sets did not generated unspecific PCR products; thus, it seems that the real-time PCR results from primer set I should be more accurate than that from primer set III. The difference might be due to the different binding efficiencies of the two forward primers to DNA templates as the reverse primers are the same. So far, quantification of amoA gene abundance using qPCR is still in a relatively early stage of development, and more reliable measurement methods must be applied to confirm the previous findings.

Compared with the results of Mosier and Francis (2008) , the AOA abundance in this study was at the high end of that in San Francisco Bay estuary sediments (from 1.4×10 4 to 3.9×10 7 copies per gram of sediment). However, in their study, high AOA abundance only occurred in low-salinity regions (0.2-9 psu). The salinities of five sampling sites in this study increased from 5.6 to 24.6 psu, while the AOA abundance kept at high level along the salinity gradient. It might due to the high ammonia concentrations at these sites. These results also provide the possible reason for the previous report of high nitrification rates in summer at the Pearl River Estuary (Dai et al. 2008) .

It was reported that AOB abundance increased with salinity (Mosier and Francis 2008) . However, in this study, the results did not show the same relationship. Overall, AOB amoA gene copy number was 100 times lower than that of AOA, suggesting that the dominant ammonia-oxidizing microorganisms in the sediment of the Pearl River Estuary were AOA.

Comparing the environmental parameters of five sites, the salinity ranged from 5.6 to 24.6 psu, which possibly was the key factor shaping the AOB community composition, although total abundance of AOB did not show significant fluctuation. The results indicated that AOB in the genus Fig. 4 Relative abundance of AOB T-RFs in the Pearl River Estuary sediments, determined by the corresponding normalized T-RFLP peak area. PCR was performed by using fluorescently labeled primer amoA-2R-Hex for AOB. PCR products were digested using AluI. The relative abundance of the selected T-RF peak (intensity above 100) was determined by calculating the ratio of its area and the total area of all selected peaks in one sample Nitrosomonas in this study might prefer low salinity, while in high salinity, the species in Nitrosospira were the dominant AOB. It was consistent with the previous findings in the Chesapeake Bay (Francis et al. 2003) , Plum Island Sound , Ythan (Scotland) estuaries (Freitag et al. 2006) , and the San Francisco Bay estuary (Mosier and Francis 2008) , where similar distribution patterns of AOB in two genera were reported. T-RFLP analysis was also applied to investigate the dynamic shift of AOA community. However, it was found that there were only three major peaks in the profiles (Fig. S3 ). The number of T-RFs was largely smaller than that of OTUs in AOA clone library. It means the resolution of the method was too low. The sequences of AOA amoA gene obtained in this study were aligned and virtually digested by over twenty restricted enzymes. But, there was not any enzyme which separated these sequences very well. It was concluded that T-RFLP was not a suitable method to analyze the difference of AOA community in this study.

Selection of primers for AOA amoA gene Junier et al. (2010) summarized all the published primer sets for AOA amoA gene amplification. The primer selection is a key step for microbial community study based on PCR. However, the specificity and sensitivity of these primers were not compared so far.

In this study, three primer sets for AOA amoA gene were selected for comparison. The results of clone libraries showed that both shared and distinct clones were recovered by using different primer sets. Generally, by using any one of three primer sets, the dominant AOA species in this sediment sample, represented by the shared six OTUs, could be recovered from the clone library. However, those distinct OTUs might be recovered only by specific primers. Specially, the sequences in cluster I were only amplified by primer sets II and III. Of course, the difference among three clone libraries was also possibly due to the low abundance of those OTUs and insufficient coverage of the clone library.

Comparison of Chao 1 among three clone libraries showed that the primer set III resulted in the highest value and significantly higher than the primer set II, indicating that the primer set III might match more AOA amoA gene sequences than other two primer sets. Theoretically, using this primer set III may recover the highest diversity of AOA amoA genes in this sediment sample, although we got the highest number of OTU from clone library by using primer set I. The significant difference was also shown when primer sets were applied in qPCR. Our previous AOA study in activated sludge already reported that two primer sets showed impact in clone library results (Zhang et al. 2009 ). These results suggested that two or more primer sets should be tried when starting the PCR for a new environment sample. It was also indicated that the primer selection was critical for clone library construction, qPCR, and other methods based on PCR.

Diversity of ammonia-oxidizing archaea and bacteria in the sediments of a hypernutrified subtropical estuary: Bahiadel Tobari

Loss of diversity of ammonia-oxidizing bacteria correlates within creasing salinity in an estuary system

Nitrification and inorganic nitrogen distribution in a large perturbed river/estuarine system: the Pearl River Estuary

Cultivation of a thermophilic ammonia oxidizing archaeon synthesizing crenarchaeol

Quantification of Nitrosomonas oligotropha-like ammonia-oxidizing bacteria and Nitrospira spp. from full-scale wastewater treatment plants by competitive PCR

Diversity of ammonia monooxygenase (amoA) genes across environmental gradients in Chesapeake Bay sediments

Ubiquity and diversity of ammonia-oxidizing archaea in water columns and sediments of the ocean

Changes in the community structure and activity of betaproteobacterial ammonia-oxidizing sediment bacteria along a freshwater-marine gradient

The population frequencies of species and the estimation of population parameters

Global patterns of marine nitrogen fixation and denitrification

Pathways of carbon assimilation and ammonia oxidation suggested by environmental genomic analyses of marine Crenarchaeota

The characteristics of nutrients and eutrophication in the Pearl River estuary, South China

Phylogenetic and functional marker genes to study ammonia-oxidizing microorganisms (AOM) in the environment

A gene encoding a membrane protein exists upstream of the amoA/amoB genes in ammonia-oxidizing bacteria; a third member of the amo operon

Isolation of an autotrophic ammonia-oxidizing marine archaeon

Analysis of ammonia-oxidizing bacteria of the β subdivision of the class proteobacteria in costal sand dunes by denaturing gradient gel electrophoresis and sequencing of PCR-amplified 16 S ribosomal DNA fragments

ARB: a software environment for sequence data

Relative abundance and diversity of ammonia-oxidizing archaea and bacteria in San Francisco Bay estuary

Diversity of ammonia monooxygenase operon in autotrophic ammonia-oxidizing bacteria

Phylogeny of all recognized species of ammonia oxidizers based on comparative 16 S rRNA and amoA sequence analysis: implications for molecular diversity surveys

The ammonia monooxygenase structural gene amoA as a functional marker: molecular fine-scale analysis of natural ammonia-oxidizing populations

Diversity and spatio-temporal distribution of ammonia-oxidizing Archaea and Bacteria in sediments of the Westerschelde estuary

Denitrification in fresh water and coastal marine ecosystems: ecological and geochemical significance

Environmental genome shotgun sequencing of the Sargasso Sea

Detection of ammonium-oxidizing bacteria of the beta-subclass of the class Proteobacteria in aquatic samples with the PCR

Prokaryotes: the unseen majority

Archaeal nitrification in the ocean

Occurrence of ammonia-oxidizing archaea (AOA) inactivated sludges of a laboratory scale reactor and two wastewater treatment plants

The authors wish to thank the Hong Kong General Research Fund (HKU7197/08E) for the financial support of this study, and Lin Ye wish to thank HKU for the postgraduate studentship. This work was partially supported by the CAS/SAFEA International Partnership Program for Creative Research Teams to PY Qian.Open Access This article is distributed under the terms of the Creative Commons Attribution Noncommercial License which permits any noncommercial use, distribution, and reproduction in any medium, provided the original author(s) and source are credited.
