Previous Article | Next Article ![]()
Applied and Environmental Microbiology, May 2009, p. 3127-3136, Vol. 75, No. 10
0099-2240/09/$08.00+0 doi:10.1128/AEM.02806-08
Copyright © 2009, American Society for Microbiology. All Rights Reserved.

Department of Biological Sciences, Microbiology, Aarhus University, Ny Munkegade, Building 1540, DK-8000 Aarhus C, Denmark,1 Department of Chemistry and Water Technology, Danish Technological Institute, Kongsvang Allé 29, DK-8000 Aarhus C, Denmark2
Received 10 December 2008/ Accepted 16 March 2009
|
|
|---|
|
|
|---|
The first key step of coupled nitrification-denitrification, the oxidation of ammonia to nitrite, is catalyzed by two groups of prokaryotes—the ammonia-oxidizing bacteria (AOB) (24) and the only recently recognized ammonia-oxidizing archaea (AOA) (22). For both groups, the gene encoding the alpha-subunit of ammonia monooxygenase (amoA) has been widely used as a functional marker to analyze their community compositions (15, 25); recent studies demonstrated the ubiquity of AOA and their predominance over AOB in a broad range of environments (32, 38). AOA, but not AOB, were also strongly enriched in the rhizosphere of the freshwater macrophyte Littorella uniflora in a mesotrophic Danish lake, suggesting that AOA were primarily responsible for increased rates of nitrification in the rhizosphere of this plant species (19). Moreover, ammonia oxidizer communities differed between rhizosphere and unvegetated sediment, indicating a plant-specific effect on AOA and AOB community composition. The objectives of this study were therefore to test whether (i) AOA generally predominate over AOB in freshwater sediments and especially in macrophyte rhizospheres and (ii) macrophytes have species-specific effects on abundance and community composition of AOA and AOB in rhizosphere sediments and on root surfaces.
To address these questions, two shallow heathland pools and five lakes in Denmark and Germany, ranging from low-pH and dystrophic sites to neutral-pH and oligotrophic and mesotrophic sites, were chosen, and three macrophyte species—Littorella uniflora, Juncus bulbosus L. (bulbous rush), and Myriophyllum alterniflorum DC. (alternate water milfoil)—were selected as model systems. These plant species differ in nitrogen nutrition, extent of radial oxygen loss, and lifestyle, presumably resulting in differential, plant species-specific effects on rhizosphere- and root-associated AOA and AOB communities. L. uniflora prefers nitrate as the nitrogen source, while J. bulbosus prefers ammonium (41, 45); oxygen release is high to moderate from the roots of L. uniflora and J. bulbosus (9, 12) but is minor from the roots of M. alterniflorum (M. Herrmann, P. Stief, and A. Schramm, unpublished results); L. uniflora and J. bulbosus remain photosynthetically active throughout the year, while only the below-ground parts of M. alterniflorum are retained during winter.
Rhizosphere sediments and roots from each plant species were sampled from three different sites per species, and unvegetated sediment was obtained from all seven sites. The comparison of samples from these different sites and compartments (rhizosphere, root surface, unvegetated sediment) allowed an evaluation of the importance of plant species relative to that of environmental conditions related to lake trophic status and pH on ammonia oxidizer communities.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Sampling sites, sample types, and chemical composition of sediment and sediment pore watera
|
Sediment chemical analysis.
Pore water samples were extracted from unvegetated sediment and rhizosphere sediment by the centrifugation of about 75 g of sediment at 1,000 rpm for 5 min in combination with filtration through glass microfiber filters (Whatman) and were analyzed for pH and NO2– plus NO3– (6), NO2– (16), and NH4+ (5) concentrations. Concentrations of free NH3 were calculated based on pH and NH4+ concentrations (pKa of NH4+/NH3 is 9.4 at 20°C). For analyses of total organic C and total N contents, sediment samples were dried at 110°C for 76 h, homogenized using a mortar, treated with H2SO3 for 64 h, and dried again at 70°C for 24 h. C/N analyses were performed on an NA 1500 nitrogen analyzer (Carlo Erba Strumentazione).
DNA extraction, PCR amplification, cloning, and sequence analysis.
DNA was extracted in triplicate from 200 mg sediment by combining enzymatic and chemical cell lysis with the FastDNA Spin kit for soil (Qbiogene Inc.) (13). For the analysis of microbial communities on root surfaces, DNA was extracted from 100 mg plant roots following the same protocol. Archaeal and bacterial amoA genes were amplified using the HotStarTaq master mix kit (Qiagen) with published protocols and primer sets Arch-AmoAF/Arch-AmoAR (14) and AmoA-1F/AmoA-2R-TC (34). Bands of the correct size were excised and purified (GenElute gel extraction kit; Sigma). Purified PCR products from triplicate DNA extracts of the triplicate sediment cores were pooled for subsequent cloning (pGEM-T cloning kit; Promega). For samples of unvegetated sediment, one clone library was constructed per site for each gene. For bacterial amoA, one combined rhizosphere sediment/root surface clone library was constructed for each plant species and site, and denaturing gradient gel electrophoresis (DGGE) (34) was used to differentiate rhizosphere and root surface populations by assigning clones to specific bands in the environmental DGGE patterns. Bands not represented by cloned inserts were excised and reamplified for the direct sequencing of purified PCR products (34). For archaeal amoA, separate clone libraries were constructed for rhizosphere sediment and root samples of each plant species and each site. In total, 16 bacterial and 25 archaeal amoA clone libraries were constructed, plasmids were extracted using the GenElute plasmid preparation kit (Sigma), and inserts were sequenced by Macrogen (South Korea). Alignments, translations to amino acids, and phylogenetic analyses were done using the ARB package (30). Phylogenetic trees of bacterial and archaeal AmoA were constructed as consensus trees of neighbor-joining and parsimony analysis (1,000 replicates). Maximum likelihood analysis was also performed but with a reduced number of sequences. Operational taxonomic units (OTUs) were defined based on a 97% identity cutoff on the amino acid level using DOTUR (43), and Good's coverage was calculated for each clone library (47).
qPCR.
Relative abundances of AOA and AOB were determined by quantitative PCR (qPCR) as described previously (19). In brief, copy numbers of 16S rRNA genes and bacterial or archaeal amoA genes were determined in triplicate for each sediment sample, using universal primers 907F/1492R (27) and amoA primers Arch-AmoAF/Arch-AmoAR (14) and AmoA-1F/AmoA-2R-TC (34) with brilliant SYBR green qPCR master mix (Stratagene) on an Mx3005P instrument (Stratagene). qPCR reactions contained 25 ng environmental DNA or 25 ng salmon sperm DNA in standard reactions. Standard curves were prepared from serial dilutions of Bacteroides fragilis ATCC 25285 genomic DNA for 16S rRNA genes and serial dilutions of plasmids containing an environmental archaeal amoA sequence (clone AS_AOA_2 [GenBank accession number EU309881]) or the amoA sequence of Nitrosospira multiformis ATCC 25196, respectively. Standard curves were linear from 107 to 102 16S rRNA gene copies and from 5 x 108 to 5 x 101 amoA gene copies; i.e., the amoA genes could not be reliably quantified when there were <50 gene copies (the limit of quantification). Detection limits were five gene copies for the functional genes. Efficiencies of the qPCR reactions were 88 to 97% for amoA and 75 to 83% for 16S rRNA genes. The specificities of PCR products were confirmed by melting curve analysis and agarose gel electrophoresis. Environmental DNA templates were spiked with standard plasmid DNA to test for the presence of PCR inhibitors. The abundances of AOA and AOB were expressed relative to the number of prokaryotic cells as estimated from 16S rRNA gene copy numbers, assuming 2.5 and one amoA gene copies per AOB and AOA, respectively (29), one 16S rRNA gene copy per AOA, and 3.6 copies of the 16S rRNA gene per average prokaryotic cell (21). Relative abundances were preferred over absolute values for this study, as they are less sensitive to differences in DNA extraction efficiencies; these can be expected when working with sample types that differ substantially in their contents of humic substances.
Statistical analysis.
To compare the clone libraries between different sites, macrophyte species, and compartments, Jaccard similarity values were calculated (44). These values consider how many OTUs are shared between clone libraries relative to the total number of OTUs in these clone libraries. Similarity matrices were turned into dendrograms to visualize the clustering of AOA and AOB communities based on their similarities in OTU composition, utilizing the programs "Clustering Calculator" (http://www2.biology.ualberta.ca/jbrzusto/cluster.php) and "Interactive Tree of Life" (http://itol.embl.de). Differences between sites in the chemical parameters and relative abundances of AOA were analyzed by one-way analysis of variance (ANOVA) and Tukey's least significant difference test (P = 0.05; XLSTAT). To test if differences in compartment, plant species, or lake trophic status had a general effect on the number of detected OTUs for AOA and AOB, samples obtained from the different sites were grouped according to compartment, plant species, or lake trophic status, and groups were compared using one-way ANOVA and Tukey's least significant difference test (P = 0.05; XLSTAT). Spearman's rank correlation coefficients were calculated to analyze if relative abundances of AOA or the numbers of detected OTUs of AOA and AOB were correlated to environmental parameters (XLSTAT). Relative abundances of AOB were excluded from correlation analysis, since the copy numbers of bacterial amoA genes usually remained below the quantification or detection limit of the qPCR, and the relative abundances given represent only maximum estimates based on the quantification limit or detection limit.
Nucleotide sequence accession numbers.
Nonredundant sequences of archaeal and bacterial amoA genes obtained in this study have been deposited in the GenBank database under accession numbers EU667645 to EU667773 (for bacterial amoA) and EU667774 to EU667992 (for archaeal amoA). Most of the sequences obtained from site A have been published previously (accession numbers EU309859 to EU309918) (19).
|
|
|---|
Community composition of AOA and AOB.
Between 8 and 24 clones were sequenced per archaeal amoA clone library to achieve coverages of 82 to 100%, resulting in a total of 306 archaeal amoA sequences. Based on 97% translated amino acid identity, they were assigned to 13 OTUs, most of which included a large number of sequences from different sites; only three OTUs were represented by single sequences. The number of detected OTUs per clone library ranged from 1 to 5. For bacterial amoA, between 16 and 90 clones per library were screened by DGGE prior to sequencing, and sequences were obtained from 1 to 33 clones per library. Additionally, 1 to 10 sequences per library were obtained directly from environmental DGGE bands not represented by clones. Including sequences derived from clones and DGGE bands, library coverages were 71 to 100%. The resulting 191 bacterial AmoA sequences were assigned to 15 OTUs, with one to six OTUs per clone library.
The number of detected OTUs of both AOA and AOB was highest in the rhizospheres of L. uniflora and M. alterniflorum at mesotrophic Lake Hampen (site A). The acidic heathland pool (site H) showed the lowest diversity, with only a single OTU for AOA AmoA and no bacterial amoA genes detected (Fig. 1). One-way ANOVA indicated that neither plant species nor compartment had a significant site-independent influence on the number of detected OTUs of AOA or AOB, while mesotrophic sites could significantly be distinguished from dystrophic and oligotrophic sites (P = 0.05). For archaeal and bacterial AmoAs, four and six OTUs, respectively, occurred exclusively in the rhizosphere or on root surfaces, while two and four OTUs, respectively, were restricted to unvegetated sediment. Out of these compartment-specific OTUs, only a single bacterial AmoA OTU occurred at more than one site.
![]() View larger version (33K): [in a new window] |
FIG. 1. Number of detected OTUs of AOA (closed circles) and AOB (open circles) and relative abundances of AOA (gray columns) and AOB (white columns) in unvegetated sediment and rhizosphere sediment. For definitions of the abbreviations for samples, see Table 1. Relative abundance values are the means ± the standard deviations of nine qPCR reactions with triplicate DNA extracts from three environmental replicates. Columns represent maximum estimates, where amoA gene copy numbers remained below the quantification limit of 50 copies (<) or the detection limit of 5 copies (<<) per reaction. Letters above columns indicate significant differences between samples.
|
![]() View larger version (24K): [in a new window] |
FIG. 2. Similarities of clone libraries of archaeal AmoA (A) and bacterial AmoA (B) among sites, plant species, and compartments, based on Jaccard values. The origins of sequences are indicated by closed circles (mesotrophic sites), open triangles (oligotrophic sites), and closed squares (dystrophic sites). For the definitions of abbreviations, see Table 1.
|
![]() ![]() View larger version (46K): [in a new window] |
FIG. 3. Phylogenetic affiliation of archaeal AmoA (A) and bacterial AmoA (B). Protein distance trees are displayed with bootstrap values (>50%) from parsimony analysis (1,000 replicates). Branchings not supported by both methods are drawn as multifurcations. The major lineages were also supported by maximum likelihood analysis with a reduced number of sequences. For all the OTUs obtained in this study, accession numbers and clone numbers for representative sequences from each clone library are shown in boldface type. For the definitions of abbreviations, see Table 1. The origins of sequences are indicated by closed circles (mesotrophic sites), open triangles (oligotrophic sites), and closed squares (dystrophic sites).
|
Phylogenetic affiliation of bacterial AmoA sequences.
Bacterial AmoA sequences affiliated with groups commonly found in freshwater habitats (24), i.e., the Nitrosospira, Nitrosospira marina/Nitrosomonas oligotropha, and Nitrosomonas europaea lineages (Fig. 3B). With few exceptions (GS_AOB_1_DGGE, CS_AOB_40, JGR_AOB_14), Nitrosomonas-like sequences were only obtained from mesotrophic sites (A, D, and F); only two sequences from the mesotrophic site A were related to Nitrosomonas europaea/Nitrosomonas eutropha, while all other Nitrosomonas-like sequences affiliated with the N. marina/N. oligotropha lineage. The majority of bacterial AmoA sequences was related to Nitrosospira-like sequences, particularly to a group of uncultured AOB retrieved from freshwater and rice paddy soils. Low-diversity AOB communities at oligotrophic sites, such as Lake Kalgård (site B) and Lake Grane Langsø (site C), consisted almost exclusively of Nitrosospira-related AOB.
Relative abundances of AOA and AOB.
AOA dominated numerically at most of the sites and were up to 3 orders of magnitude more abundant than were AOB. Unvegetated sediment of mesotrophic Lake Søby (site D) was the only sample in which no archaeal amoA genes were detected, while bacterial amoA genes were not detected in acidic heathland pool H (site H), in samples from the root surface of J. bulbosus and from the rhizosphere of L. uniflora from oligotrophic Lake Grane Langsø (C), or in samples from the root surface of M. alterniflorum from mesotrophic Lake Almind (F) (Fig. 1). For some of the other sites, bacterial amoA genes were not detected in all replicates of a certain sample type. Variation among replicate sediment cores was higher than variation among replicate DNA extracts for all three compartments, which indicated a high spatial heterogeneity of the distribution of the two groups of ammonia oxidizers.
With few exceptions, relative abundances of AOA were high (0.1 to 5.6% of all prokaryotes) in the dystrophic and oligotrophic sediments and in the rhizosphere/root surface system of L. uniflora and J. bulbosus. In contrast, archaeal amoA gene copy numbers were <50 per reaction in most mesotrophic sediments, including the rhizosphere/root surface system of M. alterniflorum (Fig. 1).
Bacterial amoA gene copy numbers were always below the limit of quantification (<50 copies per reaction) and, in many samples, even below the detection limit of the qPCR (five copies per reaction). However, bacterial amoA genes were still detectable in most samples by the standard PCR; therefore, 50 or 5 amoA gene copies per reaction were used to estimate the maximum relative abundance of AOB (Fig. 1). Based on these maximum estimates, relative abundances of AOB were highest at mesotrophic sites in association with both the rhizosphere and root surfaces of M. alterniflorum (maxima 0.005% and 0.055% of all prokaryotes, respectively) (Fig. 1).
AOA and AOB on root surfaces.
With a few exceptions, OTUs of AOA and AOB that occurred on the root surface of a given plant species at a given site could also be retrieved from the corresponding rhizosphere, with the number of detected OTUs on the root surface being equal to or lower than that detected in the rhizosphere. In particular, a decrease in the OTU numbers of root surface communities was observed for J. bulbosus at site C and for L. uniflora at site A (AOA and AOB) and C (AOA only). Relative abundances of AOA were, in general, lower on the root surfaces of L. uniflora and J. bulbosus than those in the respective rhizospheres, while there was no clear trend for M. alterniflorum or for the relative abundances of AOB (data not shown).
Relationship between abundance and community composition of AOA and AOB and environmental parameters.
Relative abundances of AOA were negatively correlated to ammonia availability and pH and were positively correlated to the sediment content of organic C and N, while there was no significant correlation with C/N ratios (Table 2). No such correlation analysis was performed for the relative abundances of AOB, since these data were only maximum estimates. For a range of 2 to 213 nmol liter–1 ammonia and pH 4.8 to 6.8, the number of detected OTUs of both AOA and AOB was strongly positively correlated to ammonia availability and pH. However, this correlation was no longer significant for AOA if samples from Lake Søby (site D) with the maximum ammonia availability and pH were included in the analysis (data not shown). The number of detected OTUs of AOA was negatively correlated to C/N ratios and sediment content of C and N (Table 2).
|
View this table: [in a new window] |
TABLE 2. Spearman's rank correlation coefficients of environmental parameters and relative abundances of AOA or numbers of detected OTUs of AOA and AOB in unvegetated sediment and rhizosphere sedimenta
|
|
|
|---|
Comparison of the compositions of AOA and AOB communities from seven lakes and three plant species in this study clearly showed that plant species-specific traits, such as nitrogen source utilization, extent of oxygen release, and lifestyle, did not select for specific ammonia oxidizers that would always be associated with the root system of a particular plant species. On the contrary, environmental conditions related to lake trophic status, such as ammonia availability, pH, and contents of humic acids, probably had a stronger effect on the conditions in the rhizosphere than did the plant species itself, resulting in a direct correspondence of especially the AOA communities to lake trophic status. AOB communities did not show any obvious clustering patterns. The lack of plant species-specific associations also indicates the absence of specific symbiotic interactions between AOA or AOB and the macrophyte species investigated.
Phylogenetic affiliation of AOA and AOB.
The affiliation of sequences obtained from freshwater environments with sequences derived from marine and estuarine habitats, wastewater treatment plants, and hot springs indicates that the phylogeny of archaeal AmoA is not primarily linked to general habitat types, as was originally proposed (14). In contrast, the results obtained in this study suggest that environmental factors related to lake trophic status, such as pH and ammonia availability, may play a critical role in the phylogenetic affiliation of freshwater AOA. While cluster A II harbored sequences retrieved from all lake types investigated in this study, affiliation with clusters A I and B was largely restricted to AmoA sequence types obtained from mesotrophic lakes. Moreover, sequences obtained from unvegetated sediment and J. bulbosus rhizosphere and root surfaces at acidic heathland pool H formed a single distinct OTU within cluster A II, suggesting a specific adaptation of the corresponding AOA phylotype to extreme site conditions of low pH, high C content of the sediment, high C/N ratios, and extremely low availability of ammonia. A clear effect of pH on AOA community composition and phylogenetic affiliation of the respective phylotypes in soils was recently reported by Nicol et al. (33), with similar patterns for both archaeal amoA and crenarchaeotal 16S rRNA genes.
Nitrosospira- and Nitrosomonas oligotropha-related AOB dominated most of the AOB communities found in this study. These organisms have previously been reported from oligotrophic freshwater environments and from neutral to acidic soils (23, 49). Affiliation with Nitrosomonas lineages other than N. oligotropha was restricted to sequences obtained from mesotrophic sites at a higher availability of ammonia, which agrees with the higher substrate requirements reported for the corresponding AOB lineages (23).
Regulating factors for AOA in freshwater environments.
AOA dominated the ammonia oxidizer communities associated with all investigated plant species from all sites, which confirms and extends previous observations of archaeal predominance in the rhizospheres of L. uniflora (19) and Oryza sativa (11). AOA can contribute significantly to ammonia oxidation in certain soils and estuarine environments (8, 33, 42); since increased nitrification rates (19, 36) and increased nitrate concentrations (found in this study) in the rhizospheres of L. uniflora and J. bulbosus coincided with an increased abundance of AOA but not AOB, it is likely that AOA also contribute to nitrification in freshwater environments, at least in macrophyte rhizospheres. If so, AOA are apparently adapted to low ammonia availability and low pH, as their relative abundance was negatively correlated with pH and ammonia availability, and AOA were absent from unvegetated sediment at site D, which has the highest pH and ammonia availability. These data that point to a less important role of AOA under mesotrophic to eutrophic conditions and at higher pHs agree well with the pH-dependent, contrasting patterns of bacterial and archaeal amoA transcript numbers reported for soils (33).
Rhizosphere-associated nitrification is stimulated by oxygen release from the roots of aquatic macrophytes (36, 40), and plant species-specific differences in the extent of oxygen release are therefore likely to influence the abundance of rhizosphere-associated ammonia oxidizers. In fact, the plant species with the highest oxygen release—L. uniflora and J. bulbosus (9, 12)—harbored the highest relative abundances of AOA, while the low oxygen availability in the rhizosphere of M. alterniflorum (Herrmann et al., unpublished) resulted in relative AOA abundances only slightly higher than those in unvegetated sediment. However, since oxygen release by the three plant species was not recorded when sediment samples for molecular analyses were taken, the observed trend cannot be directly proven by statistical means. Differences in the extent of oxygen release between plant species may also be responsible for differences in the distribution patterns of AOA within the rhizosphere/root system. Relative abundances of AOA on the root surface of L. uniflora that are lower than those in the rhizosphere may be due to the sensitivity of ammonia oxidizers to very high oxygen concentrations of up to 450 µM (12). Differences between compartments were less pronounced for J. bulbosus with its moderate oxygen release (9), and no clear trend was observed for M. alterniflorum with its low oxygen release.
The positive correlation of relative abundances of AOA to sediment content of organic C and N indicates that not only root oxygen release but also root exudation of organic carbon and nitrogen compounds may play a role for the association of AOA with the rhizosphere of freshwater macrophytes. Uptake of amino acids has been shown for marine Crenarchaeota (18, 37), and utilization of root exudates has been proposed for mesophilic Crenarchaeota in the rhizosphere or mycorrhizosphere of terrestrial plants (4, 46). The C and N contents of L. uniflora and J. bulbosus rhizospheres were increased compared with those in unvegetated sediment, indicating a substantial release of organic substances via root exudation (10, 48). In contrast, the effect of M. alterniflorum on sediment chemistry was much weaker, and earlier studies did not indicate a stimulating effect for M. alterniflorum on microbial processes in the sediment (20). The abundance pattern of AOA in line with the extent of root exudation may suggest mixotrophic or heterotrophic utilization of root exudates by freshwater AOA; such a metabolism would also provide an explanation for the high relative abundances of AOA in unvegetated sediment with high contents of organic C and N (sites C and H). With the data available, no statistical evidence can be provided for the suggested relationship between root exudation and abundance of rhizosphere-associated ammonia-oxidizing archaea. Here, further studies including a more detailed characterization of root exudation and the carbon flow in the rhizosphere, i.e., by employing techniques such as stable isotope probing, are needed to gain more insight into the potential utilization of root exudates by this group of microorganisms.
Conclusion.
The community composition of AOA in oligotrophic to mesotrophic freshwater sediments is strongly determined by factors related to lake trophic status, such as pH, ammonia availability, and sediment content of organic C and N; there was no indication for species-specific effects of freshwater macrophytes on the community composition of rhizosphere-associated AOA or AOB. Where present, AOA were more abundant than AOB; the macrophytes L. uniflora and J. bulbosus, especially, supported high numbers of AOA, possibly by releasing oxygen and organic root exudates into their rhizospheres.
We thank Britta Poulsen for excellent assistance with field and laboratory work; Bente Lomstein, Rikke Holm, and Peter Stief for help with the chemical analyses; and Kilian Stoecker for sharing his AOB amoA ARB database. Kasper Urup Kjeldsen is acknowledged for his help with clone library analysis. We thank the owners of Lake Kalgård and Lake Grane Langsø as well as the "Westfälisches Museum für Naturkunde" for permission to obtain sediment samples.
Published ahead of print on 20 March 2009. ![]()
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»