Previous Article | Next Article ![]()
Applied and Environmental Microbiology, April 2005, p. 2001-2007, Vol. 71, No. 4
0099-2240/05/$08.00+0 doi:10.1128/AEM.71.4.2001-2007.2005
Copyright © 2005, American Society for Microbiology. All Rights Reserved.
Institute of Soil Ecology, GSF-National Research Center for Environment and Health, Neuherberg,1 Department of Organic Farming and Cropping Systems, University of Kassel, Witzenhausen, Germany,2 Swiss Federal Research Station for Agroecology and Agriculture, Zürich, Switzerland3
Received 31 August 2004/ Accepted 28 October 2004
|
|
|---|
|
|
|---|
Biological denitrification is a respiratory process defined as the enzymatic, stepwise reduction of the nitrogen oxides associated with electron transport phosphorylation and evolution of the gases nitric oxide (NO), nitrous oxide (N2O), and molecular nitrogen (N2). These end products are concomitantly released, leading to a loss to the environment of fixed nitrogen (3). Furthermore, the accumulation of the intermediates NO and N2O contributes to global warming and the destruction of the stratospheric ozone layer (8). NO and N2O fluxes between soil and atmosphere result from complex and little-understood microbial activities that proceed simultaneously in soil. Thus, the underlying communities of denitrifying bacteria and their responses to environmental conditions need to be studied in detail. Analyzing the process of denitrification in a system, e.g., the rhizosphere of a grain legume, that is known to enrich soil by fixing nitrogen is, therefore, highly relevant.
Denitrifying bacteria are phylogenetically diverse. Defined as a physiological group, these facultative anaerobes can switch from oxygen to nitrogen oxides as terminal electron acceptors when kept under anoxic conditions. The reduction of nitrite to nitric oxide by nitrite reductase is the first step that distinguishes denitrifiers from nitrate-respiring bacteria, which do not reduce nitrite to gas. Nitrite reductase is one of the key enzymes in the dissimilatory denitrification process. Two structurally different nitrite reductases are found among denitrifiers (35): one contains copper (Cu-Nir) encoded by the nirK gene, and the other contains heme c and heme d1 (cd1-Nir) encoded by the nirS gene. The two genes seem to occur mutually exclusively in a given strain, but both types have been found in different strains of the same species (9). No functionally significant differences between the two nitrite reductases have been reported.
Detection of mRNAs, which typically have short half-lives, provides a strong indication of specific gene expression at the time of sampling (25). More recently, in a number of studies, workers have used reverse transcription PCR (RT-PCR) approaches to investigate gene expression in a variety of environments, including expression of nahA (30) and pmoA (7) in groundwater; nifH expression in termite guts and lake water (21, 34); rbcL expression in lake water (32); merR expression in lake water sediment (20); narG, napA, nirS, nirK, and nosZ expression in estuarine sediments (22); and expression of the Desulfovibrio [NiFe] hydrogenase gene in anaerobic reactor (29).
In the present study we investigate the diversity of transcripts of the nirK and nirS genes in the rhizospheres of three economically important grain legumes, V. faba, L. albus, and P. sativum, and the occurrence of plant rhizosphere-specific phylotypes. The presence and expression of the genes were initially checked by amplification of the nirK and nirS genes with specific primers. The amplified products were then resolved by the technique of denaturing gradient gel electrophoresis (DGGE) to screen for variations between the replicates and also to calculate the level of similarities or differences among the three legumes under study. Representative samples were further cloned to identify the differences or similarities between the legumes.
|
|
|---|
Ten rectangular boxes (60 by 45 by 40 cm) per legume were prepared with about 70 kg of soil per box. Seeds of faba beans (V. faba L., cv. Scirocco; Saaten-Union, Iserlohn, Germany), white lupin (L. albus L., cv. Amiga; Saaten-Union), and peas (P. sativum L., cv. Duel; Saatbau Linz, Linz, Austria) were inoculated with legume-specific Rhizobium inoculants (Rhizobium leguminosarum bv. viciae for V. faba and P. sativum and Rhizobium lupini for L. albus from Radicin Institute, Iserlohn, Germany) according to the manufacturer's instructions, and three seeds were sown in each box. The plants were grown at 18°C with 70% humidity and a photoperiod consisting of 16 h of light and 8 h of darkness. The plants were watered equally every alternate day. The health of the plants was monitored during this period by visual observation of the color of their leaves. Sampling was performed after fruiting (5 months from the date of sowing). Plants and soil were removed from the pots. Subsequently, excess bulk soil was removed from the roots by shaking, leaving firmly adhering soil, which was defined as rhizosphere soil. Thirty replicates (10 pots with three plants in each) for each legume were sampled. Samples were stored separately at 80°C.
Nucleic acid extraction.
Nucleic acid extraction from the rhizosphere soil was performed by using the method of coextraction of DNA and RNA described by Griffiths et al. (14) with a few modifications (2). Prior to nucleic acid extraction, water was treated overnight with 0.1% diethyl pyrocarbonate, and all solutions were prepared in diethyl pyrocarbonate-treated water. Also, the glassware was baked at 180°C for 4 h. The method principally involved bead beating and solvent extraction of the nucleic acids. To prevent the degradation of mRNA, all incubations were performed on ice. To obtain pure DNA, incubation for 10 min at 37°C with RNase A (Sigma) at a final concentration of 100 µg ml1 of water was performed. Prior to RT, DNA was removed from RNA by treatment with DNase (1 U µl1) (RNase free; Promega) according to the manufacturer's instructions.
PCR and RT-PCR amplification.
PCR and RT-PCR targeting nirK and nirS were performed with the primers nirK1F-nirK5R and nirS1F-nirS6R, respectively (5). RT was performed as described below. Two microliters of total RNA sample was added to an 18-µl RT mixture prepared according to the manufacturer's instructions (Omniscript RT Kit; Qiagen, Hilden, Germany). nirK5R and nirS6R were used in the reaction for the nirK and nirS genes, respectively. Synthesized cDNA and DNA were further used for amplification according to the following protocol. One microliter of cDNA/DNA was added to a 48.5-µl PCR mixture that consisted of 29 µl of nuclease-free water, 5 µl of 10x reaction buffer, 5 µl of 2 mM deoxyribonucleoside triphosphate mixture, 1 µl each of 10 µM primers (nirK1F-nirK5R or nirS1F-nirS6R primer for nirK or nirS, respectively), 5 µl of 3% bovine serum albumin, and 2.5 µl of dimethyl sulfoxide. For DGGE fingerprinting, a 40-bp GC clamp (23) was attached to the 5' end of the forward primer, nirK1F. The reaction involved a hot start at 95°C for 10 min, followed by the addition of 0.5 µl of cloned Pfu DNA polymerase (2.5 U µl1; Stratagene, Amsterdam, The Netherlands). The cycling parameters for 30 cycles were 94°C for 1 min, 56°C for 1 min, and 72°C for 1 min, followed by a final extension at 72°C for 10 min. DNase-treated nucleic acids, without being subjected to RT, were used as controls in PCR to check for residual DNA in RNA preparations. The amplification products were analyzed by electrophoresis on 1.5% (wt/vol) agarose gels, followed by a 15-min staining with ethidium bromide (0.5 mg l1). Bands were excised and purified by using a QIAEX II gel extraction kit (Qiagen).
DGGE of RT-PCR products.
DGGE was performed by using 6% polyacrylamide gels (ratio of acrylamide to bisacrylamide, 37:1) with a 50 to 60% denaturant. A 100% denaturant is defined as 7 M urea plus 40% formamide (1). Appropriate volumes containing about 2 µg of the purified RT-PCR products, measured by absorbance at 260 nm, were loaded. The gels were electrophoresed at 60°C at 50 V for 17 h by using a D-Gene system (Bio-Rad Laboratories, Munich, Germany). The gels were silver stained according to the protocol described by Heukeshoven and Dernick (15). Dried gels were scanned by an HP Scanjet 7400c. The DGGE profiles obtained were analyzed by clustering via the unweighted pair group method with mathematical averages (Dice coefficient of similarity) by using GelCompar II software (Applied Maths, Kortrijk, Belgium). The position tolerance was set at 1%, and background subtraction was applied. Both strong and weak bands were included in the analysis, thus taking into account the presence and absence of bands at specific positions.
Cloning and restriction fragment length polymorphism (RFLP).
Three microliters of purified amplification products was cloned into the pCR-Blunt II-TOPO vector of a Zero Blunt TOPO PCR cloning kit (Invitrogen, Karlsruhe, Germany) according to the manufacturer's instructions. Ligations were transformed into chemically competent One Shot DH5
-T1 cells provided in the kit. A total of 100 clones were randomly selected and inoculated in Luria-Bertani medium (supplemented with 50 µg of kanamycin ml1), and plasmid was isolated by using a Qiagen Plasmid Mini Kit. Purified plasmids were tested for inserts by EcoRI (MBI Fermentas, Heidelberg, Germany) digestion with 1 U in a 50-µl reaction mixture.
Inserts were amplified by using the M13 forward and reverse primers provided in the Zero Blunt TOPO PCR Cloning Kit (Invitrogen) according to the manufacturer's instructions. Amplicons (2 µl each) were then screened by RFLP by digestion with 0.3 U of the restriction enzyme MspI (MBI Fermentas) in a final volume of 20 µl at 37°C overnight according to Braker et al. (6). The digested products were separated by gel electrophoresis on 4% (wt/vol) QA-Agarose High Resolution gels (Qbiogene, Heidelberg, Germany). To assess the number of clones sufficient to encompass the bacterial diversity, the collector's curve or species abundance curve (the number of operational taxonomic units detected plotted versus the number of clones analyzed) was constructed by using the software Analytic Rarefaction (http://www.uga.edu/
strata/software/Software.html). The RFLP types obtained were defined as operational taxonomic units (OTUs).
Sequencing and sequence analysis.
Inserts of representative samples for each RFLP type were sequenced on an ABI PRISM 310 genetic analyzer (Applied Biosystems, Foster City, Calif.) according to the standard protocol. Sequences were compared by using the National Center for Biotechnology Information BLAST program and aligned by using the ClustalW program of the European Molecular Biology Laboratory. Phylogenetic analysis was performed with a neighbor-joining algorithm and distance calculation according to Jukes and Cantor (16) by using TREECON version 1.3b (28). The aniA gene from Neisseria gonorrhoeae (accession no. M97926) was used as an outgroup for phylogenetic analysis, and tree topology was evaluated by bootstrap analysis by using 100 replicates.
Nucleotide sequence accession numbers.
The clone sequences determined in this study have been submitted to the GenBank database under accession numbers AY386223 to AY386234.
|
|
|---|
![]() View larger version (57K): [in a new window] |
FIG. 1. Agarose gel for detection of PCR and RT-PCR fragments for nirK. Lane M, 100-bp markers; lanes 1 to 3, PCR products for Vicia, Lupinus, and Pisum rhizospheres, respectively; lanes 4 to 6, RT-PCR products for Vicia, Lupinus, and Pisum rhizospheres, respectively; lanes 7 to 9: PCR products of DNase-treated nucleic acids not subjected to RT for Vicia, Lupinus, and Pisum rhizospheres, respectively; lane 10, negative control.
|
![]() View larger version (120K): [in a new window] |
FIG. 2. DGGE profiles of nirK RT-PCR products obtained from the three legumes. Lanes 1 to 3, profiles of Vicia, Lupinus, and Pisum, respectively; lanes 4 to 13, profiles obtained from 10 different pots of Pisum.
|
![]() View larger version (91K): [in a new window] |
FIG. 3. RFLP patterns of nirK inserts generated by using the restriction enzyme MspI. Lane M, 100-bp markers; lanes 1 to 12, the 12 different patterns into which all the clones could be grouped.
|
![]() View larger version (13K): [in a new window] |
FIG. 4. Collector's curve obtained by plotting the number of different OTUs obtained versus the number of clones analyzed for the clone libraries of Vicia, Lupinus, and Pisum.
|
![]() View larger version (33K): [in a new window] |
FIG. 5. Relative distribution of phylotypes (RFLP patterns) in the rhizospheres of Vicia, Lupinus, and Pisum.
|
![]() View larger version (44K): [in a new window] |
FIG. 6. Phylogram for nirK based on partial gene fragments. The tree is based on distance matrix analysis and a neighbor-joining method. The GenBank accession numbers follow pattern numbers. For environmental clones NCBI accession numbers have been given. The accession numbers of reference organisms have been indicated in brackets beside the names of the organisms. The scale bar indicates the expected number of changes per sequence position. Bootstrap values greater than 50 from 100 replicate trees are reported at the nodes. The sequence of aniA from N. gonorrhoeae (accession no. M97926) served as an outgroup to root the phylogram.
|
|
|
|---|
The three legumes used for the analysis of nir transcript diversity were selected because of reports on differences in rhizodeposition among the three plants (17, 18). Studies of microbial diversity based on 16S analysis from the rhizospheres of the legumes revealed major differences among the bacterial communities of the three legume rhizospheres under investigation (26), which formed the basis for detailed analysis of specific genes in the systems.
In the present study, successful amplification of nirK from both DNA and mRNA revealed not only its presence in the gene pool but also its expression in the rhizospheres of the three legumes under investigation. In contrast, the nirS gene fragment could not be detected even with rhizosphere DNA, thus indicating the absence of denitrifiers containing nirS sequences similar to those that can be targeted by the primer pair used. The apparent predominance of either nirK or nirS in rhizosphere soils may reflect differences in denitrifier enzyme properties in the rhizosphere, where this process is strongly selected, or it may be due to traits of the organism that are not connected to denitrification (24). However, the presence of nirS genes or transcripts in very low numbers, which go undetected by PCR, cannot be ruled out. Another possible explanation for the lack of detection of nirS from the three rhizospheres might be that the nirS primers selectively amplify nirS genes closely related to those used for the development of the primers (5). Thus, they may not have detected more distantly related nirS fragments, which might have been present in the rhizosphere soils under investigation. Therefore, it is not surprising that other authors who have used the same primer pair as the present study could not amplify nirS from soil samples. Avrahami et al. (4) could not detect nirS fragments from a silt loam soil from Germany, while Priemé et al. (24) could not detect nirS from a mixed deciduous upland forest soil in Michigan. Wolsing and Priemé (31) were able to detect nirS at only one time point from soil sites receiving mineral fertilizer.
The DGGE profiles for the three legumes contained 12 to 20 bands. This number is lower than that obtained by Throbäck et al. (27) on analysis of denitrifiers in arable soil. The discrepancy could be attributed to the fact that we used a cDNA-derived profile, which is a measure of only the active denitrifiers in the present study and not the total denitrifying potential present in the community. Furthermore, a rhizosphere effect, which causes selection of only a few genotypes, might have caused a further reduction in bands compared to the number obtained by Throbäck et al. (27), wherein bulk soil was analyzed. Similar to the results reported by Throbäck et al. (27), most of the visible bands had similar intensities while only a couple of bands were comparatively more intense. Interestingly, clear differences were visible among the profiles of the three legumes investigated. Vicia and Lupinus rhizosphere samples were more similar to each other than to the Pisum profile. The fingerprint generated by the Lupinus nirK RT-PCR product consisted of the maximum number of bands, indicating its denitrifying population to be the most diverse compared to the other two legumes under investigation. In a study by Sharma et al. (26) of structural bacterial diversity based on 16S rDNA and rRNA analysis of the three rhizosphere samples, the same trend was revealed, indicating a possible correlation between the structural community and the diversity of denitrifiers. However, DGGE has the limitations resulting from the comigration of fragments and formation of heteroduplexes, which may obscure the number of truly different observable bands (27). In a related study, Mayer et al. (17, 18) reported a rhizodeposition value of 21 mg kg1 of soil in Pisum, while the value was 61 and 52 mg kg1 of soil in Vicia and Lupinus, respectively. The percentages of N derived from rhizodeposition that could be traced in microbial biomass in Vicia and Lupinus were 20 and 22%, respectively; in Pisum the value was only 8%. The similarity in the profiles of active denitrifying communities in the rhizospheres of Vicia and Lupinus could be attributed to similar N rhizodeposition values between these two legumes compared to the value in Pisum.
RFLP analysis revealed that the 266 clones (83, 92, and 91 from the libraries of Vicia, Lupinus, and Pisum, respectively) could be grouped into 12 different patterns. The Lupinus library was the most diverse, with all 12 phylotypes represented. This confirmed the result obtained by the DGGE analysis of nirK products wherein Lupinus was found to be the most diverse when the banding patterns were observed. Priemé et al. (24) reported a higher diversity of nir genes from upland and marsh soils. In upland soil, 18 different phylotypes could be observed from the analysis of 50 clones. In marsh habitat the number of phylotypes increased to 37 when 93 clones were screened. This study was based on DNA analysis from bulk soil. In the rhizosphere, unlike in bulk soil, plant roots directly influence the survival, growth, and activity of bacteria in their vicinity (10). This explains the relatively reduced number of phylotypes observed in the present study. The representative sequences for each pattern were very divergent from those described for cultured denitrifiers forming a distinct cluster. Mesorhizobium sp. (accession no. AY078254) was the nearest culturable denitrifier according to the phylogram. This finding could be explained either by the lack of nirK gene sequences in the database or as an indication that only a particular group of bacteria that has a distinct nir operon is involved in denitrification in the rhizosphere of the three legumes. In a study by Mayer et al. (18) almost identical values for the mineral nitrogen content in the rhizospheres of the three legumes under investigation were reported.
These results reflect the diversity of active denitrifiers in the rhizospheres of grain legumes. A plant-dependent rhizosphere effect on the transcripts of a gene could be shown. Furthermore, the approach used here and the quantification of gene dosages, demonstrated for nitrate reductase genes (4, 19), should significantly enhance our understanding of bacterial denitrification in the environment. However, to use these data in agricultural practices to reduce the emission of N2O back to the environment, data about the correlation between transcription and expression are necessary. Additional research about the described denitrifying bacterial communities, mainly their distribution in time and space, and possibilities for the repression of their nir gene clusters are of great importance. The results of the present study, however, are an important step toward understanding denitrification in the rhizospheres of grain legumes.
|
|
|---|
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»