Skip to main content
  • Research article
  • Open access
  • Published:

Survey and analysis of microsatellites from transcript sequences in Phytophthora species: frequency, distribution, and potential as markers for the genus

Abstract

Background

Members of the genus Phytophthora are notorious pathogens with world-wide distribution. The most devastating species include P. infestans, P. ramorum and P. sojae. In order to develop molecular methods for routinely characterizing their populations and to gain a better insight into the organization and evolution of their genomes, we used an in silico approach to survey and compare simple sequence repeats (SSRs) in transcript sequences from these three species. We compared the occurrence, relative abundance, relative density and cross-species transferability of the SSRs in these oomycetes.

Results

The number of SSRs in oomycetes transcribed sequences is low and long SSRs are rare. The in silico transferability of SSRs among the Phytophthora species was analyzed for all sets generated, and primers were selected on the basis of similarity as possible candidates for transferability to other Phytophthora species. Sequences encoding putative pathogenicity factors from all three Phytophthora species were also surveyed for presence of SSRs. However, no correlation between gene function and SSR abundance was observed. The SSR survey results, and the primer pairs designed for all SSRs from the three species, were deposited in a public database.

Conclusion

In all cases the most common SSRs were trinucleotide repeat units with low repeat numbers. A proportion (7.5%) of primers could be transferred with 90% similarity between at least two species of Phytophthora. This information represents a valuable source of molecular markers for use in population genetics, genetic mapping and strain fingerprinting studies of oomycetes, and illustrates how genomic databases can be exploited to generate data-mining filters for SSRs before experimental validation.

Background

Phytophthora spp. are notorious, world-wide pathogens because of their devastating effects on many crop species, that often result in significant economic losses. All members of the genus Phytophthora infect plants, although some exhibit a broad host range and others infect only a few species [1]. Among the most important species are P. infestans, P. sojae and the newly described species, P. ramorum [2]. P. infestans affects several host plants from the Solanaceae family including potato, tomato and a number of tropical fruits economically important for producer countries [3, 4]. Current losses in potato production worldwide reach $5 billion [5]. The soybean pathogen, P. sojae, causes root rot and damping off, and in the US alone has recently accounted for annual losses in the order of $1–2 million [1]. P. ramorum has only been described relatively recently, so that basic information about its biology and ecology is scant. This pathogen causes sudden oak death disease and is currently decimating trees and shrubs in the coastal oak forest in California, including keystone tanoak and coast live oak species, and might be expanding to other hosts such as redwoods and to other regions in North America [6].

Given the economic relevance of these pathogens, it is important to standardize high-throughput molecular methods for routinely characterizing Phytophthora populations. Genomic and bioinformatics resources have grown exponentially in recent years, generating information more rapidly than current data-processing tools can handle. Thus, the generation of new tools and the application of existing ones for exploring databases constitute a practical and inexpensive approach to elucidating biological systems. The application of bioinformatic tools to genomic databases that are already available for different species of Phytophthora could allow techniques for characterizing these pathogens to be developed rapidly.

Two of the most informative databases available for Phytophthora are the EST (Expressed Sequence Tags) public databases and the Department of Energy's Joint Genome Institute (JGI) databases in which complete genomes, predicted proteins and predicted transcripts are available for P. sojae and P. ramorum [7]. Transcript sequences constitute a rich and special source of informative molecular markers because they represent genes that are expressed in an organism. Coding sequences are generally more informative than anonymous markers because they allow for a more direct association between the molecular marker and the phenotype. Microsatellites or simple sequence repeats (SSRs) are molecular markers that consist of tandem repeats of one to six DNA base pairs. SSRs are highly versatile, PCR-based markers, usually associated with a high frequency of length polymorphism [8]. They are choice markers given that they show fairly high mutation rates and are codominant [9, 10]. They have been found in both coding and non-coding DNA sequences of all higher organisms analyzed [11, 12]. In oomycetes, SSRs have had important applications such as diagnosis and determination of mating type [13], genetic structure and disease dynamics [14, 15], and population genetics [16, 17]. Apart from their application as molecular markers, determining the abundance and density of SSRs in oomycetes may help understand whether these sequences have any functional and evolutionary significance [18].

An innovative marker system that has been developed links expressed sequence tags (ESTs) and SSRs [19, 20]. These EST-SSRs have been applied successfully in studies of genetic variation, linkage mapping, gene tagging, evolution and sequencing of several plant genomes [21]. Although they are less polymorphic than genomic SSRs, EST-SSRs tend to be more conserved, at least in the plant species in which they have been studied [22]. This characteristic makes EST-SSRs readily transferable between related organisms [17, 2326]. Therefore, SSRs from transcript sequences have considerable potential for comparative mapping studies, as well as for analyses of genetic diversity within the expressed portions of the genome in which they are located.

Although Phytophthora EST and transcript databases are publicly available, no formal analysis of SSRs in these sequences has been reported. We used an in silico approach to analyze the frequency and distribution of SSRs in transcript sequences from the oomycetes Phytophthora infestans, P. sojae and P. ramorum. Previous studies on oomycete phylogeny have suggested that these three species are monophyletic and that P. sojae and P. ramorum are more closely related to each other than to P. infestans [2729]. We also surveyed the distribution and possible patterns of SSRs in selected sequences corresponding to genes previously associated with pathogenesis and virulence. In addition, we studied in silico the transferability of these SSR-based markers between species. We generated primer pairs, where possible, for all SSRs from these three organisms. A publicly available database was generated for Phytophthora microsatellites, including primers and SSR survey results [30]. This study will serve as reference for future comparative mapping studies and for the development of strategies that take advantage of DNA sequence analyses for cross-referencing genes between species and perhaps genera.

Results

Microsatellites: motif, length and frequency

