Previous Article | Next Article 
Applied and Environmental Microbiology, June 1999, p. 2614-2621, Vol. 65, No. 6
0099-2240/99/$04.00+0
Copyright © 1999, American Society for Microbiology. All rights reserved.
Analysis of Fungal Diversity in the Wheat Rhizosphere by
Sequencing of Cloned PCR-Amplified Genes Encoding 18S rRNA and
Temperature Gradient Gel Electrophoresis
Eric
Smit,1,*
Paula
Leeflang,1
Boet
Glandorf,2
Jan
Dirk van
Elsas,3 and
Karel
Wernars1
Microbiological Laboratory for Health
Protection, National Institute of Public Health and the Environment
(RIVM), NL-3720 BA Bilthoven,1 Section
of Plant Pathology, Department of Plant Ecology and Evolutionary
Biology, Utrecht University, NL-3508 TB
Utrecht,2 and Research Institute for
Plant Protection (IPO-DLO), NL-6700
Wageningen,3 The Netherlands
Received 3 December 1998/Accepted 19 March 1999
 |
ABSTRACT |
Like bacteria, fungi play an important role in the soil ecosystem.
As only a small fraction of the fungi present in soil can be cultured,
conventional microbiological techniques yield only limited information
on the composition and dynamics of fungal communities in soil.
DNA-based methods do not depend on the culturability of microorganisms,
and therefore they offer an attractive alternative for the study of
complex fungal community structures. For this purpose, we designed
various PCR primers that allow the specific amplification of fungal
18S-ribosomal-DNA (rDNA) sequences, even in the presence of nonfungal
18S rDNA. DNA was extracted from the wheat rhizosphere, and 18S rDNA
gene banks were constructed in Escherichia coli by cloning
PCR products generated with primer pairs EF4-EF3 (1.4 kb) and EF4-fung5
(0.5 kb). Fragments of 0.5 kb from the cloned inserts were sequenced
and compared to known rDNA sequences. Sequences from all major fungal
taxa were amplified by using both primer pairs. As predicted by
computer analysis, primer pair EF4-EF3 appeared slightly biased to
amplify Basidiomycota and Zygomycota, whereas
EF4-fung5 amplified mainly Ascomycota. The 61 clones that
were sequenced matched the sequences of 24 different species in the
Ribosomal Database Project (RDP) database. Similarity values ranged
from 0.676 to 1. Temperature gradient gel electrophoresis (TGGE)
analysis of the fungal community in the wheat rhizosphere of a
microcosm experiment was carried out after amplification of total DNA
with both primer pairs. This resulted in reproducible, distinctive
fingerprints, confirming the difference in amplification specificity.
Clear banding patterns were obtained with soil and rhizosphere samples
by using both primer sets in combination. By comparing the
electrophoretic mobility of community fingerprint bands to that of the
bands obtained with separate clones, some could be tentatively
identified. While 18S-rDNA sequences do not always provide the
taxonomic resolution to identify fungal species and strains, they do
provide information on the diversity and dynamics of groups of related
species in environmental samples with sufficient resolution to produce
discrete bands which can be separated by TGGE. This combination of
18S-rDNA PCR amplification and TGGE community analysis should allow
study of the diversity, composition, and dynamics of the fungal
community in bulk soil and in the rhizosphere.
 |
INTRODUCTION |
Fungi play an important role in the
soil ecosystem as major decomposers of plant residues, releasing
nutrients that sustain and stimulate plant growth in the process. Some
fungi possess properties antagonistic towards plant pathogens
(15). In the study of soil fungi, particular attention is
given to the rhizosphere; a well-developed and diverse rhizosphere
community is thought to play a role in the suppression of pathogens
(1, 11). Knowledge of the structure and diversity of the
fungal community in the rhizosphere will lead to a better understanding
of pathogen-antagonist interactions. In this work, the development and
application of molecular techniques for assessing fungal diversity in
the rhizosphere are described. As only a certain fraction of the fungi
in soil can be cultured, molecular methods are expected to give a more realistic view of species richness and distribution. Several fungal taxa such as the saprophytic basidiomycetes and the arbuscular endomycorrhiza, which belong to the Glomales, are difficult
or impossible to isolate from soil by dilution plating (26).
Moreover, investigating fungi by estimating their numbers as CFU can be misleading, since most colonies on plates stem from fungal spores (5).
To circumvent the cultivation problem, an array of molecular
techniques, such as amplified ribosomal DNA (rDNA) sequencing, amplified rDNA restriction analysis, and temperature and denaturing gradient gel electrophoreses (TGGE and DGGE) of rDNA, has been applied
to elucidate microbial population structures in the environment (3, 8, 9, 13, 14, 17, 20, 23, 27, 29). Application of these
molecular methods has led to a tremendous increase in knowledge of
microbial ecology and has revealed the existence of formerly unknown
microorganisms (14).
Several molecular techniques have been applied to study medically
important fungi and phytopathogenic strains (22, 28). These
methods are generally not suited for soil studies, since the primers
described are not specific for fungi and some of the systems are
designed to type strains or cultivars, and therefore yield taxonomic
resolutions that are too high. The taxonomic resolution of 18S rDNA
might not always be sufficient to identify fungal species and strains.
However, amplified 18S rDNA will provide information on the fungal
diversity and dynamics of related species of fungi in soil by producing
discrete bands on TGGE gels. Among the scarce reports on this subject
is that of Kowalchuk et al. (13), who obtained interesting
results by using PCR and DGGE for the analysis of fungi associated with
Marram grass in coastal dunes. These investigators used the primers of
Sogin (24) and White et al. (30) and reported the
coamplification of plant rDNA with these primer sets (13).
However, since these primers also amplify 18S rDNA from plants, algae,
and nematodes, they are not suited for direct and specific
amplification of fungal 18S rDNA. Bock et al. (2) also
reported specificity problems when using these primers for the
detection of pathogenic fungi in clinical specimens. When TGGE banding
patterns are to be used for studying fungal communities, it is
essential to know that all bands are of fungal origin. In complex
ecosystems, such as agricultural soil, other 18S-rDNA (i.e.,
eukaryotic) sequences could obscure the fungal banding pattern.
The objective of this study was to develop primers for the specific
amplification of fungal 18S rDNA and to investigate the performance of
these primers in the characterization of fungal communities. An
important objective was to achieve an amplification range with a broad
coverage of all fungal taxa without losing specificity.
A computer-aided prediction of primer specificity was made, based on
the fungal 18S-rDNA sequences in the Ribosomal Database Project (RDP)
database, and then the range of species that could be amplified by the
primers was tested with a collection of different fungal isolates from
all major taxa. Subsequently, the primer sets were used to generate PCR
products from bulk and rhizosphere microcosm soil samples. Common
0.5-kb fragments of cloned 18S rDNA obtained from both soil samples
with both primer sets were sequenced, these sequences were matched to
those in the database, and the specificity and bias of the two primer
systems were assessed. TGGE analysis of the bulk and rhizosphere
samples was done to assess differences in fungal community structure.
 |
