INTRODUCTION

TOP ABSTRACT INTRODUCTION RESULTS DISCUSSION MATERIALS AND METHODS LITERATURE CITED

The Stramenopiles comprise a diverse group of organisms that have recently been consolidated as a result of analysis of mitochondrial and ribosomal DNA sequences. Species that were previously referred to as heterokont algae and fungal oomycetes belong to this group of related organisms (Gunderson et al., 1987; Förster et al., 1990). Stramenopiles include autotrophic and heterotrophic species that may differ enormously in their morphology and mode of life, but most of which are vegetatively diploid and possess distinctive tinsel-like flagella at some point in their life cycle. Members of this group occupy key ecological niches in marine environments. For example, the diatoms constitute the most abundant component of marine plankton, while the brown algae (e.g. Fucus, Sargassum, and Laminaria spp.) may form extensive floating or attached communities such as the kelp forests. Stramenopiles have also succeeded in terrestrial environments as plant pathogens. Oomycetes (e.g. Peronospora, Pythium, and Phytophthora spp.) include many obligate and facultative parasites that cause disease on a wide spectrum of herbaceous and woody plants throughout the world in tropical and temperate environments (Erwin and Ribeiro, 1996).

Phytophthora sojae Kaufmann and Gerdemann (syn. Phytophthora megasperma f. sp. glycinea Kuan and Erwin) is the causal agent of stem and root rot of soybean (Glycine max [L.] Merr). Extensive outbreaks of Phytophthora root rot on soybeans were first noted over 40 years ago (Kaufmann and Gerdemann, 1958; Hildebrand, 1959), and the disease remains an endemic and serious problem in most soybean-producing areas (Schmitthenner, 1985; Doupnik, 1993; Wrather et al., 1995). Chemical methods may be used to control the disease, but the most cost-effective and widespread strategy to reduce the incidence of the disease has been through the development of resistant or tolerant soybean cultivars. Genetic analysis has shown that tolerance is a quantitative trait, whereas resistance is inherited qualitatively (Schmitthenner, 1985; Ward, 1990; Kasuga et al., 1997). Single dominant host resistance (Rps) genes and corresponding dominant avirulence (Avr) genes in the parasite have been described previously (Tyler et al., 1995; Whisson et al., 1995; Gijzen et al., 1996a). Despite breeding efforts to exploit such resistance-gene-mediated protection, resistant plants eventually become susceptible to disease due to the emergence of new virulent races (Anderson and Buzzell, 1992).

P. sojae is a diploid hemibiotroph that passes through several morphological phases to complete its life cycle. Asexual, single-celled zoospores are biflagellate, motile, and chemotactic to soybean plants (Ho et al., 1967; Morris et al., 1998). Zoospores encyst and germinate on the root or hypocotyl surface, and the resulting germ tube may swell to form an appressorium-like structure at the point of penetration into host tissues (Stossel et al., 1980). Resistance is manifested by rapid induction of host defenses, phytoalexin accumulation, and hypersensitive cell death at the infection site (Ward et al., 1989; Kamoun et al., 1999b). In susceptible interactions, large, water-soaked lesions develop as the pathogen rapidly invades the host. P. sojae spreads through the intercellular matrix of the plant and forms haustoria for intimate contact with host cells. Under moist conditions, zoospores are produced in sporangia that develop at the hyphal tips of vegetative mycelium. P. sojae is homothallic, and thus may reproduce sexually by the development of oospores through self-fertilization or by outcrosses between different strains (Förster et al., 1994). Estimations of genome size in P. sojae vary from 62 to 97 Mb (Rutherford and Ward, 1985; Mao and Tyler, 1991; Voglmayr and Greilhuber, 1998).

