Previous Article | Next Article ![]()
Applied and Environmental Microbiology, June 2004, p. 3541-3551, Vol. 70, No. 6
0099-2240/04/$08.00+0 DOI: 10.1128/AEM.70.6.3541-3551.2004
Copyright © 2004, American Society for Microbiology. All Rights Reserved.
R. C. Hamelin,2 L. Bernier,3 and M. St-Arnaud1*
Institut de Recherche en Biologie Végétale, Jardin Botanique de Montréal, Montréal, Québec, Canada H1X 2B2,1 Laurentian Forestry Centre, Natural Resources Canada, Sainte-Foy, Québec, Canada G1V 4C7,2 Centre de Recherche en Biologie Forestière, Université Laval, Sainte-Foy, Québec, Canada G1K 7P43
Received 14 October 2003/ Accepted 5 February 2004
|
|
|---|
|
|
|---|
Some rhizosphere-inhabiting microbes are known to participate in detrimental interactions with the plants that lead to the development of root diseases, while others are clearly beneficial to the plants by stimulating their growth, enhancing nutrient uptake, or playing a significant role in controlling soilborne pathogens and the diseases they cause (36). This biological approach to the control of root diseases, known as biological control, relies on the idea that rhizosphere communities can be manipulated to promote natural disease suppression by mechanisms such as competition, production of antibiotics, parasitism, and induced resistance (3, 12, 39). However, very little is known about which rhizosphere microbial taxa play a significant role in controlling soilborne plant pathogens and root diseases, as well as how soilborne plant pathogens and disease development may alter rhizosphere microbial communities and the ecology of the rhizosphere under natural or controlled environmental conditions. We have clearly not yet reached the full potential concerning the use and management of disease-suppressive microbial communities to control root diseases. Therefore, more information is needed to better understand the role and implication that different rhizosphere-inhabiting microorganisms might play in biological control.
Root rot diseases, mainly caused by fungal pathogens, are common in conifer nurseries (13). The impact of these diseases on nursery seedling production is reflected by economic losses. Root rot diseases may not only cause mortality at nursery sites but may also affect survival rates after outplanting to reforestation sites. Mortality rates as high as 40% for the first year after outplanting have been reported for Picea mariana (25). An important factor is that symptoms are not always visible on infected seedlings, as the disease can remain latent and develop only after transplantation (11). The main pathogens that cause root diseases of conifers under nursery conditions belong to the fungal genera Cylindrocladium, Rhizoctonia, and Fusarium and the Oomycota genera Pythium and Phytophthora (13). To date, very few efficient biologically based control strategies have been developed to control these pathogens under conifer nursery conditions (23). A better understanding of the rhizosphere microbial ecology and diversity associated with conifer seedlings grown in tree nurseries would most certainly benefit the development of efficient biological control strategies.
Therefore, the aim of this study was to use a molecular culture-independent approach to (i) profile the main bacterial and fungal rhizosphere taxa associated with healthy black spruce (Picea mariana) seedlings and seedlings with symptoms of root rot grown in a conifer nursery and (ii) contrast the differences in terms of rhizosphere microbial community composition with respect to the health status of the seedlings.
|
|
|---|
DNA extraction.
Total genomic DNA was extracted from 250 mg of a freeze-dried soil subsample using the MoBio UltraClean Soil DNA Isolation kit (MoBio Laboratories, Inc., Solana Beach, Calif.) according to the manufacturer's protocol with the following modifications: a FastPrep apparatus (Thermo Savant, Holbrook, N.Y.) was used at speed level 4 for 15 s instead of vortexing the soil samples for initial cell disruption. The purified DNA was resuspended in 50 µl of solution S5 (MoBio Laboratories, Inc.) and stored at 20°C until PCR amplification.
Bacterial and fungal rRNA gene amplification.
PCR amplification targeting bacterial partial 16S rRNA gene was performed using the 968f and 1401r primers (24), while the NS1 and NS2 primers (37) were used for fungal partial 18S rRNA gene amplification. Five microliters of 10x PCR buffer (Qiagen, Mississauga, Ontario, Canada), 5 µl of a 5 µM concentration of each bacterial or fungal primer, 1 µl of a 10 mM concentration of each deoxyribonucleoside triphosphate, 2 µl of a 1/20 dilution of soil genomic DNA or 2 µl of double-distilled H2O (negative control), 1.25 U of Taq DNA polymerase (Qiagen), and double-distilled H2O were combined in a final volume of 50-µl PCR mixtures. Reaction mixtures were overlaid with mineral oil prior to PCR amplification. PCR amplification was performed using a MJ Research PT-100 thermal cycler and the following steps: (i) an initial denaturation step (3 min at 94°C); (ii) 35 cycles, with 1 cycle consisting of denaturation (1 min at 94°C), annealing (1 min at 62°C [primers 968f and 1401r] or 55°C [primers NS1 and NS2]), and elongation (2 min at 72°C); and (iii) a final extension step (10 min at 72°C). Products were visualized by 1.5% agarose gel electrophoresis in 1x Tris-acetate-EDTA (TAE) buffer and stained with ethidium bromide (1% [wt/vol]). Bands were excised using a sterile scalpel, and DNA was purified from gel slices by using the QIAquick gel extraction kit (Qiagen).
Shotgun cloning and construction of bacterial and fungal rRNA gene libraries.
Each replicate library (six libraries of genes from the rhizosphere of three healthy seedlings and three diseased seedlings) was produced using different soil samples. PCR amplicons obtained from each of the three samples of the rhizosphere from three healthy or diseased seedlings using the 968f-1401r or NS1-NS2 primer set (total of 12 different samples) were each ligated into pCR4-TOPO vectors and transformed into One Shot TOP10 chemically competent Escherichia coli following the manufacturer's protocol (Invitrogen, Carlsbad, Calif.). Seventy randomly picked colonies per sample were each grown overnight at 37°C in 5 ml of Luria-Bertani broth containing 100 µl of ampicillin/ml and stored in 25% glycerol at 80°C. Plasmid DNA extraction was performed using the QIAprep Spin Miniprep kit (Qiagen), and the size of the insert was determined by both EcoRI restriction digest analysis and PCR amplification using the 968f-1401r or NS1-NS2 primer set. For each cloning reaction (total of 12), 50 clones containing an insert of the expected size (approximately 433 bp using 968f-1401r or 555 bp using NS1-NS2) were used for ARDRA.
ARDRA and sequencing.
After PCR amplification of cloned rRNA genes using the 968f-1401r or NS1-NS2 primer set, 10 µl of each amplicon (600 amplicons in total) were digested for 2 h at 37°C with 2.5 U of HinfI, 2.5 U of HaeIII, and 2.5 U of RsaI (New England Biolabs, Mississauga, Ontario, Canada) in 15-µl reaction mixtures. Restriction patterns were visualized by 3% agarose gel electrophoresis in 1x TAE buffer and stained with ethidium bromide (1% [wt/vol]). Gel images were digitalized, and restriction patterns were analyzed using the Quantity One software (Bio-Rad, Mississauga, Ontario, Canada). Each different restriction pattern was defined as an operational taxonomical unit (OTU). Plasmid DNA from one of each different representative OTU was commercially sequenced (DNA Landmarks, St-Jean-sur-Richelieu, Québec, Canada) on an ABI Prism 3700 DNA analyzer (Applied Biosystems, Foster City, Calif.) using T3 and T7 primers.
Sequence analysis.
DNA sequences were edited, and consensus sequences were obtained using DNAMAN version 4.13 (Lynnon BioSoft.). All sequenced clones were analyzed for the presence of chimeras using the CHIMERA CHECK program (version 2.7) from the Ribosomal Database Project (RDP) (http://rdp.cme.msu.edu/html/). Sequences suspected of being chimeric were not included in further analysis. All remaining sequences were screened against those in the GenBank database using BLASTn (1). On the basis of BLASTn similarity searches, sequences that were suspected of not being of bacterial or fungal origin were discarded. The most homologous sequences found in the GenBank database were used to construct a multiple-sequence alignment using ClustalX (version 1.81) (33) for both the bacterial and fungal clone libraries. A neighbor-joining analysis was performed using PAUP (version 4.0b10) (31) with 1,000 bootstrap replicates. Phylogenetic trees were edited using Treeview (version 1.6.6).
Diversity indices.
Clones having the same OTU were grouped, and the frequency of each OTU was determined. A number of diversity indices were calculated for each of the bacterial and fungal libraries. These indices included three different species richness indicators: (i) the total number of OTUs per library; (ii) the Shannon diversity index, a general diversity index that considers both species richness and evenness; and (iii) the Margalef index, which measures richness by using a ratio between the number of OTUs and the ln function of the number of clones analyzed (15). Three different indices were also calculated for evenness and dominance: (i) the Shannon evenness index, which is the ratio of the Shannon diversity index to the maximum possible value that could theoretically be obtained with the observed number of OTUs; (ii) the Simpson's index, which gives the probability that two clones chosen at random will be from the same OTU; and (iii) the Berger-Parker index, which is the relative abundance of the most abundant OTU (15). To determine whether the rRNA gene sequences obtained for each library of genes from the rhizosphere of healthy and diseased seedlings were different, homologous coverage curves, Cx(D), and heterologous coverage curves, Cxy(D) were calculated by a Cramér-von Mises test statistic and compared by a Monte Carlo test procedure (28). Similarity coefficients, which reflected the proportions of shared OTUs, were calculated by performing pairwise comparisons of the libraries of genes from the rhizosphere of healthy and diseased seedlings (21).
Nucleotide sequence accession numbers.
The nucleotide sequences determined in this study have been deposited in the National Center for Biotechnology Information (NCBI) database under accession numbers AY321198 to AY321281 (bacterial clones) and AY321682 to AY321712 (fungal clones).
|
|
|---|
To confirm that the grouping of the different OTUs using 16S and 18S partial rRNA gene digestion with HinfI, HaeIII, and RsaI was reproducible and reliable for discriminating between the different bacterial and fungal clones, three replicates of three different bacterial and fungal OTUs were randomly selected and sequenced. The results revealed no important difference in sequence composition between the replicates of each of these different OTUs (data not shown), suggesting that the OTU grouping based on restriction digestion patterns was efficient for discriminating between bacterial and fungal rRNA gene clones analyzed in this study.
Of the partial bacterial 16S rRNA gene fragments amplified from the rhizosphere of healthy and diseased P. mariana seedlings and shotgun cloned, 300 different inserts (150 from healthy and 150 from diseased seedling treatments) were analyzed by ARDRA. ARDRA profiles revealed 84 different OTUs, and each OTU was sequenced. Thirty-four OTUs were found to be associated only with the rhizosphere of healthy seedlings, 26 OTUs were associated only with the rhizosphere of diseased seedlings, and 24 were found associated with the rhizosphere of both healthy and diseased seedlings. Phylogenetic analysis of the different OTUs and their closest relatives in the GenBank database was used to assign each environmental clone to a major bacterial group. The results of the sequence analysis of these OTUs, their phylogenetic affiliation, and their incidence among the rhizosphere of healthy and diseased seedlings are shown in Table 1.
|
View this table: [in a new window] |
TABLE 1. Bacterial sequences from the GenBank database with the highest similarity to each OTU and distribution of OTUs in the rhizosphere of healthy and diseased P. mariana seedlings
|
![]() View larger version (42K): [in a new window] |
FIG. 1. Unrooted neighbor-joining tree representing the phylogenetic relationships of 84 bacterial OTU sequences from the rhizosphere of healthy and diseased P. mariana seedlings to the most closely related sequences obtained from BLAST searches. The numbers at the nodes of the tree are bootstrap values for each node with 1,000 bootstrap resamplings (values below 50 are not shown). The bar represents 0.1 substitution per site.
|
|
View this table: [in a new window] |
TABLE 2. Relative abundance of clones with respect to different bacterial and fungal taxa from triplicate soil samples obtained from the rhizosphere of healthy and diseased P. mariana seedlings
|
Phylogenetic analysis of the different OTUs and their closest relatives in the GenBank database was performed to assign each environmental clone to a major fungal group. The results of the sequence analysis of these OTUs, their phylogenetic affiliation, and their incidence among the rhizosphere of healthy and diseased seedlings are shown in Table 3. Phylogenetic analysis indicated that the fungal DNAs amplified from the rhizosphere of P. mariana are taxonomically diverse and can be grouped into seven different groups belonging to the Ascomycota and Basidiomycota (Fig. 2). Healthy rhizosphere-associated fungal 18S rRNA gene clones were not associated with the same main groups and in the same proportions as the ones associated with diseased seedlings (Table 2). The majority of the clones found in the rhizosphere of healthy seedlings grouped with the Homobasidiomycetes (76%), while the ones found in the rhizosphere of diseased seedlings mainly grouped with the Homobasidiomycetes (38%) and Soradiomycetes (29%). Clones affiliated with the Eurothiomycetes and Pezizomycetes were detected from the rhizosphere of only one healthy seedling, while clones affiliated with the Dothideomycetes, mitosporic Ascomycota, and Heterobasidiomycetes were detected only from the rhizosphere of diseased seedlings.
|
View this table: [in a new window] |
TABLE 3. Fungal sequences from the GenBank database with highest similarity to each OTU and distribution of OTUs in the rhizosphere of healthy and diseased P. mariana seedlings
|
![]() View larger version (27K): [in a new window] |
FIG. 2. Unrooted neighbor-joining tree representing the phylogenetic relationships of 31 fungal OTU sequences from the rhizosphere of healthy and diseased P. mariana seedlings to the most closely related sequences obtained from BLAST searches. The numbers at the nodes of the tree indicate bootstrap values for each node with 1,000 bootstrap resamplings (values below 50 are not shown). The bar represents 0.1 substitution per site.
|
Similarity coefficients were calculated by performing pairwise comparisons of the different replicates of bacterial and fungal libraries (Table 4). The similarity coefficients obtained between healthy-healthy and diseased-diseased seedling replicates for bacterial or fungal libraries were all higher than the ones obtained for comparisons between healthy-diseased seedling treatments, indicating that all replicates from libraries of genes from the rhizosphere of healthy or diseased seedlings were more similar to each other than to different treatments.
|
View this table: [in a new window] |
TABLE 4. Pairwise comparisons of bacterial and fungal rRNA gene library compositions obtained from triplicate rhizosphere soil samples associated with healthy and diseased P. mariana seedlings
|
|
View this table: [in a new window] |
TABLE 5. Diversity indices for bacterial and fungal OTUs associated with the rhizosphere of healthy and diseased P. mariana seedlings
|
|
|
|---|
To our knowledge, the present work is the first reported culture-independent study aimed at characterizing the microbial diversity that exists under conifer nursery conditions. As a small number of samples was studied, caution has to be taken when drawing conclusions from such a limited set of data. Nevertheless, the results obtained clearly indicated that the bacterial and fungal diversity associated with the rhizosphere of P. mariana seedlings grown under nursery conditions was considerable. Different diversity and phylogenetic analyses (see Results) supported the hypothesis that the bacterial and fungal communities associated with the rhizosphere of healthy and diseased P. mariana seedlings are different. Similarity coefficients, calculated by performing pairwise comparisons to assess the variability associated with the sampling strategy, clearly indicated that the triplicate rRNA gene libraries used for healthy and diseased seedling-associated treatments were more similar within than among treatments). These analyses supported the view that our sampling was sufficient to provide a realistic portrait of the microbial diversity in this system and that the variability obtained, always inherent with any such studies, did not compromise the interpretation of the results.
Our results clearly support the hypothesis that the bacterial and fungal communities were altered in relation to the health status of P. mariana seedlings. As reflected by a higher number of different OTUs, healthy P. mariana seedlings showed a higher rhizosphere-associated bacterial diversity than diseased seedlings, while the opposite trend was observed for fungi. The high frequency of some fungal OTUs suggested that very few fungal species dominated the rhizosphere of P. mariana seedlings (mainly organisms closely related to Ossicaulis lignitalis, Thelephora sp., and Tricholoma myomyces), while for bacteria, no OTU was found at high frequency. Among the most important differences observed between the bacterial and fungal communities associated with the rhizosphere of healthy and diseased seedlings, clones related to the Nitrospirae, Verrucomicrobia, Dothideomycetes, mitosporic Ascomycota, and Heterobasidiomycetes were found associated only with some diseased seedlings, while clones related to the Deltaproteobacteria, Eurothiomycetes, and Pezizomycetes were found associated only with healthy seedlings. Although it remains difficult to interpret these differences ecologically and functionally due to the important biodiversity detected, these results, which were clearly correlated with the health status of P. mariana seedlings, certainly are of interest and may help to narrow searches using PCR primer sets targeting smaller bacterial and fungal groups.
Very few culture-independent studies have looked at the differences that may exist between the rhizosphere microbial communities associated with healthy and diseased plants and how these changes may influence disease development and/or resistance (22, 39). Yang et al. (39) detected significant changes in bacterial communities associated with healthy and Phytophthora-infected avocado roots, a conclusion also reached by McSpadden Gardener et al. (22) while investigating changes in populations of rhizosphere bacteria associated with take-all disease of wheat caused by Gaeumannomyces graminis var. tritici. Although it was hypothesized in both studies that Pseudomonas populations might have been involved in biological control against those two different fungal soilborne pathogens, clear evidence supporting this hypothesis was not found.
The bacterial and fungal sequences amplified from the rhizosphere of P. mariana belonged to different phylogenetic groups that have also been reported in other culture-independent studies as important rhizosphere-inhabiting organisms. For example, Smalla et al. (29) reported that gram-positive bacteria with a high G+C content (or Actinobacteria) were among the most dominant bacteria in the rhizosphere of strawberry (Fragaria ananassa), oilseed rape (Brassica napus), and potato (Solanum tuberosum) grown under field conditions. McCaig et al. (21) reported that the most important bacterial groups found in a grassland dominated by grass species, such as Agrostis capillaries, Festuca ovina, Trifolium repens, and Trifolium perenne, belonged to the Proteobacteria (mainly the
-group) and Actinobacteria, while Chow et al. (5) reported that the Proteobacteria (mainly the
- and ß-groups) and Acidobacteria dominated the rhizosphere of lodgepole pine (Pinus contorta) under forest conditions. The main bacterial groups identified in the present study, mainly the Proteobacteria, Acidobacteria, and Actinobacteria, were the same as those identified in other studies, even if the plant species and the ecosystems studied were completely different. It is interesting that members of different bacterial groups for which no or only a very limited number of species have been successfully cultured so far, such as the Verrucomicrobia, Acidobacteria, candidate division BD, candidate division TM7, and environmental clones, were detected at high levels (up to 44% of the bacterial rRNA gene sequences amplified) in this study. Because different molecular markers and methodological approaches have been used so far to detect and identify fungi and because very few culture-independent studies have been performed in this field (19, 27, 30), it is not yet possible to make reliable comparisons and obtain a general culture-independent portrait of the most dominant fungal groups commonly associated with the rhizosphere of plants. In the present study, Homobasidiomycetes, Sordaryomycetes, Dothideomycetes, and mitosporic Ascomycota, all belonging to the Ascomycota or Basidiomycota, were the groups most often found associated with the rhizosphere of P. mariana seedlings.
The determination of the causal agent(s) was not the objective of this study and since only the rhizosphere soil and not the root system was investigated for the presence of microorganisms, we cannot yet clearly conclude on which main root rot-causing pathogen(s) was involved. In light of the results obtained, a better hypothesis would be that organisms related to Fusarium (Nectria ochroleuca, 99% similarity) and Cylindrocladium (Cylindrocladium floridanum, 99% similarity) might have contributed to the disease symptoms clearly observed in P. mariana seedlings with symptoms of root rot, a hypothesis also supported by the fact that these organisms were found only in the rhizosphere of diseased seedling. It is interesting that other microbial members that are generally known for their beneficial impact on plants and biological control potential against soilborne plant pathogens (35), such as species related to Pseudomonas (99% similarity with Pseudomonas sp.), Bacillus (98% similarity with Bacillus sp.), and Paenibacillus (99% similarity with Paenibacillus alginolyticus), were detected only in the rhizosphere of healthy seedlings. However, because the specificity range of the PCR primers used was large and only a limited number of samples was studied, these results should mainly be considered as interesting avenues for further investigations, where selected group-specific PCR primers could be used to reveal a more precise taxonomical identification.
As most studies performed so far to identify potential biocontrol agents against given soilborne plant pathogens were conducted using isolation techniques and culture-based approaches, only a few potential rhizosphere competent biocontrol agents have been sampled (18). Hence, it appears reasonable to believe that the discovery of new rhizosphere competent biocontrol agents will most likely benefit from culture-independent studies and that we might eventually identify yet uncultured microorganisms that may play a significant role in biological control of soilborne pathogens. Better knowledge of these organisms may also promote the development of new approaches to culture them and/or manage their population.
We have not yet been able to clearly demonstrate whether the communities associated with healthy roots in fact promote disease suppression or whether their presence is simply a direct consequence of the absence of the pathogen. Nevertheless, the clear demonstration that the communities in the rhizosphere of diseased seedlings are different (as seen by the distinct rRNA gene amplification profiles obtained in this study) certainly represents a good starting point for further analyses. The next steps will involve determining the ecological role of microbial communities associated with selected healthy seedlings and clearly demonstrating the mode of action of biological control.
We are indebted to Vladimir Vujanovic for providing biological material and expertise. We thank Etienne Yergeau for technical assistance.
Present address: Department of Biology, University of Moncton, Moncton, New Brunswick, Canada E1A 3E9. ![]()
|
|
|---|
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»