MATERIALS AND METHODS |
Primer design and computation of primer pair specificity.
In
order to develop primers for the specific amplification of fungi, a
number of aligned 18S-rDNA sequences of fungi and other eukaryotes were
retrieved from the RDP database (16). Several primers were
chosen to allow the amplification of 18S-rDNA sequences from a wide
range of fungal taxa (Fig. 1): EF3
(5'-TCCTCTAAATGACCAAGTTTG-3'), EF4
(5'-GGAAGGG[G/A]TGTATTTATTAG-3'), and fung5
(5'-GTAAAAGTCCTGGTTCCCC-3'). The EF4-EF3 primer set was
chosen to amplify a major part of the 18S rDNA, and primer set
EF4-fung5 was chosen to amplify an approximately 550-bp fragment. Both
primer sets were developed for the amplification of fungal 18S rDNA
directly from the soil for cloning and sequencing. The primer sequence
NS3 was previously described (30), and a G-C tail (NS3-GC)
was added for TGGE (20). NS3-GC can be combined with EF4 to
amplify a DNA fragment that can then be analyzed by TGGE. NS3 is,
however, not specific for fungi and can only be used in a nested
amplification approach with the specific primer pairs.

View larger version (38K):
[in this window]
[in a new window]
|
FIG. 1.
Primers for amplification of fungal 18S-rDNA sequences.
Primer pair EF4-EF3 (fragment PCR1a) is specific for fungal 18S rDNA,
and primer pair EF4-fung5 (fragment PCR1b) amplifies a smaller, 550-bp
fragment specific for fungi. Primer pair EF4-NS3 (PCR2) can be used
directly on clones or in a nested approach on PCR1 products for
community analysis by TGGE. The primers which were developed in this
work are given in alignment with 18S-rDNA sequences of species from all
major fungal taxa.
|
|
To obtain a computer-aided evaluation of the amplification range of the
primers, each sequence was compared to 436 aligned eukaryotic 18S-rDNA
sequences with the Check Probe program of the RDP database. Fungal
sequences with full homology or one mismatch for each primer were
selected from the analysis listing and are shown with relation to
particular taxa in Table 1.
View this table:
[in this window]
[in a new window]
|
TABLE 1.
Computer analysis giving the percentage of database
sequences showing fewer than two mismatches in both primers of a
given pair
|
|
Fungal strains, culture conditions, and amplification range.
Various fungal species from all major taxa were used to test the primer
sets (Table 2). Lyophilized cultures were
rehydrated with a sterile 0.85% NaCl solution and applied to potato
dextrose agar plates. Plates were incubated at 28°C for several days
to 2 weeks, until sufficient hyphal growth was observed.
Fungi were also cultured from the soil. Samples of 10 g of soil
were shaken in 100 ml of a MgSO
4 solution (10 mM) for 15 min.
Serial dilutions were plated onto 2% malt extract agar containing
0.33% Solacol and 200 ppm of aureomycine (
7). Plates were
incubated
at 20°C for 1 week. Fungi with different morphologies were
then
selected and streaked onto cornmeal agar. The plates were
incubated
for 1 to 2 weeks, and some of the fungi were identified by W.
Gams at the Centraal Bureau voor Schimmelcultures, Baarn, The
Netherlands. The species that were found are listed in Table
2.
DNA was isolated by scraping hyphae from the agar surface. Cells were
disrupted, DNA was released by bead beating, and the
resulting lysate
was purified as described by Smit et al. (
23).
Prior to the examination of the amplification range of the primer sets
EF4-EF3, EF4-fung5, and EF4-NS3-GC on a collection
of fungal isolates,
fungal cell lysis and the quality of the extracted
DNA were tested by
using the general eukaryotic primer set 106-107
(
24). Only 1 strain (
Coniothyrium sporulosum) out of 27 was
poorly lysed
(results not shown). A positive signal obtained by
using these primers
indicated that the fungal cells were lysed
by the bead-beating method
and that the quality of the DNA was
sufficiently high for PCR
amplification.
Soil microcosms and DNA extraction.
To study the fungal
community in wheat rhizosphere soil, we collected soil from a small
field plot on the campus of the University of Utrecht located on the
Uithof, Utrecht, The Netherlands. This is a clay soil containing 4%
organic matter with a pH of 5.0. Samples from this soil were used for
plating culturable fungi and for setting up the microcosm experiment.
The soil was air dried and sieved, and nine small pots (diameter, 10 cm) were filled. Soil was seeded with eight seeds of Triticum
aestivum cv. Baldus per pot. Fluctuations in moisture content were
minimized by supplying water daily to keep the soil moisture content at
20%. Microcosms were incubated in a climate chamber with a light-dark
regimen of 16 and 8 h at 20 and 15°C, respectively. Microcosms
were sampled in duplicate on days 5 and 10. Bulk soil samples of 3 g were taken from root-free soil. Rhizosphere soil was obtained by
gently shaking the soil from the roots, and roots with adhering soil
were added to 50-ml polypropylene tubes with 10 ml of sterile sodium
phosphate buffer (120 mM; pH 8) and 1 g of gravel. Tubes were
vortexed for 30 s, and the buffer-soil slurry mixture was poured
into a new tube, leaving the gravel and roots behind. Total DNA was
extracted from the rhizosphere soil slurry by using a bead beater
(23). One microliter of the extract was used for PCR amplification.
The extent to which this method lyses fungi present in the soil is not
known. However, bead beating has been shown to lyse
bacterial spores
(
18), and all cultured fungal strains in this
study could be
disrupted by bead beating. Moreover, the Braun
bead beater also lyses
Saccharomyces cerevisiae cells under the
conditions
used.
PCR amplification.
Total DNA extracts from duplicate soil
samples were pooled before PCR. Primer sets EF4-EF3 and EF4-fung5 were
used for direct amplification of 18S-rDNA sequences from extracted
wheat rhizosphere DNA.
18S-rDNA clones amplified with both primer sets were also analyzed by
TGGE. For this purpose, they were reamplified from
Escherichia coli colony material. PCR was performed by using primer set
EF4-NS3-GC
to generate 18S-rDNA fragments suitable for analysis by
TGGE.
In order to generate TGGE profiles of the soil fungal community,
a nested approach was necessary. An initial PCR was done with
either
EF4-EF3 or EF4-fung5, as described above. Subsequently,
the reaction
mixture was diluted 1:500, and a second amplification
round was
performed with EF4-NS3-GC (Fig.
1). With this approach,
corresponding
fragments can be analyzed by TGGE, although different
primer sets had
been initially used for the specific amplification
of the fungal
community in the wheat rhizosphere. Differences
in banding patterns
will reflect the differences in specificity
of the initial PCR (Fig.
2). 18S-rDNA fragments were obtained
from
cloned material (see below) for TGGE analysis by the direct
use of
EF4-NS3-GC on a small amount of a specific
E. coli colony
as template source.