Despite thorough investigations into the biology of disease development, little is known concerning the pathogenic determinants of P. sojae. Molecular studies of pathogenicity and virulence of oomycetes are relatively rare compared with those on plant pathogenic fungi, bacteria, and viruses. Plant pathogenic basidiomycetes and ascomycetes, and model organisms such as Saccharomyces, Aspergillus, or Neurospora spp. are taxonomically far removed and fundamentally different in their cell wall composition, reproductive biology, and genetics from oomycete species (Judelson, 1997). Thus, expressed sequence tags (ESTs) and genome sequence information from fungal species may be of limited utility in studies of oomycetes. For these reasons, we undertook a study to profile gene expression patterns in P. sojae by analysis of ESTs. A model strain of P. sojae, P6497, was chosen for this analysis because it has been used for genome size estimations, genetic mapping of avirulence loci, and for construction of large insert bacterial artificial chromosome libraries arranged in ordered arrays. Sequence and annotation information that results from this research will also help to establish a public database to facilitate research on P. sojae and other oomycetes (Waugh et al., 2000).

We have constructed three cDNA libraries using mRNA from axenically grown mycelium, zoospores, and in planta infection sites. More than 1,000 ESTs were generated from each library to determine how patterns of gene expression change during development and pathogenesis. We show how the percentage G+C content of an EST sequence provides a measure that can help distinguish P. sojae from soybean cDNAs, and validate this with selected transcripts by hybridization analysis. Sequence data from each of the cDNA libraries were also systematically compared with EST data from Phytophthora infestans and all other non-human ESTs. Overall, our results indicate that P. sojae gene expression patterns and metabolic processes fundamentally shift during pathogenesis. We have identified transcripts encoding proteins that are likely to be involved in the recognition, attachment, penetration, and invasion of the host soybean plant.

	RESULTS

TOP ABSTRACT INTRODUCTION RESULTS DISCUSSION MATERIALS AND METHODS LITERATURE CITED

ESTs from Infected Plant Tissues Are More Likely to Produce Highly Significant Matches

Three directional cDNA libraries were constructed from in vitro-grown mycelium, from zoospores, and from P. sojae-infected soybean tissue. The average size of the cDNA inserts for each of the libraries was 0.83, 1.5, and 1.1 kb, respectively. Over 1,000 clones from each library were randomly chosen for plasmid purification and partial DNA sequencing, resulting in a total of 3,035 robust quality ESTs. Most cDNA clones were properly oriented in the vector, since only 2.7% of the ESTs possessed poly(A⁺) tracts adjacent to the T3 vector site. Thus, the vast majority of the ESTs reflect 5' end sequences of cDNA from mRNA transcripts.

After editing, an average length of 575 bp was used in database searches. Figure 1 shows the distribution of ESTs for each of the source libraries based upon deduced amino acid sequence homology to known or hypothetical proteins. From 20% to 23% of the ESTs from all three libraries returned a P value of >10². Of the remaining ESTs with P values of 10², fewer than 10% of these matched proteins of no known function. Highly significant matches occurred most frequently (45%) with sequences originating from the infected plant library. This reflects the current bias of the DNA sequence databases for plant proteins and genes relative to sequences from oomycetes and related organisms. Thus, soybean EST sequences from the infected plant library are more likely to result in highly significant match scores.

View larger version (35K): [in this window] [in a new window]   Figure 1.   Frequency distributions of ESTs from each of the source cDNA libraries according to P-value scores returned from BLASTX version 1.4.11 searches.

ESTs with P values of 10² were grouped into seven functional categories, as outlined in Figure 2, for tissue-specific expression profile comparisons. Matches to structural proteins comprised the most abundant class of ESTs in both the zoospore and mycelium libraries, whereas matches to metabolic enzymes were most numerous in the infected plant library. ESTs that matched ribosomal proteins and other factors required for protein synthesis were also highly represented in all three libraries.

View larger version (33K): [in this window] [in a new window]   Figure 2.   Classification of ESTs according to putative biological function. Protein matches resulting from the BLASTX searches were assigned to one of seven functional categories for life cycle comparisons. Shown is the percentage frequency of clones in each category for the total set of ESTs from each source library.

Two Populations of ESTs That Differ in Percentage G+C Content Are Present in the Infected Plant Library