Consensus EST databases from P. infestans and annotated transcripts from P. sojae and P. ramorum, were scanned for the presence of microsatellites, defined as short tandem repeat motifs of 1–6 bp. We also explored the existence of 7–10 bp motifs, which represent the transition between micro and minisatellites; little is known about these. A total of 84000 available EST sequences were downloaded from the Phytophthora Functional Genomics Database [31] for P. infestans, and 19276 and 16066 predicted and annotated transcripts were downloaded from the Department of Energy's Joint Genome Institute [32]. The EST set was masked for repetitive sequences, obtained from RepBase ([33]; which does not include SSR repeats), assembled and revised manually to generate consensus sequences. We used consensus EST sequences because they have the built-in advantage of eliminating redundant SSR counts, allowing us to make more precise estimates of SSR frequency. The assembly process resulted in 25965 sequences for P. infestans. Annotated transcripts were not subjected to further treatment, and along with the consensus ESTs were used for the SSR survey.

The frequency of repeat motifs in the consensus EST sequences and annotated transcripts was assessed. A first analysis was performed with only the consensus ESTs from P. infestans; more than 50% of the SSRs identified were (A)n , which is an over-representation with respect to the other mononucleotide repeats as well as to other motifs. Although mononucleotide repeats are common in genomic DNA and can be valid SSRs, most of those present in expressed sequences are the result of nucleotide additions by RNA polymerase and are not present in the genomic DNA template (e.g. poly-A tails; [34]). Therefore, the analysis of monomers was excluded and only motifs with repeats of 2 to 6 bp (for ESTs and transcripts) were included in this study. Both perfect and compound SSRs were selected with a minimum acceptable length of 12 bp for di, tri and tetra-nucleotide motifs. Only SSRs with a minimum of three repeats were included in the analyses of penta- and hexa-nucleotide repeats. The total counts, frequencies and comparisons of SSRs (2–6 bp) in each set of sequences are summarized in Tables 1, 2, 3.

Table 1 Number and distribution of SSRs in consensus EST and transcript sequences.
Table 2 Percentage, relative abundance and relative density of SSRs in the P. infestans, P. sojae, and P. ramorum sequence sets.
Table 3 Most common repeat motifs identified from perfect and compound microsatellites in the three oomycetes analyzed.

P. infestans had the most sequences analyzed but showed the lowest total SSR count and percentage of SSR-containing sequences (Table 1). However, this result has to be considered wih caution since the total size (Mb) of the sequences examined for P. infestans is approximaely half that of the other organisms. To compare the organisms more realistically, another approach was required: namely, taking the total length of each set of sequences analyzed as a reference. Thus, total relative abundance and total relative density were calculated (Table 1). Statistically significant differences (P < 0.05) were found between P. infestans and P. sojae and between P. sojae and P. ramorum when we compared the total sequence lengths (bp) contributed by SSRs with respect to the total megabases of examined sequences (relative density), but there was no difference between P. infestans and P. ramorum. The same differences were found when perfect and compound SSR proportions were compared (data not shown), demonstrating that SSR content is not in agreement with the phylogenetic distances between these organisms [6]. The numbers of SSR-containing sequences and sequences containing more than one SSR were compared only between P. sojae and P. ramorum because of the different types of source sequences (ESTs for P. infestans and transcripts for P. sojae and P. ramorum). All comparisons showed statistically significant differences (data not shown). These results might imply that the net SSR content in transcript regions of Phytophthora species could be variable, and again, not directly related to phylogenetic distance.

The total numbers of all types of microsatellite motifs are shown in Table 2 and Figures 1 and 2. All three sequence sets contained SSRs that were mainly trinucleotide repeats (> 60%), while the dinucleotide repeats represented less than 8%. This agrees with results from other eukaryotes, where trinucleotide repeats are overrepresented in coding sequences in comparison with dinucleotide repeats [35, 36]. Hexanucleotide repeats constituted the second most frequent motif in P. sojae and P. ramorum, differing markedly from P. infestans, where this motif has one of the lowest percentages. In general, statistically significant differences were not found between P. sojae and P. ramorum (proportions test, P > 0.05), showing that in contrast to the net SSR content, SSR distributions could be related to phylogenetic distances among Phytophthora species. Relative abundance and relative density allowed the similarities and differences among SSR distributions to be represented graphically (Figures 1 and 2). The high relative abundance and density values observed for tri- and hexa-nucleotides might correlate with coding region stability ([37]; see discussion).

Figure 1
figure 1

Relative abundance of SSRs among oomycete consensus ESTs and transcript sequence sets in the oomycetes analyzed. Abundance is defined as the total number of SSRs per Mb of sequence analyzed.

Figure 2
figure 2

Relative density of SSRs across consensus oomycete ESTs and transcript sequence sets. Density is defined as the total sequence length (bp) contributed by each SSR per Mb of sequence analyzed.

The fifteen most frequent motifs were analyzed for all three organisms in terms of percentage and total counts (Table 3). Most of the motif types were present, and the five most frequent motifs were the same, in all three Phytophthora species; however, more motifs were shared between P. sojae and P. ramorum (~80%) than between these two species and P. infestans (~46%). The most common trinucleotide repeats in all cases were (AGC/CGT), (ACG/CTG) and (AAG/CTT); the most common dinucleotide repeat for P. sojae and P. ramorum was (CG/GC), in complete contrast to P. infestans, where this motif was least frequent. Tetra-, penta- and hexa-nucleotide repeat motifs showed no clear trend among the three organisms. In general, these results and those from the SSR distributions, reflect the more close relation between P. sojae and P. ramorum at the sequence level.

Since the most common motifs in the three organisms analyzed were trinucleotide repeats, we attempted to identify the amino acid(s) encoded by these. Twenty sequences from each organism, containing the most common triplet in each case, were randomly selected and used for ORF analyses. The most probable open reading frame and consequently the location of the SSR in this reading frame were determined. The triplets analyzed (in their canonical forms and possible variations) were (AGC)n for P. sojae and P. ramorum and (AAG)n for P. infestans respectively. In P. sojae and P. ramorum, CAG (glutamine) was predominant at 80% and 60% respectively, followed by AGC (serine) at 10% and 35%, and finally GCA (alanine) at 10% and 5%. In P. infestans, GAG (glutamic acid) was predominant at 50%, followed by AAG (lysine) at 30% and finally CTT (leucine) at 10%. These results were further investigated for correlations between the amino acids encoded by the trinucleotide repetitions and the codon usage preferences reported for each organism in the Codon Usage Database from the Kazuka DNA Research Institute [38]. In all cases, the triplets analyzed encoded for amino acids that were normally overrepresented in the corresponding organism. The differences detected between P. infestans and P. sojae-P.ramorum emerging in our trinucleotide analysis were also corroborated by codon usage frequencies.

