Previous Article | Next Article ![]()
Applied and Environmental Microbiology, February 2006, p. 1019-1026, Vol. 72, No. 2
0099-2240/06/$08.00+0 doi:10.1128/AEM.72.2.1019-1026.2006
Copyright © 2006, American Society for Microbiology. All Rights Reserved.
Department of Ecological Microbiology, University of Bayreuth, D-95440 Bayreuth, Germany
Received 15 September 2005/ Accepted 8 November 2005
|
|
|---|
|
|
|---|
Earthworms emit N2O (27, 31) and N2 (21), and these in vivo emissions are coincident with a very abundant denitrifying population in the earthworm gut (25, 27, 31). The earthworm gut is anoxic and has a high content of water, sugars, organic and amino acids, nitrite, and ammonium, thus making it an ideal habitat for denitrifiers (22). The denitrifier population of the earthworm gut has been evaluated by cultivation methods (25, 27), and denitrifiers that occur in soil also occur in the earthworm gut (25). Analyses of 16S rRNA genes in earthworm guts and casts indicate that the microbiota in the earthworm gut is largely food derived rather than endemic (17, 18, 43). Based on these data, denitrifiers in the earthworm gut are hypothesized to be ingested soil microorganisms that are activated by the special microenvironment of the gut (22, 25). This hypothesis implies that denitrifiers in the earthworm gut are members of the soil microbial community and not endemic to the earthworm gut. If this hypothesis is valid, the species composition of denitrifiers in the earthworm gut should be equivalent to that of the soil. To test this hypothesis, the denitrifiers of soil and the earthworm gut were comparatively assessed by sequence analysis of nosZ.
|
|
|---|
Denitrifiers from earthworm gut.
Dechloromonas denitrificans (DSM 15892T), Flavobacterium denitrificans (DSM 15936T), and Pseudomonas sp. strain ED3 (GenBank accession number AJ318919) were cultured as previously described (20, 25) (The superscript T indicates a type strain; the 16S rRNA gene sequence of strain ED3 is 99.9% identical to that of Pseudomonas fluorescens [ATCC 17428]).
Extraction of nucleic acids.
DNA was isolated from soil, gut contents, and denitrifiers using the FastDNA SPIN kit (BIO 101, Carlsbad, CA) according to the manufacturer's protocol, with the following modifications to minimize shearing of large DNA fragments: bead beating was omitted, and samples (each, 1 g [fresh weight]) were suspended in 500 µl lysis buffer (150 mM NaCl, 100 mM EDTA, 10 mg lysozyme ml1, pH 8.0) and incubated for 30 min at 37°C in a thermomixer (model 5436; Eppendorf, Hamburg, Germany). Samples were supplemented with 60 µg of proteinase K, incubated 30 min at 37°C, and then subjected to three freeze-thaw cycles (80°C for 10 min and 60°C for 1 min). Chemical lysis yielded insufficient amounts of DNA from the soil of site HW; thus, bead beating (BIO 101 protocol) was used for soil from this site.
Primers and PCR conditions.
Fragments of nosZ genes were amplified using primers nosZ661F (CGG CTG GGG GCT GAC CAA) and nosZ1773R (ATR TCG ATC ARC TGB TCG TT) (38). Primers PsNosZ175F (TTC ATC AAC GAC AAG GCC) and PsNosZ1144R (CGG TGG GCA GGA AGC GGT) were designed with ARB software (http://www.arb-home.de) (46) for the selective amplification of Pseudomonas-related nosZ fragments.
Primer specificity was evaluated in silico using GenBank's BLAST tool (1). PCR conditions were tested with genomic DNA from Pseudomonas aeruginosa (DSM 6195T). Each amplification reaction was performed in a total volume of 50 µl and contained 50 mM KCl, 10 mM Tris-HCl, 1.5 mM MgCl2, 0.3 mg of bovine serum albumin ml1, 10 µl of TaqM, 200 µM nucleotides, 20 pmol of each primer, 1 µl of template DNA, and 1 U of Taq polymerase (Eppendorf, Hamburg, Germany). Touchdown PCR was performed with a T-Gradient cycler (Biometra, Göttingen, Germany) with 5 min of denaturation at 94°C, followed by 8 (10) cycles, each consisting of 94°C for 1 min, primer annealing at 60 to 56°C (63 to 58°C) for 1 min, and extension at 72°C for 3 min; this was followed by 27 cycles with annealing at 56°C (58°C) (numbers of cycles and temperatures given in parentheses refer to amplifications performed with the primer pair PsNosZ175F-PsNosZ1144R). The final PCR extension step was at 72°C for 5 min.
Cloning, screening, and sequencing.
PCR products of denitrifying earthworm gut isolates were purified by cutting out the band of expected size (1.1 kb) from a low-melting-point agarose gel (2%) with a subsequent extraction using the MiniElute gel extraction kit (QIAGEN, Hilden, Germany). PCR products of all samples were ligated into pGEM-T (Promega, Mannheim, Germany) and transformed into Escherichia coli JM109 (Promega) according to the manufacturer's protocol. For screening of the nosZ libraries, nosZ fragments were directly amplified from 1 µl of resuspended clones and electrophoresed on agarose gels (1%); those of the expected size (1.1 kb) were analyzed. The gut and soil clone libraries were screened by restriction fragment length polymorphism (RFLP); nosZ fragments were digested simultaneously with 2 U of HaeIII and RsaI (Promega) for 2 h at 37°C prior to separation on an 8% polyacrylamide gel (DCode; Bio-Rad, Richmond, CA) for 4 h at 120 V. One to 30 representative clones of each RFLP pattern per library were selected for sequencing. For the clone libraries from the denitrifying earthworm gut isolates, two clones of each library were selected for sequencing without prior RFLP analysis. The Wizard Plus Minipreps DNA Purification system (Promega, Mannheim, Germany) was used for the preparation of plasmids, and DNA was sequenced commercially at MWG Biotech (Ebersberg, Germany).
nosZ sequence analysis.
Sequence analysis was performed with the ARB software package (http://www.arb-home.de) (46). nosZ sequences and publicly available sequences overlapping in the same region were translated in silico and prealigned with the ClustalW algorithm, and the alignment was refined manually. Regions of primer binding sites were excluded from further sequence analysis. Neighbor-joining, parsimony, and maximum-likelihood algorithms were applied to sequences that were longer than 350 amino acids; a consensus tree was drawn using consistent branchings from two to three methods, while inconsistent branches were drawn as multifurcations (30). Partial sequences (<350 amino acids) were added to the consensus tree by maximum-parsimony analysis without changing the tree topology. Unless otherwise indicated, all sequence analyses were based on the in silico-translated nosZ sequences.
Analysis of nosZ genotype diversity.
Percent accepted mutations matrix (PAM)-corrected distance matrices were generated from aligned amino acid sequences with ARB (http://www.arb-home.de) (46). Partial nucleic acid sequences that shared >98% sequence similarity and sequences that had the same RFLP patterns with fully sequenced nosZ fragments were thereby represented by the fully sequenced nosZ fragments. DOTUR (http://www.plantpath.wisc.edu/fac/joh/dotur.html) (39) was applied to define genotypes, based on the amino acid sequence dissimilarity of proteins derived from translated nosZ fragments by the furthest-neighbor method. Analyses were performed with four arbitrary, predefined genotype definitions, i.e., sequences that had a maximal dissimilarity value (D) of
0.49, 3.49, 7.49, or 15.49% were combined into single nosZ genotypes. Additional information on DOTUR is available at the publisher's website (http://www.plantpath.wisc.edu/fac/joh/DOTUR/DOTURManual.pdf).
Coverage (C) (19), which is the number of the detected genotypes relative to their expected total number in a gene library, was calculated by S-Libshuff analysis as C = 1 yx1 (40), where y is the number of genotypes that occurred only once, and x is the number of clones screened. The diversity of genotypes represented in soil and gut gene libraries was analyzed by rarefaction analysis (23). Genotype richness (S), i.e., the number of different genotypes present in a gene library, was estimated by an abundance-based coverage estimator (8, 9) and Bootstrap (44), bias-corrected Chao1 (7, 9), and interpolated Jackknife (6) estimators. All calculations were done with DOTUR (39).
S-Libshuff analyses (40, 42) were performed to determine the significance of differences between gene libraries of soil and those of the earthworm gut independent of an arbitrary, predefined genotype definition. For comparison of two gene libraries, Libshuff analyses plot coverage against the evolutionary distance of sequences for one gene library (homologous coverage, CSoil or CGut = 1 ySoilxSoil1 or 1 yGutxGut1, where the indices depict the gene library; see above for definitions of x and y) and for sequences of one gene library compared to the other (heterologous coverage, CSoil/Gut or CGut/Soil = 1 ySoilxGut1 or 1 yGutxSoil1; see above). The coverage curves describe how well the sampling represents the entire library at various levels of evolutionary distance. Libshuff analysis calculates differences between homologous and heterologous coverage curves by a Cramer-von Mises-type statistic and compares them by a Monte Carlo test (42).
Accession numbers.
All nosZ sequences have been deposited in the EMBL nucleotide sequences database under accession numbers AJ703894 to AJ704214 and AJ704930 to AJ704933.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Summary of data from amino acid analyses of in silico-translated nosZ fragments from denitrifier communities in agricultural soils and earthworm gut contents
|
0.49%, respectively. The mean evolutionary distances among sequences from the field (B), garden (H), and meadow (HW) sites were 14.4% ± 9.0%, 12.1% ± 7.8%, and 7.0% ± 9.1%, respectively, indicating that the gut and soil sequences of site HW were more similar than those of sites B and H. Gene libraries from earthworm gut content and soil contained 40 to 94 and 48 to 80 nosZ-positive clones, respectively (Table 1). The number of genotypes depended on the threshold for genotype definition, i.e., it decreased with increasing evolutionary distance (0.49 to 15.49) used (Table 1). nosZ libraries of the garden (H) yielded the highest number of nosZ genotypes for soil and for gut, followed by the field (B) and the meadow (HW). More nosZ genotypes were detected in gut than in soil libraries from sites B and H. The numbers of nosZ genotypes in gut and soil libraries from site HW were similar. Overall, the numbers of different genotypes detected in gut and soil libraries were similar.
Estimation of nosZ genotype richness in gene libraries.
In most libraries, the estimated number of genotypes was similar to the number of genotypes found (Table 1). Estimated genotype richness was highest for the gut and soil libraries of the garden site (H), followed by the field site (B) and the meadow site (HW) (Table 1). Estimated genotype richness depended on the estimator used but was mostly higher for the gut libraries than for the soil libraries from sites B and H. The estimated genotype richness for the soil library of site HW was mostly higher than that of the corresponding gut library (Table 1). However, at evolutionary distances of 7.49 and 15.49%, differences were only marginal for all sites.
Rarefaction analysis.
At evolutionary distance levels of D
0.49%,
3.49%,
7.49%, and
15.49%, most rarefaction curves of gut libraries were higher than those of corresponding soil libraries (Fig. 1A to G and I). Rarefaction curves of the gut and soil libraries from sites H and HW were essentially identical at D
15.49% (Fig. 1H) and D
3.49 (Fig. 1J), respectively. Rarefaction curves at D
7.49 (Fig. 1K) and D
15.49% (Fig. 1L) for libraries of site HW were higher for soil than for gut. However, differences between gut and soil rarefaction curves at all distance levels were not significant; confidence intervals overlapped. Confidence intervals of gut and soil curves tended to separate for site B only (Fig. 1A to D). The collective rarefaction analyses (Fig. 1) indicated that the nosZ diversities detected in soil and gut were similar.
![]() View larger version (62K): [in a new window] |
FIG. 1. Rarefaction curves with cloned nosZ fragments from the earthworm gut () and surrounding soil ( ) for various evolutionary distances (D). Error bars represent 95% confidence intervals.
|
Phylogenetic analysis.
nosZ phylogenetic trees that were generated from truncated (1.1-kb) and full-length (approximately 2-kb), in silico-translated nosZ sequences of pure cultures retrieved from the EMBL database (http://www.ebi.ac.uk) had congruent topologies (data not shown), indicating that the 1.1-kb nosZ fragments were sufficient for phylogenetic analyses. Thus, all sequences were truncated to 1.1 kb for the following analyses of the newly retrieved nosZ fragments.
Translated nosZ sequences from soil and the earthworm gut did not cluster with marine sequences, but some sequences clustered with translated nosZ sequences derived from a Michigan soil (45; see Fig. S2 in the supplemental material), indicating that terrestrial denitrifier populations differed from those found in marine habitats. The analysis revealed 14 phylogenetically distinct nosZ clusters, 13 of which contained both gut- and soil-derived sequences (Fig.2, GS clusters I to XIII). The mean evolutionary distance of translated nosZ sequences within a cluster was 6.3% ± 5.5%. Some of the gut-derived sequences were nearly identical (99 to 100% sequence similarity) (Fig. 3) to soil-derived sequences. Only one cluster (the G cluster) contained a single gut-derived sequence but no soil-derived sequence. nosZ sequences that occurred universally in all three of the gut libraries also occurred in at least one of the soil libraries (e.g., GS clusters II and IV).
![]() View larger version (35K): [in a new window] |
FIG.2. Unrooted maximum-likelihood tree of representative nosZ fragments (approximately 350 amino acids). Branching that was not reproduced with neighbor-joining and parsimony algorithms is drawn as multifurcations. Short sequences of 100 to 260 amino acids are termed "partial" and were added to the tree without changing the topology. The nomenclature for the sequences obtained in this study is as follows: soil or gut, retrieved from soil or the earthworm gut; identification code; and study site H (garden), HW (meadow), or B (field) (see Materials and Methods). Sequences generated with the new primers PsNosZ175F and PsNosZ1144R are indicated by the sequence identifier "Ps." Accession numbers are provided. G, gut; S, soil. The bar represents an estimated sequence dissimilarity of 10%. Sequences of some other studies (37, 38) were not considered, due to the small overlap (76 amino acids) with fragments from this study.
|
![]() View larger version (15K): [in a new window] |
FIG. 3. Pairwise comparisons of nosZ () and in silico-translated protein ( ) sequences from gut and soil. Only the highest values obtained within a cluster are shown. The similarity of a translated nosZ protein sequence is plotted as 1 D (as a percentage), where D is the PAM-corrected evolutionary distance.
|
|
View larger version (34K): [in a new window] |
FIG. 4. PCR amplification of Pseudomonas-related nosZ from site B with serial dilutions of template DNA. Lanes: 1, marker (1.0 and 0.75 kb); 2, negative control (water); 3, positive control (Pseudomonas aeruginosa); 4 to 6, soil DNA at a dilution of 100, 101, and 102, respectively; 7 to 9, gut DNA at a dilution of 100, 101, and 102, respectively.
|
Comparison of 16S rRNA gene and nosZ similarity.
The comparative phylogeny of nosZ and the 16S rRNA gene of nosZ-containing denitrifiers was assessed to determine if nosZ could be used for the identification of denitrifying species. Pairwise similarities of sequences retrieved from denitrifying isolates from this study and public databases ranged from 37 to 99% and from 60 to 100% for nosZ and 16S rRNA genes, respectively, indicating that the 16S rRNA is more conserved than nosZ (Fig. 5). Denitrifiers that shared <48% nosZ similarity (or 68% protein similarity of translated nosZ) always shared <97% 16S rRNA gene similarity. Thus, nosZ sequences with a similarity of <48% (or 68% protein similarity of translated nosZ) to any nosZ of cultured denitrifiers indicate a hitherto-unknown denitrifying species.
![]() View larger version (26K): [in a new window] |
FIG. 5. Correlation of 16S rRNA gene similarity with nosZ (A) and in silico-translated nosZ sequence similarity (B) of pure cultures. The similarity of the translated nosZ protein sequence is plotted as 1 D (as a percentage), where D is the PAM-corrected evolutionary distance. The solid lines indicate 16S rRNA gene threshold values for species delineation, and the dotted lines indicate threshold values below which sequences are indicative of novel denitrifier species. The dashed lines represent the expression 1 Q, where Q is the 90% quantile of pairwise sequence comparisons with a 16S rRNA similarity of >97%.
|
|
|
|---|
The absence of an endemic, earthworm gut-specific denitrifier population is in agreement with previous studies. The passage of ingested matter through the earthworm takes <20 h (2), presumably too short a time to sustain a gut-specific population in the gut lumen; in fact, such a population would require an unrealistically high growth rate. The earthworm gut has relatively little compartmentalization (16), and a gut-specific microbiota might be restricted to the gut wall. However, cell densities of earthworm gut wall-associated microbes are low (26, 49), their association with the gut wall appears to be opportunistic (18, 29, 35, 43), and the gut wall contributes only marginally to the production of N2O (25) and N2 (21). No evidence for an endemic, gut-specific microbiota has been obtained by cultivation or 16S rRNA gene-based analyses (18, 29, 35, 43).
Based on the above considerations, a quantitatively significant, gut-specific denitrifier population in the gut lumen or at the gut wall is unlikely. Thus, the proposal that the origin of the denitrifiers in the earthworm gut is ingested soil would seem to be a given. However, in general, the detected nosZ diversity of gut exceeded that of soil. For example, for site H, up to 22 different nosZ genotypes were detected in soil compared to 39 different nosZ genotypes in the earthworm gut (Table 1, site H). Coverage tended to be higher for soil libraries than for gut libraries (Table 1) and thus cannot account for this discrepancy. Selective feeding on microbe-rich particles, e.g., rhizosphere soil (15), and the proposed activation during gut passage might explain the slightly higher nosZ diversity in the gut libraries of sites B and H. The activation of microorganisms increases cell volume and genetic material (33). Active cells might also have increased lysis efficiency, which in comparison to inactive microorganisms would make them easier to detect by molecular methods.
Although pseudomonads are common in soil and can be readily cultured from the earthworm gut (25, 34), Pseudomonas-related nosZ sequences were only detected in the soil samples of the field site (B) but not in the corresponding gut samples when the general primer pair nosZ662F-nosZ1772R was used (Fig. 2). Although Pseudomonas-related nosZ fragments were obtained from both soil and gut samples from site B with newly designed Pseudomonas-specific primers PsNosZ175F and PsNosZ1144R (Fig. 2), semiquantitative PCR indicated that numbers of Pseudomonas-related nosZ targets in soil DNA extract were higher than in gut DNA extracts (Fig. 4) and, consequently, that the cell numbers of Pseudomonas-related denitrifiers were higher in soil than in the earthworm gut. Gammaproteobacteria, of which the genus Pseudomonas is a member, are abundant in the food source of Lumbricus but are not detected in anterior parts of the gut (41), and the number of Pseudomonas-related species decreases during gut passage (10, 47). Soil pseudomonads might be selectively digested during passage through the earthworm (10, 41), and cells with a high volume (e.g., cells that are highly active) are more likely to be mechanically disrupted in the gizzard of earthworms than are smaller cells (e.g., cells that are dormant or inactive) (16). The assumed digestion of soil pseudomonads during gut passage is consistent with the low detected Pseudomonas-related nosZ diversity in the gut compared to that of soil.
Libshuff analyses revealed no significant differences between nosZ gene libraries from gut and soil for the meadow site (HW) (see Table S1 in the supplemental material), while gut and soil libraries from the field (B) and garden (H) sites were only significantly different at a nosZ evolutionary distance of <12 and 5%, respectively (see Fig. S1 in the supplemental material). However, even a 12% evolutionary distance is not a major qualitative difference for nosZ (Fig. 5) and other functional genes such as the dissimilatory sulfite reductase genes (dsrAB), the ammonia monooxygenase gene (amoA), and the particulate methane monooxygenase gene (pmoA). PAM-corrected amino acid sequence dissimilarity of nosZ, dsrAB, amoA, and pmoA from organisms of the same genus ranges from 0 to 32% (mean, 8%), 2 to 26% (mean, 12%), 0 to 13% (mean, 6%), to 0 to 2% (mean 1%), respectively (data not shown), indicating that nosZ is the least-conserved marker of these four genes, apparently subjected to horizontal gene transfer. The mean PAM-corrected dissimilarity of dsrAB, amoA, and pmoA from organisms of the same genus equals 9%, and only ammonium oxidizers that show more than 15% amoA dissimilarity are always different species (28). The equivalent value for nosZ would be 32%, although 12% might be a more realistic threshold (Fig. 5). Consequently, a nosZ dissimilarity of 5 or 12% does not prove the occurrence of different denitrifying species; in addition, the obvious existence of extensive horizontal gene transfer leaves some uncertainty about the identity of the nosZ-carrying organisms. Nevertheless, the collective nosZ analyses (Table 1; Fig. 2, 3, and 5) demonstrated that the nosZ gene pool of the denitrifying populations in gut and soil had no significant differences. However, there is still some uncertainty about whether the denitrifier species composition in gut and soil is also similar.
In conclusion, the present study lends support to the hypothesis that the high denitrification potential of the earthworm gut (22, 25, 27, 31) and the concomitant in vivo emission of both N2O (27) and N2 (21) are not due to an earthworm gut-specific denitrifier population but due to the activation of ingested soil denitrifiers during gut passage.
We thank Achim Schmalenberger and Steffen Kolb for provision of dsrAB and pmoA alignments, respectively, and Erik Malzahn for technical assistance.
Supplemental material for this article may be found at http://aem.asm.org/. ![]()
Present address: Department of Microbiology, University of Aarhus, 8000 Aarhus C, Denmark. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»