Figure 3 shows a comparison of the percentage G+C content of ESTs from each of the cDNA libraries described in this study, and 1,300 soybean ESTs derived from a seed coat cDNA library. The average G+C content of soybean ESTs was 46%, whereas P. sojae zoospore and mycelium ESTs clustered around a mean of 58% G+C. The graphs indicate that the percentage G+C content of ESTs from soybean and P. sojae produce distinct, but overlapping, normal distribution curves. When ESTs from the P. sojae-infected soybean tissue are similarly plotted, two separate peaks of percentage G+C content are evident, with each peak closely corresponding to the mean percentage G+C content of the uninfected soybean and axenically grown P. sojae ESTs. This bimodal distribution of ESTs based on percentage G+C content can be used to determine the proportion of P. sojae and soybean transcripts within the infected plant library. By fitting normal distribution curves to each of the two different populations of ESTs present in the infected plant library, and comparing the area covered by each of these curves, we estimate that 60% to 70% of the EST sequences from the infected plant library are derived from P. sojae cDNA transcripts.

View larger version (22K): [in this window] [in a new window]   Figure 3.   Distribution of ESTs on the basis of percentage G+C content. The ESTs were grouped in 1% increments of percentage G+C content and the resulting distributions plotted. Source cDNA library is indicated for each histogram.

Most ESTs Occur as Singletons

Nucleotide sequences were clustered by similarity in each of the libraries to deduce tissue-specific redundancy. Figure 4 shows that the majority of ESTs from all three libraries occurred as singletons. Thus, 680/1,002 (68%) mycelium, 691/1,031 (67%) zoospore, and 674/1,002 (67%) infected plant ESTs did not assemble into contiguous sequences with other ESTs when compared for similarity. Total redundancy was also determined by pooling EST sequences from all three libraries. By this analysis 1,894/3,035 (62%) of the ESTs are singletons. The total EST assembly further enabled us to positively identify many P. sojae transcripts in the infected plant library. Sequences from the mycelium or zoospore library that match ESTs from the infected plant library most likely originate from P. sojae mRNAs.

View larger version (30K): [in this window] [in a new window]   Figure 4.   Assessment of EST redundancy in each of the three cDNA libraries by contiguous sequence assembly analysis. Contigs that consist of one sequence are considered singletons, whereas contigs comprised of two or more sequences are classified as redundant ESTs. Shown is a plot of the frequency of singleton and redundant ESTs for each of the source libraries.

Highly Represented Transcripts Vary among the Three cDNA Libraries

Table I shows a comparison of the most abundantly represented transcripts (excluding ribosomal proteins) from each library. This includes all contiguous sequences of four or more ESTs that resulted from the assembly analysis. Transcripts encoding ribosomal proteins were abundantly represented in all three libraries at levels up to 9-fold (0.9%) redundancy. However, other redundant ESTs were generally specific for a particular library. From the infected soybean library, we conclude that five of the eight redundant ESTs shown in Table I are P. sojae transcripts, and the remaining three represent soybean cDNAs. This was determined by hybridization analysis, or was predicted from the representation, percentage G+C content, and BLAST similarity match of the ESTs comprising the sequence contig. The three abundant soybean transcripts correspond to two different pathogenesis-related proteins, and to a soybean chalcone synthase. Abundant P. sojae transcripts from the infected soybean library include enzymes of intermediary metabolism and a mitochondrial ATP synthase subunit. In contrast, mycelium and especially zoospores showed very different patterns of gene expression. Transcripts matching structural proteins were more prevalent in the mycelium and zoospore ESTs. Almost none (1/12) of the abundant zoospore transcripts shown in Table I were represented in the mycelium or infected plant library.

Comparative Analyses Indicate Varied Patterns of Similarity to Other ESTs

Raw scores from the best TBLASTX match found in each of two targets are plotted in Figure 5, for each of the three cDNA libraries. Several patterns indicate that sequences from P. sojae resemble sequences obtained from P. infestans mycelium to varying degrees. In general, some P. sojae sequences match one target but have no corresponding matches in the other. Those sequences that match both targets tend to resemble the P. infestans sequences more closely than the non-human target, as a majority of points lie below the identity function (dotted line). This pattern is least prominent in the infected soybean library, likely indicating sequences originating from soybean. This is consistent with the observation that highly significant matches occurred most frequently in sequences originating from the infected plant library. The comparative results also indicate cases in which the best score was relatively high in one target but low in the other. Table II summarizes the best examples of the difference between scores from either target being the greatest but neither score being zero, with the best match against the "nr" amino acid target.