Repeats containing motifs between 7 and 14 bp were scarce in the three oomycetes transcriptome analyzed, accounting for less than 1% of the total number of SSRs. They did not merit further attention in our study because of their low abundance. Regarding SSR lengths, we found that di-, tetra- and penta-nucleotide motif types did not exceed 30 bp (Table 4), while tri and hexanucleotides were clearly longer in all organisms. Although P. sojae had the longest SSR for di-, tetra-, penta- and hexa-nucleotide repeat types, longer than reported for other organisms (Karaoglu et al., 2004), only a very few of these lengths can be considered "long microsatellites" (> 15 repetitions [39]). The motif sequences of the longest SSRs are not shared among organisms, indicating that this factor is independent in each species. Thus, our results suggest that long SSRs are absent from these consensus EST or transcripts sequences. The number of repeats found in SSR loci ranged from 3 to 33, from 3 to 32 and from 3 to 15 in P. infestans, P. sojae and P. ramorum, respectively. Most SSR loci showed seven repeats or less (96% P. infestans, 98% P. sojae and 99%P. ramorum), with a repeat number of four being the most common in all species (Fig. 3).

Table 4 The longest SSR motifs found in EST consensus and transcript sequences for the three oomycetes analyzed.
Figure 3
figure 3

Frequency of repeat unit numbers of SSRs from P. infestans (black),P. sojae (grey) and P. ramorum (white).

SSRs in pathogenicity factors

In total, 136, 318 and 171 sequences corresponding to pathogenicity factors or annotated as putatively involved in pathogenicity were selected for P. infestans, P. sojae and P. ramorum, respectively, and their SSR distributions were characterized. They included enzymes such as cell wall degrading enzymes (cutinases, glucanses, polygalacturonases, pectate lyases and cellulases), elicitins and avirulence homolog proteins with conserved RXLR motifs, and other secreted proteins potentially related to pathogenicity [7, 40]. SSRs were also surveyed in two additional sequence sets corresponding to constitutively expressed genes: ribosomal and housekeeping genes such as actin, cytochrome P450-like protein and NADH hydrogenase. The results showed that the percentage of SSR-containing sequences in the pathogenicity factors was not significantly different from that in the ribosomal and housekeeping genes (Table 5). The scarcity of SSRs in the putative pathogenicity factors suggested that SSR length variation has no significant influence on the mutation rate in these sequences.

Table 5 Distribution of SSRs in three functional categories of genes

Primer design for EST-SSRs and databases

Sequences flanking microsatellites from each of the three organisms were used to develop primer pairs using Primer3 software [41]. For this first survey, the parameters were not stringent since a high number of sequences were analyzed. We produced primer pairs for 61.44%, 87.47% and 88.97% of the SSRs from P. infestans, P. sojae and P. ramorum respectively; not all the SSRs were located in positions suitable for optimum design. Three different primer pairs were generated for each SSR and deposited in the developed database [30] with their sequences, the consensus EST or transcript sequence (with the original sequence ID) and a brief description of the amplifiable SSR. These primers can be used to amplify the corresponding SSR region for diverse applications, so they constitute a publicly available resource for future research. Of the three primer pairs, only the first (the best score) was used to analyze transferability.

in silico analysis of transferability

Primer pairs for each SSR locus were assayed in silico for cross-transferability. Primers for each organism were aligned against the consensus ESTs or transcript databases from the other two organisms. Three criteria were established to filter the comparisons: (i) high degree of similarity (> 90%) between primer and aligned sequence; (ii) primers aligning with only one sequence containing an SSR; (iii) a hypothetical PCR product size longer than 100 bp. Of all the primers designed (8739), 7.5% appeared transferable between at least two species. Not surprisingly, most of the virtually transferable primers were found between P. sojae and P. ramorum, since these are the two most closely related species. However, 84 (~1 %) primers were found to be virtually transferable among all three species (Table 6).

Table 6 In silico determination of potential cross-transferability of the SSRs between Phytophthora species

Discussion

The present study was designed to create microsatellite databases for P. infestans, P. sojae and P. ramorum, taking advantage of publicly available sequences for these organisms. In the case of P. infestans, EST sequences were first assembled and then explored for SSRs. For P. ramorum and P. sojae, annotated transcripts from the available whole genome sequences were used directly to search for SSRs and to analyze the distribution and organization of SSRs in the transcribed regions of these organisms. Approximately 6–22% of the sequences contained SSRs, which shows higher frequencies than previously reported for plant ESTs [25, 37, 42], fungal endophytes [43] or other higher eukaryotes [8, 44]. However, these differences might reflect the different criteria used to select the SSRs. Repeat sequences of at least 12 bp for di-, tri- and tetra-nucleotides, and three or more repeated units for penta- and hexa-nucleotides, were chosen for this study. These lengths were used because they have been considered the minimum acceptable microsatellite lengths [34] and are efficient thresholds for detecting high levels of polymorphism [22]. Our results showed that SSR lengths are very restricted in the coding regions of oomycete genomes; approximately 99% of all the SSRs analyzed were shorter than 30 bp. Only a few SSRs had higher numbers of repeat units, perhaps because of their location in well-conserved regions of the genome. Strong evolutionary and functional constraints limit the expansion of microsatellite repeats in expressed regions of the genome [45, 46], because longer repeats have higher mutation rates and could therefore be less stable [47, 48]. Short microsatellites are probably generated by random mutations and then expanded by DNA polymerase slippage. Thus, the base composition of a sequence that seeds the evolution of repeats is expected to influence microsatellite density [49, 50]. Therefore, the similarities between SSR motifs within oomycetes and between oomycetes and fungi may indicate that specific common sequence composition contributes to the evolution of SSRs.

