Empresa Brasileira de Pesquisa Agropecuária, Embrapa Agrobiologia, Seropédica, Rio de Janeiro, Brazil,1 Center for Terrestrial Ecology, Netherlands Institute of Ecology, 6666 ZG Heteren,2 Microbiological Laboratory for Health Protection, National Institute for Public Health and the Environment, NL-3720 BA Bilthoven, The Netherlands3
Received 4 September 2003/ Accepted 19 November 2003
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
rRNA genes have become the most widely used targets for detection of AMF in environmental samples (13). Several PCR-based strategies targeting rRNA genes have recently been developed to detect AMF in DNA extracted from roots, soil, or spores (23, 29, 30, 56). Such strategies have provided new insights into AMF diversity by circumventing the need for trap cultures and morphological identifications, which can be highly biased, time-consuming, and inaccurate. Despite these advances, the operational taxonomic units obtained in most of these works can only be identified precisely to genus level or above. Thus, little progress has been made in species characterization and identification per se, which are still strongly dependent on morphological analysis and the investigator's level of expertise. Few studies have actually used the rRNA genes to identify species of AMF (38), with most analyses being limited to the detection of defined species of interest (32, 56).
Molecular analyses have revealed that a single AMF isolate or even individual spores may contain substantial heterogeneity among rRNA gene copies (2, 12, 31, 33, 34, 40; for a recent review, see reference 41), which may be unevenly distributed in the heterokaryotic nuclei of AMF spores (31, 54). Intraspecific rRNA heterogeneity seems to be a common phenomenon in AMF as well as in other groups of organisms, such as bacteria (1, 37), plants (11), insects (52), and crustaceans (18). However, little progress has been made in the interpretation of this heterogeneity. Such heterogeneity may lead to overestimations of the number of species when interpreting clone libraries of rRNA recovered from the environment (14). However, if the heterogeneity is consistent within a species, intraisolate heterogeneity might be used as an advantage to generate species-specific rDNA fingerprints for AMF detection and identification.
The genus Gigaspora represents an ecologically and economically (4, 42) important group within the Glomeromycota, and numerous studies have been dedicated to identify species within this genus and to study their ecology. The taxonomy of the genus Gigaspora has recently been revised by morphological (5), fatty acid methyl ester (6), molecular (3), and combined (32) approaches. Among the eight Gigaspora species described to date, Bentivenga and Morton (5) considered five to be valid species based on spore morphology: G. albida, G. decipiens, G. gigantea, G. margarita, and G. rosea. Two species were considered synonymous with previously described Gigaspora species (G. candida = G. rosea, and G. ramisporophora = G. margarita), and one species, G. tuberculata, was considered synonymous with Scutellospora persica. However, there are few useful morphological characters for Gigaspora species determination, and character ranges such as spore size and color often overlap between species (5). Bago et al. (3) used molecular signatures within the V9 region of the 18S rRNA gene as diagnostic characters for Gigaspora spp. identification. These authors were able to identify three distinct groups among the currently recognized species of Gigaspora: the Gigaspora rosea group (G. rosea and G. albida), Gigaspora margarita group (G. margarita and G. decipiens), and Gigaspora gigantea. Their analysis also revealed intraspore heterogeneity, which caused ambiguities in the signatures they found. Lanfranco et al. (32) continued the molecular characterization of selected Gigaspora species and proposed a number of species-specific primer sets, and Yokoyama et al. (58) developed primers based on satellite fragments to identify specific isolates of G. margarita. Thus, while important strides have been made recently, rapid and reliable methods to assess Gigaspora diversity are still lacking. Furthermore, little is known about the establishment, distribution, diversity, and competitiveness of Gigaspora spp. in the field (4, 42).
PCR-denaturing gradient gel electrophoresis (DGGE) was initially developed to study mutations. Nowadays, it has become one of the most applied culture-independent techniques to study the community structure of microorganisms (35). Separation in DGGE is based on differences in sequence composition that affect the melting behavior of the amplicons, causing a decrease in the electrophoretic mobility of a partially melted DNA molecule in a polyacrylamide gel containing a linearly increasing gradient of DNA denaturants (for more information, see reference 35). Recently, Kowalchuk et al. (30) successfully applied PCR-DGGE to study the community structure of AMF associated with Ammophila arenaria.
This study had three main goals: (i) to develop a rapid and reliable method to characterize and identify Gigaspora species, based upon PCR-DGGE analysis of variable regions of the 18S rRNA gene; (ii) to assess the level of intraspore and intraisolate 18S rDNA heterogeneity within the genus Gigaspora and evaluate PCR-DGGE as a method for studying this phenomenon; and (iii) to test the application of a Gigasporaceae-specific PCR-DGGE strategy for the assessment of Gigaspora diversity in environmental samples.
| MATERIALS AND METHODS |
|---|
|
|
|---|
|
Soil samples.
Samples were taken from an 8-year-old grassland field in a cattle farm in Brazil. The farm, Agropecuária Lopes, was located in Santo Antônio, Goiás State (16°28'00"S, 49°17'00"W, at 823 m above sea level), Brazil. The grassland was dominated by Brachiaria decumbens, which had replaced the native vegetation in the Cerrado (savannah) biome.
Intact soil cores (7.5 cm in diameter; 8.0 cm deep) were collected with polyvinyl chloride cylinders. The soil, a clayey dark red oxissol, had a pH of 5.5 (soil/water ratio, 1:2.5 [vol/vol]). It contained 0, 2.4, and 0.7 cmol of charge per kg (dry weight) of soil (cmolc) of Al, Ca, and Mg, respectively, dm-3 in a 1 N KCl extraction and 0.16 and 1 mg of P and K, respectively, dm-3 in a Mehlich I extraction. The samples were transported to Embrapa Agrobiologia, Seropédica, Rio de Janeiro, Brazil, and used to establish 10 trap cultures to assess the diversity of AMF. Brachiaria decumbens was used as the host plant. The trap cultures were sent to the Netherlands for further analysis. Two of these trap cultures were selected: one contained large numbers (sample A) and the other contained small numbers (sample B) of Gigaspora spores (Table 2). Spore identification and counting as well as DNA extraction were performed with three replicates of 30 g of soil inoculum each. The procedure for morphological identification of spores is described in de Souza et al. (16).
|
Control experiments with multiple target and nontarget AMF species.
In order to evaluate the effect of different target and nontarget species combinations on the detection limits and reproducibility of PCR-DGGE banding patterns, we combined DNA from strains Scutellospora heterogama CNPAB2 and G. margarita CNPAB16 at different ratios. The ratios used were Scutellospora heterogama CNPAB2 to G. margarita CNPAB16 at 1:1; 1:5; 1:10; 1:25; 1:50, and 1:100 and G. margarita CNPAB16 to S. heterogama CNPAB2 at the same ratios. In addition, we also combined nontarget DNA obtained from Glomus clarum CNPAB5 in ratios ranging up to 100:1 with G. margarita CNPAB16. Three replicas were performed for each combination.
Greenhouse experiment.
To ensure that the PCR-DGGE approach was sensitive enough to detect multiple species of Gigaspora and Scutellospora in root samples, we performed a greenhouse experiment. Clover plants (Trifolium pratense) were inoculated or not with a mixture of soil inoculum containing G. margarita (CNPAB16), S. gregaria (CNPAB7), S. heterogama (CNPAB2), and S. reticulata (CNPAB11). To ensure nodulation, the clover plants were also inoculated with soil filtrate containing rhizobia collected from clover field plots. The plants were grown in plastic cone pots containing 250 ml of a mixture of clay soil and sand (1:1, vol/vol) at pH 6.2 (soil/water ratio, 1:2.5). The pots were fertilized intermittently with 1/10th-strength nutrient solution (25) without N and P. After 2 months of growth, the pot contents were harvested. The soil was carefully removed from the roots with tap water, and the root system was cleaned by ultrasonication (60 W; B-2200 E1; Bransonic) twice for 3 min each in autoclaved water, followed by a final wash with autoclaved water. The roots from each pot were divided into subsamples for either DNA extraction or assessment of colonization rate by the method of Giovannetti and Mosse (21).
DNA extraction from spores, roots, and soil samples.
Spores were removed from the freezer (-80°C) and crushed with a micropestle (Treff AG, Degersheim, Switzerland) in 40 µl of 10 mM Tris-HCl buffer, pH 8.0, with 10 µl of 20% Chelex 100 (Bio-Rad Laboratories, Hercules, Calif.), for single spores. The same procedure was used with multiple spores except that the reagent concentrations were 80 and 40 µl for buffer and Chelex, respectively. The tubes were then incubated at 95°C for 10 min, chilled on ice, and centrifuged at 10,000 x g for 2 min. The supernatant was carefully transferred to a new tube and stored at -20°C until use. Samples from multiple spore isolations were treated with RNase before being stored. DNA extractions from trap plants with bulk soil and root material were performed with the UltraClean soil DNA isolation kit according to the manufacturer's instructions (MoBio Laboratories, Solana Beach, Calif.). Prior to DNA extraction, the samples (10 g of soil or 2 g of root) were homogenized and ground under liquid N2 with a mortar and pestle. A subsample of 0.5 g of bulk soil or root material was used for each DNA extraction. After extraction, the soil- and root-derived DNA was purified once more with the Wizard DNA purification kit (Promega, Madison, Wis.) as described by the manufacturer. For the greenhouse experiment, we extracted DNA from plant roots with the protocol described by Edwards et al. (17) with 50 mg of liquid nitrogen-powdered roots.
Nested PCR conditions for amplification from spore DNA.
The DNA from spores was first amplified with the forward primer NS1 in combination with reverse primer ITS4, covering the region from the beginning of the 18S rRNA gene through the 5' end of the 25S rRNA gene (57). Primer positions are given in Fig. 1, and primer sequences, references, and PCR conditions are provided in Table 3. These reactions were performed in a final volume of 15 µl, with 5 µl of template DNA. The PCR mixture was composed of 200 µM each of the four deoxynucleoside triphosphates, 1.5 µM MgCl2, a 0.4 µM concentration of each primer, and 1 U of Expand high-fidelity DNA polymerase (Roche Diagnostics, Nederland B.V., Almere, The Netherlands) according to the manufacturer's recommended buffer conditions. All reactions were performed in a PTC200 thermal cycler (MJ Research; Waltham, Mass.).
|
|
Gigasporaceae-specific PCR conditions for the analysis of environmental samples.
To obtain Gigasporaceae-specific products from soil, roots, and spores from trap cultures and the greenhouse experiment, the first step of the nested PCR combined primers FM6 (this study) and GIGA5.8R (39) (see Fig. 1 and Table 3 for primer positions and sequences). PCR mix was prepared as described above, and the DNA extracted from soil and root samples was diluted 1:50 to 1:100 and used as template. The product of this first PCR amplification was diluted 500- to 1,000-fold, depending on product concentration, and used as the template for a second PCR with the primer pair NS7 (57) with GC clamp in combination with primer F1Ra (this study) as described above.
DGGE analysis.
All DGGE analyses were performed with the D-Gene system (Bio-Rad), with gradients of 25 to 40% and 32 to 42% denaturant and running conditions of 75 V for 16 h for the NS31 and GC/AM1 and 95 V for 17 h for the NS7 and GC/F1Ra primer combinations, respectively. Gels were run in 0.5x TAE (Tris-acetate-EDTA) buffer at a constant temperature of 60°C. Gels were stained for 20 min in MilliQ water containing 0.5 mg of ethidium bromide liter-1 and destained twice for 15 min in MilliQ water prior to UV transillumination. Gel images were digitally captured with the ImaGo system (B & L, Maarssen, The Netherlands).
DGGE banding patterns were assessed by cluster analysis with a Jaccard similarity coefficient, and similarities between profiles were depicted as a dendrogram constructed by the unweighted pair group method with arithmetic average (UPGMA) within the Bionumerics program, version 2 (Applied Maths, Kortrijk, Belgium). The banding pattern of G. rosea FL105 was used as a marker to standardize different gels.
Tests for reproducibility.
All the PCR amplifications and DGGE analyses were performed with three independent single-spore DNA isolates and compared with multiple-spore isolates. This was done to detect potential artifacts due to spore contamination (43) and to evaluate the reproducibility of the method. For soil and root DNA, three subsamples were compared for each trap culture.
Recovery of DNA from DGGE gels.
The most prominent bands obtained in the DGGE profiles of the trap culture spores were sequenced. The middle portion of selected DGGE band was excised, and approximately 60 mg of acrylamide gel material per band was transferred to a 0.5-ml microcentrifuge tube containing 40 µl of MilliQ water and frozen at -80°C for 1 h. Subsequently, the gel material was crushed with a plastic pellet mix (Treff AG), and the tubes were incubated at 37°C for 3 h. After centrifugation at 11,000 x g for 60 s, the supernatant was transferred to a new tube, and 1 µl of it was used as the template for subsequent PCR-DGGE analysis to check band position and purity. This procedure was repeated until a single sharp band was detected. After that, PCR was performed with the same primer pair used in the DGGE analysis without the GC clamp, and the product was prepared for sequencing analysis.
Cloning of AMF rDNA.
In order to obtain clones for different variants of the ribosomal genes present (ribotypes) in one species (spore), amplicons were obtained from DNA extracted from individual spores after PCR amplification with the primer pair NS1 and ITS4 (57) as described previously. PCR products were purified with the High Pure PCR product purification kit (Boehringer, Mannheim, Germany). Purified PCR products were then cloned into the pGEM-T easy vector, with Escherichia coli strain JM109 used for transformation, according to the procedure given by the manufacturer (Promega Benelux, Leiden, The Netherlands). The clones obtained were cultured and, after plasmid extraction by the Wizard Plus SV miniprep DNA purification system (Promega, Benelux), used as templates for PCR (see below).
Clone selection with DGGE and sequence analysis.
Plasmids containing an insert obtained from different Gigaspora isolates were used as the templates for reactions with the forward primer NS7-GC in combination with the reverse primer F1Ra as described above. DGGE screening was performed as described above, and the PCR products obtained from the original isolates were used as reference to select clones that corresponded to each of the ribotypes detected in each isolate examined.
Prior to sequencing, DNA templates were purified with Qiaquick purification columns. Sequencing reactions were performed with the Perkin Elmer Biosystems Big Dye terminator sequence reaction kit (Perkin Elmer, Foster City, Calif.) and run on a Perkin Elmer 3700 capillary sequencer at the National Institute for Public Health and the Environment (Bilthoven, The Netherlands).
Sequence alignments.
Sequences recovered from the GenBank/EMBL database or generated in this work were first aligned with Clustal-X (53), and then the alignment was improved by manual inspection. Phylogenetic analysis was conducted with the parsimony method in PAUP* version 4.0 Beta 10 (51).
Nucleotide sequence accession numbers.
The sequences and alignment generated in this study were deposited in the EMBL database under accession numbers AJ539236 to AJ539305 and alignment number ALIGN_000606.
| RESULTS |
|---|
|
|
|---|
DGGE profiles targeting the V9 region of the 18S rRNA gene.
PCR-DGGE analysis of the V9 region of the 18S rRNA gene could differentiate all Gigaspora species based on the type materials used, including those (G. candida and G. ramisporophora) declared to be invalid by morphological analysis (Fig. 2). Almost all isolates yielded multiple bands, indicating the presence of intraspore variation between ribotypes in this region in the strains examined. To test the consistency of the PCR-DGGE profile for a given isolate, we examined several single- and multiple-spore DNA isolations per isolate. No between-spore variations could be detected for any of the 48 isolates tested with the exception of G. albida CL151 and G. margarita UFLA36. The former produced two very similar patterns (CL151a and CL151b in Fig. 3), but in the CL151b type, one of bands were absent and the intensity of the lower band was higher than that in the CL151a type. The latter accession produced two very different banding patterns (UFLA36-T1 and UFLA36-T2 in Fig. 3). These two strains may each actually represent two coisolated populations (see the Discussion). No difference was observed between different cultures of the same accession strain maintained in different laboratories (data not shown).
|
|
The banding patterns for strains of the same species were generally highly similar and clustered together (Fig. 3). Isolate-specific bands could be identified for some strains of the species G. albida, G. gigantea, and G. margarita but generally not for G. decipiens and G. rosea, although accession G. rosea MA457C appeared to contain a number of unique ribotypes (Fig. 3).
PCR-DGGE detected potential misidentifications by revealing some patterns that did not cluster with their respective type material. For instance, within the species G. albida, accession UFLA24 clustered with the G. ramisporophora type material (CNPAB22) and accession INVAM927 produced a pattern identical to that of the large group of G. rosea isolates. Similarly, G. rosea BR235 grouped together with G. albida isolates; it was the blind sample sent by INVAM's curator. This result confirms the supposed misidentification of this accession based upon morphological characteristics (J. B. Morton, personal communication).
Although G. candida BEG17 and G. ramisporophora CNPAB22 are both considered to represent nonvalid species names based upon morphological evaluation (5), they could be clearly distinguished from the other species examined. Interestingly, the pattern of G. candida BEG17 was more similar to that of the G. albida isolates than to those of G. rosea, the species to which it was previously assigned (Fig. 3). G. ramisporophora CNPAB22 generated a banding pattern that was more similar to those of G. rosea (Fig. 3) than those of G. margarita, the species to which G. ramisporophora was assigned based upon morphological characteristics. A robust molecular characterization of the species G. candida and G. ramisporophora will require the analysis of additional isolates. Unfortunately, no other well-defined isolates of G. candida are available at this time, to the best of our knowledge, and very few well-defined isolates of G. ramisporophora are available.
Sequence analysis of PCR-DGGE banding patterns.
To gain further insight into the nature of the intra- and interspecies heterogeneity detected, sequence information was obtained for the region analyzed by PCR-DGGE for each of the ribotypes observed within representative isolates of each species (Table 4). Sequence analysis confirmed the identity of bands that displayed the same migratory behavior. Furthermore, sequence and phylogenetic analyses (Table 4) also confirmed a number of the relationships depicted in the dendrogram analysis of the PCR-DGGE patterns. For instance, G. margarita and G. decipiens are closely related and distinct from the group formed by G. albida, G. candida, G. ramisporophora, and G. rosea and that of G. gigantea. In addition, G. gigantea UFLA872 is quite distinct from all other accession strains and contains a DNA signature, CGCGTG, that has been reported to occur in Scutellospora species (3). With reference to the characterization of putative invalid species, sequence analysis confirmed ribotype overlap between G. candida BEG17 and G. albida BR607A, including a shared DNA signature, TAGGTT, which is distinct from that of G. rosea (3). The ribotype sequences obtained for G. ramisporophora were also more similar to those of G. albida, G. candida, and G. rosea than to those of G. margarita despite apparent similarity in the morphology of G. ramisporophora and G. margarita spores. Those results do not support the reclassification based on morphological analysis by which G. candida was reclassified as being synonymous with G. rosea and G. ramisporophora was reclassified as being synonymous with G. margarita (6).
|
(ii) Detection of Gigasporaceae species in greenhouse experiments and environmental samples.
In a controlled greenhouse experiment, clover plants were inoculated with four AMF species, one species of Gigaspora and three of Scutellospora. Despite a colonization level of less than 20%, as determined by microscopic inspection, AMF-specific products could be easily detected with a nested PCR and DGGE approach. All four AMF species could be detected, although secondary bands had to be used to determine the presence of the two Scutellospora species, since the two species used presented a prominent DGGE band in the same position (Fig. 4).
|
|
| DISCUSSION |
|---|
|
|
|---|
Due to the overlap of some ribotypes across species boundaries (e.g., G. albida, G. candida, G. ramisporophora, and G. rosea or G. decipiens and G. margarita; see Fig. 3 and Table 4), it became clear that the entire pattern of 18S rDNA types needed to be used for species characterization and identification, as opposed to just a single sequence. Thus, in contrast to other molecular approaches, such as PCR followed by cloning and sequencing, in which intraspore rDNA heterogeneity impairs interpretation, the PCR-DGGE approach uses this heterogeneity to advantage for the generation of highly reproducible isolate- or species-specific patterns. Furthermore, in comparison to some other molecular methods, PCR-DGGE is easy, rapid, and inexpensive to perform.
We used two different nested PCR strategies. The first nested PCR strategy uses family-specific primers in the first PCR. With that strategy, it was possible to analyze DNA extracted from environmental samples with small amounts of material extracted directly from spores, roots, or soil without the problem of nonspecific amplification. However, despite the family-specific nature of the nested PCR strategy and the presence of both Gigaspora and Scutellospora species in the mixed spores, roots, and soils analyzed, the latter genus was not detected. However, in a control experiment, we could detect the simultaneous presence of all four of the AMF species introduced into the root system (Fig. 4). PCR-DGGE strategies will generally fail to detect some minority populations (i.e., <1% of the total target) (10, 35), and our experiments showed that a total spore volume of 10% of the total was typically necessary to ensure the detection of a minority target species.
Furthermore, the DGGE signals detected are based on the relative number of template molecules of a given species in a sample. Thus, even assuming equal amplification efficiencies, it remains difficult to translate the levels of DGGE signal detected for certain AMF species, as the ratio of 18S rDNA target to fungal biomass and the ecological importance of this ratio are not known at present.
The second nested PCR strategy used universal primers that, in our case, amplified the whole 18S rRNA gene plus the internal transcribed spacer (ITS) region. That fragment was cloned, followed by PCR-DGGE selection targeting variable regions in the cloned fragment. That strategy proved to be an excellent approach to study inter- and intraspecies heterogeneity in the 18S rDNA. The principles of that strategy can be applied to study polymorphism in any other gene found in AMF. The application of a similar approach with faster-evolving genes might allow better discrimination between isolates and also might help to shed some light on evolutionary issues such as genetic drift and recombination in these ancient asexual organisms (19).
It is interesting that the DGGE system used was highly efficient in detecting sequence variation within the V9 region of the 18S rDNA gene. Within the 344 bp (without GC clamp) of sequence targeted, 303 bp were constant, 17 bp were parsimoniously uninformative, and 24 bp were parsimoniously informative (Table 4). The V9 18S rDNA fragments analyzed by PCR-DGGE are within the optimal size range for DGGE analysis, i.e., below 500 bp. Most of the mutations found were transitions (37 of 41), which allow discrimination by DGGE, as they always cause a change in the melting temperature. However, transversions were detected in four of the parsimoniously informative positions (positions 274, 275, 279, and 299; see Table 4), and some of them did not affect the melting temperature of the amplicons; consequently, they might not be detected if they were the only mutation present in the amplicon. Despite having the ability to discriminate most single-base-pair differences within 18S rDNA V9 region fragment analyzed here, any DGGE analysis will be limited by the amount of heterogeneity, type of mutation, and fragment size of the target region (36). This was exemplified by our analyses of the V3-V4 region, which contained less variation than the V9 region and also used a longer fragment (550 bp without the GC clamp). Interestingly, the V3-V4 region could be used to discriminate between Glomus species (30). Although the PCR-DGGE banding patterns were highly reproducible within a given isolate across various DNA isolations, exact banding patterns are dependent on the electrophoresis conditions used (data not shown), and analytical consistency and the use of type strains (isolates) as markers are critical when comparing samples.
The similarity between the sequences from the V9 region analyzed was high both between and within species (range, 98.6 to 100%; Table 4). Of the variable positions described, only 24 were phylogenetically informative, hampering robust phylogenetic analysis. As such, the comparison of V9 18S rDNA fingerprints (Fig. 3) provided a more reliable method of comparison than tree construction based on the sequence variations described in Table 4. With proper primer design, PCR-DGGE strategies as implemented here are also ideal for determining specific ribotypes of other AMF genera or other loci in experiments designed to address AMF reproduction and evolution (19, 41), as well as similar issues with respect to intraspecies heterogeneity among rDNA copies in bacteria (1, 14, 37). It should be noted that our analysis of single spores does not address the homokaryotic or heterokaryotic nature of AMF nuclei (20, 54), although the sensitivity of PCR should permit PCR-DGGE of single AMF nuclei.
Comparison between single- and multiple-spore DNA isolations and geographically diverse isolates.
The PCR-DGGE banding patterns of all single-spore DNA isolations tested for the same isolate were identical for 46 of 48 accession strains tested (Fig. 3). Thus, at least to the level of detection afforded by the system used here, a single spore appeared to contain the full range of variation of ribotypes present in an entire spore population. The two exceptions were both isolates recovered by trap cultures. Trap cultures are the most common way to isolate AMF from field samples. However, if more than one morphologically closely related species are coisolated in the same culture, further discrimination by spore morphology is difficult. One, G. margarita isolate, UFLA36, produced spores with two very different DGGE patterns that clustered apart from each other; one type (UFLA36-T1) clustered in the G. margarita-G. decipiens cluster, and the other type (UFLA36-T2) clustered with G. ramisporophora and G. rosea. In the case of G. albida CL151, the two spore types were similar, differing only by the absence of one band in one spore type and the relative intensity of one of the other bands (Fig. 3). In those cultures, further purification by single-spore pot culture followed by spore identification with PCR-DGGE can be used to purify and distinguish those populations.
Interestingly, some geographically distinct isolates were more similar than some isolates that were found at the same site. Within the species G. gigantea, for instance, strains NC110A and NC150 are rather different single-spore cultures recovered from the same site (7), whereas isolates CUT, MA453A, and MN453A-7 had identical rDNA patterns despite being isolated from disparate locations. Similarly, G. albida CL151 type a from the United States and BR601 from Brazil were identical (Fig. 3), while other sympatric populations showed more diversity. The presence of strain-specific bands should allow the tracking of specific AMF populations in studies dedicated to unraveling the ecological significance of such sympatric populations. The band intensities observed for different geographic isolates differed in some cases. This result suggests that the proportion of the different ribotypes may differ between isolates of the same species, but more quantitative methods, such as introduction of an internal standard in DGGE experiments (10), will be necessary to address this question.
PCR-DGGE characterization versus other schemes applied previously.
The lack of discrete and diagnostic characters for species identification within the genus Gigaspora has been a major obstacle to ecological studies of this genus. The use of PCR-DGGE as applied in this study clearly offers a higher level of discrimination and reliability than previous methods used to address this issue. For instance, based on spore morphology, G. candida was considered synonymous with G. rosea and G. ramisporophora was considered synonymous with G. margarita (5). However, our results not only do not support this reclassification, they show that the so-called invalid species are actually less related to their supposed synonymous species than to other species, based on the rDNA marker. We could also distinguish between species that were grouped together on the basis of the molecular signatures proposed by Bago et al. (3). Our detection of high degrees of heterogeneity within the V9 region of the 18S rRNA gene, which encompasses the sequence stretch examined by Bago et al. (3), also explains the ambiguities found within the signature sequences that they defined. The mixed PCR products recovered from an isolate could be resolved by PCR-DGGE (Table 4) but produced ambiguities at heterogeneous positions upon direct sequencing as preformed by Bago et al. (3).
Our results suggest that G. albida INVAM927 should be reassigned as G. rosea (Fig. 3). G. albida INVAM927 was one of the strains used by Bentivenga and Morton (5) as type material to redescribe this species. In contrast, the material we used for this isolate came from the INVAM in Florida and not from West Virginia University. Unfortunately, a comparison between the material examined in this study and the original material used by Bentivenga and Morton (5) is no longer possible, as this accession strain has been lost from the INVAM collection (J. B. Morton, personal communication). The accession strain G. albida UFLA24 must also be reassigned as G. ramisporophora on similar grounds. One of the strains identified as G. rosea obtained from INVAM was actually G. albida, based on information sent by the curator of the collection. This isolate was correctly identified by our PCR-DGGE analysis (Fig. 3) as strain BR235, confirming the morphological analysis.
Despite being characterized as G. gigantea, accession UFLA872 presented a PCR-DGGE pattern and a DNA sequence that did not match those of any other isolate examined (Fig. 3). Spores of G. gigantea UFLA872 have the same size range as expected for G. gigantea and also exhibit a cytoplasm color typical of G. gigantea. This feature is considered a unique identifying characteristic of this species (46). Other authors have also described conflicts between morphological and molecular identifications of AMF. Bago et al. (3) suggested the reassignment of isolate DAOM194757, morphologically identified as G. margarita, to the G. rosea group. Lanfranco et al. (32) confirmed this result and also suggested reassignment of isolate E29 (G. margarita based on morphology) to the G. rosea group on similar grounds.
Implications of 18S rDNA heterogeneity for ecological and evolutionary studies.
Heterogeneity between rRNA markers within a species or single individual is a phenomenon that has been described for a wide range of organisms (1, 11, 14, 18, 31, 37, 41, 52). In our PCR-DGGE approach, we have not only used this heterogeneity to characterize and identify species, but have also combined this with a cloning and screening strategy to tease apart the relationships between the 18S rDNA variants detected in a single spore (species). Knowledge of this intraspecific variation is fundamental for proper assignment of operational taxonomic units for molecular sequence identification and phylogenetic analysis (14). It is possible that cloning and sequencing will underestimate the number of species if the clones chosen for analysis happen to harbor those ribotypes held in common between the different species present in a sample. It is more likely, however, that such studies will overestimate species diversity (14, 40). Our results show that some different morphospecies share certain ribotypes but also contain species-specific variants. Given the asexual nature of AMF (31), the mechanisms for the establishment and maintenance of such patterns of sequence diversity remain to be discovered.
Although rRNA genes have been highly useful markers for the phylogenetic study of microorganisms, including Gigaspora (44), data gained from this marker alone clearly cannot resolve all phylogenetic and evolutionary issues for these organisms. The analysis of other molecular markers is an urgent issue for reconstructing phylogenetic relationships within this genus and for other AMF. Thus, although we were able to detect and utilize interesting patterns of rDNA heterogeneity within the Gigasporaceae for detection and identification purposes, further studies are necessary to explain these patterns and understand how they fit into the scheme of AMF life history and evolution.
The ability to assess Gigaspora diversity directly in environmental samples opens new possibilities for studying the ecology of this group under field conditions without the need for trap cultures. Recently, G. margarita was found in Europe, and its occurrence seems to be affected by the tillage system used (28). Some Gigaspora species are known to harbor an endosymbiont of a proposed new bacterial genus, candidatus Glomeribacter gigasporarum (8, 9), and the characterization of fungal and bacterial partners might clarify the evolutionary and ecological aspects of that symbiosis. PCR-DGGE targeting the V9 region of the 18S rDNA provides a fast and reliable method to identify Gigaspora species, to assess Gigasporaceae diversity in field conditions, and to characterize inter- and intraspecies rDNA heterogeneity.
| ACKNOWLEDGMENTS |
|---|
We thank the following researchers for providing access to germ plasm: L. Abbott, G. Bécard, V. Bianciotto, D. D. Douds, V. Gianninazzi-Person (curator, European Bank of Glomeromycota), R. Herrera-Peraza, L. C. Maia, J. Morton (curator, INVAM), M. Saito, J. O. Siqueira, and C. Walker. We thank the anonymous reviewers for helpful comments.
| FOOTNOTES |
|---|
Publication number 3268 of the NIOO-KNAW. ![]()
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||