View larger version (28K): [in this window] [in a new window]   Figure 5.   Comparative analysis of ESTs derived from each of the source cDNA libraries using TBLASTX, version 2.0.6. P. sojae and soybean EST sequences generated in this study were compared with two distinct data sets, one consisting of P. infestans mycelium ESTs and another, larger data set of non-human ESTs from all taxa, as described in "Materials and Methods." Raw scores from the best match to each target are shown from infected soybean (top), mycelium (middle), and zoospore (bottom) libraries. The black circles indicate the 10 sequences from each library that the had the greatest difference between scores, but where neither score was zero, as shown in Table II.

	DISCUSSION

TOP ABSTRACT INTRODUCTION RESULTS DISCUSSION MATERIALS AND METHODS LITERATURE CITED

We analyzed gene expression patterns in P. sojae by comparing populations of ESTs generated from three different cDNA libraries. We chose to harvest diseased soybean tissue for library construction 48 h after challenge to increase the probability of cloning P. sojae transcripts involved in pathogenesis.

Determining the percentage G+C content of these partially resolved cDNA sequences is a useful indicator in the analysis of ESTs from the infected plant library, since soybean and P. sojae transcripts are different in their average G+C content. This measure provides a means to estimate the total proportion of soybean and P. sojae cDNA transcripts present in the mixed-source library. Individual soybean and P. sojae ESTs may also be distinguished from one another by comparison with other sequence databases. As sequence information accumulates and ESTs approach saturation coverage for soybean and P. sojae, determining the origin of a particular EST from a mixed library may simply be accomplished by a comparative representation analysis.

The high representation (60%-70%) of P. sojae cDNAs present in the infected plant library was surprising, despite the fact that this pathogen is aggressive and fast growing. P. sojae is a considered a hemibiotroph, and therefore it first establishes itself in host tissues as a biotroph but then switches to a more necrotrophic type of growth, rapidly invading and killing host cells. The high proportion of P. sojae ESTs from the infected library may result from this pathogenic strategy, since by 48-h post infection, P. sojae had ramified throughout the hypocotyl tissues and caused large, water-soaked lesions. Further EST analyses at different stages of the infection process and of other host-pathogen interactions would help to put these results in context.

By considering the spectrum and abundance of ESTs represented in the infected library, the present study provides much information about the mechanisms of pathogenesis and of the host defense response. P. sojae genes that appear to be differentially expressed in planta are predominantly associated with intermediary metabolism. The abundant ESTs that match formate dehydrogenase, alcohol dehydrogenase, and glycolytic enzymes such as 2-phosphoglycerate dehydratase suggest that P. sojae relies on glycolysis and mixed alcohol/formic acid fermentation for substrate catabolism and energy generation during growth in host tissues. The oxidative deamination of Glu to -ketoglutarate also seems to be an important route in the degradation and assimilation of amino acid carbon skeletons during pathogenesis, since at least six different Glu dehydrogenase transcripts were represented by 10 ESTs from the infected library. This enzyme plays a key role in nitrogen metabolism in many organisms, because free ammonia is released as a product of catalysis (Garnier et al., 1997).

The high proportion of infection site ESTs that encode enzymes of intermediary metabolism (most of which appear to be P. sojae transcripts) indicates that rapid growth and invasion of host tissues puts a massive demand on central metabolic processes to furnish simple precursors and ATP. Thus, critical steps within these pathways could provide new targets to control the growth of P. sojae and reduce its virulence on soybean plants. In this regard, many other ESTs from the infected library were also potentially significant. For example, several different glucanases and proteinases were also identified from the infection site and zoospore libraries, as shown in Table III. Secreted hydrolytic enzymes are important components that aid in the processes of physical ingress and nutrient solubilization, and thus may constitute quantitative factors that contribute to the overall virulence of the pathogen. Furthermore, it has been well documented that Phytophthora spp. lack a complete sterol synthesis pathway and require an exogenous source of these lipids for normal growth and development (Nes, 1987). A putative progesterone receptor and an estradiol dehydrogenase homolog are thus noteworthy for their possible role in sterol metabolism.