View larger version (92K):
[in this window]
[in a new window]
|
FIG. 2.
TGGE of amplified 18S-rDNA fragments representing the
fungal community in bulk and rhizophere soil samples of the microcosm
experiment. Banding patterns were obtained by mixing the duplicate
samples before PCR and adding PCR products from both primer pairs
(EF4-EF3 and EF4-fung5) the same lane. Lane 1, day 5 bulk soil; lane 2, day 5 rhizosphere soil; lane 3, day 10 bulk soil; and lane 4, day 10 rhizosphere soil. Numbered bands are explained in the text.
|
|
PCR mixtures for all three primer sets consisted of 5 µl of 10× PCR
buffer (Boehringer), 1.7 mM MgCl
2, 200 µM concentrations
of each deoxynucleoside triphosphate, 300 nM concentrations of
each
primer, 2.6 U of Expand enzyme mix (Boehringer), 1 µl of
1:10-diluted
template DNA, and sterile Millipore water to a final
volume of 50 µl.
The following thermocycling pattern was used:
94°C for 3 min (1 cycle); 94°C for 1 min, 48°C for 1 min, and
72°C for 3 min (40 cycles); and 72°C for 10 min (1
cycle).
Cloning of EF4-EF3 and EF4-fung5 PCR products.
Only PCR
products from the pooled, day 10 rhizosphere samples were cloned and
sequenced. PCR products from the soil microbial community were
separated on a 1% agarose-Tris-borate-EDTA gel. Bands were excised,
and the DNA was purified by centrifuging the gel pieces for 15 min at
16,000 × g in a Wizard column without resin. The
flowthrough containing the DNA was collected and used without further
purification. DNA fragments were then ligated into the pGEM-T vector,
which has 3' T overhangs to facilitate the cloning of PCR products
(Promega, Madison, Wis.). Ligation mixes were used to transform
ultracompetent E. coli XL1-Blue cells (Stratagene,
Cambridge, United Kingdom) according to the manufacturer's instructions. White colonies were selected, cultured in 2 ml of Luria
broth, and stored at
70°C. Subsequently, a quick clone screening
procedure based on PCR amplification was performed directly from colony
material. Clones obtained from the rhizosphere soil sample with both
primer sets and containing inserts of the correct size were selected
for sequencing.
TGGE analysis of the fungal rhizosphere community and the
18S-rDNA clones.
The products amplified from the fungal community
of the microcosm experiment and from the clones were analyzed by TGGE.
For TGGE analysis, a Diagen electrophoresis unit was used according to
the instructions of the manufacturer. The polyacrylamide gel was
composed of 8% (wt/vol) acrylamide, 0.21% (wt/vol) bisacrylamide, 8 M
urea, 20 mM MOPS (morpholinepropanesulfonic acid), 1 mM EDTA (pH 8),
and 20% (vol/vol) formamide. The gel was polymerized onto a gel
support film (Gelbond PAG; FMC). Electrophoresis was performed at a
constant of 110 V for 17 h with a temperature gradient from 36 to
44°C.
Gels were silver stained with the PlusOne DNA silver staining kit
(Pharmacia Biotech) in a Hoefer automated gel stainer (Pharmacia
Biotech) according to the standard DNA-staining protocol with
increased
solution volumes (total volume, 187.5 ml), a developing
time of 8 min,
an extra washing step with H
2O, and a
Na
2CO
3 fixing
step.
Stained gels were air dried, and the images were analyzed with Bioprint
(version 96.11) and Biogene (version 96.15) software
packages (Vilber
Lourmat, Marne-la-Vallee, France). Banding patterns
were analyzed by
cluster analysis using the unweighted pair group
method, and genetic
similarity was calculated according to the
method of Nei and Li
(
6a) [
a = 2
nxy/(
nx+
ny)] by using Biogene
software (Vilbert
Lourmat).
Partial sequencing of the cloned fungal 18S-rDNA fragments.
A 0.5-kb fragment of each clone obtained from wheat rhizosphere soil
was sequenced in both directions on an ABI377 DNA sequencer by cycle
sequencing using the dye terminator system (Eurogentec, Seraing,
Belgium). EF4-fung5 clones were sequenced with the SP6 and T7 primers.
The corresponding region of EF4-EF3 clones, which contains a much
longer insert, was sequenced with primers EF4 and FS7
(5'-GCTTTGAACACTCTAATT-3'). Consensus sequences were
obtained by using the DNASIS software package (version 2.5; Hitachi
Software Engineering America, Ltd., San Francisco, Calif.). Sequences
were then manually edited and analyzed with the Sequence Match (version 2.7) program of the RDP database (containing 10,723 unaligned ribosomal
sequences) and the Blast 2 Advanced Blast program at the National
Center for Biotechnology Information (NCBI) (containing 351,020 sequences). Fungal species sequences in the database most similar to
those of the clones are listed in Tables
3 and 4. To
check for occurrences of presumptively chimeric sequences, all clones
were analyzed by using the Chimera Check program of the RDP database
(version 2.7). Since it is impossible to unambiguously determine
whether a sequence is chimeric, all clones were designated either N
(not likely to be chimeric) or P (possibly chimeric) (Tables 3 and 4).
Nucleotide sequence accession numbers.
All nucleotide
sequences were submitted to NCBI and assigned accession no. AF095648 to
AF095686 and AF096351 to AF096369.
 |
