This paper presents a DNA macroarray for the identification of the majority of
Pythium species. This relatively new technique is emerging as a powerful molecular tool for multiplex detection of plant-pathogenic fungi. The ability to detect, identify, and differentiate
Pythium spp. by PCR-based macroarray technology, without the need for cultivation, has great potential for enhancing the throughput of detection and diagnostic procedures, adding a new tool for the detection and identification of
Pythium species in a wide range of studies. To date, most DNA macroarrays have been designed to detect and identify only a very limited number of species in a single reaction mixture. Lievens et al. (
29) developed a DNA array for the detection of
Fusarium oxysporum f. sp.
lycopersici,
Verticillium albo-atrum, or
V. dahliae, while Sholberg et al. (
46) reported on an array for monitoring five bacterial and fungal pathogens of apple diseases. DNA array technology can, in principle, be used to detect an unlimited number of different pathogens in a single assay mixture (
33,
52). The advantage of this is that it combines nucleic acid amplification with the unlimited screening capacity of DNA arrays, resulting in high sensitivity, specificity, and throughput capacity (
30). A DNA array for simultaneous detection of more than 40 different plant-pathogenic soilborne fungi from several genera and 10 bacterial pathogens that frequently occur in greenhouse crops was developed (
30). The DNA macroarray designed and tested in this study is one of the first to cover, almost completely, all of the known species of a genus which comprises pathogens of economic importance in field and greenhouse crops, biological control agents, and the causal agent of pythiosis in animals and humans.
In the present study, an oligonucleotide array was designed based on sequence heterogeneity within the ITS region of the rRNA gene cistron. The ITS is useful for detection and identification of species of oomycetes (
5,
26,
32,
33). The ITS rRNA gene of
Pythium species is flanked by regions that are highly conserved within the oomycetes, allowing the same primers to be used to amplify and label template DNA across a wide range of species and samples. In this study, there was a comprehensive database of sequences from the large rRNA gene subunit in oomycetes (
25,
40,
41), allowing the design of a reverse PCR primer with high specificity for this group. This primer specificity minimizes the amount and number of nonspecific amplicons produced, improving the efficiency of the PCR. For field or environmental samples, digoxigenin-labeled PCR products were generated with these oomycete primers, eliminating the need for preliminary isolation and identification of the target
Pythium species.
Closely related species within a sample were differentiated by their respective hybridization patterns on the DNA array. In most cases, the species targeted in this study exhibited specific hybridization with their corresponding oligonucleotides spotted on the array under the single-stringency conditions used. More than 92% of the spotted oligonucleotides exhibited positive and specific hybridization signals with their target species of
Pythium, indicating very high specificity of the sequences. This allowed reliable identification and differentiation of most of the plant-pathogenic species of
Pythium, including differentiation of very closely related species (e.g.,
P. aphanidermatum and
P. deliense) and genetic groups that eventually may be separated into multiple species. For example, two oligonucleotides (ulti-40 and ulti-41) were designed to detect
P. ultimum at the species level, while four were spotted to differentiate the two different intraspecific groups (for var.
ultimum, ult1-42 and ult1-43; for var.
sporangiiferum, ult2-45 and ult2-46). Amplicons from both groups gave positive hybridization signals with ult40 and ult41 but reacted only with oligonucleotides designed to be specific for each subgroup. The differentiation of these two groups with specific oligonucleotides supports suggestions by Klassen et al. (
20) and Lévesque et al. (
26) that there are at least two molecular groups within
P. ultimum. However, it appears to be a coincidence that the ex-type of
P. ultimum var.
sporangiiferum represents a less common genetic group than
P. ultimum var.
ultimum (
26).
Group- and species-specific oligonucleotides on DNA arrays have been reported previously (
29,
46), and their importance for differentiating clusters of closely related species is demonstrated in this study. Hybridization patterns of group- and species-specific oligonucleotides were consistent for some phylogenetic clusters within
Pythium that they were designed to mirror and revealed similarities between species with different names but identical or nearly identical ITS sequences. DNA arrays designed and reported to date (
26,
29,
30,
46) focused on simple detection of the target plant pathogen, with no further information at a higher taxonomic level than the species level. In this study, oligonucleotides were designed that reacted with all species having specific morphological structures such as filamentous (fil51 and fil52) or globose (IP110) sporangia. Forty-two species with filamentous sporangia exhibited positive reactions with oligonucleotide fil52, correlating this phylogenetic characteristic with a simple hybridization signal on a membrane. All of these species are classified in clade B in the recent study of
Pythium phylogeny (
25), which consists almost entirely of species with inflated and noninflated filamentous sporangia. Each species without a putative synonym(s) had a unique hybridization pattern made of a combination of cluster- and species-specific oligonucleotides. Some genetic subgroups also had unique hybridization patterns (e.g.,
P. ultimum and
P. irregulare).
Oligonucleotides pyril236 (location E19) and pyril237 (location F19), designed for
P. pyrilobum, and oligonucleotide scler679 (location I16), derived from
P. scleroteichum, cross-hybridized with some nontarget species of
Pythium. The results of a comprehensive BLAST search indicated that some of the cross-hybridizations could be explained by some degree of sequence match within the ITS region of the nontarget
Pythium species. In previously published studies, the species-specific oligonucleotide Val1, designed for the detection of
Verticillium albo-atrum, cross-hybridized with
V. tricorpus (
29) while oligonucleotide EA-H3d, designed for
Erwinia pyrifoliae, cross-reacted with a strain (G-5) of
E. amylovora (
46). Lievens et al. (
29) could not explain this discrepancy by the number of mismatches between the
V. tricorpus amplicons and the spotted oligonucleotide. Our results also did not always show hybridization consistent with known mismatches, suggesting that other factors could influence specificity when only a few mismatched nucleotides are involved. For example, work with the cystic fibrosis gene has shown that a single base mismatch is enough to design a mutation-specific oligonucleotide for a reverse dot blot assay (
19) and for differentiating tRNA species (
38) and strains of
Citrus tristeza virus (
6) by standard hybridization methods. In addition, PCR assays, enzyme-linked immunosorbent assays, and fluorescence assays (
36,
37) have been reported to distinguish several strains of
C. tristeza virus with only one or two nucleotide differences. Fessehaie et al. (
15) were able to differentiate closely related species or subspecies of
Clavibacter michiganensis with a mismatch of two nucleotides in oligonucleotide design. Even though discriminatory specificity is greatest when mismatches are in the center of the oligonucleotide (
18), the base composition of the entire immobilized fragment could greatly affect its hybridizing potential and stability. Lévesque et al. (
26) managed to eliminate some cross-reactions by increasing the stringency conditions. With the extensive number of oligonucleotides that are specific for closely related species in our array, this is not a practical option. Two separate hybridizations would be needed for each sample because the signal from some of the specific oligonucleotides would decrease or disappear when there is a perfect match at higher stringency. We believe that further advances in methods and software to find oligonucleotides, design them, and test them in silico may alleviate many of the cross-reaction problems.
The reasons why seven oligonucleotides (4%) did not show detectable gray scale signals with any
Pythium species, including the targeted species, are unknown. This could be partly attributed to the base composition of the oligonucleotides and conditions that are too stringent for them. Lievens et al. (
29) reported similar results when amplicons generated from
Verticillium nubilum DNA did not hybridize to any of the oligonucleotides designed for
Verticillium detection. The oligonucleotides might be situated within strong secondary structures that do not open up under the hybridization conditions used. Better prediction models for hybridizations that include the potential secondary structures of the template should reduce this problem.
The number of GenBank entries of
Pythium species isolates increased substantially subsequent to the initial design of the oligonucleotide probes that were used in this work, and this provided an additional opportunity to investigate oligonucleotide specificity and intraspecific variation. A BLAST analysis was done for all of our oligonucleotides against the 513 complete ITS sequences available in GenBank to assess specificity and intraspecific variations among the different isolates within a species. The collection of GenBank sequences represented several strains of each of the
Pythium species implicated in major and economically important plant diseases. The BLAST analyses of our oligonucleotides against the 500 entries of known species, synonyms, and newly described species within their respective clades confirmed the results obtained after hybridization of the oligonucleotides with the labeled amplicons of the ex-type or representative strains. In addition to using BLAST analyses to investigate specificity and intraspecificity, hybridization patterns of 50 isolates that included 12 major pathogenic
Pythium species demonstrated that isolates of a given species displayed no pattern variation or only minor quantitative differences in hybridization signals. All sequences of each species from GenBank matched some specific oligonucleotides, except for very few recently described species. Several of the newly described species were putative synonyms (
25), explaining why these entries perfectly matched oligonucleotides spotted on the DNA array to detect already existing specific
Pythium species or clades. The BLAST results also showed that the species-specific oligonucleotides on the array would not react with any recently described species that were not part of our array and that are not putative synonyms of known species. Most of these unequivocally new species matched their clade- or cluster-specific oligonucleotides. Therefore, as planned, a new species might be suspected to be present in a sample if a clade-specific oligonucleotide reacts but none of the species-specific oligonucleotides within that clade is positive. In cases like these, group oligonucleotides at a finer resolution would be necessary to highlight a potential new species or genotype, something we managed to incorporate for some groups. The only indisputably new species that we were unable to differentiate by BLAST analysis from its closest relative is
P. rostratifingens from
P. rostratum. We had a few different isolates of
P. rostratum, some of which ended up being described as the new species
P. rostratifingens (
10). Not surprisingly, what we originally designed as a species-specific oligonucleotide to cover all of the variations we observed within
P. rostratum became cluster specific for these two related species.
The array was used for direct detection and identification of some
Pythium species present in soil samples. Direct processing of the samples and hybridization of labeled amplicons to the array of oligonucleotides detected the presence of 13 species. These results were supported by results obtained by agar plating, root baiting, and partial subcloning. The only exception was
P. nunn, which could not be isolated by traditional methods. However, it should be pointed out that
P. nunn was isolated repeatedly in a nearby field (unpublished data). Results of this study demonstrate the feasibility and reliability of using DNA arrays to study
Pythium communities in the environment. Root baiting isolated the fewest isolates, confirming the work of Arcate et al. (
3). They observed increased diversity in communities through direct DNA extraction and hypothesized that this might be due to oospore populations of oomycetes that have gone undetected by baiting procedures. It is possible that some species detected by DNA array analysis may be less responsive to baits or may be less competitive under bait incubation conditions. The choice of baiting material can also influence the
Pythium species that can be detected. This inherent variation in results with the type of bait material and culture medium used is problematic. It is also clear from our results that
Pythium species with no corresponding oligonucleotide(s) on the array would not be detected. For example,
P. longandrum, a newly described species with no species-specific oligonucleotide on our array, was isolated by soil dilution plating and its identity was confirmed by sequencing but not on the array. BLAST analysis also did not show homology to any species-specific oligonucleotides but revealed some degree of similarity to group-specific oligonucleotides multi90 and marsi99, which were designed for the
P. multisporum and
P. marsipium subgroups, respectively. However, this was a situation where other species from this clade were present in the sample, namely,
P. minus and
P. acrogynum, which reacted with the cluster-specific oligonucleotides, preventing the detection of a possible new species. The reaction of clade-specific oligonucleotides without any species-specific oligonucleotide reactions would indicate either a new species or a new genotype of an existing species. This soil sample analysis demonstrated the specificity of the spotted oligonucleotides but also a shortcoming of the DNA array with new species. However, as new species are described and sequences are made available, the array could be updated to detect these species.
Molecular methods have become increasingly important to specifically detect pathogens, and different regions of the genome can be targeted to obtain the desired specificity (
28). However, there remain limitations that can hamper accurate pathogen detection. The lack of adequate sequence information, as well as finding DNA sequences that are shared by all members of a given species or genus, may be challenging. Currently, there is a trend to build up DNA sequence databases for identification known as “DNA barcodes” (
17) or microcodes (
47). From the sequence barcodes, species-, clade-, or genus-specific oligonucleotides can be identified (
12,
44). These oligonucleotide barcodes could form the basis of a high-density DNA array for environmental monitoring of species from several different genera.
Overall, the results reported here show that the oligonucleotide array designed and tested could be useful for the identification and detection of most species of
Pythium. It is possible that new genotypes of the species studied will be found and necessitate some modification of the current array. However, we are convinced that the hybridization patterns observed for this array are consistent with known genotypes, which probably encompass the majority of the pathogenic or economically relevant variation in this genus. This array could become a tool for rapid detection of
Pythium species and for determining population diversity in environmental samples. The hybridization signal on DNA arrays appears to be proportional to the quantity of target DNA for a given species (
30). Therefore, future work aimed at standardizing and quantifying hybridization signals could be instrumental in enhancing the direct estimation or enumeration of different species of
Pythium in soil samples in a single assay. Simplifying detection and enumeration of soilborne pathogens such as
Pythium could have a significant impact on soil ecological studies.