The ESTs from the infection site also offer insight into host plant responses, since an estimated 300 to 400 of these sequences originate from soybean cDNA transcripts. Matches to pathogenesis-related and other defense-related proteins (e.g. peroxidase and polyphenol oxidase) were prevalent among these ESTs, as were matches to enzymes in the phenylpropanoid pathway (e.g. cinnamyl alcohol dehydrogenase, chalcone synthase/reductase, isoflavone reductase). These enzymes produce antimicrobial compounds such as phytoalexins and quinones, and are also involved in the polymerization of phenolic compounds to impede pathogen spread. The isoflavonoid glyceollins are the main phytoalexins produced by soybean via a metabolic pathway of at least 11 different enzyme-catalyzed steps from the precursor coumaroyl-CoA. Molecular characterization of the enzymes and corresponding genes of many of these steps remain to be accomplished. Nonetheless, at least five of the 11 steps in glyceollin synthesis were represented by EST matches, including a cDNA encoding a P450 enzyme (4-1B-HA) that was subsequently characterized and shown to catalyze the aryl migration reaction of isoflavonoid biosynthesis (Steele et al., 1999). Clearly, the infection site ESTs are a rich source of genes involved in host defense responses and offer good opportunities for identifying new enzymes and proteins that participate in these processes.

The collection of zoospore ESTs will also be useful for studies on infection processes and pathogenicity, since zoospores are considered to be the primary infectious propagules of most Phytophthora spp. (Hardham et al., 1991; Erwin and Ribeiro, 1996). Highly represented in the zoospore ESTs were matches to cell structural proteins, including transcripts that may be important for adhesion. Also abundant was a transcript encoding a protein similar to the Bremia lactucae HAM34 protein (Judelson and Michelmore, 1990). The HAM34 mRNA is an abundant transcript in B. lactucea spores, and promoter sequences from this gene have been used in the construction of oomycete transformation vectors (Judelson et al., 1993). Other zoospore ESTs, such as matches to myotubulin and kinesin motor protein, may be involved in motility as components of the prominent tinsel and whiplash flagella.

Mycelium ESTs were more likely to match sequences from the infected plant library than were zoospore ESTs, and were therefore useful to identify many P. sojae transcripts from the mixed infection site. Among the most abundant transcripts in the mycelium library was a match to RIC1, a putative membrane protein from P. infestans that is also similar to stress-induced proteins from yeast (van West et al., 1999). A mycelium superoxide dismutase transcript also present in the infected-plant library is potentially significant because H₂O₂ and other reactive oxygen species generated at infection sites are part of the plant defense response.

The EST data generated in this study represent an important complement to ongoing mapping studies that focus on the genetic localization of avirulence genes by molecular cloning strategies (Tyler et al., 1995; Whisson et al., 1995; Gijzen et al., 1996a). Linkages found between avirulence genes, and sequencing of contiguous stretches of P. sojae genomic DNA have raised the possibility of pathogenicity islands in P. sojae. Thus, positioning ESTs onto physical maps could also provide useful information about gene organization and function in this species.

P. sojae ESTs could also be used for developing new transformation methods. Transient expression assays have shown that vectors constructed using promoters derived from ascomycetes and basidiomycetes were nonfunctional in P. sojae (Judelson et al., 1993). Current P. sojae transformation vectors rely on promoter sequences from Bremia lactucae that may be less than optimal for expression in P. sojae. Promoters from tissue-specific or constitutively present ESTs identified in this study can be isolated from genomic DNA and used for construction of Phytophthora specific transformation vectors.

In summary, the comparative EST-based study presented here represents a resource for genetic and biochemical studies on P. sojae and soybean. We have successfully produced an inventory of cDNA clones and corresponding sequences that will help unravel the underlying mechanisms defining virulence, pathogenicity, and host specificity of P. sojae and defense responses of the host plant.

	MATERIALS AND METHODS

TOP ABSTRACT INTRODUCTION RESULTS DISCUSSION MATERIALS AND METHODS LITERATURE CITED

Culture and Growth Conditions