RESULTS AND DISCUSSION |
Primer design, computation of specificity, and analysis of
amplification range.
The heterogeneity of fungal taxa and the
homology of fungi to other eukaryotes made it unfeasible to develop a
single specific primer set with both a satisfactory fungal
amplification range and a sufficiently low 18S homology to other
eukaryotes. The target sequences of the primers shown in Fig. 1
demonstrate that the EF4-EF3 primer set exhibited only an occasional
mismatch with some fungi. From the computer comparison of EF4-EF3 to
the RDP database, in which one mismatch was allowed for each primer, it was predicted that 78% of all fungal 18S-rDNA sequences in the aligned
database would be amplified (Table 1). This is only an estimation,
since it is not known how many mismatches can be tolerated under the
given PCR conditions. Primer set EF4-EF3 exhibits a relatively high
coverage of the Basidiomycota (95 to 100%) compared to
EF4-fung5. The latter primer pair shows a relatively high coverage of
the Ascomycota and a lower coverage of the
Basidiomycota. Neither set of primers has full coverage of
the Zygomycota; EF4-EF3 will theoretically amplify 75% of
the species of this taxon, while EF4-fung5 will amplify only 40%.
According to the computer analysis, both sets of primers for the
analysis of fungal communities should adequately amplify most species
of the different fungal taxa, except for the
Chytridiomycota, which are regarded as having little importance for soil ecosystems.
The results in Table
2 demonstrate that all fungi studied, except
C. sporulosum,
Candida lipolytica, and
Aureobasidium pullulans,
could be amplified by using primer
set EF4-EF3. Although
A. pullulans and
C. sporulosum could be amplified with EF4-fung5 and EF4-NS3-GC,
little difference was found between the amplification ranges of
EF4-EF3
and EF4-fung5 for the
Ascomycota. Quite a difference can
be
observed for the
Zygomycota, since EF4-fung5 amplified only
two of the four species amplified by EF4-EF3. The four species
of
Basidiomycota tested could be amplified by all three primer
sets, except for
Ustilago avenae, which failed to be
amplified
by EF4-fung5. Although the results of the computer prediction
(Table
1) are based on the assumption that only one mismatch
per primer
is allowed, the amplification coverage of the isolates
(Table
2) seems
to show the same trend. However, the amplification
of 18S rDNA from
fungal isolates does not necessarily predict
the performance of the
primer sets in environmental samples which
harbor complex
communities.
Analysis of the EF4-fung5 clones from wheat rhizosphere soil.
About 70 clones of the 18S-rDNA sequences amplified with EF4-fung5
obtained from the wheat rhizosphere sample were checked for the
presence of an insert of the expected size (approximately 0.5 kb) by
reamplifying the fragments directly from colony material. Thirty-nine
E. coli clones with the inserts of the correct size were
further analyzed by sequencing.
Sequences were analyzed by using the Blast 2.0 (advanced) program at
the NCBI and by applying the Sequence Match program of
the RDP. Both
programs yielded identical results for almost all
clones. In Table
3,
the species of fungus with 18S-rDNA sequences
most similar to the clone
are listed. In cases where different
results were obtained, both
alignments were manually checked and
the most likely species were
chosen. Among the 39 sequences, 17
different species were identified.
Some nonidentical cloned sequences
matched the sequences of identical
species. This is probably due
to the incompleteness of the database,
although PCR-related amplification
errors cannot be excluded. As shown
in Table
3, most clones detected
by PCR with EF4-fung5 belong to the
euascomycetes of the
Ascomycota.
No
Basidiomycota
were detected, although the computer-based analysis
and isolate tests
(Table
2) indicate that EF4-fung5 could amplify
some members of the
basidiomycetes. Besides PCR bias, cloning
bias could also play a role.
One sequence of the
Zygomycota was
detected (UUF64), and a
sequence of the
Chytridiomycota was found,
although the
computer analysis predicts that only 14% of the species
of this taxon
can be amplified. Only one clone, UUF31, completely
matched a sequence
in the database (
Eupenicillium javanicum).
All other clones
exhibited similarities ranging from 0.676 to
0.978 to various database
sequences. Clones with a similarity
to a fungal species below 0.95 can
be expected to have originated
from a fungal species only
phylogenetically related to the given
name. This suggests that most of
the species names are only very
approximate indications of the
identities of the cloned sequences.
Clones with low similarities
probably originated from species
distantly related to those present in
the database and whose 18S-rDNA
sequences have not been determined, or
they are from fungi yet
to be isolated. For instance, several clones
matched
Neocosmospora vasinfecta, which is not common in the
soil of temporate regions
(
7a). Since
N. vasinfecta has 18S sequences similar to those
of
Fusarium
solani, the sequences detected could represent
Fusarium species from which no 18S sequences are present in the database.
Nevertheless, it is of importance for judging the primer performance
that genera matching clones of such organisms as
Cladosporium,
Eupenicillium,
Exophiala, and
Pleospora and the potential plant
pathogens
Alternaria and
Verticillium are common
in soil (
4,
6,
10,
11,
19). Both
Spizellomyces-
and
Mycosphaerella-like
sequences (as in clones UUF54 and
UUF68) were detected by Kowalchuk
et al. (
13) in association
with Marram grass
roots.
Analysis of the EF4-EF3 clones from wheat rhizosphere soil.
The clones obtained by PCR with EF4-EF3 are depicted in Table 4. Of 50 clones, 26 were found to have inserts of the expected size, and out of
these 26 clones, only 22 could be sequenced with the use of the
internal primer FS7. Nine different species were identified among the
twenty-two clones. This primer set obviously is biased to amplifying
Basidiomycota, since most of the species that were detected
belong to this taxon. Clones resembling the yeast-like symbionts of the
genus Symbiotaphrina were detected by using both primer
sets. The relatively low similarities suggest that the cloned 18S rDNA
could represent distant euascomycete relatives of these symbionts
(21). A large number of clone sequences matched the fungus
Mortierella polycephala. Obviously these sequences were
related only to those of M. polycephala since this species is not a common soil inhabitant, and it might be more likely that the
sequences represent fungi such as Mortierella alpina or
Mortierella elongata, whose 18S sequences are not present in
the database (7a). A thorough phylogenetic analysis of the
clones with sequences of related species could give more information;
this is, however, beyond the scope of the present report.
Examining the fungal diversity of the wheat rhizosphere with both
primer systems revealed that 24 different species were detected
among
61 clones. Clone diversity is much greater since no clones
were found
with similar sequences, and the number of species in
the database
limits the potential matches that can be found. The
high sequence
diversity and the broad taxonomic range of the clones
are indicative of
the importance of using both primer sets in
a complementary fashion in
studies of rhizosphere fungal
communities.
TGGE analysis of the fungal community and clones from a wheat
rhizosphere.
TGGE analysis of amplified fungal 18S-rDNA fragments
from wheat rhizosphere samples was done by using a nested PCR approach in two amplification rounds, the first with EF4-EF3 or EF4-fung5 and
the second with EF4-NS3-GC. Analysis of the fungal community profiles
from soil samples of the microcosm experiment by using either EF4-EF3
or EF4-fung5 in the first selective amplification round produced
different banding patterns when the different primer sets were used
(not shown). This result supported the biases of both primer systems
which were found by using the clone libraries. This confirms that, in
order to visualize as many bands as possible, both primer systems
should be used in combination to study fungal diversity in soil.
Repeated analyses of the TGGE patterns of duplicate soil samples, with
either primer pair, showed little variation, which suggests good
reproducibility of the DNA extraction procedure, the PCR amplification,
and the TGGE analysis.
The soil samples from the microcosm experiment could be adequately
analyzed by mixing the PCR products of both primer systems
before
electrophoresis. TGGE analysis of mixed, duplicate samples
revealed
differences between both the day 5 and 10 fungal communities
and the
bulk and rhizosphere samples (Fig.
2 and
3). Banding patterns
of the bulk and
rhizosphere samples at days 5 and 10 revealed
considerable differences:
in the bulk soil sample at day 5, 24
bands were present, while in the
rhizosphere sample, 20 bands
could be visualized. On day 10, 26 bands
could be produced from
the bulk sample, and 25 bands from the
rhizosphere sample could
be visualized. A dendrogram representing the
differences between
the profiles indicates 60% similarity between the
bulk and rhizosphere
samples at day 5 and 75% at day 10 and a
similarity of 55% between
day 5 and 10 profiles (Fig.
3). Results
suggest that both the
presence of plant roots and the time of sampling
of the microcosms
affects fungal diversity and composition in soil.
Considering
the number of bands, diversity seems to be somewhat lower
in the
rhizosphere than in bulk soil.