Statistically, there were no quantitative differences in the distributions of SSR motifs between P. sojae and P. ramorum. However, marked differences were found when these two species were compared with P. infestans. This is particularly interesting because the frequency distribution of SSRs does not follow a genus-wide pattern; if it is strongly species-dependent, it could indicate evolutionary events specific to the organisms compared. In addition, SSRs derived from exonic regions, which are more conserved than genomic SSRs, might consequently show only minor differences in distribution among related species. Thus, differences in genomic organization could explain the SSR distributions observed in the organisms examined, reflecting the phylogenetic distances between them [29, 51]. More sequence data on ESTs/transcripts will soon become available for Stramenopiles, as P. infestans, P. capsici, Pythium ultimum, and Hyaloperonospora parasitica are or will be sequenced in the near future that will provide further resolution on the evolution of SSR motifs in Oomycete coding regions.

The most common SSRs comprised trinucleotide repeat units with low repeat numbers. A wide variety of repeat motifs were represented at high percentages in these trinucleotide arrays. The abundance of repeat motifs differed slightly, especially between P. infestans and the other two species; (AGC)n, (ACG)n and (AGG)n were the most abundant triplets in all organisms but their abundances differed among species. This finding was expected, since EST or transcript-derived microsatellites are likely to be conserved in frequency, abundance and distribution across closely-related species [52]. A database search of all possible trinucleotide repeat motifs (>20 bp) showed that (AGG)n, (AAT)n and (ATC)n are relatively common in fungi, but (ACG)n and (CCG)n are relatively rare [53]. Differences in abundance and density among trimeric repeats could be explained by species-specific cellular factors that interact with the motifs and play an important role in generating the repeats [18]. Among plant species, the abundances of different repeat motifs in EST-derived SSRs vary greatly. However, trinucleotide units with low repeat numbers are common features of EST-SSRs [25, 37, 42, 54]. Our results suggest that Phytophthora spp. might have a set of common motifs, as is the case with fungi, whereas the motifs in oomycetes may vary widely, as they do in plants.

High dinucleotide repeat abundances in whole genomes have been reported for fungi, Drosophila,Caenorhabditis elegans and a subset of plant genomes [36, 51, 54, 55]. Dinucleotide repeats have been characterized as being the most important SSRs because of their higher mutation rates [51]. This suggests an explanation for their high abundance in genomic regions and low abundance in coding regions, which must be conserved to maintain functionality. On the other hand, many studies have reported that trinucleotides are most abundant in coding regions of higher eukaryotic genomes [46, 56, 57]. Previous studies have shown that trinucleotide repeats predominate in plant EST libraries, supporting our observations on oomycetes [12, 58]. Among all SSRs, expansions or deletions in coding regions can be tolerated for tri- and hexa-nucleotides, which do not perturb reading frames [36]. This could explain why the longest SSRs in all the organisms examined belong to these motif categories, although the three-dimensional structure of a protein translated from a sequence containing such SSRs would not necessarily be unchanged. Interestingly, our analyses showed that glutamine is the amino acid encoded by the most highly abundant repeat in P. sojae and P. ramorum: (CAG)n . This agrees with previous reports on Drosophila, C. elegans and yeast [36], in which the same glutamine-coding triplet was the most common in tandem triplet repeats. In contrast, glutamic acid was encoded by the most common repeat in P. infestans: (AAG)n . Like glutamine, this is also one of the most common tandemly repeated amino acids in the aforementioned organisms. Katti et al. [59] analyzed all the protein sequences from the SWISS-PROT database for single amino acid repeats, tandem oligo-peptide repeats, and periodically conserved amino acids; the results showed that repeats of glutamine, serine, glutamic acid, glycine and alanine seem to be fairly well tolerated in many proteins. Four of these five amino acids were found in our analysis and have also been reported in recent studies of complete genome coding sequences [36]. Such triplet repeat patterns in ORFs of oomycetes and other organisms could reflect functional selection of amino acid reiterations in the encoded proteins. Whole genome analyses have shown that repeat stretches of small/hydrophilic amino acids are more frequent in proteins [59]. The expansion of codons encoding such amino acids might be better tolerated than the expansion of hydrophobic amino acid stretches because they probably would not change the three-dimensional protein conformation as drastically. Therefore, nucleotide composition might strongly affect the structures and functions of encoded proteins, and it could be a determining force in the selection of SSRs in coding sequences.

We also explored the possibility that valuable information for mapping and for diversity or population structure studies could be obtained from these data. For this reason, microsatellite markers were evaluated for in silico cross-transferability among all three species. Many studies on organisms from different kingdoms have led to the development of markers amplifiable across species [34, 44, 54], even among Phytophthora species [17]. To our knowledge, this possibility has never been assessed in silico for any organism (including oomycetes). SSRs designed from EST/transcripts sequences are especially valuable owing to their genome location, which implies constraints on length, motif, abundance and flanking regions. The last of these is of particular interest for characterizing interspecies genetic diversity using common primers. Therefore, we designed primers for all SSR sets observed by species and identified those that were potentially cross-transferable. This approach was similar to a virtual PCR and provided candidate SSRs for experimental screening. The results revealed a small proportion of primers that could theoretically be used as transferable molecular markers (using a similarity criterion of 100%). This was an unexpected result because EST/transcript types of sequence are highly conserved and there is a close phylogenetic relationship among the Phytophthora spp. [27]. However, if the similarity criterion for alignment was less stringent (90%), 653 primer pairs were virtually transferable between at least two species. This number represents a remarkably large, diverse reservoir of markers that is potentially useful for diversity and population studies. With far fewer primers, an acceptable level of polymorphism has been found in previous studies of different organisms, where the primers are used to amplify a DNA region that has undergone SSR expansion in one lineage but not in related ones [8, 54, 60]. The transferability results might suggest that the organization of oomycete genomes shows marked variations. However, a more plausible explanation could be that the evolutionary distance between the selected taxa is greater than in other systems or organisms in which SSRs have been shown to be transferable (e.g. plants). Experimental validation of these hypothetically transferable SSRs and their polymorphism is on-going. This will verify the potential effectiveness of this in silico tool for finding transferable molecular markers for evaluating intra- and inter-specific diversity, instead of spending resources to validate them.