Phytophthora sojae strain P6497 is a race 2 phenotype originally isolated in Mississippi and obtained from the Phytophthora culture collection at the University of California, Riverside (Förster et al., 1994). Working stocks of the organism were routinely maintained on vegetable juice agar at 25°C in the dark. Sporangial development was induced by repeatedly flooding 5- to 7-d-old mycelium colonies with sterile distilled water. Zoospores were collected by centrifugation at 2,000g (Ward et al., 1979). To obtain axenically prepared mycelium, P6497 was grown for 7 d on synthetic agar media selective for Phytophthora growth (Hoitink and Schmitthenner, 1969). Mycelium discs cut from the growing edge of each colony were transferred into flasks containing 50 mL of liquid synthetic medium and grown for 3 weeks at 25°C in the dark. Liquid cultures were vacuum decanted onto filter paper, and mycelium tissue was collected, frozen in liquid N₂, lyophilized, and stored at 80°C.

Soybean (Glycine max [L.] Merr) cv Harosoy seeds were from the collection at Agriculture and Agri-Food Canada. This cultivar possesses the Rps7 gene for resistance to P. sojae, but it is susceptible to infection by P6497 and most other strains of P. sojae. Etiolated seedlings of cv Harosoy were placed in trays and each hypocotyl inoculated with 10 drops of a zoospore suspension, with each 10-µL drop containing approximately 10³ zoospores (Ward et al., 1979; Gijzen et al., 1996b). Challenged seedlings were incubated for 48 h at 25°C in the dark. Infected, water-soaked tissue was excised, frozen in liquid N₂, lyophilized, and stored at 80°C.

Nucleic Acid Extraction

Freshly collected zoospores, lyophilized mycelium, or infected plant tissues were pulverized to a fine powder in liquid N₂ using a mortar and pestle. Total RNA from zoospores and mycelium was extracted in a solution of phenol-guanidine isothiocyanate (TRIZOL, Life Technologies/Gibco-BRL, Cleveland) according to instructions provided by the manufacturer. Total RNA was isolated from inoculated soybean hypocotyls following the procedure of Wang and Vodkin (1994), because these tissues contained interfering phenolic compounds. P. sojae and soybean genomic DNA was prepared according to the method of Dellaporta et al. (1983) with modifications. Following phenol-chloroform extraction, genomic DNA was incubated with RNAse A for 30 min at 37°C prior to isopropanol precipitation. To remove residual protein and polysaccharide complexes from soybean genomic DNA, samples were further purified by repeated chloroform-isoamyl alcohol (24:1, v/v) extractions in the presence of 0.8 M NaCl and 0.1 volume of 10% (w/v) cetyl-trimethyl-ammonium bromide in 0.7 M NaCl.

cDNA Library Construction

Poly(A⁺) RNA was purified from total RNA by oligo(dT)-cellulose chromatography. The synthesis of cDNA, size selection, addition of linkers, insertional ligation, and packaging into  vector (ZAP Express, Stratagene, La Jolla, CA) followed the manufacturer's instructions and did not deviate significantly from standard methods (Sambrook et al., 1989). The total primary titer of each library in recombinant plaque forming units was: mycelium, 1.1 × 10⁶; zoospore, 7.1 × 10⁶; and infected plant, 9.0 × 10⁵. After amplification, samples of each of the cDNA libraries were used to subclone inserts by mass excision for conversion of  to phagemid vector. The resulting phagemid libraries were plated at low density on Luria Bertani agar plates containing kanamycin (25 mg L¹). Over 1,000 randomly selected bacterial colonies from each of the three cDNA libraries were cultured for plasmid isolation and long-term storage in microtiter plates. Identification codes for ESTs were derived from plate and well numbers and from the cDNA library source (mycelium, MY; zoospore, ZO; cv Harosoy-infected tissue, HA). Sequence data from a seed coat cDNA library (Gijzen, 1997) was used to determine the percentage G+C content of soybean ESTs.

Nucleotide Sequencing

