Previous Article | Next Article ![]()
Applied and Environmental Microbiology, February 2003, p. 1004-1012, Vol. 69, No. 2
0099-2240/03/$08.00+0 DOI: 10.1128/AEM.69.2.1004-1012.2003
Copyright © 2003, American Society for Microbiology. All Rights Reserved.
UMR INRA/Université de Bourgogne, Microbiologie et Géochimie du Sol, INRA-CMSE, 21065 Dijon Cedex, France
Received 8 August 2002/ Accepted 18 November 2002
|
|
|---|
|
|
|---|
Zones of high denitrifying activity occur predominantly in specific soil microsites (21) and in "activation sites," such as rhizospheric soil (15). Indeed, the proportion of fluorescent pseudomonads able to reduce nitrates appeared to be significantly higher in the rhizosphere than in the bulk soil. Even more, this proportion increases gradually and significantly in the vicinity of root (4). This observation has been made in the rhizosphere of different plant species cultivated in different soils (5). These results indicate that bacteria able to dissimilate nitrogen are selected in the rhizosphere. Furthermore, the nitrate reductase encoded by narG was recently shown to be involved not only in the rhizospheric competitiveness of Pseudomonas fluorescens C7R12 but more generally in its saprophytic competence in soil environments (17).
Altogether, these data suggest that enzymatic activities involved in denitrification would be involved in the saprophytic competence of fluorescent pseudomonads in soil environments. The aim of the present study was to compare the genetic diversity of the narG and nosZ genes to that of 16S ribosomal DNA (rDNA), commonly used as an evolution marker, within a large collection of fluorescent pseudomonads previously isolated from different bulk and rhizospheric soils (13). These two genes were chosen because the enzymes that they encode, nitrate reductase and nitrous oxide reductase, catalyze the first and last steps of denitrification, respectively. Furthermore, the nitrate reductase is known to be involved not only in denitrification but also in another major respiratory process, the dissimilatory reduction of nitrate into ammonia. In addition, the ability of the strains studied to reduce nitrate and/or nitrous oxide was assessed.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Phenotypic and genotypic diversity of the fluorescent pseudomonad strains analyzed in this study
|
Test for nitrogen dissimilation.
For preparation of inocula, strains were cultured in 10 ml of Luria broth (LB) medium (16). After 24 h of incubation, bacterial cells were collected by centrifugation at 8,000 x g for 20 min. The cells were resuspended in LB medium supplemented with KNO3 (1 mM) to obtain an absorbance at 580 nm of 2. For each strain, 2 ml of culture was aseptically introduced into 20-ml Vacutainer (Becton Dickinson) test tubes with a syringe though a rubber stopper. Nitrate reducers were revealed by the presence of NO2- in the culture, determined colorimetrically with the Griess Illosway reagent and by the absence of significant amounts of N2O when grown in tubes in which the atmosphere had been evacuated and replaced with a 90:10 He-C2H2 mixture to ensure anaerobiosis and avoid any reduction of N2O to N2 (34). Denitrifiers able to reduce N2O were considered total denitrifiers and determined on the basis of the absence of N2O accumulation when grown in tubes in which the atmosphere had been evacuated and replaced with helium to ensure anaerobiosis. Nitrous oxide concentration was measured with an MTI high-speed micro-gas chromatograph equipped with a catharometer detector (SRA Instruments). Strains unable to dissimilate NO3- in anaerobiosis were considered nondissimilators.
PCR amplification.
Amino acid sequences derived from the narG genes of P. aeruginosa PAO1 and P. fluorescens CFBP 5759 and AK15 (accession nos. Y15252, AF197465, and U71398, respectively) were aligned, and conserved regions were used for the design by back-translation of the PCR primers narGf (5'-GA[C/T]ATGCA[C/T]CC[A/C/G/T]TT-3') and narGr (5'-A[C/T]CCA [A/G]TC[A/G]TT[A/G]TC-3'). Amino acid sequences derived from the nosZ genes of P. aeruginosa PAO1, P. fluorescens CFBP 5759, P. denitrificans, Rhizobium meliloti, Bradyrhizobium japonicum, and Ralstonia eutropha (accession nos. X65277, AF197468, X74792, U47133, AJ002531, and X65278, respectively) were aligned, and conserved regions were used for the design by back-translation of the PCR primers nosZf (5'-AACGACAAG[G/A/T][C/T]CAA-3') and nosZr (5'-A[G/T][G/C]GC[A/G]TGGCAGAA-3'). The primers narGr plus nosZf or nosZr were used to amplify a fragment of narG (1,008 bp) and a fragment of nosZ (1,433 bp), respectively. The 16S rRNA gene (rDNA) was amplified with primer pair fD1 and rD1 (31).
PCR amplifications were carried out with cell suspensions obtained as described previously (14). PCR amplifications of narG and nosZ were performed in a total volume of 100 µl, by mixing 5-µl aliquots of cell suspension with 0.5 µM each primer, 50 µM each dATP, dCTP, dGTP, and dTTP, 4 U of Taq DNA polymerase (Q-BIOgene, Illkirsch, France), and PCR buffer (10 mM Tris-HCl [pH 9.0 at 25°C], 50 mM KCl, 1.5 mM MgCl2, 0.1% Triton X-100, 0.2 mg of bovine serum albumin per ml). The PCR mix for 16S rDNA was similar to those for narG and nosZ except for the primer (0.2 µM) and the deoxynucleotide triphosphate (40 µM each) concentrations. A negative control for amplification was prepared with 5 µl of H2O instead of 5 µl of cell suspension.
Amplifications were conducted in a thermal cycler GeneAmp PCR System 9600 (Perkin Elmer Applied Biosystems, Foster City, Calif.) with an initial denaturation step (3 min at 95°C) followed by 35 cycles of 1 min at 94°C, 1 min at 50°C for narG and 55°C for nosZ and 16S rDNA, 1 min for narG and 2 min for nosZ and 16S rDNA at 72°C, and a final extension of 3 min at 72°C. Aliquots (5 µl) of each PCR product were analyzed by electrophoresis in 0.8% agarose gels stained with ethidium bromide and photographed under UV illumination with Polaroid type 665 positive-negative film.
Extraction and digestion of total genomic DNA.
For each strain, extraction and purification of the total genomic DNA were performed with previously described methods (2). Aliquots of 1 µg of total genomic DNA were digested overnight in 10 µl with 30 U of EcoRI (Roche Diagnostic, Meylan, France). This enzyme was chosen because there is only one site of attack in the narG sequence and no site in the nosZ sequence of P. fluorescens CFBP 5759, both used as probes for hybridization.
Hybridization of total genomic DNA and PCR products.
narG and nosZ probes were prepared from PCR products of strain P. fluorescens CFBP 5759 obtained as described above. narG and nosZ PCR products, corresponding to 1,008-bp and 1,433-bp DNA fragments, respectively, were purified after electrophoresis on a 0.8% agarose gel with a Qiaex II gel extraction kit (Qiagen S.A., Courtaboeuf, France). The purified fragments were analyzed in a 0.8% agarose gel and quantified by comparison with the molecular mass markers Smart Ladder (Eurogentec, Seraing, Belgium). One hundred nanograms of each fragment was labeled for 20 h with digoxigenin with a DIG-High Prime kit (Roche Diagnostic), which theoretically produces 1,500 ng of labeled DNA.
The restriction fragments of total genomic DNA were separated by electrophoresis in 0.8% agarose gels in TAE buffer (40 mM Tris-HCl [pH 7.9], 4 mM sodium acetate, and 1 mM EDTA). DNA was transferred from the gel to a positively charged nylon membrane (Pall Gelman Laboratory, Ann Arbor, Mich.) with the Vacugene XL apparatus (Pharmacia Biotech, Orsay, France) according to the manufacturer's instructions. The DNA was fixed to the membrane by baking at 120°C for 30 min. Southern blots were hybridized successively with the narG and nosZ probes. Hybridization was performed according to the manufacturer's instructions with a probe concentration of 25 ng per ml of hybridization solution. Hybridized fragments were revealed with the DIG luminescent detection kit (Roche Diagnostic).
The PCR products of narG and nosZ were hybridized with the corresponding probes defined above to confirm that these products correspond to the targeted sequences. Aliquots of 1 µl of narG PCR product or 5 µl of nosZ PCR product were analyzed by electrophoresis in TAE buffer. PCR products were then transferred from the gel to a positively charged nylon membrane (Pall Gelman Laboratory, Ann Arbor, Mich.) and hybridized as described above.
Restriction fragment length polymorphism assay of narG, nosZ, and 16S rDNA.
Discriminant restriction enzymes for narG and nosZ were selected from comparison of the sequences previously used to design the PCR primers. Aliquots (5 µl) of PCR products were digested overnight with 5 U of the following restriction endonucleases: AvaII, BamHI, HinfI, NciI, and NlaIII for narG; AvaII, EcoRI, HincII, HindIII, HinfI, and NarI for nosZ; and the 13 endonucleases described by Laguerre et al. (12) for 16S rDNA. The restriction fragments were separated by electrophoresis at 100 V in TAE buffer with 4% Nusieve 3:1 agarose (FMC, Rockland, Maine) on gels containing 0.4 µg of ethidium bromide per ml. Gels were photographed under UV illumination with Polaroid type 665 positive-negative film.
For each gene, each strain was assigned a composite type defined by the combination of the patterns obtained with the different restriction endonucleases. The computer program NTSYS (27) was used to estimate the relationships among narG, nosZ, and 16S rDNA types from the fragments obtained with the different restriction endonucleases. The narG, nosZ, and 16S rDNA types were compared by pairwise coefficients of similarity calculated by the Jaccard similarity coefficient (28). For each gene, the matrix of similarities was displayed as a dendrogram with the UPGMA (unweighted pair group method arithmetic) algorithm (28). For each possible combination, narG/16S rDNA, nosZ/16S rDNA, and narG/nosZ, the level of correlation between the genes was determined with pairwise similarities between strains.
|
|
|---|
The 28 strains isolated from soil environments, including the two reference strains AK15 and CFBP 5759, were either total denitrifiers (57.1%), nitrate reducers (14.3%), or nondissimilators (28.6%). Strains of P. fluorescens were total denitrifiers (81.8%) or nondissimilators (18.2%), whereas strains of P. putida were nitrate reducers (37.5%) or nondissimilators (62.5%). All P. lini strains were total denitrifiers, and the two strains of P. chlororaphis were nitrate reducers.
Detection of narG and nosZ.
Primers narGf and narGr permitted the amplification of a single DNA fragment of about 1,000 bp for 19 of the 32 strains analyzed (Table 1). A positive hybridization was recorded for all the PCR products. Additionally, hybridization of the genomic DNA with this probe gave a positive signal for the 19 strains.
Primers nosZf and nosZr permitted the amplification of a single DNA fragment of about 1,450 bp for 19 of the 32 strains analyzed. A positive hybridization was recorded for all the PCR products. Additionally, hybridization of the genomic DNA with this probe gave a positive signal in the 19 strains. Despite optimization of PCR conditions (data not shown), the quantities of PCR products for both P. chlororaphis strains were too low to allow restriction fragment length polymorphism analysis of their nosZ genes but high enough to hybridize with the nosZ probe.
Analysis of restriction fragment length polymorphisms of narG and nosZ.
The PCR products of the narG genes were digested with each of the five restriction enzymes AvaII, BamHI, HinfI, NciI, and NlaIII for the 19 strains that gave a positive amplification with the narGf and narGr primers. Depending on the restriction enzyme, two to nine distinct restriction patterns were resolved. A total of 69 different restriction fragments were obtained, all of them being polymorphic. The 19 strains analyzed with the five restriction enzymes gave 14 different combinations of patterns representing 14 narG types (Table 1). The level of similarity among the narG types ranged from 0.05 to 0.79, indicating a high polymorphism of narG among the strains of fluorescent pseudomonads tested (Fig. 1).
![]() View larger version (22K): [in a new window] |
FIG. 1. Dendrogram showing the relationships between the fluorescent Pseudomonas strains based on variation in restriction sites in the narG gene. The pairwise coefficients of similarity (Jaccard) were clustered with the UPGMA algorithm of NTSYS-pc. Numbers 1 through 14 indicate the 14 narG types.
|
The PCR products of the nosZ genes were digested with each of the six restriction enzymes AvaII, EcoRI, HincII, HindIII, HinfI, and NarI for the 17 strains that gave a clear positive amplification with the nosZf and nosZr primers. Depending on the restriction enzyme, three to seven distinct restriction patterns were resolved. A total of 57 different restriction fragments were obtained, all of them being polymorphic. The 17 strains analyzed with the six restriction enzymes gave 11 different combinations of patterns representing nine nosZ types (Table 1). The level of similarity among the nosZ types ranged from 0.04 to 0.88, indicating a high polymorphism of nosZ among the strains of fluorescent pseudomonads tested (Fig. 2).
![]() View larger version (20K): [in a new window] |
FIG. 2. Dendrogram showing the relationships between the fluorescent Pseudomonas strains based on variation in restriction sites in the nosZ gene. The pairwise coefficients of similarity (Jaccard) were clustered with the UPGMA algorithm of NTSYS-pc. Numbers 1 to 11 indicate the 11 nosZ types.
|
Analysis of restriction fragment length polymorphisms of 16S rDNA.
The combined patterns obtained with 13 restriction endonucleases gave 15 distinct ribotypes among the 32 strains analyzed (Fig. 3). The ribotypes of most reference and type strains were published previously (12). The level of similarity among the ribotypes ranged from 0.77 to 0.99, indicating a low polymorphism of 16S rDNA among the strains of fluorescent pseudomonads tested. The 16S rDNA type Aer, including the P. aeruginosa type strain (CFBP 2466T), was the most distant from the others and was thus used to root the dendrogram.
![]() View larger version (27K): [in a new window] |
FIG. 3. Dendrogram showing the relationships between the 16S rDNA types. The pairwise coefficients of similarity (Jaccard) were clustered with the UPGMA algorithm of NTSYS-pc. P. lini CFBP 5732 and CFBP 5733 were identified by the double C3/C2 16S rDNA type described previously (13). Bacterial strains underlined were shown to harbor both the narG and nosZ genes; bacterial strains not underlined were shown to harbor neither the narG gene nor the nosZ gene.
|
Relationships among polymorphisms of narG, nosZ, and 16S rDNA.
For each combination of gene types, narG/16S rDNA (Fig. 4A), nosZ/16S rDNA (Fig. 4B), and nosZ/narG (Fig. 4C), the level of correlation was determined by pairwise similarities between strains.
![]() View larger version (14K): [in a new window] |
FIG. 4. Correlation between pairwise coefficients of similarity (Jaccard) of the narG gene and 16S rDNA (A), the nosZ gene and 16S rDNA (B), and the narG and nosZ genes (C). Equation of linear regressions 1 and 2 were y = 1.1492x + 0.2117 and y = 0.2059x + 0.0368, respectively. Linear regressions 1 and 2 were calculated with the pairwise coefficients of similarity represented by triangles and diamonds, respectively.
|
Comparison of the two functional genes, narG and nosZ, allowed us to identify two different groups among the pairwise comparisons (Fig. 4C). For each of these two groups, a linear regression between the levels of similarity of the two genes in the different strains was calculated. The slope values, obtained after calculation of the equations of the two regressions, were equal to 1.15 (regression 1) and 0.21 (regression 2), and their correlation coefficients were equal to 0.87 and 0.45, respectively. All species having narG and nosZ genes were represented in these two groups.
|
|
|---|
The ability of the strains studied to reduce nitrate and nitrous oxides was assessed. A large majority of the strains analyzed (71.4%) expressed the ability to reduce nitrate, and among them, 78.2% were able to reduce nitrate into dinitrogen. The ability to only reduce nitrate in nitrite was present at a low frequency within our collection of soil-borne pseudomonads (28.6%). Overall, the ability to express the dissimilatory activities was well related to the presence of the corresponding genes (narG and nosZ). Most of the strains harbored both genes, a few of them harbored only the narG gene, and none of them harbored only the nosZ gene. Our experimental procedures did not enable us to check the absence of nitrous oxide reduction by strains missing the nosZ gene. Nitrate reductase activity was recorded in some strains missing the narG gene (P. putida CFBP 5744, CFBP 5745, and CFBP 5746), suggesting the presence of another enzyme performing nitrate reduction, such as the periplasmic nitrate reductase previously described in other fluorescent pseudomonads (3, 24). In contrast, the two strains of P. chlororaphis were not able to reduce nitrous into dinitrogen despite a positive hybridization signal of total DNA and PCR products with the nosZ probe. However, the quantity of these PCR products was very low, suggesting either that the primers were not specific enough, that other sequences than nosZ were amplified, or that the corresponding genes were not functional.
The presence in the bacteria of the dissimilatory narG and nosZ genes and the diversity of these genes were analyzed in relation to the ribotypes of the bacterial strains. 16S rDNA is commonly used for diversity analyses of bacterial populations (19), especially of fluorescent pseudomonads (11, 12, 13). 16S rDNA is also considered a genetic marker for phylogeny studies in bacteria (32) and has been applied by several authors when studying the phylogeny of functional genes in bacteria (10). In contrast to functional genes, which are able to accumulate mutations without being impaired (e.g., variation of the third position of the codon), the level of acceptable variation in the 16S rRNA gene (not lethal) is low due to the functionality in this gene, which is strongly related to its structure. Consequently, the 16S rRNA gene has an evolutionary rate significantly lower than that of functional genes (0.02 and 0.7 to 0.8 per 106 years, respectively) (18).
The distribution of strains among the different ribotypes was clearly related to the presence of the narG and nosZ genes. Indeed, the narG and nosZ genes were never detected in the strains belonging to the ribotypes included in the first cluster of the 16S rDNA dendrogram. In contrast, both narG and nosZ were detected in all strains belonging to the ribotypes included in the second cluster of the 16S rDNA dendrogram. This observation suggests that the strains with narG and nosZ genes have an evolutionary history different from that of those harboring neither of the dissimilative genes.
Polymorphism analysis of the two denitrifying genes indicated clearly the wide range of similarities within the bacterial strains tested. Indeed, the similarity of the narG and nosZ genes varied from 0.05 to 0.79 and from 0.04 to 0.88, respectively. As expected, this range of similarities was a lot lower than that of the 16S rRNA gene. This observation suggests that the evolution rates of ribosomal and functional genes differ. This is in agreement with the estimations made by Ochman and Wilson (18). No correlation was found between the similarity indexes of 16S rDNA and denitrifying genes (narG or nosZ). In contrast, Philippot (22) reported a congruence between the 16S rDNA and narG genes when analyzing denitrifying genes of strains belonging to different genera in the major microbial groups. The lack of congruence between 16S rDNA and the dissimilatory genes recorded in the present study could be related to the strains analyzed, which belong to species which are closely related (fluorescent Pseudomonas spp., oxidase positive) and not well defined (1), with ribotypes showing a high level of similarity (12).
Since narG and nosZ genes are known to be organized in two distinct operons (22, 35) and to encode enzymes contributing to nitrogen dissimilation, a similar rate of evolution of these two functional genes could be expected. However, the pairwise comparison of similarity indexes of narG and nosZ in the different strains leads to the identification of two different groups of strains. These groups were distributed along two linear regressions showing different slopes (1.15 and 0.21). The slope value close to l in the first linear regression indicates that the similarity indexes of the narG and nosZ genes harbored by the corresponding strains varied within the same range. This observation suggests that these two genes have followed a similar evolution and therefore have a similar molecular clock. In the second group of strains, the slope value of the linear regression was lower than 1, indicating that the range of diversity of the nosZ gene is higher than that of the narG gene. This suggests that the two genes had different evolutionary rates.
Altogether, our data suggest that in the strains belonging to the first group (slope equal to 1.15), the narG and nosZ genes were present together in more ancestral strains than in the strains belonging to the second group (slope equal to 0.21). According to Petri and Imhoff (22) and to Zumft (35), the respiratory nitrate reductase gene had a long evolutionary history and already played a key function in energy metabolism during preoxic times. To our knowledge, bacteria harboring only the nosZ gene have never been described. These observations, taken together, support the hypothesis that the acquisition by bacteria of nitrous oxide reduction is probably a more recent event than that of the ability to reduce nitrate. This would suggest that strains harboring both genes were favored in soil environments.
This hypothesis is supported by the significantly higher frequency among soil-borne pseudomonads of strains which are total denitrifiers (57.1%) compared to those which are nitrate reducers (14.3%), this difference having already been stressed in our previous work (4). This hypothesis is also supported by the observations recently made in our group indicating that the fluorescent pseudomonads showing the best rhizosphere competence are total denitrifiers (S. Delorme, T. Corberand, Y. Dessaux, L. Gardan, J.-M. Meyer, J. Raaijmakers, and P. Lemanceau, unpublished data). Altogether, these observations suggest that the ability to reduce nitrate to molecular nitrogen confers a selective advantage to fluorescent pseudomonads in soil environments.
This work was partly supported by the INCO-DC Program ERBIC18CT970180 and the Conseil Régional de Bourgogne.
|
|
|---|
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»