Evolutionary forces keep pathogen lineages and their hosts in an arms race, as each evolves new strategies for attack and defense [61]. The success of this arms race over time depends on the mechanisms by which pathogenicity factors evolve. Microsatellite mutation rates, ranging from 10-6 to 10-2 per generation, are higher than base substitution rates [10], and it is thus reasonable to assume that their presence in pathogenicity factors may contribute to variability in the relevant sequences [62]. In this study, we surveyed SSR abundance and distribution in putative pathogenicity factor sequences. The results were compared with SSR abundance and distribution in ribosomal and housekeeping sequences, which are presumed to be less variable because their fundamental functions in the organism are conserved. However, our results did not show higher SSR abundance in the pathogenicity factors. Differences from reference sequences (ribosomal and housekeeping) were not statistically significant (P > 0.05). This suggests that SSRs probably do not make a clear contribution to functional changes in pathogenicity factors, which allow P. infestans, P. sojae and P. ramorum to cause disease, but our study does not rule out the possibility that SSR polymorphisms may be present in pathogenicity genes.

Conclusion

Databases for P. sojae, P. infestans and P. ramorum were inexpensive sources of SSRs. In silico investigation of SSRs is a significant step towards understanding the biological functions and nature of these important portions of the DNA. SSR markers are potentially useful tools for identifying Phytophthora species and for assessing genetic variation among populations. Evaluation of intraspecific genetic variation in Phytophthora concomitant with equivalent studies of host genetic variation could allow the relationships between host genotype variability, Phytophthora genotype and the selection pressures acting between them to be established. An important direct application of this study is as a resource for future population research.

Methods

Sources of ESTs and annotated transcript sequences

EST sequences from P. infestans (81974) were downloaded from the Phytophthora Functional Genomics Database [31]. The annotated transcript sequences available for P. sojae (19276) and P. ramorum (16066) were downloaded from the Department of Energy's Joint Genome Institute [32]. These redundant ESTs and transcript sequence sets were obtained in FASTA format. In addition, we downloaded 136, 318, and 171 sequences corresponding to putative or validated pathogenicity factors from P. infestans (National Center for Biotechnology Information [63]), P. sojae (Virginia Bioinformatics Institute [64]) and P. ramorum (DOE Joint Genome Institute [32]). From the same databases, ribosomal and housekeeping gene sequences were obtained for each of the organisms studied.

Masking and assembly of the ESTs sequences from P. infestans

UniVec build # 3.2 was obtained from the NCBI web site FTP directory [65] and the latest repetitions database RepBase [33] EST were masked using the cross-match utility. PHRAP version 0.990329 was used to cluster the ESTs and to generate consensus sequences. Parameters phrap-miniscore 100 -minmatch 50 were used to generate consensus. The total consensus sequence count was 25965.

Scanning of non-redundant EST and transcript sequences from P. infestans, P. sojae and P. ramorum for SSR survey

Non-redundant EST and transcript sequences were scanned using a local version of the Microsatellite Identification Tool (MISA) available from the Plant Genome Resources Center (PGRC) [66]. This program searches for both perfect and compound SSRs with 2 to 6 nucleotides in the basic repeat unit. It records repeat numbers and SSR locations inside the EST sequence and deposits these results in an output file. SSR redundancy was minimized by counting only a single match when there was more than one record for the same SSR locus. Minimum SSR length was determined by the number of repeats, which were (2/6) (3/4) (4/3) (5/3) (6/3); the first digit refers to the SSR repeat type and the second to the minimum number of repeat units (Table 1). For example, (2/6) means that a motif consisting of 2 bp and should have at least 6 repeats to be considered a microsatellite in our analysis. Pathogenicity factors, ribosomal genes and housekeeping genes were also scanned for SSRs using the same parameters (Table 6).

Total SSR numbers were normalized by calculating relative abundance. This was calculated as the number of SSRs per Mb of sequence analyzed so that the three sequence sets of different sizes could be compared. The relative density (bp/Mb) of each set of ESTs sequences was calculated by dividing the number of base pairs in the sequences (bp) contributed by each SSR by the total length of sequences examined (Mb); see Figures 1 and 2

Primer design and cross-transferability

P3_in.pl and P3_out.pl perl scripts, which complement the MISA program, were adapted and used for automated selection and transfer of SSR-containing sequences from the database to the Primer3 program. Parameters used for the Primer3 program were as follows: optimal Tm of 57°C with a minimum and maximum of 50°C and 64°C respectively, and 50–75% GC content. The probability of dimer or hair-loop formation was low. The size of the PCR product is expected to be between 50 and 300 bp and no secondary structure.

To search for cross-transferability between paired P. infestans-P. sojae, P. infestans- P. ramorum and P. ramorum- P. sojae, each designed primer pair from each species was compared against all the SSR-containing sequences from the other species. The EMBOSS 3.0 Primersearch utility was used for this purpose. Cross-transferability among the three species was assessed manually.

In order to track all the data generated during this research, a MySQL relational database was designed and implemented. This database can be accessed at the URL listed in the references section [30]

Statistical analysis

A proportion test was performed to evaluate differences between variables [67]. All statistical tests were performed using Statistix 8.0.

Abbreviations

EST:

Expressed Sequence Tag

SSR:

Simple Sequence Repeat

EST-SSR:

EST-derived SSR