Plasmid DNA was purified from Escherichia coli cultures by alkaline lysis, vacuum filtration, and anion-exchange chromatography using a high-throughput, 96-well format system (Qiagen, Mississauga, Ontario, Canada). Automated cycle sequencing of DNA was carried out using T3 primer and dye-labeled terminators, and products were resolved by electrophoresis through acrylamide gels (model 377, Applied Biosystems, Foster City, CA). All EST data are publicly available through the National Center for Genome Resources (NCGR), coordinating site of the Phytophthora Genome Initiative Project (http://www.ncgr.org/pgi/).

Sequence Data Analysis

Raw DNA sequence data were edited to remove vector sequences and poor quality data, and the percentage G+C content was determined using a computer program (LASERGENE software, DNASTAR, Inc., Madison, WI). Computer-processed sequences were checked manually, compared with electropherograms, and further edited if necessary to improve the quality and reliability of the data. The six deduced amino acid translations of the partial cDNAs were searched against a non-redundant database available through the National Center for Biotechnology Information (http://www.ncbi.nlm.nih.gov/) using the BLASTX 1.4.11 program for comparative identification (Altschul et al., 1990). The probability (P) value was rounded to an order of magnitude and used to classify ESTs as highly significant (P < 10¹⁹), moderately significant (P values from 10⁵-10¹⁹), or weakly significant (P values from 10²-10⁴) matches, where P represents the smallest sum of probability in a Poisson distribution (Adams et al., 1991). ESTs that returned P values greater than 10² were deemed not statistically significant matches.

Assembly of ESTs into Sequence Contigs

The edited EST sequences were arranged into nucleotide-matched clusters using an alignment program (Seqman, LASERGENE, DNASTAR, Inc.) installed in a personal computer (Dell OptiPlex GX1, Pentium II 350 MHz processor, 128 MB RAM) to determine the frequency of sampling redundancy. A pre-pass function optimizes assembly order by constructing a matrix of 16-mer content for each sequence and grouping candidate overlapping ESTs with the greatest number of 16-mers in common. From this data set, sequences were processed into contigs based on the strategy of pair-wise alignment. This approach involves implementation of the Martinez algorithm (Martinez, 1983), which assesses regions of a perfect match between paired sequences, followed by the Needleman-Wunch method, which aligns regions between these matched sites (Needleman and Wunch, 1970). Default parameters (including a 75% minimum match threshold) were used to construct the alignments for each of the databases derived from EST sequences from mycelium, zoospores, or infected plant tissue. For the final assembly, EST sequences from all three libraries were assembled using a minimum match threshold of 65%.

Comparative TBLASTX Analyses

To investigate further the similarity of P. sojae ESTs to other sequences, we performed a comparative analysis against two distinct sets of target sequences: one consisting of 1,490 Phytophthora infestans ESTs from a mycelium library (Kamoun et al., 1999a), and another that comprised 1,429,415 EST sequences from all taxa, excluding humans. Human sequences were excluded because these data comprise more than one-half (55%) of all ESTs and are of limited use for the purposes of a comparative analysis across taxa. Both target sets are accessible publicly at http://seqsim.ncgr.org. We used the TBLASTX algorithm, version 2.0.6 (Altschul et al., 1997), to compare all six reading frames of both query and target nucleotide sequences as amino acid translations. This program returns Expect (E) values in place of P values. Matches with E values < 10⁶ were compared as raw scores from both target sets. Raw scores were the basis for comparison rather than E values, because E values are scaled according to target size, but the two targets differ in size by three orders of magnitude. In this approach, a discrepancy between scores from either target indicates that the P. sojae sequence resembles a sequence from one target more closely than from the other. For the 10 sequences from each library that had the greatest difference between scores, we ran a BLASTX search against the "nr" amino acid library to infer the function of these sequences.

Hybridization Analysis of ESTs to P. sojae and Soybean Genomic DNA

For Southern analysis, 2 µg of P. sojae or 30 µg of soybean genomic DNA was digested with EcoRI endonuclease. Restriction fragments were separated by electrophoresis on a 0.7% (w/v) agarose gel, transferred onto nylon filters, and fixed by UV cross-linkage. The membranes were pre-hybridized in 0.25 M Na₂HPO₄ (pH 7.2), 1% (w/v) bovine serum albumin (BSA), 1 mM EDTA, and 7% (w/v) SDS at 65°C for 5 h. Probes were ³²P-labeled by random primer extension according to the manufacturer's instructions (Pharmacia, Uppsala). After an 18-h hybridization at 65°C, the filters were washed in 20 mM Na₂HPO₄ (pH 7.2), 1 mM EDTA, 1% (w/v) SDS at 68°C for four 20-min incubations, followed by three 5-min washes in the same buffer at room temperature. The blots were exposed to Kodak X-AR film (Eastman Kodak, Rochester, NY) with an intensifying screen at 70°C.