View larger version (13K):
[in this window]
[in a new window]
|
FIG. 3.
A dendrogram based on the presence or absence of bands
(Fig. 2) was constructed to represent the percentages of genetic
identity between the profiles of the microcosm samples of bulk and
rhizosphere (rhizo) soils.
|
|
By comparing the single-clone bands with those of the community
profiles, the presumptive identities of some of the bands
in the
community profile were obtained (Fig.
2). It should be
stated that this
method provides only circumstantial evidence
of identity. For more
exact identification, bands should be excised
and sequenced. Some
bands, such as band 1, which might be UUF10,
band 3, which could be
UUF30-31, and band 5, of unknown identity
(Fig.
2), appear in all
sample profiles (Fig.
2). Clones UUF30
and UUF31 have high similarities
to
E. javanicum, which is common
in both soil and
rhizosphere (
6). Other bands appear in specific
samples,
such as band 2, which is of high intensity in day 10
samples, and band
4, which matches clone UUF103, which is mainly
found in day 10 samples.
Band 7, of unknown origin, is found only
in day 5 samples. The
group of bands around the position of band
6 matched some bands of the
clones with
Mortierella-like sequences,
which appeared to
have different mobilities (not
shown).
Primer performance and resolution.
The results of this study
show that primer sets enabled us to detect sequences of all major
fungal taxa, except the hemiascomycetes (yeasts) and archaeascomycetes
in wheat rhizosphere soil samples. Analysis of the fungal community by
TGGE allows the study of the differences between samples and the
dynamics of certain species. Comparison of cultured isolates from soil
to 18S-rDNA sequences was not feasible, since this study did not
thoroughly isolate fungi. Nevertheless, any overlap in species
identification due to the use of both methods is expected to be small,
since plate count techniques favor the isolation of fast-growing,
low-substrate-specific, spore-producing fungi (25), while
molecular techniques probably favor numerically dominant fungi with
relatively high amounts of vegetative mycelium. Moreover, the sequences
should be efficiently amplified and cloned in order to be detected.
Fusaria are often found in soil and especially in the wheat rhizosphere
by plate techniques (11, 12). One striking observation of
this study is that no Fusarium species were identified. This
is in line with the work of Kowalchuk et al. (13), who
studied Marram grass roots and also did not detect fusaria with their
PCR method. However, no matches can be found, since a screening of the
database showed that no 18S-rDNA sequences from Fusarium
species were present (only 28S and internal transcribed spacer
sequences were present in the database). This clearly illustrates that
a shortcoming of the use of 18S rDNA to detect and identify fungi lies
in the incompleteness of existing databases and the taxonomic
resolution of the sequence. For certain fungi, either the ITS region or
the 28S rDNA gives a higher resolution which enables discrimination between closely related species. At the moment, a major disadvantage of
using these sequences to develop primers for all fungi in soil is the
relatively low number of sequences present in the 28S-rDNA database.
However, for the amplification of specific fungal groups, such as
Fusarium, development of primers for 28S-rDNA sequences will
probably be much more suitable and feasible.
The combination of cloning and sequencing with whole-community TGGE
analysis, although still in its infancy for fungi, has
been shown to be
a powerful technique for elucidating the differences
and dynamics of
the fungal community in the
rhizosphere.
 |
ACKNOWLEDGMENTS |
This investigation was performed for the Dutch Ministry of
Housing, Spatial Planning and Environment, Directorate General for
Environmental Protection.
We are indebted to K. Smalla and coworkers for their help and
instructions with the TGGE technique. We also thank W. Gams and G. Kowalchuk for critically reading the manuscript and for their valuable advice.
 |
FOOTNOTES |
*
Corresponding author. Mailing address: Microbiological
Laboratory for Health Protection, National Institute of Public Health and the Environment, P.O. Box 1, NL-3720 BA Bilthoven, The Netherlands. Phone: 31 30 2743924. Fax: 31 30 2744434. E-mail:
Eric.Smit{at}rivm.nl.
 |