References

  1. Erwin DC, Ribeiro OK: Phytophthora diseases worldwide. 1996, St. Paul , APS Press

    Google Scholar 

  2. Werres S, Marwitz R, Man In´T Veld WA, e Cock AWAM, Bonants PJM, De Weerdt M, Themann K, Ilieva E, Baayen RP: Phytophthora ramorum sp. Nov., a new pathogen on rhododendron and viburnum. Mycol Res. 2001, 105: 1155-1165.

    Article  CAS  Google Scholar 

  3. Fry WE, Goodwin SB: Resurgence of the irish potato famine fungus. BioScience. 1997, 47 (6): 363-371. 10.2307/1313151.

    Article  Google Scholar 

  4. Grunwald NJ, Flier WG: The biology of Phytophthora infestans at its center of origin. Annu Rev Phytopathol. 2005, 43: 10.1-10.2. 10.1146/annurev.phyto.43.040204.135906.

    Article  Google Scholar 

  5. Duncan JM: Phytophthora - an abiding threat to our crops. Microbiology Today. 1999, 26: 114-116.

    Google Scholar 

  6. Rizzo DM, Garbelotto M, Hansen EM: Phytophthora ramorum: integrative research and management of an emerging pathogen in California and Oregon forests. Annu Rev Phytopathol. 2005, 43: 309-335. 10.1146/annurev.phyto.42.040803.140418.

    Article  PubMed  Google Scholar 

  7. Tyler BM, Tripathy S, Zhang X, Dehal P, Jiang RHY, Aerts A, Arredondo F, Baxter L, Damasceno CMB, Dickerman A, Dorrance AE, Dou D, Dubchak I, Garbelotto M, Gijzen M, Gordon S, Govers F, Grunwald N, Huang W, Ivors K, Jones RW, Kamoun S, Krampis K, Lamour K, Lee MK, McDonald WH, Medina M, Meijer HJG, Nordberg E, MacLean DJ, Ospina-Giraldo MD, Morris P, Phuntumart V, Putnam N, Rash S, Rose JKC, Sakihama Y, Salamov AA, Savidor A, Scheuring C, Smith B, Sobral BWS, Terry A, Torto-Alalibo T, Win J, Xu Z, Zhang H, Grigoriev I, Rokhsar D, Boore J: Phytophthora genome sequences uncover evolutionary origins and mechanisms of pathogenesis. Science. 2006, In Press:

    Google Scholar 

  8. Zhan A, Bao ZM, Wang XL, Hu JJ: Microsatellite markers derived from bay scallop Argopecten irradians expressed sequence tags. Fisheries Science. 2005, 71: 1341-1346. 10.1111/j.1444-2906.2005.01100.x.

    Article  CAS  Google Scholar 

  9. Schlotterer C: The evolution of molecular markers--just a matter of fashion?. Nat Rev Genet. 2004, 5 (1): 63-69. 10.1038/nrg1249.

    Article  PubMed  Google Scholar 

  10. Schlotterer C: Evolutionary dynamics of microsatellite DNA. Chromosoma. 2000, 109 (6): 365-371.

    Article  PubMed  CAS  Google Scholar 

  11. Tautz D, Renz M: Simple sequences are ubiquitous repetitive components of eukaryotic genomes. Nucleic Acids Res. 1984, 12 (10): 4127-4138.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  12. Gupta PK, Balyan HS, Sharma PC, Ramesh B: Microsatellites in plants: a new class of molecular markers. Current Science. 1996, 70: 45-54.

    CAS  Google Scholar 

  13. Prospero S, Black JA, Winton ML: Isolation and characterization of microsatellite markers in Phytophthora ramorum, the causal agent of sudden oak death. Molecular Ecology. 2004, 4: 672-674. 10.1111/j.1471-8286.2004.00778.x.

    Article  CAS  Google Scholar 

  14. Gobbin D, Pertot I, Gessler C: Identification of microsatellites markers for Plasmopara viticola and establishment of high-throughput method for SSR analysis. European journal of Plant Pathology. 2003, 109: 153-164. 10.1023/A:1022565405974.

    Article  CAS  Google Scholar 

  15. Knapova G, Gisi U: Phenotypic and genotypic structure of Phytophthora infestans populations on potato and tomato in France and Switzerland. Plant Pathology. 2002, 51: 641-653. 10.1046/j.1365-3059.2002.00750.x.

    Article  Google Scholar 

  16. Ivors K, Garbelotto M, Vries ID, Ruyter-Spira C, Te Hekkert B, Rosenzweig N, Bonants P: Microsatellite markers identify three lineages of Phytophthora ramorum in US nurseries, yet single lineages in US forest and European nursery populations. Mol Ecol. 2006, 15 (6): 1493-1505. 10.1111/j.1365-294X.2006.02864.x.

    Article  PubMed  CAS  Google Scholar 

  17. Lees AK, Wattier R, Shaw DS, Sullivan L, Williams NA, Cooke DEL: Novel microsatellite markers for the analysis of Phytophthora infestans populations. Plant Pathology. 2006, 55: 311-319. 10.1111/j.1365-3059.2006.01359.x.

    Article  CAS  Google Scholar 

  18. Toth G, Gaspari Z, Jurka J: Microsatellites in different eukaryotic genomes: survey and analysis. Genome Res. 2000, 10 (7): 967-981. 10.1101/gr.10.7.967.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  19. Kantety RV, La Rota M, Matthews DE, Sorrels ME: Data mining for simple sequence repeats in expressed sequence tags from barley, maize, rice, sorghum and wheat. Plant Molecular Biology. 2002, 48: 501-510. 10.1023/A:1014875206165.

    Article  PubMed  CAS  Google Scholar 

  20. Smulders MJM, Bredemeijer G, Rus-Koretekaas W, Arens P: Use of short-microsatellites from database sequences to generate polymorphism among Lycopersicon species. Theor Appl Genet. 1997, 97: 264-272. 10.1007/s001220050409.

    Article  Google Scholar 

  21. Qureshi SN, Saha S, Kantety RV, Jenkins JN: EST-SSR: a new class of genetic markers in cotton. Journal of Cotton Science. 2004, 8: 112-123.

    CAS  Google Scholar 

  22. Coulibaly I, Gharbi K, Danzmann RG, Yao J, Rexroad CE: Characterization and comparison of microsatellites derived from repeat-enriched libraries and expressed sequence tags. Anim Genet. 2005, 36 (4): 309-315. 10.1111/j.1365-2052.2005.01305.x.

    Article  PubMed  CAS  Google Scholar 

  23. Herron BJ, Silva GH, Flaherty L: Putative assignment of ESTs to the genetic map by use of the SSLP database. Mamm Genome. 1998, 9 (12): 1072-1074. 10.1007/s003359900929.

    Article  PubMed  CAS  Google Scholar 

  24. Eujayl I, Sorrels ME, Baum M, Wolters P, Powell W: Isolation of EST-derived microsatellite markers for genotyping the A and B genomes of wheat. Theor Appl Genet. 2002, 104: 399-407. 10.1007/s001220100738.

    Article  PubMed  CAS  Google Scholar 

  25. Scott KD, Eggler P, Seaton G, Rossetto M, Ablett EM, Lee LS, Henry RJ: Analysis of SSR-derived from grape ESTs . Theor Appl Genet. 2000, 100: 723-726. 10.1007/s001220051344.

    Article  CAS  Google Scholar 

  26. Sorrells ME: The evolution of comparative plant genetics: Columbia, MO.Edited by: Gustafson JP. 1998, ,

    Google Scholar 

  27. Cooke DEL, Drenth A, Duncan JM, Wagels G, Brasier CM: A Molecular Phylogeny of Phytophthora and Related Oomycetes. Fungal Genetics and Biology. 2000, 30: 17-32. 10.1006/fgbi.2000.1202.

    Article  PubMed  CAS  Google Scholar 

  28. Torto-Alalibo T, Tian M, Gajendran K, Waugh ME, van West P, Kamoun S: Expressed sequences tags from the oomycete fish pathogen Saprolegnia parasitica reveal putative virulence factors. BMC Microbiology. 2005, 5: 46-10.1186/1471-2180-5-46.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Kroon LPNM, Bakker FT, van den Bosch GBM, Bonants PJM, Flier WG: Phylogenetic analysis of Phytophthora species based on mitochondrial and nuclear DNA sequences. Fungal Genetics and Biology. 2004, 41: 766-782. 10.1016/j.fgb.2004.03.007.

    Article  PubMed  CAS  Google Scholar 

  30. Phytophthora Microsatellites Database [http://bioinf.ibun.unal.edu.co/phytossr].

  31. Phytophthora Functional Genomics Database [http://www.pfgd.org/].

  32. Department of Energy's Joint Genome Institute [http://www.jgi.doe.gov/].

  33. RepBase [http://www.girinst.org].

  34. La Rota M, Kantety RV, Yu JK, Sorrells ME: Nonrandom distribution and frequencies of genomic and EST-derived microsatellite markers in rice, wheat, and barley. BMC Genomics. 2005, 6: 23-10.1186/1471-2164-6-23.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Prasad MD, Muthulakshmi M, Madhu M, Archak S, Mita K, Nagaraju J: Survey and analysis of microsatellites in the silkworm, Bombyx mori: frequency, distribution, mutations, marker potential and their conservation in heterologous species. Genetics. 2005, 169 (1): 197-214. 10.1534/genetics.104.031005.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  36. Katti MV, Ranjekar PK, Gupta VS: Differential distribution of simple sequence repeats in eukaryotic genome sequences. Mol Biol Evol. 2001, 18 (7): 1161-1167.

    Article  PubMed  CAS  Google Scholar 

  37. Cardle L, Ramsay L, Milbourne D, Macaulay M, Marshall D, Waugh R: Computational and experimental characterization of physically clustered simple sequence repeats in plants. Genetics. 2000, 156 (2): 847-854.

    PubMed  CAS  PubMed Central  Google Scholar 

  38. Kazuka DNA Research Institute [http://www.kazusa.or.jp/codon].

  39. Harr B, Schlotterer C: Long microsatellite alleles in Drosophila melanogaster have a downward mutation bias and short persistence times, which cause their genome-wide underrepresentation. Genetics. 2000, 155 (3): 1213-1220.

    PubMed  CAS  PubMed Central  Google Scholar 

  40. Kamoun S: A Catalogue of the Effector Secretome of Plant Pathogenic Oomycetes. Annu Rev Phytopathol. 2006

    Google Scholar 

  41. Rozen S, Skaletsky H: Primer3 on the WWW for general users and for biologist programmers. Methods Mol Biol. 2000, 132: 365-386.

    PubMed  CAS  Google Scholar 

  42. Cho YG, Ishii T, Temnykh S, Chen X, Lipovich L, McCouch SR, Park WD, Ayres N, Cartinhour S: Diversity of microsatellites derived from genomic libraries and GenBank sequences in rice (Oryza sativa L.). Theor Appl Genet. 2000, 100: 713-722. 10.1007/s001220051343.

    Article  CAS  Google Scholar 

  43. van Zijll de Jong E, Guthridge KM, Spangenberg GC, Forster JW: Development and characterization of EST-derived simple sequence repeat (SSR) markers for pasture grass endophytes. Genome. 2003, 46 (2): 277-290. 10.1139/g03-001.

    Article  PubMed  CAS  Google Scholar 

  44. Yue GH, Ho MY, Orban L, Komen J: Microsatellites within genes and ESTs of common carp and their applicability in silver crucian carp. Aquaculture. 2004, 234: 85-98. 10.1016/j.aquaculture.2003.12.021.

    Article  CAS  Google Scholar 

  45. Dokholyan NV, Buldyrev SV, Havlin S, Stanley HE: Distributions of dimeric tandem repeats in non-coding and coding DNA sequences. J Theor Biol. 2000, 202 (4): 273-282. 10.1006/jtbi.1999.1052.

    Article  PubMed  CAS  Google Scholar 

  46. Metzgar D, Bytof J, Wills C: Selection against frameshift mutations limits microsatellite expansion in coding DNA. Genome Res. 2000, 10 (1): 72-80.

    PubMed  CAS  PubMed Central  Google Scholar 

  47. Wierdl M, Dominska M, Petes TD: Microsatellite instability in yeast: dependence on the length of the microsatellite. Genetics. 1997, 146 (3): 769-779.

    PubMed  CAS  PubMed Central  Google Scholar 

  48. Kruglyak S, Durrett RT, Schug MD, Aquadro CF: Equilibrium distributions of microsatellite repeat length resulting from a balance between slippage events and point mutations. Proc Natl Acad Sci U S A. 1998, 95 (18): 10774-10778. 10.1073/pnas.95.18.10774.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  49. Kruglyak S, Durrett RT, Schug MD, Aquadro CF: Distribution and abundance of microsatellites in the yeast genome can be explained by a balance between slippage events and point mutations. Molecular Biology and Evolution. 2000, 17: 1210-1219.

    Article  PubMed  CAS  Google Scholar 

  50. Bachtrog D, Weiss S, Zangerl B, Brem G, Schlotterer C: Distribution of dinucleotide microsatellites in the Drosophila melanogaster genome. Mol Biol Evol. 1999, 16 (5): 602-610.

    Article  PubMed  CAS  Google Scholar 

  51. Karaoglu H, Lee CMY, Meyer W: Survey of simple sequence repeats in completed fungal genomes. Molecular Biology and Evolution. 2004, 22: 639-649. 10.1093/molbev/msi057.

    Article  PubMed  Google Scholar 

  52. O'Brien SJ, Womack JE, Lyons LA, Moore KJ, Jenkins NA, Copeland NG: Anchored reference loci for comparative genome mapping in mammals. Nat Genet. 1993, 3 (2): 103-112. 10.1038/ng0293-103.

    Article  PubMed  Google Scholar 

  53. Groppe K, Sanders I, Wiemken A, Boller T: A microsatellite marker for studying the ecology and diversity of fungal endophytes (Epichloe spp.) in grasses. Appl Environ Microbiol. 1995, 61 (11): 3943-3949.

    PubMed  CAS  PubMed Central  Google Scholar 

  54. Cordeiro GM, Casu R, McIntyre CL, Manners JM, Henry RJ: Microsatellite markers from sugarcane (Saccharum spp.) ESTs cross transferable to erianthus and sorghum. Plant Sci. 2001, 160 (6): 1115-1123. 10.1016/S0168-9452(01)00365-X.

    Article  PubMed  CAS  Google Scholar 

  55. Dieringer D, Schlotterer C: Two distinct modes of microsatellite mutation processes: evidence from the complete genomic sequences of nine species. Genome Research. 2003, 13: 2242-2251. 10.1101/gr.1416703.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  56. Borstnik B, Pumpernik D: Tandem repeats in protein coding regions of primate genes. Genome Res. 2002, 12 (6): 909-915. 10.1101/gr.138802.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  57. Subramanian S, Madgula VM, George R, Mishra RK, Pandit MW, Kumar CS, Singh L: Triplet repeats in human genome: distribution and their association with genes and other genomic regions. Bioinformatics. 2003, 19 (5): 549-552. 10.1093/bioinformatics/btg029.

    Article  PubMed  CAS  Google Scholar 

  58. Powell W, Machray GS, Provan J: Polymorphism revealed by simple sequence repeats. Trends In Plant Science. 1996, 7: 215-222.

    Article  Google Scholar 

  59. Katti MV, Sami-Subbu R, Ranjekar PK, Gupta VS: Amino acid repeat patterns in protein sequences: their diversity and structural-functional implications. Protein Sci. 2000, 9 (6): 1203-1209.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  60. Woodhead M, Russell J, Squirrell J, Hollingsworth PM, Cardle L, Ramsay L, Gibby M, Powell W: Development of EST-SSRs from the Alpine lady-fern, Athyrium distentifolium. Molecular Ecology Notes. 2003, 3: 287-290. 10.1046/j.1471-8286.2003.00427.x.

    Article  CAS  Google Scholar 

  61. Hines PJ, Marx J: The endless race between plant and pathogen. Science. 2001, 292 (5525): 2269-10.1126/science.292.5525.2269.

    Article  PubMed  CAS  Google Scholar 

  62. Hood DW, Deadman ME, Jennings MP, Bisercic M, Fleischmann RD, Venter JC, Moxon ER: DNA repeats identify novel virulence genes in Haemophilus influenzae. Proc Natl Acad Sci U S A. 1996, 93 (20): 11121-11125. 10.1073/pnas.93.20.11121.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  63. National Center for Biotechnology Information [http://www.ncbi.nih.gov].

  64. Virginia Bioinformatics Institute [http://phytophthora.vbi.vt.edu].

  65. UniVec [http://ftp.ncbi.nih.gov/pub/UniVec].

  66. Plant Genome Resources Center (PGRC) [http://pgrc.ipk-gatersleben.de].

  67. Zar J: Biostatistical analysis . 1999, Upper Saddle River, NJ , Prentice Hall Press, 4th

    Google Scholar 

Download references

Acknowledgements

We would like to thank Sophien Kamoun, Valérie Verdier, Angela Chaparro and Christine Smart for critical reading of the manuscript. Funding was provided to DG by the Facultad de Ciencias of Universidad de los Andes.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Silvia Restrepo.

Additional information

Authors' contributions

DPG performed the data analysis and wrote the manuscript. AMP carried out all the programming activities and designed the database. LMQ, AJB and EB contributed to the conception and design of the study. AJB contributed to the writing of the manuscript. SR, EB and NJG directed and oriented the project and revised the manuscript. All authors read and approved the final manuscript.

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Garnica, D.P., Pinzón, A.M., Quesada-Ocampo, L.M. et al. Survey and analysis of microsatellites from transcript sequences in Phytophthora species: frequency, distribution, and potential as markers for the genus. BMC Genomics 7, 245 (2006). https://doi.org/10.1186/1471-2164-7-245

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2164-7-245

Keywords