REFERENCES |
| 1.
|
Alabouvette, C.
1990.
Biological control of Fusarium wilt pathogens in suppressive soils, p. 27-34.
In
D. Hornby (ed.), Biological control of soil-borne pathogens. CAB International, Wallingford, United Kingdom.
|
| 2.
|
Bock, M.,
M. Maiwald,
R. Kappe, and H. Näher.
1994.
Polymerase chain reaction-based detection of dermatophyte DNA with a fungus-specific primer system.
Mycoses
37:79-84[Medline].
|
| 3.
|
Borneman, J.,
P. W. Skroch,
K. M. O'Sullivan,
J. A. Palus,
N. G. Rumjanek,
J. L. Jansen,
J. Nienhuis, and E. W. Triplett.
1996.
Molecular microbial diversity of an agricultural soil in Wisconsin.
Appl. Environ. Microbiol.
62:1935-1943[Abstract].
|
| 4.
|
Christensen, M.
1981.
Species diversity and dominance in fungal communities, p. 201-232.
In
D. T. Wicklow, and G. C. Carroll (ed.), The fungal community. Marcel Dekker, Inc., New York, N.Y.
|
| 5.
|
De Ley, F. A. A. M., and J. M. Lynch.
1997.
Functional diversity of the rhizosphere, p. 38-43.
In
A. Ogoshi, K. Kobayashi, Y. Homma, F. Kadoma, N. Kondo, and S. Akino (ed.), Plant growth-promoting rhizobacteria. Proceedings of the Fourth International Workshop on plant-growth promoting rhizobacteria. OECD, Paris, France.
|
| 6.
|
Domsch, K. H.,
W. Gams, and T.-H. Anderson.
1980.
Compendium of soil fungi, vol. 1.
Academic Press, London, England.
|
| 6a.
|
Ferguson, A.
1980.
Biochemical systematics and evolution.
Blackie, Glasgow, United Kingdom.
|
| 7.
|
Gams, W., and W. Van Laar.
1982.
The use of solacol (validamycin) as a growth retardant in the isolation of soil fungi.
Neth. J. Plant Pathol.
88:39-45.
|
| 7a.
| Gams, W. Personal communication.
|
| 8.
|
Heuer, H.,
M. Krsek,
P. Baker,
K. Smalla, and E. M. H. Wellington.
1997.
Analysis of actinomycete communities by specific amplification of genes encoding 16S rRNA and gel-electrophoretic separation in denaturing gradients.
Appl. Environ. Microbiol.
63:3233-3241[Abstract].
|
| 9.
|
Hiorns, W. D.,
R. C. Hastings,
I. M. Head,
A. J. McCarthy,
J. R. Saunders,
R. W. Pickup, and G. H. Hall.
1995.
Amplification of 16S ribosomal RNA genes of autotrophic ammonia-oxidizing bacteria demonstrates the ubiquity of Nitrospiras in the environment.
Microbiology
141:2793-2800[Abstract/Free Full Text].
|
| 10.
|
Hurek, T.,
B. Wagner, and B. Reinhold-Hurek.
1997.
Identification of N2-fixing plant- and fungus-associated Azoarcus species by PCR-based genomic fingerprints.
Appl. Environ. Microbiol.
63:4331-4339[Abstract].
|
| 11.
|
Jarosik, V.,
E. Kovacikova, and H. Maslowska.
1996.
The influence of planting location, plant growth stage and cultivars on microflora of winter wheat roots.
Microbiol. Res.
151:177-182.
|
| 12.
|
Khanna, R.,
S. Chandra, and K. K. Khanna.
1993.
Rhizosphere microflora of Triticale.
Biol. Mem.
19:111-121.
|
| 13.
|
Kowalchuk, G. A.,
S. Gerards, and J. W. Woldendorp.
1997.
Detection and characterization of fungal infections of Ammophila arenaria (Marram grass) roots by denaturing gradient gel electrophoresis of specifically amplified 18S rDNA.
Appl. Environ. Microbiol.
63:3858-3865[Abstract].
|
| 14.
|
Liesack, W., and E. Stackebrandt.
1992.
Occurrence of novel groups of the domain Bacteria as revealed by analysis of genetic material isolated from an Australian terrestrial environment.
J. Bacteriol.
174:5072-5078[Abstract/Free Full Text].
|
| 15.
|
Lumsden, R. D.
1981.
Ecology of mycoparasitism, p. 295-328.
In
D. T. Wicklow, and G. C. Carroll (ed.), The fungal community. Marcel Dekker, Inc., New York, N.Y.
|
| 16.
|
Maidak, B. L.,
G. J. Olsen,
N. Larson,
R. Overbeek,
M. J. McCaughey, and C. R. Woese.
1997.
The RDP (Ribosomal Database Project).
Nucleic Acids Res.
25:106-111.
|
| 17.
|
Martinez-Murcia, A. J.,
S. G. Acinas, and F. Rodriguez-Valera.
1995.
Evaluation of prokaryotic diversity by restrictase digestion of 16S rDNA directly amplified from hypersaline environments.
FEMS Microbiol. Ecol.
17:247-256.
|
| 18.
|
Moré, M. I.,
J. B. Herrick,
M. C. Silva,
W. C. Ghiorse, and E. L. Madsen.
1994.
Quantitative cell lysis of indigenous microorganisms and rapid extraction of microbial DNA from sediment.
Appl. Environ. Microbiol.
60:1572-1580[Abstract/Free Full Text].
|
| 19.
|
Mueller-Dombois, D.
1981.
Ecological measurements and microbial populations, p. 173-184.
In
D. T. Wicklow, and G. C. Carroll (ed.), The fungal community. Marcel Dekker, Inc., New York, N.Y.
|
| 20.
|
Muyzer, G.,
E. C. De Waal, and A. G. Uitterlinden.
1993.
Profiling of complex microbial populations by denaturing gradient gel electrophoresis analysis of polymerase chain reaction-amplified genes coding for 16S rRNA.
Appl. Environ. Microbiol.
59:695-700[Abstract/Free Full Text].
|
| 21.
|
Noda, H., and K. Kodama.
1996.
Phylogenetic position of yeastlike endosymbionts of Anobiid beetles.
Appl. Environ. Microbiol.
62:162-167[Abstract].
|
| 22.
|
Rand, K. H.,
H. Houck, and M. Wolf.
1994.
Detection of candidemia by polymerase chain reaction.
Mol. Cell. Probes
8:215-222[Medline].
|
| 23.
|
Smit, E.,
P. Leeflang, and K. Wernars.
1997.
Detection of shifts in microbial community structure and diversity in soil caused by copper contamination using amplified ribosomal DNA restriction analysis.
FEMS Microbiol. Ecol.
23:249-261.
|
| 24.
|
Sogin, M. L.
1990.
Amplification of ribosomal RNA genes for molecular evolution studies, p. 307-320.
In
M. A. Innes, D. H. Gelfand, J. J. Sninsky, and T. J. White (ed.), PCR protocols. Academic Press, San Diego, Calif.
|
| 25.
|
States, J. S.
1981.
Useful criteria in the description of fungal communities, p. 185-200.
In
D. T. Wicklow, and G. C. Carroll (ed.), The fungal community. Marcel Dekker, Inc., New York, N.Y.
|
| 26.
|
Thorn, R. G.,
C. A. Reddy,
D. Harris, and E. A. Paul.
1996.
Isolation of saprophytic basidiomycetes from soil.
Appl. Environ. Microbiol.
62:4288-4292[Abstract].
|
| 27.
|
Torsvik, V.,
R. Sørheim, and J. Goksør.
1996.
Total bacterial diversity in soil and sediment communities a review.
J. Ind. Microbiol.
17:170-178.
|
| 28.
|
Walsh, T. J.,
A. Francesconi,
M. Kasai, and S. J. Chanock.
1995.
PCR and single-strand conformational polymorphism for recognition of medically important opportunistic fungi.
J. Clin. Microbiol.
33:3216-3220[Abstract].
|
| 29.
|
Weidner, S.,
W. Arnold, and A. Pühler.
1996.
Diversity of uncultured microorganisms associated with the seagrass Halophila stipulacea estimated by restriction fragment length polymorphism analysis of PCR-amplified 16S rRNA genes.
Appl. Environ. Microbiol.
62:766-771[Abstract].
|
| 30.
|
White, T. J.,
T. Burns,
S. Lee, and J. Taylor.
1990.
Amplification and direct sequencing of fungal ribosomal RNA genes for phyologenetics, p. 315-322.
In
M. A. Innes, D. H. Gelfand, J. J. Sninsky, and T. J. White (ed.), PCR protocols. Academic Press, San Diego, Calif.
|
Applied and Environmental Microbiology, June 1999, p. 2614-2621, Vol. 65, No. 6
0099-2240/99/$04.00+0
Copyright © 1999, American Society for Microbiology. All rights reserved.
This article has been cited by other articles:
-
Freeman, K. R., Martin, A. P., Karki, D., Lynch, R. C., Mitter, M. S., Meyer, A. F., Longcore, J. E., Simmons, D. R., Schmidt, S. K.
(2009). Evidence that chytrids dominate fungal communities in high-elevation soils. Proc. Natl. Acad. Sci. USA
106: 18315-18320
[Abstract]
[Full Text]
-
Le Calvez, T., Burgaud, G., Mahe, S., Barbier, G., Vandenkoornhuyse, P.
(2009). Fungal Diversity in Deep-Sea Hydrothermal Ecosystems. Appl. Environ. Microbiol.
75: 6415-6421
[Abstract]
[Full Text]
-
Petrosino, J. F., Highlander, S., Luna, R. A., Gibbs, R. A., Versalovic, J.
(2009). Metagenomic Pyrosequencing and Microbial Identification. Clin. Chem.
55: 856-866
[Abstract]
[Full Text]
-
Gao, Z., Li, B., Zheng, C., Wang, G.
(2008). Molecular Detection of Fungal Communities in the Hawaiian Marine Sponges Suberites zeteki and Mycale armata. Appl. Environ. Microbiol.
74: 6091-6101
[Abstract]
[Full Text]
-
Oliver, K. L., Hamelin, R. C., Hintz, W. E.
(2008). Effects of Transgenic Hybrid Aspen Overexpressing Polyphenol Oxidase on Rhizosphere Diversity. Appl. Environ. Microbiol.
74: 5340-5348
[Abstract]
[Full Text]
-
Pitkaranta, M., Meklin, T., Hyvarinen, A., Paulin, L., Auvinen, P., Nevalainen, A., Rintala, H.
(2008). Analysis of Fungal Flora in Indoor Dust by Ribosomal DNA Sequence Analysis, Quantitative PCR, and Culture. Appl. Environ. Microbiol.
74: 233-244
[Abstract]
[Full Text]
-
Fierer, N., Breitbart, M., Nulton, J., Salamon, P., Lozupone, C., Jones, R., Robeson, M., Edwards, R. A., Felts, B., Rayhawk, S., Knight, R., Rohwer, F., Jackson, R. B.
(2007). Metagenomic and Small-Subunit rRNA Analyses Reveal the Genetic Diversity of Bacteria, Archaea, Fungi, and Viruses in Soil. Appl. Environ. Microbiol.
73: 7059-7066
[Abstract]
[Full Text]
-
Mansfield, M.A., Kuldau, G.A.
(2007). Microbiological and molecular determination of mycobiota in fresh and ensiled maize silage. Mycologia
99: 269-278
[Abstract]
[Full Text]
-
Hunter, P. J., Petch, G. M., Calvo-Bado, L. A., Pettitt, T. R., Parsons, N. R., Morgan, J. A. W., Whipps, J. M.
(2006). Differences in Microbial Activity and Microbial Populations of Peat Associated with Suppression of Damping-Off Disease Caused by Pythium sylvaticum.. Appl. Environ. Microbiol.
72: 6452-6460
[Abstract]
[Full Text]
-
Griffiths, L. J., Anyim, M., Doffman, S. R., Wilks, M., Millar, M. R., Agrawal, S. G.
(2006). Comparison of DNA extraction methods for Aspergillus fumigatus using real-time PCR.. J Med Microbiol
55: 1187-1191
[Abstract]
[Full Text]
-
Paulino, L. C., Tseng, C.-H., Strober, B. E., Blaser, M. J.
(2006). Molecular analysis of fungal microbiota in samples from healthy human skin and psoriatic lesions.. J. Clin. Microbiol.
44: 2933-2941
[Abstract]
[Full Text]
-
Jumpponen, A., Johnson, L. C.
(2005). Can rDNA analyses of diverse fungal communities in soil and roots detect effects of environmental manipulations--a case study from tallgrass prairie.. Mycologia
97: 1177-1194
[Abstract]
[Full Text]
-
Luo, Q., Krumholz, L. R., Najar, F. Z., Peacock, A. D., Roe, B. A., White, D. C., Elshahed, M. S.
(2005). Diversity of the Microeukaryotic Community in Sulfide-Rich Zodletone Spring (Oklahoma). Appl. Environ. Microbiol.
71: 6175-6184
[Abstract]
[Full Text]
-
Berg, G., Zachow, C., Lottmann, J., Gotz, M., Costa, R., Smalla, K.
(2005). Impact of Plant Species and Site on Rhizosphere-Associated Fungi Antagonistic to Verticillium dahliae Kleb.. Appl. Environ. Microbiol.
71: 4203-4213
[Abstract]
[Full Text]
-
Jorgensen, H. B., Johansson, T., Canback, B., Hedlund, K., Tunlid, A.
(2005). Selective foraging of fungi by collembolans in soil. Biol Lett
1: 243-246
[Abstract]
[Full Text]
-
Anastasi, A., Varese, G. C., Filipello Marchisio, V.
(2005). Isolation and identification of fungal communities in compost and vermicompost. Mycologia
97: 33-44
[Abstract]
[Full Text]
-
Filion, M., Hamelin, R. C., Bernier, L., St-Arnaud, M.
(2004). Molecular Profiling of Rhizosphere Microbial Communities Associated with Healthy and Diseased Black Spruce (Picea mariana) Seedlings Grown in a Nursery. Appl. Environ. Microbiol.
70: 3541-3551
[Abstract]
[Full Text]
-
Green, S. J., Freeman, S., Hadar, Y., Minz, D.
(2004). Molecular tools for isolate and community studies of Pyrenomycete fungi. Mycologia
96: 439-451
[Abstract]
[Full Text]
-
Seghers, D., Wittebolle, L., Top, E. M., Verstraete, W., Siciliano, S. D.
(2004). Impact of Agricultural Practices on the Zea mays L. Endophytic Community. Appl. Environ. Microbiol.
70: 1475-1482
[Abstract]
[Full Text]
-
Gomes, N. C. M., Fagbola, O., Costa, R., Rumjanek, N. G., Buchner, A., Mendona-Hagler, L., Smalla, K.
(2003). Dynamics of Fungal Communities in Bulk and Maize Rhizosphere Soil in the Tropics. Appl. Environ. Microbiol.
69: 3758-3766
[Abstract]
[Full Text]
-
Viebahn, M., Glandorf, D. C. M., Ouwens, T. W. M., Smit, E., Leeflang, P., Wernars, K., Thomashow, L. S., van Loon, L. C., Bakker, P. A. H. M.
(2003). Repeated Introduction of Genetically Modified Pseudomonas putida WCS358r without Intensified Effects on the Indigenous Microflora of Field-Grown Wheat. Appl. Environ. Microbiol.
69: 3110-3118
[Abstract]
[Full Text]
-
Jumpponen, A., Newsham, K. K., Neises, D. J.
(2003). Filamentous ascomycetes inhabiting the rhizoid environment of the liverwort Cephaloziella varians in Antarctica are assessed by direct PCR and cloning. Mycologia
95: 457-466
[Abstract]
[Full Text]
-
Landeweert, R., Leeflang, P., Kuyper, T. W., Hoffland, E., Rosling, A., Wernars, K., Smit, E.
(2003). Molecular Identification of Ectomycorrhizal Mycelium in Soil Horizons. Appl. Environ. Microbiol.
69: 327-333
[Abstract]
[Full Text]
-
Valinsky, L., Della Vedova, G., Jiang, T., Borneman, J.
(2002). Oligonucleotide Fingerprinting of rRNA Genes for Analysis of Fungal Community Composition. Appl. Environ. Microbiol.
68: 5999-6004
[Abstract]
[Full Text]
-
Klamer, M., Roberts, M. S., Levine, L. H., Drake, B. G., Garland, J. L.
(2002). Influence of Elevated CO2 on the Fungal Community in a Coastal Scrub Oak Forest Soil Investigated with Terminal-Restriction Fragment Length Polymorphism Analysis. Appl. Environ. Microbiol.
68: 4370-4376
[Abstract]
[Full Text]
-
Ranjard, L., Poly, F., Lata, J.-C., Mougel, C., Thioulouse, J., Nazaret, S.
(2001). Characterization of Bacterial and Fungal Soil Communities by Automated Ribosomal Intergenic Spacer Analysis Fingerprints: Biological and Methodological Variability. Appl. Environ. Microbiol.
67: 4479-4487
[Abstract]
[Full Text]
-
Glandorf, D. C. M., Verheggen, P., Jansen, T., Jorritsma, J.-W., Smit, E., Leeflang, P., Wernars, K., Thomashow, L. S., Laureijs, E., Thomas-Oates, J. E., Bakker, P. A. H. M., van Loon, L. C.
(2001). Effect of Genetically Modified Pseudomonas putida WCS358r on the Fungal Rhizosphere Microflora of Field-Grown Wheat. Appl. Environ. Microbiol.
67: 3371-3378
[Abstract]
[Full Text]
-
Röling, W. F. M., Kerler, J., Braster, M., Apriyantono, A., Stam, H., van Verseveld, H. W.
(2001). Microorganisms with a Taste for Vanilla: Microbial Ecology of Traditional Indonesian Vanilla Curing. Appl. Environ. Microbiol.
67: 1995-2003
[Abstract]
[Full Text]
-
Whipps, J. M.
(2001). Microbial interactions and biocontrol in the rhizosphere. J Exp Bot
52: 487-511
[Abstract]
[Full Text]
-
Borneman, J., Hartin, R. J.
(2000). PCR Primers That Amplify Fungal rRNA Genes from Environmental Samples. Appl. Environ. Microbiol.
66: 4356-4360
[Abstract]
[Full Text]
-
Peters, S., Koschinsky, S., Schwieger, F., Tebbe, C. C.
(2000). Succession of Microbial Communities during Hot Composting as Detected by PCR-Single-Strand-Conformation Polymorphism-Based Genetic Profiles of Small-Subunit rRNA Genes. Appl. Environ. Microbiol.
66: 930-936
[Abstract]
[Full Text]